base/diagSoftAbsMetric.h

#ifndef _BETA_DIAGSOFTABSMETRIC_

#include <Eigen/Dense>

#include "dynamMetric.h"


class diagSoftAbsMetric: public dynamMetric
{
    
    public:
        
        explicit diagSoftAbsMetric(int dim);
        virtual ~diagSoftAbsMetric() {};
        
        //                   Accessors                  //
        
        double T();
        double tau();
        
        void bounceP(const VectorXd& normal);
        
        //                   Mutators                   //
        
        void sampleP(Random& random);
        
        void setSoftAbsAlpha(double a) { mSoftAbsAlpha = a; }
        
        //              Auxiliary Functions             //
        
        void checkEvolution(const double epsilon = 1e-6);
        
        void displayState();
        
    protected:
        
        double mSoftAbsAlpha;   
        
        VectorXd mDiagH;        
        MatrixXd mGradDiagH;    
        
        VectorXd mLambda;       
        VectorXd mGradHelper;   
        
        virtual void fComputeDiagH() = 0;
        
        virtual void fComputeGradDiagH() = 0;
        
        void fComputeMetric();
        void fPrepareSpatialGradients() { fComputeGradDiagH(); gradV(); }

        VectorXd& dTaudp();
        VectorXd& dTaudq();
        VectorXd& dPhidq();
    
};

#define _BETA_DIAGSOFTABSMETRIC_
#endif