Usage

Installation

To use pydiffsol, install using pip:

(.venv) $ pip install pydiffsol

Basic Usage

import pydiffsol as ds
import numpy as np

ode = ds.Ode(
   """
   in { r = 1.0 }
   k { 1.0 }
   u { 0.1 }
   F { r * u * (1.0 - u / k) }
   """,
   matrix_type=ds.nalgebra_dense,
)

# Solve up to t = 0.4, overriding r input param = 2.0
params = np.array([2.0])
solution = ode.solve(params, 0.4)
print(solution.ys, solution.ts)

# Above defaults to bdf. Try esdirk34 instead
ode.ode_solver = ds.esdirk34
solution = ode.solve(params, 0.4)
print(solution.ys, solution.ts)