EMM_DemagnetizedOperator.h

Go to the documentation of this file.

#ifndef EMM_DemagnetizedOperator_H
#define EMM_DemagnetizedOperator_H


#include "EMM_MagneticExcitationLinearOperator.h"

#include "EMM_MultiScaleGrid.h"

#include "MATH_MultiLevelsToeplitzMatrix.h"

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

    

  
    
private:
    tString mMultiScaleGridName;
    
    //levels number
    tUInteger mLevelsNumber;
    tInteger mZonalLevelsNumber;
    tBoolean mIsLevelsNumberAutoComputed;
    
    // ASSOCIATION
    SP::MATH_MultiLevelsToeplitzMatrix mToeplitzMatrix;
    SP::EMM_MultiScaleGrid mMultiScaleGrid;

   
protected: 

    // METHODS
  
    // CONSTRUCTORS 
  
    EMM_DemagnetizedOperator(void);
  
    // DESTRUCTORS 
    
    virtual ~EMM_DemagnetizedOperator(void);   
  
  
public:
   
    
    
    // SET methods
    // ===========

    inline void setMultiScaleGridName(const tString& name) {
        mMultiScaleGridName=name;
    }
    void setStorage(const tString& storage);

    inline void setLevelsNumber(const tUInteger& l) {
        mLevelsNumber=l;
    }
    inline void setZonalLevelsNumber(const tInteger& l) {
        mZonalLevelsNumber=l;
    }

    inline void setIsLevelsNumberAutoComputed(const tBoolean& v) {
        mIsLevelsNumberAutoComputed=v;
    }
    
    // GET methods
    // ===========
    inline EMM_MultiScaleGrid& getMultiScaleGrid() {
        if (mMultiScaleGrid.get()==null) {
            throw EMM_Exception("emicrom/operators/demagnetized",
                                "EMM_DemagnetizedOperator::getMultiScaleGrid()",
                                "Fatal Error: no multi grid set");
        }
        return *mMultiScaleGrid.get();
    }
    
    
    inline const tString& getMultiScaleGridName() const {
        return mMultiScaleGridName;
    }
    inline const tBoolean& isLevelsNumberAutoComputed() const {
        return mIsLevelsNumberAutoComputed;
    }
    inline const tUInteger& getLevelsNumber() const {
        return mLevelsNumber;
    }
    
    virtual tULLInt getMemorySize() const {
        return EMM_MagneticExcitationLinearOperator::getMemorySize()+mToeplitzMatrix->getMemorySize();
    }

   

    //DATA Methods
    //=============

    virtual tBoolean discretize(const EMM_LandauLifschitzSystem& system);


    inline const MATH_MultiLevelsToeplitzMatrix& getToeplitzMatrix() const {
        return *mToeplitzMatrix.get();
    }
 


    //MAGNETIC Excitation methods
    //============================
    
    virtual void computeMagneticExcitationField(const tUIndex& nCells,
                                                const tDimension& dim,
                                                const tReal *sigmaM,
                                                tReal *H) const;

public:
    virtual tString toString() const;
        
};

#endif