search

SKopecz / PositiveIntegrators.jl

A Julia library of positivity-preserving time integration methods
https://skopecz.github.io/PositiveIntegrators.jl/
MIT License
13 stars 4 forks source link

use LinearSolve.jl for in-place `MPRK22` and set default alg to `LUFactorization()` #60

Closed ranocha closed 7 months ago

ranocha commented 7 months ago

Since the default algorithm choice of LinearSolve.jl has quite some overhead (see https://github.com/SKopecz/PositiveIntegrators.jl/issues/59), I set the default algorithm to LUFactorization(). This avoids setting up all the other caches and improves performance.

I didn't specify the default choice. This allows us to change it later in any non-breaking release.