EmbeddedRKModule.h

Go to the documentation of this file.

 
#pragma once
 
#include "modules/ExplicitModule.h"
#include <jem/base/Class.h>
 
class EmbeddedRKModule : public ExplicitModule
{
public:
  JEM_DECLARE_CLASS(EmbeddedRKModule, ExplicitModule);
 
  static const char *TYPE_NAME; 
 
  explicit EmbeddedRKModule(const String &name = "EmbeddedRK");
 
  virtual Status init(const Properties &conf,
                      const Properties &props,
                      const Properties &globdat) override;
 
  virtual void solve(const Properties &info,
                     const Properties &globdat) override;
 
  static Ref<Module> makeNew(const String &name,
                             const Properties &conf,
                             const Properties &props,
                             const Properties &globdat);
 
  static void declare();
 
  virtual ~EmbeddedRKModule();
 
private:
  void initODE23_();
 
  void initODE45_();
 
  void correctDisp_(const Vector &uncorrected, const Vector &delta);
 
  double bernoulliCoeff_(const idx_t i);
 
  Matrix adjoint_(const Matrix &point,
                  const Matrix &about,
                  const idx_t iterate);
 
  void invDerivExpMap_(const Matrix &res,
                       const Matrix &point,
                       const Matrix &about);
 
private:
  String kind_;     
  idx_t order_;     
  idx_t butchSize_; 
  Matrix a_;        
  Vector b_;        
  bool fsal_;       
  Vector c_;        
};