AutoAMRCore.h

#ifndef gridripper_amr1d_AutoAMRCore_h
#define gridripper_amr1d_AutoAMRCore_h

#include "AMRCore.h"
#include <gridripper/math/expr/Evaluator.h>
#include <string>
#ifdef HAVE_GSTREAM
#include <gstream/gstream.h>
#else
#include <fstream>
#endif

namespace gridripper { namespace amr1d {

using namespace std;
using namespace gridripper::math::expr;

class AutoAMRCore: public AMRCore
{
public:
    class OutFileExistsException: public gridripper::Exception
    {
    public:
        OutFileExistsException(const Exception& other):
            gridripper::Exception(other) { }

        OutFileExistsException(const string& fname):
            gridripper::Exception(fname + ": file exists",
                                  "AutoAMRCore::AutoAMRCore") { }
    };
private:
    Evaluator::Function* dtwriteFunc;
    GReal_t dtwriteValue;

    int blankLines;

    int fileNo;
    int nfile;
    string outputFile;
    int maxBDataSize;

    bool outputAppended;

    bool endless;

    int initialTimeStepCount;

    GReal_t prevSaveTime;

    int prevSaveTimeStep;

    GReal_t theLastTime;

    int msecsPerDisplay;

    int nextTimeStepToDisplay;

    double calcTimeWAvg;

    long calcTimeTot;

    long waitTimeTot;

    string dataFormat;
    string energyFormatStr;

    unsigned startTime;
    int maxRuntime;

    bool binary;

#ifdef HAVE_GSTREAM
    ogstream outputFileStream;
#else
    ofstream outputFileStream;
#endif

public:
    AutoAMRCore(PDE* pde, InitCond& ini, bool constrforced, Integrator* gi,
                Regridder* regridder, const string& dtdx,
                int maxi, int npabovemax,
                int maxl, int r,
                AMError* amerr, double maxerr,
                int errchkf, int regridf,
                int bufzoner,
                int fms, string f, int fileno, int bdatasize, int blines,
                const string& dtwrite, bool endless,
                GReal_t tmax, int maxruntime, 
                const string& datafmt, const string& energyfmt, bool enable)
            throw(InitCond::Exception&, IllegalArgumentException&,
                  OutFileExistsException&);

    virtual ~AutoAMRCore();

    AMRCore* cloneAMRCore() const;

    long getCalcTime() const {
        return calcTimeTot;
    }

    long getWaitTime() const {
        return waitTimeTot;
    }

    bool isOutputAppended() const {
        return outputAppended;
    }

    bool isTimeSliceSaved(int tsc) const;

    void output();

protected:
    bool isLastTime(const Grid& g) const;

    virtual void coarseStepped(GReal_t t);
};

} } // namespace gridripper::amr1d

#endif /* gridripper_amr1d_AutoAMRCore_h */