Jason Rader
Jason Rader
Fixes an edge case where `dfx.LinearInterpolation` will return a `nan` value when input arrays have a repeated leading constant value. ie. the following example will return a `nan` value in...
Changes: - Add the Störmer-Verlet method - Add `all_symplectic_solvers` to `test/helpers`. - Change `test_semi_implicit_euler` to a parameterised test `test_symplectic_solvers` which includes Störmer-Verlet - Add an order test `test_symplectic_ode_order` similar to...
GMRES often fails to converge when restart is less than the problem dimension, but will return **without indicating failure.** This failure (likely due to stagnation) seems to depend on the...
I mentioned here https://github.com/JuliaLinearAlgebra/IterativeSolvers.jl/issues/338 that GMRES was failing silently from stagnation. Turns out BiCGStab is also failing silently from (some form of) breakdown. Not so surprising that this fails given...
GMRES often fails to converge when restart is less than the problem dimension, but will return **without indicating failure.** This failure (likely due to stagnation) seems to depend on the...
When solving against matrices of the form $A + \lambda I$ with $\lambda$ set to the max dtype value, QR solves return `nan` instead of zero. This does not occur...
Vmapping across `y0` with any method using `AbstractGaussNewton` throws a TypeError. MWE ```python import jax import jax.numpy as jnp import optimistix as optx def rosenbrock(x, args): del args term1 =...