Updated SPRKStep and LSRKStep to only request Lagrange interpolation …

…(not to construct the Lagrange module directly).
LLNL · Jul 16, 2024 · b45e563 · b45e563
1 parent 1943331
commit b45e563
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 10 deletions.
diff --git a/src/arkode/arkode_lsrkstep.c b/src/arkode/arkode_lsrkstep.c
@@ -147,15 +147,8 @@ void* LSRKStepCreate(ARKRhsFn fe, ARKRhsFn fi, sunrealtype t0, N_Vector y0,
     return (NULL);
   }
 
-  /* Specify Interpolation Type */
-  retval = ARKodeSetInterpolantType(ark_mem, ARK_INTERP_LAGRANGE);
-  if (retval != ARK_SUCCESS)
-  {
-    arkProcessError(ark_mem, retval, __LINE__, __func__, __FILE__,
-                    "Unable to specify interpolation type");
-    ARKodeFree((void**)&ark_mem);
-    return (NULL);
-  }
+  /* Specify preferred interpolation type */
+  ark_mem->interp_type = ARK_INTERP_LAGRANGE;
 
   return ((void*)ark_mem);
 }

diff --git a/src/arkode/arkode_sprkstep.c b/src/arkode/arkode_sprkstep.c
@@ -151,7 +151,7 @@ void* SPRKStepCreate(ARKRhsFn f1, ARKRhsFn f2, sunrealtype t0, N_Vector y0,
 
   /* SPRKStep uses Lagrange interpolation by default, since Hermite is
      less compatible with these methods. */
-  ARKodeSetInterpolantType(ark_mem, ARK_INTERP_LAGRANGE);
+  ark_mem->interp_type = ARK_INTERP_LAGRANGE;
 
   /* Initialize main ARKODE infrastructure */
   retval = arkInit(ark_mem, t0, y0, FIRST_INIT);