Hunch.h

#ifndef gridripper_phys_gr_fixmp_kerrhiggs_Hunch_h
#define gridripper_phys_gr_fixmp_kerrhiggs_Hunch_h


#include <gridripper/Parameters.h>
#include <gridripper/amr1d/FuncInitCond.h>
#include <gridripper/amr1d/PDE.h>
#include <gridripper/amr1d/FuncInitCond.h>


namespace gridripper { namespace phys { namespace gr { namespace fixmp {
namespace kerrhiggs {


using namespace gridripper;
using namespace gridripper::util;
using namespace gridripper::amr1d;
using namespace std;


class Hunch: public FuncInitCond
{
private:
    int maxOrder;
    int arraySize;

    bool selfInteraction;
    GReal_t Phi0;
    // The field is f=(Phi-Phi0)*r if selfinteraction is on, 
    // otherwise f=Phi*r.

    bool gaussianHunch;
    GReal_t amplitudeF, centerF, widthF, rampF;
    int lF, mF;
    GReal_t amplitudeFt, centerFt, widthFt, rampFt;
    int lFt, mFt;
    GReal_t omega;

    GReal_t M;
    GReal_t a;
    GReal_t e;
    bool isMinkowski;

    GReal_t minrh;
    GReal_t maxrh;

    GReal_t t0;

    string message;

public:
    Hunch(string& args, const Parameters* p, const PDE& pde)
          throw(InitCond::Exception&, IllegalArgumentException&);
    ~Hunch();
    GReal_t getTime() const;

    void init(PDE* pde, GReal_t rh, FieldWrapper& v)
        throw(InitCond::Exception&);

    string getMessage() const;
};


} } } } } // namespace gridripper::phys::gr::fixmp::kerrhiggs


#endif /* gridripper_phys_gr_fixmp_kerrhiggs_Hunch_h */