*src_ypred,
*src_fpred,
jcurPtr_tmp);
-# if DEAL_II_SUNDIALS_VERSION_GTE(2, 0, 0)
- *jcurPtr = jcurPtr_tmp ? SUNTRUE : SUNFALSE;
-# else
- *jcurPtr = jcurPtr_tmp ? TRUE : FALSE;
-# endif
+ *jcurPtr = jcurPtr_tmp ? SUNTRUE : SUNFALSE;
return err;
}
int
solve_with_jacobian_callback(ARKodeMem arkode_mem,
N_Vector b,
-# if DEAL_II_SUNDIALS_VERSION_LT(3, 0, 0)
- N_Vector,
-# endif
- N_Vector ycur,
- N_Vector fcur)
+ N_Vector ycur,
+ N_Vector fcur)
{
Assert(arkode_mem->ark_user_data != nullptr, ExcInternalError());
ARKode<VectorType> &solver =
template <typename VectorType>
int
- solve_with_mass_matrix_callback(ARKodeMem arkode_mem,
-# if DEAL_II_SUNDIALS_VERSION_LT(3, 0, 0)
- N_Vector b,
- N_Vector
-# else
- N_Vector b
-# endif
- )
+ solve_with_mass_matrix_callback(ARKodeMem arkode_mem, N_Vector b)
{
Assert(arkode_mem->ark_user_data != nullptr, ExcInternalError());
ARKode<VectorType> &solver =
if (setup_jacobian)
{
ARKode_mem->ark_lsetup = setup_jacobian_callback<VectorType>;
-# if DEAL_II_SUNDIALS_VERSION_LT(3, 0, 0)
- ARKode_mem->ark_setupNonNull = true;
-# endif
}
}
else
if (setup_mass)
{
ARKode_mem->ark_msetup = setup_mass_matrix_callback<VectorType>;
-# if DEAL_II_SUNDIALS_VERSION_LT(3, 0, 0)
- ARKode_mem->ark_MassSetupNonNull = true;
-# endif
}
}