EMM_GaussLegendreRelaxation.h

Go to the documentation of this file.

#ifndef EMM_GaussLegendreRelaxation_H
#define EMM_GaussLegendreRelaxation_H

#include "MATH_GaussLegendreIntegration.h"

#include "EMM_LandauLifschitzRelaxation.h"


DEFINE_SPTR(EMM_GaussLegendreRelaxation);
DEFINE_SVPTR(EMM_GaussLegendreRelaxation);
class EMM_GaussLegendreRelaxation : public  EMM_LandauLifschitzRelaxation {
  
    SP_OBJECT(EMM_GaussLegendreRelaxation);
    // ATTRIBUTES
  
public:
   
    
private:

    //GL integration 
    SP::MATH_GaussLegendreIntegration mIntegrator;

    //local field for GL integration
    SP::EMM_RealField mMs;
    SP::EMM_RealField mFs;

   
    
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    EMM_GaussLegendreRelaxation(void);
  
  
  
    // DESTRUCTORS 
  
    
    virtual ~EMM_GaussLegendreRelaxation(void);   
  
  
public:

    inline static SP::EMM_GaussLegendreRelaxation New() {
        SP::EMM_GaussLegendreRelaxation sp(new EMM_GaussLegendreRelaxation(),EMM_GaussLegendreRelaxation::Delete());
        sp->setThis(sp);
        return sp;
    }

    virtual void toDoAfterThisSetting();
    

      
    
    // SET methods 
    //============
  
    // GET methods
    // ============
public:
    inline const EMM_RealField& getIntermediateMagnetizationField()  const {
        return *mMs.get();
    }
    inline const EMM_RealField& getIntermediateMagnetizationTimeDerivativeField()  const {
        return *mFs.get();
    }

    inline EMM_RealField& getIntermediateMagnetizationField()  {
        return *mMs.get();
    }
   
    inline EMM_RealField& getIntermediateMagnetizationTimeDerivativeField()   {
        
        return *mFs.get();
    }
    
  
    
public:
 
    virtual tBoolean discretize();

    
    virtual tBoolean integrateMagnetizationFieldAtTime(const tReal& t,
                                                       const tFlag& order,
                                                       const EMM_RealArray& sigma,
                                                       const EMM_RealField& dM_dt0,
                                                       const EMM_RealField& M0,
                                                       EMM_RealField& M);
    
   

 
private:                                             
    tBoolean computeMagnetizationFieldAtTimeWithGL1Interpolation(const tReal& dt,
                                                                 const EMM_RealArray& sigma,
                                                                 const EMM_RealField& dM_dt0,
                                                                 const EMM_RealField& M0,
                                                                 EMM_RealField& Mt);
    
    tBoolean computeMagnetizationFieldAtTimeWithGLnInterpolation(const tReal& dt,
                                                                 const EMM_RealArray& sigma,
                                                                 const EMM_RealField& dM_dt0,
                                                                 const EMM_RealField& M0,
                                                                 EMM_RealField& Mt);
    
};

#endif