ODERelaxInit.h

#ifndef gridripper_amr1d_initcond_ODERelaxInit_h
#define gridripper_amr1d_initcond_ODERelaxInit_h

#include <gridripper/amr1d/FuncInitCond.h>
#include <gridripper/util/GrrProperties.h>
#include <gridripper/Parameters.h>
#include <gridripper/odesolver/ODE.h>

namespace gridripper { namespace amr1d { namespace initcond {

using namespace gridripper;
using namespace gridripper::util;
using namespace gridripper::amr1d;
using namespace gridripper::odesolver;

class ODERelaxInit: public FuncInitCond
{
private:
    GrrProperties properties;

    tvalarray<GReal_t> data;

    GReal_t minX;

    GReal_t maxX;

    int maxi;

    ODE* ode;

    string message;

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

public:
    string getMessage() {
        return message;
    }

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

    void endInit();

private:
    void init(const PDE* pde) throw(InitCond::Exception&);
};

} } } // namespace gridripper::amr1d::initcond

#endif /* gridripper_amr1d_initcond_ODERelaxInit_h */