EMM_MagneticExcitationOperator.h

Go to the documentation of this file.

#ifndef EMM_MagneticExcitationOperator_H
#define EMM_MagneticExcitationOperator_H

#include "EMM_Operator.h"

DEFINE_SPTR(EMM_MagneticExcitationOperator);
DEFINE_SVPTR(EMM_MagneticExcitationOperator);
class EMM_MagneticExcitationOperator : public EMM_Operator{
  
    SP_OBJECT(EMM_MagneticExcitationOperator);
    // ATTRIBUTES
    
public:
    
    
    
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    EMM_MagneticExcitationOperator(void): EMM_Operator() {
    }
  
  
  
    // DESTRUCTORS 
  
    
    virtual ~EMM_MagneticExcitationOperator(void) {
    }
  
  
    // SET methods 
 
   
public:
   
    
    //Magnetic excitation field methods
    //=================================
    
    virtual void  computeMagneticExcitationField(const EMM_RealArray& sigma,const EMM_RealField& M,EMM_RealField& H) const=0;
    
    virtual tBoolean computeMagneticExcitationFieldGradient(const EMM_RealArray& sigma,
                                                            const EMM_RealField& M,
                                                            const EMM_RealField& D,
                                                            EMM_RealField& gradH) const=0;
    
  


    //ENERGY Methods
    //==============
    
    inline tReal computeMagneticExcitationFieldAndEnergy(const tReal& t,
                                                         const EMM_RealArray& sigma,
                                                         const EMM_RealField& Mt,
                                                         EMM_RealField& Ht) const {
        computeMagneticExcitationField(sigma,Mt,Ht);
        return computeEnergyWithMagneticExcitation(t,sigma,Mt,Ht);
    }
    

protected:
  
    virtual tReal computeEnergyWithMagneticExcitation(const tReal& t,
                                                      const EMM_RealArray& sigma,
                                                      const EMM_RealField& Mt,
                                                      const EMM_RealField& Ht) const=0;
public:

    virtual tReal computeEnergy(const tReal& t,
                                const EMM_RealArray& sigma,
                                const EMM_RealField& Mt,
                                EMM_RealField& W) const=0;
    
    
    
  
};

#endif