SM_StratonovichSystem.h

Go to the documentation of this file.

#ifndef SM_StratonovichSystem_H
#define SM_StratonovichSystem_H
 
#include "SM_TemplatedSystem.h"
 
class SM_StratonovichSystem : public SM_TemplatedSystem<SM_StratonovichSystem> {
   
    //attributes
private :
 
    tReal mK;//the stratonovich constant for LL function
    
  
    
public:
    // CONSTRUCTORS 
    SM_StratonovichSystem(void) {
        
        
    }
    
    // DESTRUCTORS 
    virtual ~SM_StratonovichSystem(void) {
    }
 
 
public : 
   
    inline static CORE_UniquePointer<SM_StratonovichSystem> New() {
        CORE_UniquePointer<SM_StratonovichSystem> p=CORE_UniquePointer<SM_StratonovichSystem>(new SM_StratonovichSystem(),SM_StratonovichSystem::Delete());
        p->getLandauLifschitzFunction().setIsMuNormalized(true);
        return p;
    }
 
    //MEMORY
    
    virtual tMemSize getMemorySize() const {
        return sizeof(*this)+getContentsMemorySize();
    }
    
 
public:
  
  
 
    inline const tReal& getK() const {
        return mK;
    }
    
public:
    virtual void discretize() override {
        SM_TemplatedSystem<SM_StratonovichSystem>::discretize();
        const SM_LandauLifschitzFunction &ll=this->getLandauLifschitzFunction();
        const tReal& alpha=ll.getAlpha();
        const tReal& beta=ll.getBeta();
        mK=alpha*alpha*(1+beta*beta);
        
    }
    
    inline tBoolean computeTemplatedMuAtNextTimeStep(const tReal& dt,const tReal& epsilon_t,
                                                 const tIndex& nParticles,const tDimension& dim,
                                                 const tReal* mu_t,
                                                 tReal *mu_tpdt)  const {
        tBoolean succeeds =true;
        tIndex i,N=nParticles*dim;
        const tReal *mu_t_i=mu_t;
        tReal *mu_tpdt_i=mu_tpdt;
        for (i=0;i<N;i++) {
            //update mu(t+dt)
            (*mu_tpdt_i)+=(1-mK*epsilon_t*epsilon_t*dt)*(*mu_t_i);
 
            //verify the convergence of the algorithm
            succeeds=succeeds && (fabs(*mu_tpdt_i)<10);
            
            //next element
            mu_t_i++;
            mu_tpdt_i++;
        }
 
        return succeeds;
        
    }
    
};
 
 
#endif