18 releases

new 0.2.3 Dec 7, 2024
0.2.2 Nov 8, 2024
0.2.0 Aug 31, 2024
0.1.9 Jul 10, 2024
0.1.1 Mar 19, 2024

#279 in Programming languages

Download history 2/week @ 2024-08-23 158/week @ 2024-08-30 5/week @ 2024-09-06 206/week @ 2024-09-13 154/week @ 2024-09-20 36/week @ 2024-09-27 3/week @ 2024-10-04 29/week @ 2024-10-11 48/week @ 2024-10-18 4/week @ 2024-10-25 25/week @ 2024-11-01 243/week @ 2024-11-08 18/week @ 2024-11-15 12/week @ 2024-11-22 26/week @ 2024-11-29

301 downloads per month
Used in diffsol

Custom license and GPL-3.0-or-later

11MB
149K SLoC

C++ 69K SLoC // 0.1% comments LLVM 64K SLoC // 0.4% comments Rust 12K SLoC // 0.0% comments C 3K SLoC // 0.2% comments Bazel 729 SLoC // 0.0% comments Forge Config 253 SLoC // 0.9% comments Julia 145 SLoC // 0.1% comments FORTRAN Modern 113 SLoC // 0.2% comments Python 82 SLoC // 0.3% comments Shell 75 SLoC // 0.0% comments Pest 50 SLoC

DiffSL

CI build status badge

A compiler for a domain-specific language for ordinary differential equations (ODEs) of the following form:

$$ M(t) \frac{d\mathbf{u}}{dt} = F(\mathbf{u}, t) $$

As an example, the following code defines a classic DAE testcase, the Robertson (1966) problem, which models the kinetics of an autocatalytic reaction, given by the following set of equations:

$$ \begin{align} \frac{dx}{dt} &= -0.04x + 10^4 y z \ \frac{dy}{dt} &= 0.04x - 10^4 y z - 3 \cdot 10^7 y^2 \ 0 &= x + y + z - 1 \end{align} $$

The DiffSL code for this problem is as follows:

in = [k1, k2, k3]
k1 { 0.04 }
k2 { 10000 }
k3 { 30000000 }
u_i {
  x = 1,
  y = 0,
  z = 0,
}
dudt_i {
  dxdt = 1,
  dydt = 0,
  dzdt = 0,
}
M_i {
  dxdt,
  dydt,
  0,
}
F_i {
  -k1 * x + k2 * y * z,
  k1 * x - k2 * y * z - k3 * y * y,
  1 - x - y - z,
}
out_i {
  x,
  y,
  z,
}

DiffSL Language Features

See the DiffSL language documentation for a full description.

  • Tensor types:
    • Scalars (double precision floating point numbers)
    • Vectors (1D arrays of scalars)
    • N-dimensional tensor of scalars
    • Sparse/dense/diagonal tensors
  • Tensor operations:
    • Elementwise operations
    • Broadcasting
    • Tensor contractions/matmul/translation etc via index notation

Usage

Generally the easiest way to make use of DiffSL is via an ode solver that supports the language, for example the diffsol library. Please see the diffsol documentation and consult the DiffSL language documentation for more information.

If you are writing your own ode solver and want to make use of the DiffSL compiler, please either get in touch by opening an issue, contacting the author or by looking at the diffsol source code.

Dependencies

You will need to install the LLVM project. The easiest way to install this is to use the package manager for your operating system. For example, on Ubuntu you can install these with the following command:

sudo apt-get install llvm

Installation

You can install DiffSL using cargo. You will need to indicate the llvm version you have installed using a feature flag. For example, for llvm 14:

cargo add diffsl --features llvm14-0

Other versions of llvm are also supported given by the features llvm13-0, llvm14-0, llvm15-0, llvm16-0, llvm17-0.

Dependencies

~13–22MB
~345K SLoC