EMMG_SLSDXPeriodicMultiScale.h

Go to the documentation of this file.

#ifndef EMMG_SLSDXPeriodicMultiScale_H
#define EMMG_SLSDXPeriodicMultiScale_H


#include "EMMG_SLPeriodicMultiScale.h"

DEFINE_SPTR(EMMG_SLSDXPeriodicMultiScale);
class EMMG_SLSDXPeriodicMultiScale : public virtual EMMG_SLPeriodicMultiScale   {
  
    SP_OBJECT(EMMG_SLSDXPeriodicMultiScale);
    // ATTRIBUTES
  
private:

    EMMG_RealField mHz;
    EMMG_RealField mYz;
    SV::EMMG_RealField mLevelMagnetisms;
    tBoolean mIsHighPrecision;
    
 
    
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    EMMG_SLSDXPeriodicMultiScale(void);
  
    // DESTRUCTORS 
    
    virtual ~EMMG_SLSDXPeriodicMultiScale(void);   
  
  
public:
    inline static SP::EMMG_SLSDXPeriodicMultiScale New() {
        SP::EMMG_SLSDXPeriodicMultiScale p(new EMMG_SLSDXPeriodicMultiScale(),
                                     EMMG_SLSDXPeriodicMultiScale::Delete());
        p->setThis(p);
        return p;
    };
  

    // SET methods 
public:
  
    // GET

    inline void setIsHighPrecision(const tBoolean& h) {
        mIsHighPrecision=h;
    }

   


public :
  
    void computeMultiGridExcitationField(const tUIndex& nCells,
                                         const tDimension& dim,
                                         const tReal* M,
                                         tReal * H); 

protected:
    virtual void computeZonalDemagnetizedFieldAndNextLevelMagnetizationField(const tUInteger& twoPowerL,
                                                                             const tUIndex& nCells,
                                                                             const tDimension& dim,
                                                                             const tReal* Ml,
                                                                             tReal *Mz,
                                                                             tReal *Mlp1,
                                                                             tReal *H0) const {
    }
    virtual void computeZonalCenteredDemagnetizedFieldFromLevel(const tUInteger& twoPowerL,
                                                                  const tUIndex& nCells,
                                                                  const tDimension& dim,
                                                                  tReal* Ml,
                                                                  tReal *Mz,
                                                                  tReal *H0) const {
      }
  
    
    
private:


    void computeLevelDemagnetizedField(const tUIndex& nCells,const tDimension& dim,
                                       const int& l,
                                       const tReal *Hl,tReal *Yl) const;
    
    void computeHighLevelDemagnetizedField(const tUInteger &nx,const tUInteger& ny,const tUInteger& nz,
                                           const tUIndex& nCells,const tDimension& dim,
                                           const int& l,
                                           tReal *Hz,tReal *Yz,
                                           const tReal* Hl,tReal* Yl) const;
    
    void addLevelPeriodicDemagnetizedField(const tUInteger &nx,const tUInteger& ny,const tUInteger& nz,
                                           const tUIndex& nCells,const tDimension& dim,
                                           const tBoolean isMeshPeriodic[3],
                                           const int& l,
                                           tReal* Hz,tReal* Yz,
                                           const tReal* Hl_1,tReal* Yl) const;
        
    
    void addLevelContributionToSubLevel(const tUInteger& nx,const tUInteger& ny,const tUInteger& nz,
                                        const tUIndex& nCells,const tDimension& dim,
                                        const int& l,const tReal *Yl,tReal* Yl_1) const;
        
    void addLevelContributionToLevel0(const tUInteger& nx,const tUInteger& ny,const tUInteger& nz,
                                      const tUIndex& nCells,const tDimension& dim,
                                      const int& l,const tReal* Yl,tReal* Ydem) const;
    
    
};

#endif