Sundials/ML is an interface to the Sundials suite of numerical solvers.
Full documentation is online.
Sundials is a collection of six numerical solvers: CVODE, CVODES, IDA, IDAS, ARKODE, and KINSOL. This interface provides access to all features of the underlying library except those related to GPUs and external libraries which have no existing OCaml support (see "Known Limitations" below).
Sundials/ML is more than just a simple low-level wrapper around a subset of the functions and data structures of Sundials. It exploits the type system of OCaml to document the underlying library and prevent some utilization errors. Additional dynamic checks detect and signal other problems in user code. The module system of OCaml is exploited to structure the interface and increase its convenience. Exceptions are used for detecting and handling error conditions. Sundials/ML allows users to define nvectors, linear solvers, and nonlinear solvers in OCaml. It provides access to and reimplementations of nvector and matrix routines. The OCaml interface is completely documented with references back to the Sundials library and documentation.
Many programs just use a small subset of Sundials functions, e.g., a handful of CVODE functions with serial nvectors, so Sundials/ML tries to isolate OCaml programs from differences in the underlying library with support back to Sundials 2.7.0. The idea is to program in OCaml to the most up-to-date Sundials/ML interface (which evolves as Sundials does) but to be able to install on platforms where an older version of Sundials is installed. Functions are provided to query the version and features of the underlying library and exceptions are raised if use of an unsupported feature is attempted.
See our detailed notes, or try:
- Download Sundials, extract, and install it:
cmake -Wno-dev ../sundials-5.8.0, optionally adding:
- Either install from OPAM:
opam install sundialsml, or
- Download Sundials/ML, extract, and install it:
- Start coding!
Our goal is to provide access from OCaml to as much of Sundials as we can. We do not, however, expose features for which there is no corresponding OCaml library (e.g., some of the nvector modules) or whose implementation would be overly complicated (e.g., in terms of types or mixed C/OCaml data structures) or inefficient.
- Unsupported nvector modules: CUDA, Hypre ParVector, PETSC, RAJA, Trilinos, HIP, and SYCLC. (If you provide the required OCaml interface and help with testing, we'll add then vector).
- No support for CUDA (including SUNMATRIX_CUSPARSE and SUNLinearSolver_cuSolverSp_batchQR linear solver).
- No support for the MAGMA linear algebra library (if you provide the required OCaml interface, we'll add support).
- No support for fused kernels.
- No support for the SUNLinSol_KLUGetSymbolic, SUNLinSol_KLUGetNumeric, and SUNLinSol_KLUGetCommon functions (patches welcome).
- No support for SUNMatrix_SLUNRloc and SUNLinearSolver_SuperLUDIST (if you provide a solid OCaml interface to SuperLU_DIST, we'll add support for this matrix type and linear solver).
- No support for the XBraid library or the three associated NVector operations: N_VBufSize, N_VBufPack, and N_VBufUnpack.
- No support for the N_VGetDeviceArrayPointer operation.