BData.h

#ifndef gridripper_amr1d_initcond_BData_h
#define gridripper_amr1d_initcond_BData_h

#include <gridripper/io/BDataReader.h>
#include <gridripper/io/IOException.h>
#include <gridripper/amr1d/InitCond.h>
#include <gridripper/util/GrrProperties.h>
#include <gridripper/Parameters.h>
#include <gridripper/amr1d/interpolation/DataInterpolation.h>
#include <gridripper/lang/IntArray.h>

namespace gridripper { namespace amr1d {

class PDE;

namespace initcond {

using namespace gridripper;
using namespace gridripper::util;
using namespace gridripper::io;
using namespace gridripper::amr1d::interpolation;

class BData : public InitCond
{

private:
    static const GReal_t EPSILON=1e-10;
    bool debug;
    BDataReader* bData;
    DataInterpolation* dint;
    DataInterpolation* dint2;
    std::string inputFile;
    tvalarray<GReal_t>* energyLosses;
    int xMultiplier;
    GReal_t theTime0;
    GReal_t theTime;
    tvalarray<int>* kxArray;
    tvalarray<GReal_t> theData;
    tvalarray<IntArray>* timeStepCounts;
    tvalarray<IntArray>* lastRegridTimes;
    std::string message;

public:
    BData(const std::string &args, const Parameters* p, const PDE& pde)
          throw(IOException&);

    ~BData();

    bool isOnInitialHypersurface() const;

    bool isAppended(const std::string &f) const;

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

    GReal_t getTime() const;

    std::string getMessage() const;

private:
    int init0(tvalarray< tvector<Grid*> >& grids, PDE* pde, GReal_t dt, int r,
              int maxlevel);

    static int gridSpacingAt(const tvalarray<int>& kxarr, int i);

    int createInitialSubgrids(tvalarray< tvector<Grid*> >& grids, Grid* g,
                              PDE* pde, GReal_t dt, int r,
                              int finestl, int maxlevel, int dkx, int ikx0);

    void init(Grid& g, int ikx, int dkx, PDE* pde);

};


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

#endif /* gridripper_amr1d_initcond_BData_h */