EMM_MagneticExcitationLinearOperator.h

Go to the documentation of this file.

#ifndef EMM_MagneticExcitationLinearOperator_H
#define EMM_MagneticExcitationLinearOperator_H

#include "EMM_MagneticExcitationOperator.h"


DEFINE_SPTR(EMM_MagneticExcitationLinearOperator);
class EMM_MagneticExcitationLinearOperator : public virtual EMM_MagneticExcitationOperator {
  
    SP_OBJECT(EMM_MagneticExcitationLinearOperator);
    // ATTRIBUTES

  
private:
  
    
    
public:
   

    // ASSOCIATION
 
private:
   
 
protected: 

    // METHODS
  
    // CONSTRUCTORS 
  
    EMM_MagneticExcitationLinearOperator(void) {
    }
  
    // DESTRUCTORS 
    
    virtual ~EMM_MagneticExcitationLinearOperator(void) {
    }
  
  
public:
   
    virtual tBoolean isAffine() const {
        return true;
    };
    virtual tBoolean isGradientComputationable() const {
        return true;
    };


    //FIELDS Methods
    //==============
    virtual tBoolean resetToInitialState(const EMM_LandauLifschitzSystem& system) {
        return true;
    }
    
    virtual tBoolean computeFieldsAtTime(const tReal& t,
                                         const tFlag& order,
                                         const EMM_RealArray& sigma,
                                         const EMM_RealField& dM_dt0,
                                         const EMM_RealField& M0) {
        return true;
    }
    
   
    virtual tBoolean updateAtNextTimeStep(const tReal& dt,
                                          const EMM_RealArray& sigma,
                                          const EMM_RealField& Mt) {
        return true;
    }

    //Magnetic excitation field method
    //=================================
    
    virtual tBoolean computeMagneticExcitationFieldGradient(const EMM_RealArray& sigma,
                                                            const EMM_RealField& M,
                                                            const EMM_RealField& D,
                                                            EMM_RealField& gradH) const {
        computeMagneticExcitationField(sigma,D,gradH);
        return true;
    }
    


    //ENERGY methods
    //==============

protected:
    virtual tReal computeEnergyWithMagneticExcitation(const tReal& t,
                                                      const EMM_RealArray& sigma,
                                                      const EMM_RealField& Mt,
                                                      const EMM_RealField& Ht) const {
        return -getElementVolume()*Ht.dot(sigma,Mt);
    }

    
public:

    virtual tReal computeEnergy(const tReal& t,
                                const EMM_RealArray& sigma,
                                const EMM_RealField& Mt,
                                EMM_RealField& W) const {
        return EMM_MagneticExcitationOperator::computeMagneticExcitationFieldAndEnergy(t,sigma,Mt,W);
    }
    
     
  
};

#endif