EMMG_DemagnetizedOperator.h

Go to the documentation of this file.

#ifndef EMMG_DemagnetizedOperator_H
#define EMMG_DemagnetizedOperator_H


#include "EMM_DemagnetizedOperator.h"



DEFINE_SPTR(EMMG_DemagnetizedOperator);
class EMMG_DemagnetizedOperator : public virtual EMM_DemagnetizedOperator {
    
    SP_OBJECT(EMMG_DemagnetizedOperator);
    // ATTRIBUTES
    
    
private:

    

  
    
public:
   
   
    
    
protected: 

    // METHODS
  
    // CONSTRUCTORS 
  
    EMMG_DemagnetizedOperator(void) {
    }
  
    // DESTRUCTORS 
    
    virtual ~EMMG_DemagnetizedOperator(void) {
    }
  
  
public:
   
    inline static SP::EMMG_DemagnetizedOperator New() {
        SP::EMMG_DemagnetizedOperator p(new EMMG_DemagnetizedOperator(),
                                          EMMG_DemagnetizedOperator::Delete());
        p->setThis(p);
        return p;
    };
    
    // SET methods
    // ===========

  
    
    // GET methods
    // ===========
  


 

public:
    virtual void computeMagneticExcitationField(const EMM_RealArray& sigma,
                                                const EMM_RealField& M,
                                                EMM_RealField& H) const {  
        // size of the operator
        tUIndex nCells=M.getSize();
        
        //set the size of the magnetic excitation field
        H.setDimension(M.getDimension());
        H.setSize(nCells);
        
        //no computation
        if (nCells==0) return;
        
        //verify the type of H
        const EMMG_RealField *gM=dynamic_cast<const EMMG_RealField *>(&M);
        if (gM==null) throw EMM_Exception("generic","EMMG_DemagnetizedOperator:: computeMagneticExcitationField(..)",
                                          "M is an incompatible real field");
        
        //verify the type of H
        EMMG_RealField *gH=dynamic_cast<EMMG_RealField *>(&H);
        if (gH==null) throw EMM_Exception("generic","EMMG_DemagnetizedOperator:: computeMagneticExcitationField(..)",
                                          "H is an incompatible real field");

       
        
        //H=Ad.(sigma.M)
        if (sigma.getSize()>1) {
          
            H.initField(sigma,M);
            EMM_DemagnetizedOperator::computeMagneticExcitationField(nCells,M.getDimension(),
                                                                     &(*gH)[0],&(*gH)[0]);
        } else {
            EMM_DemagnetizedOperator::computeMagneticExcitationField(nCells,M.getDimension(),
                                                                     &(*gM)[0],&(*gH)[0]);
        }
        
    }
    

        
};

#endif