I94-checkpointing

src/error.rs

@@ -74,6 +74,8 @@ pub enum OdeSolverError {
     #[error("atol must have length 1 or equal to the number of states")]
     AtolLengthMismatch,
     #[error("t_eval must be increasing and all values must be greater than or equal to the current time")]
+    StateProblemMismatch,
+    #[error("State is not consistent with the problem equations")]
     InvalidTEval,
     #[error("Sundials error: {0}")]
     SundialsError(String),

src/lib.rs

@@ -153,10 +153,11 @@ use nonlinear_solver::{
 };
 use ode_solver::jacobian_update::JacobianUpdate;
 pub use ode_solver::{
-    bdf::Bdf, builder::OdeBuilder, equations::OdeEquations, equations::OdeSolverEquations,
-    method::OdeSolverMethod, method::OdeSolverState, method::OdeSolverStopReason,
-    problem::OdeSolverProblem, sdirk::Sdirk, sens_equations::SensEquations,
-    sens_equations::SensInit, sens_equations::SensRhs, tableau::Tableau,
+    bdf::Bdf, bdf_state::BdfState, builder::OdeBuilder, equations::OdeEquations,
+    equations::OdeSolverEquations, method::OdeSolverMethod, method::OdeSolverStopReason,
+    problem::OdeSolverProblem, sdirk::Sdirk, sdirk_state::SdirkState,
+    sens_equations::SensEquations, sens_equations::SensInit, sens_equations::SensRhs,
+    state::OdeSolverState, tableau::Tableau,
 };
 pub use op::{
     closure::Closure, constant_closure::ConstantClosure, linear_closure::LinearClosure,