EMM_GradGaussLegendreRelaxation.h

Go to the documentation of this file.

#ifndef EMM_GradGaussLegendreRelaxation_H
#define EMM_GradGaussLegendreRelaxation_H

#include "EMM_GaussLegendreRelaxation.h"


DEFINE_SPTR(EMM_GradGaussLegendreRelaxation);
DEFINE_SVPTR(EMM_GradGaussLegendreRelaxation);
class EMM_GradGaussLegendreRelaxation : public  EMM_GaussLegendreRelaxation {
    
    SP_OBJECT(EMM_GradGaussLegendreRelaxation);
    // ATTRIBUTES
    
public:
   
    
private:
 
    
    //the landau lifschitz grad function F in the direction F = mMs 
    SP::EMM_RealField mGradF_F;
    //the landau lifschitz grad function H in the direction F  
    SP::EMM_RealField mGradH_F;
    
    
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    EMM_GradGaussLegendreRelaxation(void);
  
  
  
    // DESTRUCTORS 
  
    
    virtual ~EMM_GradGaussLegendreRelaxation(void);   
  
  
public:

    inline static SP::EMM_GradGaussLegendreRelaxation New() {
        SP::EMM_GradGaussLegendreRelaxation sp(new EMM_GradGaussLegendreRelaxation(),EMM_GradGaussLegendreRelaxation::Delete());
        sp->setThis(sp);
        return sp;
    }
    
    virtual void toDoAfterThisSetting();
    
    // SET methods 
    //============
  
    virtual tBoolean addOperator(SP::EMM_MagneticExcitationOperator op) {
        if (op.get()==null) return false;
        if (!op->isGradientComputationable()) return false;
        return EMM_GaussLegendreRelaxation::addOperator(op);
    }
   
    // GET methods
    // ============
  
public:
  
    inline const EMM_RealField& getMagneticExcitationFieldGradient()   const {
        return *mGradH_F.get();
    }
    
     
    inline const EMM_RealField& getMagnetizationFieldTimeDerivativeGradient()   const {
        return *mGradF_F.get();
    }
 
  
    virtual tBoolean discretize();

    virtual void computeLandauLifschitzFields(const EMM_RealField& M,
                                               const EMM_RealField& H,
                                              EMM_RealField& F);

public:
    inline void computeMagnetizationFieldTimeDerivativeGradient(const EMM_RealField& M,
                                                                const EMM_RealField& H,
                                                                const EMM_RealField& D,
                                                                const EMM_RealField& GradH,
                                                                EMM_RealField& GradF) const {
        computeGradFunction(M,H,D,GradH,GradF);
    }

    
    void computeMagneticExcitationFieldGradient(const EMM_RealField& M,
                                                const EMM_RealField& D,
                                                EMM_RealField& gradH) const;
    
    tReal computeEnergyTimeDerivativesAtTime(const tReal& t,const tFlag& order,tReal& dE_dt,tReal& d2E_dt2);

    tReal computeEnergyTimeDerivatives(tReal& dE_dt,tReal& d2E_dt2) const ;
   

 
 
    
};

#endif