FuncInitCond.h

#ifndef gridripper_amr1d_FuncInitCond_h
#define gridripper_amr1d_FuncInitCond_h

#include "InitCond.h"
#include "FieldWrapper.h"
#include <gridripper/Parameters.h>

namespace gridripper { namespace amr1d {

class FuncInitCond: public InitCond
{
private:
    bool numericalConstraintSettingNeeded;
    int mostRefinedGridMaxI;
    GReal_t mostRefinedGridDeltaX;

    int interpolatedLevel;

protected:
    FuncInitCond(string& args, const Parameters* p)
        throw(InitCond::Exception&);

public:
    int init(tvalarray< tvector<Grid*> >& grids, PDE* pde, GReal_t dt, int r,
             int maxlevel, int bufzone, AMError* amerr, double maxerr)
        throw(InitCond::Exception&);

    virtual GReal_t getTime() const {
        return 0.0;
    }

    int getInterpolatedLevel() const {
        return interpolatedLevel;
    }

    void setInterpolatedLevel(int l);

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

    int getMostRefinedGridMaxI() const {
        return mostRefinedGridMaxI;
    }

    GReal_t getMostRefinedGridDeltaX() const {
        return mostRefinedGridDeltaX;
    }

    virtual bool isNumericalConstraintSettingNeeded() const {
        return numericalConstraintSettingNeeded;
    }

    void setNumericalConstraintSettingNeeded(bool v) {
        numericalConstraintSettingNeeded = v;
    }

    virtual void endInit() {
    }

private:
    void parseArgs(const Parameters*, string& args)
        throw(IllegalArgumentException&);

    void init(Grid* g, PDE* pde) throw(InitCond::Exception&);

    int createInitialSubgrids(tvalarray< tvector<Grid*> >& grids, Grid* g,
                              PDE* pde, GReal_t dt, int r,
                              int finestl, int maxlevel, int bufzone,
                              AMError* amerr, double maxerr)
        throw(InitCond::Exception&);
};

} } // namespace gridripper::amr1d

#endif /* gridripper_amr1d_FuncInitCond_h */