EMMG_SLRPPeriodicMultiScale.h

Go to the documentation of this file.

#ifndef EMMG_SLRPPeriodicMultiScale_H
#define EMMG_SLRPPeriodicMultiScale_H

#include "EMMG_SLPeriodicMultiScale.h"

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


    SV::EMMG_RealField mLevelMagnetisms;
    SV::EMMG_RealField mZoneMagnetisms;
    SV::EMMG_RealField mNearZoneMagnetisms;

 
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    EMMG_SLRPPeriodicMultiScale(void);
  
    // DESTRUCTORS 
    
    virtual ~EMMG_SLRPPeriodicMultiScale(void);   
  
  
public:
    inline static SP::EMMG_SLRPPeriodicMultiScale New() {
        SP::EMMG_SLRPPeriodicMultiScale p(new EMMG_SLRPPeriodicMultiScale(),
                                     EMMG_SLRPPeriodicMultiScale::Delete());
        p->setThis(p);
        return p;
    };
  

    // SET methods 
public:
   
    


public :
  
    virtual 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 computeNearMagnetism(const tUInteger& nx,
                              const tUInteger& ny,
                              const tUInteger& nz,
                              const tDimension& dim,
                              const tBoolean isMeshPeriodic[3],
                              const tReal *M,
                              SV::EMMG_RealField& nearMs) const;
    
    void computeZoneMagnetism(const tUInteger& nx,
                              const tUInteger& ny,
                              const tUInteger& nz,
                              const EMMG_RealField& M,
                              SV::EMMG_RealField& zoneMs) const;

    void sumContributions(const tUInteger& nx,const tUInteger& ny,const tUInteger& nz,const tDimension& dim,
                          const tBoolean isMeshPeriodic[3],
                          const SV::EMMG_RealField& levelHs,
                          const SV::EMMG_RealField& zoneHs,
                          const SV::EMMG_RealField& nearHs,
                          tReal *fineMagnetism) const;

  

};

#endif