base/denseFisherMetric.h

#ifndef _BETA_DENSEMETRIC_

#include <Eigen/Cholesky>

#include "dynamMetric.h"


class denseFisherMetric: public dynamMetric
{
    
    public:
        
        explicit denseFisherMetric(int dim);
        virtual ~denseFisherMetric() {};
        
        //                   Accessors                  //
        
        double T();
        double tau();

        void bounceP(const VectorXd& normal);
    
        MatrixXd& G() { return mG; }
    
        LLT<MatrixXd>& GL() { return mGL; }
        
        //                   Mutators                   //
        
        void sampleP(Random& random);
        
        //              Auxiliary Functions             //
        
        void checkEvolution(const double epsilon = 1e-6);
        
        void displayState();
        
        void prepareEvolution() { fComputeCholeskyG(); }
        
    protected:
    
        MatrixXd mG;          
        MatrixXd mGradG;      
        LLT<MatrixXd> mGL;    
        VectorXd mC;          

        void fComputeMetric() {  fComputeCholeskyG(); }
    
        virtual void fComputeG() = 0;
    
        virtual void fComputeGradG(int i) = 0;
    
        void fComputeCholeskyG();
    
        VectorXd& dTaudp();
        VectorXd& dTaudq();
        VectorXd& dPhidq();
    
};
 
#define _BETA_DENSEFISHERMETRIC_
#endif