Update code to match pybamm

pybamm-team · Jan 6, 2025 · 07c27cc · 07c27cc
1 parent c5b4e73
commit 07c27cc
Show file tree

Hide file tree

Showing 3 changed files with 5 additions and 0 deletions.
diff --git a/src/pybammsolvers/idaklu_source/IDAKLUSolverOpenMP.inl b/src/pybammsolvers/idaklu_source/IDAKLUSolverOpenMP.inl
@@ -192,6 +192,9 @@ void IDAKLUSolverOpenMP<ExprSet>::SetSolverOptions() {
   // Initial step size
   CheckErrors(IDASetInitStep(ida_mem, solver_opts.dt_init));
 
+  // Minimum absolute step size
+  CheckErrors(IDASetMinStep(ida_mem, solver_opts.dt_min));
+
   // Maximum absolute step size
   CheckErrors(IDASetMaxStep(ida_mem, solver_opts.dt_max));
 

diff --git a/src/pybammsolvers/idaklu_source/Options.cpp b/src/pybammsolvers/idaklu_source/Options.cpp
@@ -142,6 +142,7 @@ SolverOptions::SolverOptions(py::dict &py_opts)
       max_order_bdf(py_opts["max_order_bdf"].cast<int>()),
       max_num_steps(py_opts["max_num_steps"].cast<int>()),
       dt_init(RCONST(py_opts["dt_init"].cast<double>())),
+      dt_min(RCONST(py_opts["dt_min"].cast<double>())),
       dt_max(RCONST(py_opts["dt_max"].cast<double>())),
       max_error_test_failures(py_opts["max_error_test_failures"].cast<int>()),
       max_nonlinear_iterations(py_opts["max_nonlinear_iterations"].cast<int>()),

diff --git a/src/pybammsolvers/idaklu_source/Options.hpp b/src/pybammsolvers/idaklu_source/Options.hpp
@@ -31,6 +31,7 @@ struct SolverOptions {
   int max_order_bdf;
   int max_num_steps;
   double dt_init;
+  double dt_min;
   double dt_max;
   int max_error_test_failures;
   int max_nonlinear_iterations;