RungeKutta.h

Go to the documentation of this file.

#ifndef KINETX_RUNGEKUTTA_H_
#define KINETX_RUNGEKUTTA_H_

#include "Kinetx/Integrator/Integrator.h"
#include "Kinetx/Network.h"

namespace Scine {
namespace Kinetx {
class RungeKutta : public IntegratorBase {
 public:
  RungeKutta(Network& net);

  void propagate(Eigen::VectorXd& concentrations, Eigen::VectorXd& yFlux, Eigen::VectorXd& rFlux,
                 Eigen::VectorXd& rForwardFlux, Eigen::VectorXd& rBackwardFlux, double& t, double& dt) const;

 protected:
  virtual void propagateY(Eigen::VectorXd& concentrations, double& t, double& dt) const = 0;
};

} /* namespace Kinetx */
} /* namespace Scine */

#endif // KINETX_RUNGEKUTTA_H_