Disk.h

#ifndef gridripper_phys_gr_fixmp_kerrkgaxhyp_Disk_h
#define gridripper_phys_gr_fixmp_kerrkgaxhyp_Disk_h


#include <gridripper/Parameters.h>
#include <gridripper/amr1d/PDE.h>
#include <gridripper/amr1d/FuncInitCond.h>
#include <gridripper/math.h>
#include <sstream>


namespace gridripper { namespace phys { namespace gr { namespace fixmp {
namespace kerrkgaxhyp {


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


class Disk: public FuncInitCond
{
private:
    int maxOrder;
    int maxOrderInit;
    int arraySize;

    bool selfInteraction;
    GReal_t phi0;
    // The field is f=(phi-phi0)*r if selfinteraction is on, 
    // otherwise f=phi*r.

    GReal_t amplitudeF, amplitudeFt;
    GReal_t centerxF, centerxFt;
    GReal_t widthxF, widthxFt;
    GReal_t rampxF, rampxFt;
    GReal_t centerzF, centerzFt;
    GReal_t widthzF, widthzFt;
    GReal_t rampzF, rampzFt;
    GReal_t omegaInit, mInit;
    int fieldDiscret;
    GReal_t normLimit;
    GReal_t tailErrorTolerance;

    GReal_t M;
    GReal_t a;
    GReal_t e;
    bool isMinkowski;

    GReal_t minR;
    GReal_t maxR;

    GReal_t T0;

    enum norm_T { None, L2, C0, Bound };
    enum norm_T Norm;
#ifndef HAVE_AGM
    tvector<GReal_t> theta;
    tvector<GReal_t> phi;
    tvector< tvector<GComplex_t> > field;
    tvector< tvector<GComplex_t> > fieldt;
    tvector< tvector<GReal_t> > fieldBuff;
    tvector< tvector<GReal_t> > fieldTBuff;
    tvector< tvector<GReal_t> > fieldBuffDiff;
    tvector< tvector<GReal_t> > fieldtBuffDiff;
    tvector<GReal_t> intdphiF;
    tvector<GReal_t> intdphiFt;
    tvector<GReal_t> intdphiDiffF;
    tvector<GReal_t> intdphiDiffFt;
#endif

    GReal_t t, r;

    unsigned int tstart, tnow;

    GReal_t maxNormF;
    GReal_t maxNormFt;
    GReal_t maxNormDiffF;
    GReal_t maxNormDiffFt;
    
    ostringstream message;

public:
    Disk(string& args, const Parameters* p, const PDE& pde)
         throw(InitCond::Exception&, IllegalArgumentException&);
    ~Disk();
    //GReal_t getTime() const;

    void init(PDE* pde, GReal_t rh, FieldWrapper& v)
        throw(InitCond::Exception&);

    GComplex_t funcF(const GReal_t theta, const GReal_t phi) const;
    GComplex_t funcFt(const GReal_t theta, const GReal_t phi) const;

    //string getMessage() const;

};


} } } } } // namespace gridripper::phys::gr::fixmp::kerrkgaxhyp


#endif /* gridripper_phys_gr_fixmp_kerrkgaxhyp_Disk_h */