EMM_LandauLifschitzODE_RK.h

Go to the documentation of this file.

#ifndef EMM_LandauLifschitzODE_RK_H
#define EMM_LandauLifschitzODE_RK_H


#include "EMM_LandauLifschitzODE.h"

DEFINE_SPTR(EMM_LandauLifschitzODE_RK);
class EMM_LandauLifschitzODE_RK : public  EMM_LandauLifschitzODE {
  
    SP_OBJECT(EMM_LandauLifschitzODE_RK);
    // ATTRIBUTES
  
public:
   
 
   
   
private:
   

   
   
    //fields
    SV::EMM_RealField mFs;
    
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    EMM_LandauLifschitzODE_RK(void);
  
  
  
    // DESTRUCTORS 
  
    
    virtual ~EMM_LandauLifschitzODE_RK(void);   
  
  
public:

    //NEW constructor
protected:
   
    
    //DATA methods
    //=============
public:


  
    void setOrder(const tUSInt& n);
  
    inline tUSInt getOrder() const {
        return mFs.getSize();
    }

    inline const SV::EMM_RealField& getIntermediateFields()  const {
        return mFs;
    }

    virtual const tReal* getACoefficients() const=0;
    virtual const tReal* getBCoefficients() const=0;

    virtual const tReal* getCTCoefficients() const=0;

    virtual const tReal* getCHCoefficients() const=0;
    
    
    //MAGNETIZATION field methods
    //===========================
public:
    virtual tReal integrateMagnetizationFieldAtTime(tReal& dt,
                                                    const tFlag& errorNorm,
                                                    const tReal& epsilon,
                                                    tUSInt& nEvaluations,
                                                    const EMM_RealField& dM_dt0,
                                                    const EMM_RealField& M0,
                                                    EMM_RealField& Mt);


    


};

#endif