base/diagOuterSoftAbsMetric.h

#ifndef _BETA_DIAGOUTERSOFTABSMETRIC_

#include <Eigen/Dense>

#include "dynamMetric.h"


class diagOuterSoftAbsMetric: public dynamMetric
{
    
    public:
        
        explicit diagOuterSoftAbsMetric(int dim);
        virtual ~diagOuterSoftAbsMetric() {};
        
        //                   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;   
    
        MatrixXd mH;            
        
        VectorXd mLambda;       
        VectorXd mGradHelper;   

        virtual void fComputeH() = 0;
        
        void fComputeMetric();
        void fPrepareSpatialGradients() { fComputeH(); gradV(); }
        
        VectorXd& dTaudp();
        VectorXd& dTaudq();
        VectorXd& dPhidq();
    
};

#define _BETA_DIAGOUTERSOFTABSMETRIC_
#endif