ODEintSolver.h

/* Copyright 2017-2021 Institute for Automation of Complex Power Systems,
 *                     EONERC, RWTH Aachen University
 *
 * This Source Code Form is subject to the terms of the Mozilla Public
 * License, v. 2.0. If a copy of the MPL was not distributed with this
 * file, You can obtain one at https://mozilla.org/MPL/2.0/.
 *********************************************************************************/
 
#pragma once
 
#include <iostream>
#include <list>
#include <vector>
 
#include <dpsim-models/Solver/ODEintInterface.h>
#include <dpsim-models/SystemTopology.h>
#include <dpsim/Solver.h>
 
#include <boost/numeric/odeint/integrate/integrate_const.hpp> //ODEInt Integrator with constant time step
#include <boost/numeric/odeint/stepper/runge_kutta4.hpp> //ODEInt Runge-Kutta stepper
 
namespace DPsim {
class ODEintSolver : public Solver {
 
protected:
  CPS::ODEintInterface::Ptr mComponent;
  Real mTimestep;
  int ProbDim;
  boost::numeric::odeint::runge_kutta4<std::vector<Real>> stepper;
  std::vector<std::vector<Real>> solution;
  std::vector<Real> times;
  std::vector<CPS::ODEintInterface::stateFnc> system;
 
private:
  inline static ODEintSolver *self = nullptr;
  static void StateSpaceWrapper(const std::vector<double> &y,
                                std::vector<double> ydot, double t);
 
public:
  std::vector<Real> curSolution;
  ODEintSolver(String name, CPS::ODEintInterface::Ptr comp, Real dt, Real t0);

  Real step(Real time);

  ~ODEintSolver();
};
 
} // namespace DPsim