Hunch.h

#ifndef gridripper_phys_Hunch_h
#define gridripper_phys_Hunch_h

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

namespace gridripper { namespace phys {

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

class Hunch: public FuncInitCond
{
private:
    tvector<int> componentIndices;

protected:
    GReal_t x0;
    GReal_t deltaX;
    GReal_t c;
    GReal_t d;

private:
    void init(string& args, const Parameters* p, const PDE& pde,
         const string& x0name, const string& deltaXname)
            throw(IllegalArgumentException&);

protected:
    Hunch(string& args, const Parameters* p, const PDE& pde,
          const string& x0name, const string& deltaXname)
            throw(InitCond::Exception&, IllegalArgumentException&);

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

    GReal_t getX0() const {
        return x0;
    }

    virtual GReal_t function(GReal_t x) const;

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

} } // namespace gridripper::phys;

#endif /* gridripper_phys_Hunch_h */