Disk.h

#ifndef gridripper_phys_gr_fixmp_kerrhiggs_Disk_h
#define gridripper_phys_gr_fixmp_kerrhiggs_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 kerrhiggs {


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;
    int mInit;
    int fieldDiscret;
    GReal_t normLimit;
    GReal_t tailErrorTolerance;

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

    GReal_t minrh;
    GReal_t maxrh;

    GReal_t t0;

    string normName;
    enum norm_t { None, L2, C0, Bound };
    enum norm_t Norm;
#ifdef HAVE_AGM
    GReal_t agmTolerance;
#else
    tvector<GReal_t> theta;
    tvector<GReal_t> phih;
    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;

    mutable 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 phih) const;
    GComplex_t funcFt(const GReal_t theta, const GReal_t phih) const;

    string getMessage() const;

};


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


#endif /* gridripper_phys_gr_fixmp_kerrhiggs_Disk_h */