EMM_Operator.h

Go to the documentation of this file.

#ifndef EMM_Operator_H
#define EMM_Operator_H

#include "EMM_Object.h"


#include "EMM_RealField.h"
#include "EMM_Grid3D.h"

class EMM_LandauLifschitzSystem;

DEFINE_SPTR(EMM_Operator);
DEFINE_SVPTR(EMM_Operator);
class EMM_Operator : public EMM_Object {
  
    SP_OBJECT(EMM_Operator);
    // ATTRIBUTES
    
public:

   
private:
    tReal mElementVolume;
    tBoolean mIsCubic;
    
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    EMM_Operator(void);
  
  
  
    // DESTRUCTORS 
  
    
    virtual ~EMM_Operator(void);   
  
  
    // SET methods 
 
protected:
   
    
    // GET methods
    
public:
    virtual tULLInt getMemorySize() const {
        return 0;
    }


    
    virtual tString getName() const {
        tString cn=getClassName();
        //find the first '_' from first occurence of 'EMM'
        tUIndex i0=cn.find("EMM");
        i0=cn.find("_",i0);
        //find the Operator name
        tUIndex i1=cn.rfind("Operator");
        //return the string betwen '_' & 'Operator'
        return cn.substr(i0+1,i1-i0-1);
    }
    virtual tBoolean isAffine() const=0;
    
    virtual tBoolean isGradientComputationable() const=0;

    inline const tBoolean& isCubicVolume() const {
        return mIsCubic;
    }
  
    inline const tReal& getElementVolume() const {
        return mElementVolume;
    }

    
public:
    
    //DATA methods
    //=============
    virtual void adimensionize(const tReal& Le,const tReal& Ms,const tReal& T,const tReal& L) {
    }
    
    virtual tBoolean discretize(const EMM_LandauLifschitzSystem& system);
    
    virtual tBoolean resetToInitialState(const EMM_LandauLifschitzSystem& system)=0;
    
    
    
    // FIELDS Methods
    //================
   
    virtual  tUSInt getDataFieldsNumber() const {
        return 0;
    }
    
    virtual tBoolean getDataFieldSpace(const tUSInt& index,tString& dataName,tFlag& supportType,tString& dFieldType,tUIndex& n,tDimension& dim) const {
        dataName="";
        supportType=EMM_Grid3D::ELEMENT;
        n=0;
        dim=1;
        dFieldType="";
        return false;
    }
    virtual tBoolean getDataField(const tUSInt& index,tString& dataName,tUIndex& n,tDimension& dim,const float *& values) const {
        dataName="";
        n=0;
        dim=1;
        return false;
    }
    virtual tBoolean getDataField(const tUSInt& index,tString& dataName,tUIndex& n,tDimension& dim,const double *& values) const {
        dataName="";
        n=0;
        dim=1;
        return false;
    }
    virtual tBoolean getDataField(const tUSInt& index,tString& dataName,tUIndex& n,tDimension& dim,const long double *& values) const {
        dataName="";
        n=0;
        dim=1;
        return false;
    }

 
   

public:

    virtual tBoolean backup(const tString& prefix,const tString& suffix,const tString& ext) const {
        return true;
    }
    virtual tBoolean restore(const EMM_LandauLifschitzSystem& system,
                             const tString& prefix,const tString& suffix,const tString& ext) {
        return true;
    }
    
    virtual tBoolean computeFieldsAtTime(const tReal& t,
                                         const tFlag& order,
                                         const EMM_RealArray& sigma,
                                         const EMM_RealField& dM_dt0,
                                         const EMM_RealField& M0)=0;
    
  
    virtual tBoolean updateAtNextTimeStep(const tReal& dt,
                                          const EMM_RealArray& sigma,
                                          const EMM_RealField& Mt)=0;
    

    
};

#endif