Hunch.h

#ifndef gridripper_phys_gr_fixmp_kerrkgaxhyp_Hunch_h
#define gridripper_phys_gr_fixmp_kerrkgaxhyp_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 kerrkgaxhyp {


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


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

   

    GReal_t amplitudeF, centerF, widthF, rampF;
    int lF, mF;
    GReal_t amplitudeFt, centerFt, widthFt, rampFt;
    int lFt, mFt;
    int Ym,absYm;
    int parity;
    GReal_t omega;

    GReal_t M;
    GReal_t A;
    //GReal_t e;
    bool isMinkowski;

    GReal_t minR;
    GReal_t maxR;

    GReal_t T0;

    string message;
    
    GReal_t *f, *fT;

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 R, FieldWrapper& v)
        throw(InitCond::Exception&);

    string getMessage() const;
};


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


#endif /* gridripper_phys_gr_fixmp_kerrkgaxhyp_Hunch_h */