julia> sol = solve(ode, KenCarp4(), saveat=range(first(tspan), stop=last(tspan), length=200), abstol=1.0e-7, reltol=1.0e-7)
┌ Warning: Split ODE problem do not work with factorization linear solvers. Bug detailed in https://github.com/SciML/OrdinaryDiffEq.jl/pull/1643. Defaulting to linsolve=KrylovJL()
└ @ OrdinaryDiffEq ~/.julia/packages/OrdinaryDiffEq/DGMT0/src/alg_utils.jl:186
I got the message
Xref https://github.com/SciML/OrdinaryDiffEq.jl/pull/1643