diff --git a/docs/html/md_pages_rom_galerkin_default_explicit.html b/docs/html/md_pages_rom_galerkin_default_explicit.html
index 2849402..230bf2d 100644
--- a/docs/html/md_pages_rom_galerkin_default_explicit.html
+++ b/docs/html/md_pages_rom_galerkin_default_explicit.html
@@ -152,8 +152,7 @@
const scalar_t dt = 1.;
const int num_steps = 3;
Observer obs;
- pressio::ode::advance_n_steps_and_observe(
- problem.stepper(), romState, 0., dt, num_steps, obs);
+ pressio::ode::advance_n_steps_and_observe(problem, romState, 0., dt, num_steps, obs);
pressio::log::finalize();
}
diff --git a/docs/html/md_pages_rom_galerkin_default_explicit_custom_types.html b/docs/html/md_pages_rom_galerkin_default_explicit_custom_types.html
index a740211..a885387 100644
--- a/docs/html/md_pages_rom_galerkin_default_explicit_custom_types.html
+++ b/docs/html/md_pages_rom_galerkin_default_explicit_custom_types.html
@@ -326,10 +326,10 @@
const scalar_t dt = 1.;
const int num_steps = 3;
Observer obs;
- pressio::ode::advance_n_steps_and_observe(
- problem.stepper(), romState, 0., dt, num_steps, obs);
+ pressio::ode::advance_n_steps_and_observe(problem, romState, 0., dt, num_steps, obs);
- pressio::log::finalize();
+ pressio::log::finalize();
+}
diff --git a/docs/html/md_pages_rom_galerkin_default_implicit.html b/docs/html/md_pages_rom_galerkin_default_implicit.html
index e54b7d8..0335bcd 100644
--- a/docs/html/md_pages_rom_galerkin_default_implicit.html
+++ b/docs/html/md_pages_rom_galerkin_default_implicit.html
@@ -157,16 +157,15 @@
auto problem = pressio::rom::galerkin::create_default_implicit_problem
(odescheme, fomSystem, decoder, romState, fomReferenceState);
using problem_t = decltype(problem);
- auto & stepperObj = problem.stepper();
using galerkin_jacobian_t = typename pressio::Traits<problem_t>::galerkin_jacobian_type;
using lin_solver_t = pls::Solver<pls::iterative::LSCG, galerkin_jacobian_t>;
lin_solver_t linearSolverObj;
- auto nonLinSolver = pnonls::create_newton_raphson(stepperObj, romState, linearSolverObj);
+ auto nonLinSolver = pnonls::create_newton_raphson(problem, romState, linearSolverObj);
nonLinSolver.setMaxIterations(1);
scalar_t dt = 2.;
- pressio::ode::advance_n_steps(stepperObj, romState, 0.0, dt, 1, nonLinSolver);
+ pressio::ode::advance_n_steps(problem, romState, 0.0, dt, 1, nonLinSolver);
std::cout << romState << std::endl;
pressio::log::finalize();
diff --git a/docs/latex/md_pages_rom_galerkin_default_explicit.tex b/docs/latex/md_pages_rom_galerkin_default_explicit.tex
index 5ccece2..96299c1 100644
--- a/docs/latex/md_pages_rom_galerkin_default_explicit.tex
+++ b/docs/latex/md_pages_rom_galerkin_default_explicit.tex
@@ -75,8 +75,7 @@
\DoxyCodeLine{ \textcolor{keyword}{const} scalar\_t dt = 1.;}
\DoxyCodeLine{ \textcolor{keyword}{const} \textcolor{keywordtype}{int} num\_steps = 3;}
\DoxyCodeLine{ Observer obs;}
-\DoxyCodeLine{ pressio::ode::advance\_n\_steps\_and\_observe(}
-\DoxyCodeLine{ problem.stepper(), romState, 0., dt, num\_steps, obs);}
+\DoxyCodeLine{ pressio::ode::advance\_n\_steps\_and\_observe(problem, romState, 0., dt, num\_steps, obs);}
\DoxyCodeLine{}
\DoxyCodeLine{ pressio::log::finalize();}
\DoxyCodeLine{\}}
diff --git a/docs/latex/md_pages_rom_galerkin_default_explicit_custom_types.tex b/docs/latex/md_pages_rom_galerkin_default_explicit_custom_types.tex
index 26d50f3..64d7cc7 100644
--- a/docs/latex/md_pages_rom_galerkin_default_explicit_custom_types.tex
+++ b/docs/latex/md_pages_rom_galerkin_default_explicit_custom_types.tex
@@ -264,10 +264,10 @@
\DoxyCodeLine{ \textcolor{keyword}{const} scalar\_t dt = 1.;}
\DoxyCodeLine{ \textcolor{keyword}{const} \textcolor{keywordtype}{int} num\_steps = 3;}
\DoxyCodeLine{ Observer obs;}
-\DoxyCodeLine{ pressio::ode::advance\_n\_steps\_and\_observe(}
-\DoxyCodeLine{ problem.stepper(), romState, 0., dt, num\_steps, obs);}
+\DoxyCodeLine{ pressio::ode::advance\_n\_steps\_and\_observe(problem, romState, 0., dt, num\_steps, obs);}
\DoxyCodeLine{}
\DoxyCodeLine{ pressio::log::finalize();}
+\DoxyCodeLine{\}}
\end{DoxyCode}
\ No newline at end of file
diff --git a/docs/latex/md_pages_rom_galerkin_default_implicit.tex b/docs/latex/md_pages_rom_galerkin_default_implicit.tex
index af6addc..669f481 100644
--- a/docs/latex/md_pages_rom_galerkin_default_implicit.tex
+++ b/docs/latex/md_pages_rom_galerkin_default_implicit.tex
@@ -80,16 +80,15 @@
\DoxyCodeLine{ \textcolor{keyword}{auto} problem = pressio::rom::galerkin::create\_default\_implicit\_problem}
\DoxyCodeLine{ (odescheme, fomSystem, decoder, romState, fomReferenceState);}
\DoxyCodeLine{ \textcolor{keyword}{using} problem\_t = decltype(problem);}
-\DoxyCodeLine{ \textcolor{keyword}{auto} \& stepperObj = problem.stepper();}
\DoxyCodeLine{}
\DoxyCodeLine{ \textcolor{keyword}{using} galerkin\_jacobian\_t = \textcolor{keyword}{typename} pressio::Traits::galerkin\_jacobian\_type;}
\DoxyCodeLine{ \textcolor{keyword}{using} lin\_solver\_t = pls::Solver;}
\DoxyCodeLine{ lin\_solver\_t linearSolverObj;}
-\DoxyCodeLine{ \textcolor{keyword}{auto} nonLinSolver = pnonls::create\_newton\_raphson(stepperObj, romState, linearSolverObj);}
+\DoxyCodeLine{ \textcolor{keyword}{auto} nonLinSolver = pnonls::create\_newton\_raphson(problem, romState, linearSolverObj);}
\DoxyCodeLine{ nonLinSolver.setMaxIterations(1);}
\DoxyCodeLine{}
\DoxyCodeLine{ scalar\_t dt = 2.;}
-\DoxyCodeLine{ pressio::ode::advance\_n\_steps(stepperObj, romState, 0.0, dt, 1, nonLinSolver);}
+\DoxyCodeLine{ pressio::ode::advance\_n\_steps(problem, romState, 0.0, dt, 1, nonLinSolver);}
\DoxyCodeLine{ std::cout << romState << std::endl;}
\DoxyCodeLine{}
\DoxyCodeLine{ pressio::log::finalize();}
diff --git a/docs/m.math.cache b/docs/m.math.cache
index ffd62a1..3dc9df3 100644
Binary files a/docs/m.math.cache and b/docs/m.math.cache differ
diff --git a/docs/xml/md_pages_rom_galerkin_default_explicit.xml b/docs/xml/md_pages_rom_galerkin_default_explicit.xml
index 18fac89..0f71a26 100644
--- a/docs/xml/md_pages_rom_galerkin_default_explicit.xml
+++ b/docs/xml/md_pages_rom_galerkin_default_explicit.xml
@@ -79,8 +79,7 @@
constscalar_tdt=1.;
constintnum_steps=3;
Observerobs;
-pressio::ode::advance_n_steps_and_observe(
-problem.stepper(),romState,0.,dt,num_steps,obs);
+pressio::ode::advance_n_steps_and_observe(problem,romState,0.,dt,num_steps,obs);
pressio::log::finalize();
}
diff --git a/docs/xml/md_pages_rom_galerkin_default_explicit_custom_types.xml b/docs/xml/md_pages_rom_galerkin_default_explicit_custom_types.xml
index aa29a45..c78b78c 100644
--- a/docs/xml/md_pages_rom_galerkin_default_explicit_custom_types.xml
+++ b/docs/xml/md_pages_rom_galerkin_default_explicit_custom_types.xml
@@ -267,10 +267,10 @@
constscalar_tdt=1.;
constintnum_steps=3;
Observerobs;
-pressio::ode::advance_n_steps_and_observe(
-problem.stepper(),romState,0.,dt,num_steps,obs);
+pressio::ode::advance_n_steps_and_observe(problem,romState,0.,dt,num_steps,obs);
pressio::log::finalize();
+}
diff --git a/docs/xml/md_pages_rom_galerkin_default_implicit.xml b/docs/xml/md_pages_rom_galerkin_default_implicit.xml
index 23e9e8e..87a3ec3 100644
--- a/docs/xml/md_pages_rom_galerkin_default_implicit.xml
+++ b/docs/xml/md_pages_rom_galerkin_default_implicit.xml
@@ -84,16 +84,15 @@
autoproblem=pressio::rom::galerkin::create_default_implicit_problem
(odescheme,fomSystem,decoder,romState,fomReferenceState);
usingproblem_t=decltype(problem);
-auto&stepperObj=problem.stepper();
usinggalerkin_jacobian_t=typenamepressio::Traits<problem_t>::galerkin_jacobian_type;
usinglin_solver_t=pls::Solver<pls::iterative::LSCG,galerkin_jacobian_t>;
lin_solver_tlinearSolverObj;
-autononLinSolver=pnonls::create_newton_raphson(stepperObj,romState,linearSolverObj);
+autononLinSolver=pnonls::create_newton_raphson(problem,romState,linearSolverObj);
nonLinSolver.setMaxIterations(1);
scalar_tdt=2.;
-pressio::ode::advance_n_steps(stepperObj,romState,0.0,dt,1,nonLinSolver);
+pressio::ode::advance_n_steps(problem,romState,0.0,dt,1,nonLinSolver);
std::cout<<romState<<std::endl;
pressio::log::finalize();
diff --git a/tpls/pressio b/tpls/pressio
index 049b7be..5bee44a 160000
--- a/tpls/pressio
+++ b/tpls/pressio
@@ -1 +1 @@
-Subproject commit 049b7bec173d635f8d35ad40ae56e1d7559263e0
+Subproject commit 5bee44a1cb9d2fc117d89da4a401bff8d0918605
diff --git a/tutorials/rom_default_explicit_galerkin_custom_types.cc b/tutorials/rom_default_explicit_galerkin_custom_types.cc
index 83ba53c..4c69e42 100644
--- a/tutorials/rom_default_explicit_galerkin_custom_types.cc
+++ b/tutorials/rom_default_explicit_galerkin_custom_types.cc
@@ -251,8 +251,7 @@ int main(int argc, char *argv[])
const scalar_t dt = 1.;
const int num_steps = 3;
Observer obs;
- pressio::ode::advance_n_steps_and_observe(
- problem.stepper(), romState, 0., dt, num_steps, obs);
+ pressio::ode::advance_n_steps_and_observe(problem, romState, 0., dt, num_steps, obs);
pressio::log::finalize();
}
diff --git a/tutorials/rom_default_explicit_galerkin_eigen.cc b/tutorials/rom_default_explicit_galerkin_eigen.cc
index 2134c4b..98380c0 100644
--- a/tutorials/rom_default_explicit_galerkin_eigen.cc
+++ b/tutorials/rom_default_explicit_galerkin_eigen.cc
@@ -116,8 +116,7 @@ int main(int argc, char *argv[])
const scalar_t dt = 1.;
const int num_steps = 3;
Observer obs;
- pressio::ode::advance_n_steps_and_observe(
- problem.stepper(), romState, 0., dt, num_steps, obs);
+ pressio::ode::advance_n_steps_and_observe(problem, romState, 0., dt, num_steps, obs);
pressio::log::finalize();
}
diff --git a/tutorials/rom_default_implicit_galerkin_eigen.cc b/tutorials/rom_default_implicit_galerkin_eigen.cc
index 7f79964..eadfec2 100644
--- a/tutorials/rom_default_implicit_galerkin_eigen.cc
+++ b/tutorials/rom_default_implicit_galerkin_eigen.cc
@@ -121,16 +121,15 @@ int main(int argc, char *argv[])
auto problem = pressio::rom::galerkin::create_default_implicit_problem
(odescheme, fomSystem, decoder, romState, fomReferenceState);
using problem_t = decltype(problem);
- auto & stepperObj = problem.stepper();
using galerkin_jacobian_t = typename pressio::Traits::galerkin_jacobian_type;
using lin_solver_t = pls::Solver;
lin_solver_t linearSolverObj;
- auto nonLinSolver = pnonls::create_newton_raphson(stepperObj, romState, linearSolverObj);
+ auto nonLinSolver = pnonls::create_newton_raphson(problem, romState, linearSolverObj);
nonLinSolver.setMaxIterations(1);
scalar_t dt = 2.;
- pressio::ode::advance_n_steps(stepperObj, romState, 0.0, dt, 1, nonLinSolver);
+ pressio::ode::advance_n_steps(problem, romState, 0.0, dt, 1, nonLinSolver);
std::cout << romState << std::endl;
pressio::log::finalize();