JuliaParallel · ranocha · Apr 25, 2026 · Apr 26, 2026 · Apr 26, 2026 · Apr 26, 2026
diff --git a/CITATION.bib b/CITATION.bib
@@ -2,12 +2,12 @@ @software{petsc_jl
   author       = {Kaus, Boris and
                   Shah, Viral B. and
                   Kozdon, Jeremy E. and
-                  Churavy, Valentin and 
+                  Churavy, Valentin and
                   Ranocha, Hendrik and
                   Schnetter, Erik and
                   Byrne, Simon},
-  title        = {PETSc.jl: Julia bindings for PETSc},
-  year         = 2026,
+  title        = {{PETSc.jl}: {J}ulia bindings for {PETSc}},
+  year         = {2026},
   publisher    = {Zenodo},
   doi          = {10.5281/zenodo.18274810},
   url          = {https://github.com/JuliaParallel/PETSc.jl}

diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "PETSc"
 uuid = "ace2c81b-2b5f-4b1e-a30d-d662738edfe0"
 version = "0.4.10"
-authors = ["Boris Kaus <kaus@uni-mainz.de>", "Viral B. Shah <virals@gmail.com>", "Valentin Churavy <v.churavy@gmail.com>", "Erik Schnetter <eschnetter@perimeterinstitute.ca>", "Jeremy E. Kozdon <jeremy@kozdon.net>", "Simon Byrne <simonbyrne@gmail.com>"]
+authors = ["Boris Kaus <kaus@uni-mainz.de>", "Viral B. Shah <virals@gmail.com>", "Jeremy E. Kozdon <jeremy@kozdon.net>", "Valentin Churavy <v.churavy@gmail.com>", "Hendrik Ranocha <hendrik.ranocha@uni-mainz.de>", "Erik Schnetter <eschnetter@perimeterinstitute.ca>", "Simon Byrne <simonbyrne@gmail.com>"]
 
 [deps]
 ForwardDiff = "f6369f11-7733-5829-9624-2563aa707210"
@@ -19,9 +19,11 @@ UnicodePlots = "b8865327-cd53-5732-bb35-84acbb429228"
 
 [weakdeps]
 CUDA = "052768ef-5323-5732-b1bb-66c8b64840ba"
+SciMLBase = "0bca4576-84f4-4d90-8ffe-ffa030f20462"
 
 [extensions]
 PETScCUDAExt = "CUDA"
+PETScSciMLExt = "SciMLBase"
 
 [compat]
 CUDA = "5, 6"
@@ -34,6 +36,7 @@ OffsetArrays = "1.0"
 PETSc_jll = "3.22"
 Pkg = "^1.10"
 Preferences = "1"
+SciMLBase = "2"
 SparseArrays = "1.10"
 Statistics = "^1.10"
 UnicodePlots = "3.0"
@@ -46,9 +49,10 @@ KernelAbstractions = "63c18a36-062a-441e-b654-da1e3ab1ce7c"
 Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
 Printf = "de0858da-6303-5e67-8744-51eddeeeb8d7"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
+SciMLBase = "0bca4576-84f4-4d90-8ffe-ffa030f20462"
 SparseDiffTools = "47a9eef4-7e08-11e9-0b38-333d64bd3804"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 UnicodePlots = "b8865327-cd53-5732-bb35-84acbb429228"
 
 [targets]
-test = ["ForwardDiff", "UnicodePlots", "Test", "Plots", "SparseDiffTools", "Printf", "Random", "CairoMakie", "KernelAbstractions"]
+test = ["CairoMakie", "ForwardDiff", "KernelAbstractions", "Plots", "Printf", "Random", "SciMLBase", "SparseDiffTools", "Test", "UnicodePlots"]
diff --git a/README.md b/README.md
@@ -6,17 +6,21 @@
 [![DOI](https://zenodo.org/badge/38933145.svg)](https://doi.org/10.5281/zenodo.18274809)
 
 
-`PETSc.jl` provides an interface to the Portable, Extensible Toolkit for Scientific Computation ([PETSc](https://petsc.org)) library, allowing the combination of Julia features (such as automatic differentiation) with the PETSc's infrastructure, including linear, nonlinear, and optimization solvers, timesteppers, domain management (DM), and more, in a distributed-memory (MPI) environment. 
+PETSc.jl provides an interface to the Portable, Extensible Toolkit for Scientific Computation ([PETSc](https://petsc.org)) library, allowing the combination of Julia features (such as automatic differentiation) with the PETSc's infrastructure, including linear, nonlinear, and optimization solvers, timesteppers, domain management (DM), and more, in a distributed-memory (MPI) environment.
 
-This package comprises two main components:
+This package comprises three main components:
 
 1. An automatically generated, low-level interface for large parts of the PETSc API (see `PETSc.LibPETSc`).
-2. A curated, high-level, more Julianic interface for selected functionality.
+2. A curated, high-level, more Julian interface for selected functionality.
+3. A package extension based on [SciMLBase.jl](https://github.com/SciML/SciMLBase.jl) that allows solving problems such as `ODEProblem`s with PETSc's algorithms.
 
-The low-level interface covers nearly the entire PETSc API, but may be awkward to work with and likely requires previous experience with PETSc to use effectively. The high level interface is designed to be more familiar and convenient for Julia users, and allows, for example, to set matrix entries with `A[1,2] = 3.0`, rather than having to call `LibPETSc.MatSetValue`. It, however, exposes only a small portion of the functionality of the underlying library. 
+The low-level interface covers nearly the entire PETSc API, but may be awkward to work with and likely requires previous experience with PETSc to use effectively.
+The high-level interface is designed to be more familiar and convenient for Julia users, and allows, for example, to set matrix entries with `A[1,2] = 3.0`, rather than having to call `LibPETSc.MatSetValue`.
+However, it exposes only a small portion of the functionality of the underlying library.
+The SciML package extension is work in progress and currently supports only (in-place) `ODEProblem`s.
 
 ## Installation
-This package can be added with the julia command:
+This package can be added with the Julia command:
 ```julia
 julia>]add PETSc
 ```
@@ -45,10 +49,29 @@ julia>using PETSc
 julia>[PETSc.petsclibs...]
 ```
 
-## Windows users 
-The package currently does not work on windows, mainly because `MicrosoftMPI_jll` does not function when used along with the precompiled version used in `PETSc_jll`. Windows users are therefore advised to install the [Windows Subsystem for Linux](https://en.wikipedia.org/wiki/Windows_Subsystem_for_Linux) (WSL) and run PETSc through there. 
+## Windows users
+The package currently does not work on Windows, mainly because `MicrosoftMPI_jll` does not function when used along with the precompiled version used in `PETSc_jll`.
+Windows users are therefore advised to install the [Windows Subsystem for Linux](https://en.wikipedia.org/wiki/Windows_Subsystem_for_Linux) (WSL) and run PETSc through there.
 
 ## Getting started
-Have a look at the [documentation](https://juliaparallel.org/PETSc.jl/stable/), at the [examples](./examples/) directory or at the tests in the [test](./test) directory. We do keep the tests up to date, so that is a good starting point. 
+Have a look at the [documentation](https://juliaparallel.org/PETSc.jl/stable/), at the [examples](./examples/) directory or at the tests in the [test](./test) directory.
+We do keep the tests up to date, so that is a good starting point.
 
-cause Note that we do not have tests in place for the whole library at this stage. The best supported parts are `DMDA`, `DMStag`, `DMPlex`, `KSP`, `SNES`, `Vec`, and `Mat`. Other DM types (DMForest, DMNetwork, DMSwarm) do not yet have a high-level interface; users will have to rely on the low-level `LibPETSc` interface for those.
+Note that we do not have tests in place for the whole library at this stage. The best supported parts are `DMDA`, `DMStag`, `DMPlex`, `KSP`, `SNES`, `Vec`, and `Mat`. Other DM types (DMForest, DMNetwork, DMSwarm) do not yet have a high-level interface; users will have to rely on the low-level `LibPETSc` interface for those.
+
+## SciML integration
+
+PETSc.jl ships a package extension (`PETScSciMLExt`) that lets you solve in-place `ODEProblem`s with PETSc's TS time integrators through the standard SciML interface.
+The extension activates automatically when SciMLBase.jl is loaded, e.g., when other ODE solver packages from the SciML ecosystem are used.
+
+```julia
+using PETSc, SciMLBase
+
+f!(du, u, p, t) = (du[1] = -u[1]; nothing)
+prob = ODEProblem(f!, [1.0], (0.0, 1.0))
+
+sol = solve(prob, PETSc.TSRK("3bs"))
+sol = solve(prob, PETSc.TSImplicit("bdf", ["-snes_fd"]); dt = 1e-3)
+```
+
+See the [SciML Integration](https://juliaparallel.github.io/PETSc.jl/dev/man/sciml/) documentation page for the full API, all supported algorithm types, keyword arguments, and current limitations.
diff --git a/docs/Project.toml b/docs/Project.toml
@@ -2,3 +2,4 @@
 Documenter = "e30172f5-a6a5-5a46-863b-614d45cd2de4"
 MPI = "da04e1cc-30fd-572f-bb4f-1f8673147195"
 PETSc = "ace2c81b-2b5f-4b1e-a30d-d662738edfe0"
+SciMLBase = "0bca4576-84f4-4d90-8ffe-ffa030f20462"
diff --git a/docs/make.jl b/docs/make.jl
@@ -4,7 +4,13 @@ makedocs(;
     modules=[PETSc],
     sitename="PETSc.jl",
     checkdocs=:exports,  # Only check exported functions, skip LibPETSc internals
-    warnonly=true,  # Warn but don't error for any documentation issues
+    # Warn (don't error) for the documentation issues the large auto-wrapped
+    # LibPETSc API trips, but let `@example` / `@repl` block failures error the
+    # build so the executed code examples (e.g. man/sciml.md) stay functional.
+    warnonly=[
+        :autodocs_block, :cross_references, :docs_block, :eval_block,
+        :footnote, :linkcheck, :linkcheck_remotes, :meta_block, :missing_docs,
+    ],
     format=Documenter.HTML(;
         prettyurls=get(ENV, "CI", "false") == "true",
         size_threshold_warn = nothing,  # Disable size warnings for large low-level API pages
@@ -14,6 +20,7 @@ makedocs(;
         "Home" => "index.md",
         "Installation" => "man/installation.md",
         "Getting Started" => "man/getting_started.md",
+        "SciML Integration" => "man/sciml.md",
         "High-level interface" => Any[
             "Vec" =>  "man/vec.md",
             "Mat" =>  "man/mat.md",