ODEShoot.h

#ifndef gridripper_odesolver_shoot_ODEShoot_h
#define gridripper_odesolver_shoot_ODEShoot_h

#include <gridripper/config.h>
#include <gridripper/lang/IllegalArgumentException.h>
#include <gridripper/amr1d/FieldWrapper.h>
#include <gridripper/util/GrrProperties.h>
#include <gridripper/tvector.h>
#include <string>

namespace gridripper { namespace odesolver {

class ODE;

namespace shoot {

class ODEIntegrator;

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

class ODEShoot
{
private:
    ODE* theODE;
    GReal_t deltaX;
    int subSteps;
    ODEIntegrator* odeIntegrator;
    string dataFormat;

public:
    ODEShoot(const GrrProperties& props)
        throw(IllegalArgumentException&);

    ODEShoot(const GrrProperties& props, GReal_t dx, int subSteps)
        throw(IllegalArgumentException&);

    ODEShoot(ODE* ode, GReal_t dx, int dn, ODEIntegrator* oi)
            throw(IllegalArgumentException&);

public:
    ODE& getODE() const {
        return *theODE;
    }

    GReal_t getDeltaX() const {
        return deltaX;
    }

    void solve(tvector<GReal_t>& sol) const;

    string getDataFormat() const {
        return dataFormat;
    }

private:
    void solve(tvector<GReal_t>& sol,
               GReal_t x0, const FieldWrapper& f0,
               GReal_t x1, const FieldWrapper& f1) const;

    void solve(tvector<GReal_t>& sol, GReal_t x0, const FieldWrapper& fw0)
        const;

    void init(ODE* ode, GReal_t dx, int dn, ODEIntegrator* oi)
            throw(IllegalArgumentException&);
};

} } } // namespace gridripper::odesolver::shoot

#endif /* gridripper_odesolver_shoot_ODEShoot_h */