EMM_StaticMagnetostrictionOperator.h

Go to the documentation of this file.

#ifndef EMM_StaticMagnetostrictionOperator_H
#define EMM_StaticMagnetostrictionOperator_H


#include "EMM_MagneticExcitationLinearOperator.h"

#include "EMM_2PackedSymmetricTensors.h"

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

public:
   

    // ASSOCIATION
 
private:
    
   

    //constraint
    SP::EMM_2PackedSymmetricTensors mMagneticStressTensors;

    EMM_2PackedSymmetricTensors  mStressTensor;
    
protected: 

    // METHODS
  
    // CONSTRUCTORS 
  
    EMM_StaticMagnetostrictionOperator(void);
  
    // DESTRUCTORS 
    
    virtual ~EMM_StaticMagnetostrictionOperator(void);
  
    virtual tULLInt getMemorySize() const {
        return EMM_MagneticExcitationLinearOperator::getMemorySize()+mMagneticStressTensors->getMemorySize();
        
    }
public:
    /* \brief set the values of the stress 2-tensor of size 3x3
     * @param n: the size of the values either 1, (value on {0,0}), 3 (diagonal value), 6 (packed form), 9 (full form)
     * @param v: values of the tensor \f$ v = \left[ {\begin{array}{ccc} v[0] & v[1] & v[2] \\ v[1] & v[3]  & v[4] \\ v[2] & v[4] & v[5] \end{array} } \right]  \f$
     * 
     */
    inline void setStressTensor(const tDimension& n,const tReal* v) {
        mStressTensor.init(0,0.0);
        if ((v==null) || (n==0)) return;

        //n=1,3,6
        mStressTensor(0,0,0)=*v;

        //n=3
        if (n==3) {
            mStressTensor(0,1,1)=v[1];
            mStressTensor(0,2,2)=v[2];
        }
        //n=6
        if (n==6) {
            const tReal *Vi=v;
            mStressTensor(0,0,0)=*Vi;Vi++;
            mStressTensor(0,0,1)=*Vi;Vi++;
            mStressTensor(0,0,2)=*Vi;Vi++;
            mStressTensor(0,1,1)=*Vi;Vi++;
            mStressTensor(0,1,2)=*Vi;Vi++;
            mStressTensor(0,2,2)=*Vi;
           
        }
        //n=9
        if (n==9) {
            const tReal *Vi=v;
            mStressTensor(0,0,0)=*Vi;Vi++;
            mStressTensor(0,0,1)=*Vi;Vi++;
            mStressTensor(0,0,2)=*Vi;Vi+=2;
            mStressTensor(0,1,1)=*Vi;Vi++;
            mStressTensor(0,1,2)=*Vi;Vi+=3;
            mStressTensor(0,2,2)=*Vi;
            
        }
        
        
    }
    virtual tBoolean discretize(const EMM_LandauLifschitzSystem& system);

     //Magnetic field methods
    //======================

    virtual void computeMagneticExcitationField(const tUIndex & nCells,
                                                const tDimension& dim,
                                                const EMM_RealArray& sigma,
                                                const tReal  *M,
                                                tReal *H) const;
       
};
#endif