EMM_OptimalTimeStep.h

Go to the documentation of this file.

#ifndef EMM_OptimalTimeStep_H
#define EMM_OptimalTimeStep_H

#include "EMM_RegularTimeStep.h"


#include "EMM_RealField.h"
#include "MATH_P4.h"

#include "EMM_LandauLifschitzSystem.h"
class EMM_GradGaussLegendreRelaxation;

DEFINE_SPTR(EMM_OptimalTimeStep);
class EMM_OptimalTimeStep : public EMM_RegularTimeStep {
  
    SP_OBJECT(EMM_OptimalTimeStep);
    // ATTRIBUTES
    
public:
    
    
   
private:

    mutable MATH_P4 mP4;
  
    //working array
    tReal mT[4];
    SP::EMM_RealField mGradH_gradF;
    
 
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    EMM_OptimalTimeStep(void);
  
  
  
    // DESTRUCTORS 
  
    
    virtual ~EMM_OptimalTimeStep(void);   
  
  
private:

public:
    // New share pointer constructor
    static inline SP::EMM_OptimalTimeStep New() {
        SP::EMM_OptimalTimeStep p(new EMM_OptimalTimeStep(),EMM_OptimalTimeStep::Delete());
        p->setThis(p);
        return p;
    }
    
    // SET methods 
  
    
    // GET methods 
   
public:
    virtual tBoolean isSystemCompatible(SPC::EMM_LandauLifschitzSystem system) const {
        return EMM_RegularTimeStep::isSystemCompatible(system) && system->isAffine();
    }
    
    virtual tBoolean computeNewTimeStep(tReal& dt,tUInteger& nEvaluations);
    
public:
    // OTHERS methods

    

    virtual tString toString() const {
        return EMM_RegularTimeStep::toString();
    }

   

private:
    tBoolean computeOptimalTimeStep(EMM_GradGaussLegendreRelaxation& GSystem,tReal& dt);
};

#endif