EMMG_UniaxialAnisotropyOperator.h

Go to the documentation of this file.

#ifndef EMMG_UniaxialAnisotropyOperator_H
#define EMMG_UniaxialAnisotropyOperator_H


#include "EMM_UniaxialAnisotropyOperator.h"


DEFINE_SPTR(EMMG_UniaxialAnisotropyOperator);
class EMMG_UniaxialAnisotropyOperator : public virtual EMM_UniaxialAnisotropyOperator {
    
    SP_OBJECT(EMMG_UniaxialAnisotropyOperator);
    // ATTRIBUTES
    
    
private:

    
public:
   

    // ASSOCIATION
 
private:
    
  
    
protected: 

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

protected:
    
   

    
public:
  
    void computeMagneticExcitationField(const EMM_RealArray& sigma,
                                        const EMM_RealField& M,
                                        EMM_RealField& H) const {
        
        // size of the operator
        tUIndex nCells=M.getSize();
        const tDimension &dim=M.getDimension();
        
        
        //set the H size
        H.setDimension(dim);
        H.setSize(nCells);
        
        //no computation
        if (nCells==0) return;
        
        //verify the type of M
        const EMMG_RealField *gM=dynamic_cast<const EMMG_RealField *>(&M);
        if (gM==null) throw EMM_Exception("generic","EMMG_PlanarAnisotropyOperator:: 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","EMM_PlanarAnisotropyOperator:: computeMagneticExcitationField(..)",
                                          "H is an incompatible real field");
        
        
        
        EMM_UniaxialAnisotropyOperator::computeMagneticExcitationField(nCells,dim,sigma,&(*gM)[0],&(*gH)[0]);
        
        
    }
private:
  
  
public:
   
};

#endif