Adaptive time-stepping for transient solver using SUNDIALS

cmake/ExternalMFEM.cmake

@@ -317,7 +317,7 @@ Intel C++ compiler for MUMPS and STRUMPACK dependencies")
     list(APPEND MFEM_OPTIONS
       "-DSUNDIALS_DIR=${CMAKE_INSTALL_PREFIX}"
       "-DSUNDIALS_OPT=-I${CMAKE_INSTALL_PREFIX}/include"
-      "-DSUNDIALS_LIB=-Wl,-rpath,${CMAKE_INSTALL_PREFIX}/lib64 -L${CMAKE_INSTALL_PREFIX}/lib64 -lsundials_arkode, -lsundials_cvode -lsundials_nvecserial -lsundials_kinsol -lsundials_nvecparhyp -lsundials_nvecparallel"
+      "-DSUNDIALS_LIB=-Wl,-rpath,${CMAKE_INSTALL_PREFIX}/lib64 -L${CMAKE_INSTALL_PREFIX}/lib64 -lsundials_arkode, -lsundials_cvode -lsundials_kinsol -lsundials_nvecparhyp -lsundials_nvecparallel"
       "-DSUNDIALS_REQUIRED_PACKAGES=${SUNDIALS_REQUIRED_PACKAGES}"
     )
     if(NOT "${SUNDIALS_REQUIRED_LIBRARIES}" STREQUAL "")

cmake/ExternalSUNDIALS.cmake

@@ -6,7 +6,7 @@
 #
 
 # Force build order
-set(SUNDIALS_DEPENDENCIES )
+set(SUNDIALS_DEPENDENCIES)
 
 set(SUNDIALS_OPTIONS ${PALACE_SUPERBUILD_DEFAULT_ARGS})
 list(APPEND SUNDIALS_OPTIONS
@@ -42,18 +42,11 @@ endif()
 
 if(PALACE_WITH_CUDA)
   list(APPEND SUNDIALS_OPTIONS
+    "-DCMAKE_CUDA_ARCHITECTURES=${CMAKE_CUDA_ARCHITECTURES}"
     "-DENABLE_CUDA=ON"
   )
 endif()
 
-if(PALACE_WITH_MAGMA)
-  list(APPEND SUNDIALS_OPTIONS
-    "-DENABLE_MAGMA=ON"
-    "-DMAGMA_DIR=${MAGMA_DIR}"
-  )
-endif()
-
-
 string(REPLACE ";" "; " SUNDIALS_OPTIONS_PRINT "${SUNDIALS_OPTIONS}")
 message(STATUS "SUNDIALS_OPTIONS: ${SUNDIALS_OPTIONS_PRINT}")
 

docs/src/config/solver.md

@@ -226,8 +226,7 @@ the second-order system of differential equations. The available options are:
     ODE system for the electric field with adaptive time-stepping. This option is only available when *Palace* has been [built with SUNDIALS support](../install.md#Configuration-options).
   - `"CVODE"` :  SUNDIALS CVODE implicit multistep method scheme applied to the first-order
     ODE system for the electric field with adaptive time-stepping. This option is only available when *Palace* has been [built with SUNDIALS support](../install.md#Configuration-options).
-  - `"RungeKutta"` : SUNDIALS ARKode implicit Runge-Kutta scheme applied to the first-order
-    ODE system for the electric field with fixed time-stepping. This option is only available when *Palace* has been [built with SUNDIALS support](../install.md#Configuration-options).
+  - `"RungeKutta"` : Two stage, singly diagonal implicit Runge-Kutta (SDIRK) method. Second order and L-stable.
   - `"Default"` :  Use the default `"GeneralizedAlpha"` time integration scheme.
 
 `"Excitation" [None]` :  Controls the time dependence of the source excitation. The
@@ -262,8 +261,9 @@ disk for [visualization with ParaView](../guide/postprocessing.md#Visualization)
 saved in the `paraview/` directory under the directory specified by
 [`config["Problem"]["Output"]`](problem.md#config%5B%22Problem%22%5D).
 
-`"Order" [2]` :  Order of the Runge-Kutta integrators or maximum order of the multistep method.
-Only relevant when `"Type"` is `"ARKODE"`, `"Runge-Kutta"` or `"CVODE"`.
+`"Order" [0]` :  Order of the adaptive Runge-Kutta integrators or maximum order of the multistep method.
+Only relevant when `"Type"` is `"ARKODE"` or `"CVODE"`. A value less than 1 defaults to 3 for `"ARKODE"`
+and 2 for `"CVODE"`.
 
 `"RelTol" [1e-4]` :  Relative tolerance used in adaptive time-stepping schemes. Only relevant
 when `"Type"` is `"CVODE"` or `"ARKODE"` .

palace/drivers/transientsolver.cpp

@@ -102,7 +102,7 @@ TransientSolver::Solve(const std::vector<std::unique_ptr<Mesh>> &mesh) const
     {
       Mpi::Print("\n");
       t += delta_t;
-      time_op.Init(delta_t);  // Initial conditions
+      time_op.Init();  // Initial conditions
     }
     else
     {