ode/include/nk/ode/solvers.hpp

Go to the documentation of this file.

// toolboxes/ode/include/nk/ode/solvers.hpp
//
// Initial-value problem (IVP) solvers for `y' = f(t, y)`.
//
// References:
//   • Hairer, Nørsett, Wanner, "Solving Ordinary Differential Equations I"
//     (Springer, 2nd ed. 1993) — DOPRI5 coefficients and step-size strategy.
//   • Dormand & Prince, "A family of embedded Runge-Kutta formulae",
//     J. Comput. Appl. Math. 6(1), 1980.
//   • Shampine & Reichelt, "The MATLAB ODE Suite", SIAM J. Sci. Comput.
//     18(1), 1997.
 
#pragma once
 
#include <memory_resource>
#include <nk/value/fn_handle.hpp>
#include <nk/value/value.hpp>
 
#include <tuple>
 
namespace nk::ode {

 

std::tuple<Value, Value>
ode45(FnHandle rhs, const Value &tspan, const Value &y0,
      const Value &opts, std::pmr::memory_resource *mr = nullptr);

std::tuple<Value, Value>
ode23(FnHandle rhs, const Value &tspan, const Value &y0,
      const Value &opts, std::pmr::memory_resource *mr = nullptr);
 

} // namespace nk::ode