base/dynamMetric.h

#ifndef _BETA_DYNAMMETRIC_

#include "baseHamiltonian.h"


class dynamMetric: public baseHamiltonian
{
    
    public:
        
        explicit dynamMetric(int dim);
        virtual ~dynamMetric() {};
    
        //                   Accessors                  //
        
        virtual double tau() = 0;
    
        double phi() { return 0.5 * mLogDetMetric + V(); }
        
        //                   Mutators                   //
        
        void setMaxNumFixedPoint(int n) { mMaxNumFixedPoint = n; }
    
        void setFixedPointThreshold(double t) { mFixedPointThreshold = t; }
        
        void evolveQ(const double epsilon);
        void beginEvolveP(const double epsilon);
        void finishEvolveP(const double epsilon);  
    
        //              Auxiliary Functions             //

        virtual void prepareEvolution() { fComputeMetric(); fPrepareSpatialGradients(); }
    
    protected:

        int mMaxNumFixedPoint; 
    
        double mFixedPointThreshold;
    
        VectorXd mAuxVector;
    
        VectorXd mInit; 
    
        VectorXd mDelta;
    
        double mLogDetMetric;

        virtual void fComputeMetric() = 0;
    
        virtual void fUpdateP() {};
    
        virtual void fPrepareSpatialGradients() { gradV(); }
        
        virtual void fHatT(const double epsilon);
    
        virtual void fHatTau(const double epsilon, const int numFixedPoint);
    
        void fHatPhi(const double epsilon);
    
        virtual VectorXd& dTaudp() = 0;
    
        virtual VectorXd& dTaudq() = 0;
    
        virtual VectorXd& dPhidq() = 0;
    
};

#define _BETA_DYNAMMETRIC_
#endif