#1477 only call IDAGetSens if calculating sensitivities

martinjrobins · martinjrobins · commit b3a30914caa0 · 2021-06-14T21:39:40.000+01:00
diff --git a/pybamm/solvers/c_solvers/idaklu.cpp b/pybamm/solvers/c_solvers/idaklu.cpp
@@ -445,7 +445,9 @@ Solution solve(np_array t_np, np_array y0_np, np_array yp0_np,
 
     if (retval == IDA_TSTOP_RETURN || retval == IDA_SUCCESS || retval == IDA_ROOT_RETURN)
     {
-      IDAGetSens(ida_mem, &tret, yyS);
+      if (number_of_parameters > 0) {
+        IDAGetSens(ida_mem, &tret, yyS);
+      }
 
       t_return[t_i] = tret;
       for (int j = 0; j < number_of_states; j++)

Original file line number	Diff line number	Diff line change
`@@ -445,7 +445,9 @@ Solution solve(np_array t_np, np_array y0_np, np_array yp0_np,`
`445`	`445`
`446`	`446`	`if (retval == IDA_TSTOP_RETURN \|\| retval == IDA_SUCCESS \|\| retval == IDA_ROOT_RETURN)`
`447`	`447`	`{`
`448`		`- IDAGetSens(ida_mem, &tret, yyS);`
	`448`	`+ if (number_of_parameters > 0) {`
	`449`	`+ IDAGetSens(ida_mem, &tret, yyS);`
	`450`	`+ }`
`449`	`451`
`450`	`452`	`t_return[t_i] = tret;`
`451`	`453`	`for (int j = 0; j < number_of_states; j++)`