AIBECS.jl: A tool for exploring global marine biogeochemical cycles

AIBECS v0.17.0

Diff since v0.16.7

Now opens up to all solvers from NonlinearSolve.jl and LinearSolve.jl.

Breaking changes

• SciML ecosystem opt-in via a new AIBECSNonlinearSolveExt (NonlinearSolve + ADTypes + DI + SparseConnectivityTracer + SparseMatrixColorings as weak deps):

  • AIBECS.nonlinearproblem(prob) — SteadyStateProblem → NonlinearProblem with sparse jac_prototype (sidesteps NonlinearSolve's iip-promotion of 3-arg f(u, p, t=0)).
  • AIBECS.recommended_nlalg() — NewtonRaphson + UMFPACK with sparse-AD fallback.
  • AIBECS.default_sparse_ad() — overridable.

• CTKAlg overhaul. No longer welded to factorize/\:

  • CTKAlg(linsolve = …) plugs in any LinearSolve algorithm (UMFPACK, KLU, MKL Pardiso, …); the inner loop mutates a LinearSolve.LinearCache so Shamanskii's lazy refresh rides the cache's isfresh toggle.
  • Convergence delegates to NonlinearSolveBase termination modes; abstol/reltol/termination_condition replace the old nrm/τstop kwargs.
  • LinearSolve promoted to direct dep; NonlinearSolveBase added.

• Existing tutorials and SteadyStateProblem usage continue to work modulo the nrm/τstop → abstol/reltol/termination_condition rename on CTKAlg.

Merged pull requests:

• Open up to all SciML nonlinear/linear solvers (#122) (@briochemc)