base/outerSoftAbsMetric.h

#ifndef _BETA_OUTERSOFTABSMETRIC_

#include "dynamMetric.h"


class outerSoftAbsMetric: public dynamMetric
{
    
    public:
        
        explicit outerSoftAbsMetric(int dim);
        virtual ~outerSoftAbsMetric() {};
        
        //                   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 prepareEvolution() { gradV(); fComputeH();}
        
    protected:
        
        double mSoftAbsAlpha;        
    
        MatrixXd mH;                 

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


};

#define _BETA_OUTERSOFTABSMETRIC_
#endif