EMM_MinimalExchangeOperator.h

Go to the documentation of this file.

#ifndef EMM_MinimalExchangeOperator_H
#define EMM_MinimalExchangeOperator_H

#include "EMM_ExchangeOperator.h"

#include "EMM_Grid3D.h"
#include "EMM_LandauLifschitzSystem.h"
#include "EMM_MatterField.h"

DEFINE_SPTR(EMM_MinimalExchangeOperator);
class EMM_MinimalExchangeOperator : public virtual EMM_ExchangeOperator {
  
    SP_OBJECT(EMM_MinimalExchangeOperator);
    // ATTRIBUTES

  
private:
  
   
    
public:
   

    // ASSOCIATION
 
private:

    CORE_UIndexMorseArray mNeighbors;
    
    EMM_RealArray mA;
    
    tReal *mL2;
    
    //leading dimension of the work array
    tUCInt mLDW;
    
    //work array
    tReal *mWork;

    
protected: 

    // METHODS
  
    // CONSTRUCTORS 
  
    EMM_MinimalExchangeOperator(void);
  
    // DESTRUCTORS 
    
    virtual ~EMM_MinimalExchangeOperator(void);   
  
  
public:
    
    
    //DATA Methods
    //==============
    
    virtual tBoolean discretize(const EMM_LandauLifschitzSystem& system);


public:
    virtual tULLInt getMemorySize() const {
        return EMM_ExchangeOperator::getMemorySize()+mNeighbors.getMemorySize()+mA.getMemorySize()+mLDW*sizeof(tReal);
            
    }

protected:
    
    inline const CORE_UIndexMorseArray& getNeighbors() const {
        return mNeighbors;
    }
    inline CORE_UIndexMorseArray& getNeighbors()  {
        return mNeighbors;
    }

    inline const EMM_RealArray& getExchangeField() const {
        return mA;
    }
    inline EMM_RealArray& getExchangeField()  {
        return mA;
    }

  
   
    inline const tReal* getFaceAdimensionizedLength() const {
        return mL2;
    }
   
protected:
    inline void clearWorkingArray() {
        if (mWork!=null) delete[] mWork;
        mWork=null;
    }



    //MAGNETIC Excitation field methods
    //=================================
    
protected:
    
    
    virtual void computeMagneticExcitationField(const tUIndex& nCells,
                                                const tDimension& dim,
                                                const EMM_RealArray& sigma,
                                                const tReal* M,
                                                tReal* H) const;
        

};

#endif