wHunch.h

#ifndef gridripper_phys_minkowski_ymh_moncrief_wHunch_h
#define gridripper_phys_minkowski_ymh_moncrief_wHunch_h

#include <gridripper/Parameters.h>
#include <gridripper/amr1d/initcond/ODERelaxInit.h>

namespace gridripper { namespace phys { namespace minkowski { namespace ymh {
namespace moncrief {

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

class wHunch: public gridripper::amr1d::initcond::ODERelaxInit
{
private: int I_h;
    int I_hR;
    int I_w;
    int I_wT;
    int I_wR;
    GReal_t a;
    GReal_t b;
    GReal_t c;
    GReal_t d;
    GReal_t kappa;
    GReal_t lambda;

public:
    wHunch(string& args, const Parameters* p, const PDE& pde)
           throw(InitCond::Exception&, IllegalArgumentException&);

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

protected:
    bool isNumericalConstraintSettingNeeded() {
        return false;
    }
};

} } } } } // namespace gridripper::phys::minkowski::ymh::moncrief;

#endif /* gridripper_phys_minkowski_ymh_moncrief_wHunch_h */