Hunch.h

#ifndef gridripper_phys_gr_fixmp_higgsmp_Hunch_h
#define gridripper_phys_gr_fixmp_higgsmp_Hunch_h


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


namespace gridripper { namespace phys { namespace gr { namespace fixmp {
namespace higgsmp {


using namespace gridripper;
using namespace gridripper::util;
using namespace gridripper::metric;
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.

    GReal_t amplitudeF, centerF, widthF, rampF;
    int lF, mF;
    GReal_t amplitudeFtau, centerFtau, widthFtau, rampFtau;
    int lFtau, mFtau;

    Metric *metric;
    MetricValues g;

    GReal_t tau0;

    GReal_t kappa;

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

    string getMessage() const;
};


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


#endif /* gridripper_phys_gr_fixmp_higgsmp_Hunch_h */