ODE: wrong nullptr passed as argument

kokkos · Jun 12, 2024 · 6a9d85e · 6a9d85e
1 parent 55f44a7
commit 6a9d85e
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/ode/src/KokkosODE_BDF.hpp b/ode/src/KokkosODE_BDF.hpp
@@ -103,6 +103,7 @@ struct BDF {
 
     const double dt = (t_end - t_start) / num_steps;
     double t        = t_start;
+    int count = 0;
 
     // Load y0 into y_vecs(:, 0)
     for (int eqIdx = 0; eqIdx < ode.neqs; ++eqIdx) {
@@ -118,7 +119,7 @@ struct BDF {
     KokkosODE::Experimental::ODE_params params(table.order - 1);
     for (int stepIdx = 0; stepIdx < init_steps; ++stepIdx) {
       KokkosODE::Experimental::RungeKutta<RK_type::RKF45>::Solve(
-          ode, params, t, t + dt, y0, y, update, kstack, nullptr);
+          ode, params, t, t + dt, y0, y, update, kstack, &count);
 
       for (int eqIdx = 0; eqIdx < ode.neqs; ++eqIdx) {
         y_vecs(eqIdx, stepIdx + 1) = y(eqIdx);