CompactKScODE.h

#ifndef gridripper_phys_gr_dynss_ekg_kodama_CompactKScODE_h
#define gridripper_phys_gr_dynss_ekg_kodama_CompactKScODE_h

#include <gridripper/odesolver/ODE.h>
#include <gridripper/amr1d/FieldComponents.h>
#include <gridripper/lang/IllegalArgumentException.h>

namespace gridripper {
class Parameters;
}

namespace gridripper { namespace phys { namespace gr { namespace dynss {
namespace ekg { namespace kodama {

using namespace gridripper;
using namespace gridripper::odesolver;
using namespace std;

class CompactKScODE: public ODE
{
private:
    // Field indices:
    enum {
//      I_Phi = 0,
//      I_Pi  = 1,
        I_B = 0,
        I_psi = 1,
        NUM_COMPS = 2};

    static const string COMPNAMES[NUM_COMPS];

    GReal_t rhomax;

    GReal_t a;
    GReal_t b;
    GReal_t c;
    GReal_t d;
    GReal_t mass;
    GReal_t psi0;

public:
    class Field: public FieldComponents<GReal_t>
    {
    public:
        Field(): FieldComponents<GReal_t>(NUM_COMPS) { }
        Field(const Field& f): FieldComponents<GReal_t>(f) { }
        FieldWrapper* cloneFieldWrapper() const {
            return new Field(*this);
        }
    };

    CompactKScODE(const Parameters* p) throw(IllegalArgumentException&);

    FieldArray getInitialValuesForShooting() const;

    tvalarray<GReal_t> getInitialX() const {
        return tvalarray<GReal_t>(0.0, 1);
    }

    GReal_t getMinX() const {
        return 0;
    }

    GReal_t getMaxX() const {
        return rhomax;
    }

    void initPDE(GReal_t rho, const GReal_t* fode, GReal_t* fpde);

    void eval(const GReal_t* F, int offset, GReal_t rho, GReal_t* dF);

private:
    GReal_t calc_phi(GReal_t rho) const;
};

} } } } } } // namespace gridripper::phys::gr::dynss::ekg::kodama;

#endif /* gridripper_phys_gr_dynss_ekg_kodama_CompactKScODE_h */