MilneDeviceModule.h

Go to the documentation of this file.

 
#pragma once
 
#include "modules/ExplicitModule.h"
 
class MilneDeviceModule : public ExplicitModule
{
public:
  JEM_DECLARE_CLASS(MilneDeviceModule, ExplicitModule);
 
  static const char *TYPE_NAME; 
 
  explicit MilneDeviceModule(const String &name = "MilneDevice");
 
  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 ~MilneDeviceModule();
 
protected:
  Vector updForce(const Vector &fint,
                  const Vector &fext,
                  const Properties &globdat);
 
private:
 
  inline void AMupdate_(const Vector &delta_y, const Vector &f_pre,
                        const Vector &f_cur) const;
 
  inline void AMupdate_(const Vector &delta_y, const Vector &f_pre) const;
};
 
//=======================================================================
//   Inline implementations
//=======================================================================
 
inline void MilneDeviceModule::AMupdate_(const Vector &delta_y,
                                         const Vector &f_pre,
                                         const Vector &f_cur) const
{
  delta_y = dtime_ / 2. * (1. * f_pre + 1. * f_cur);
}
 
inline void MilneDeviceModule::AMupdate_(const Vector &delta_y,
                                         const Vector &f_pre) const
{
  delta_y = dtime_ * f_pre;
}