Ok, so that’s a limitation we can’t get around…
Regarding sparse Jacobians, I read through the Automatic Sparsity Detection section of “Solving Large Stiff Equations · DifferentialEquations.jl” as much as I could, but I got lost when it started talking about brusselators, and getting arrays of random numbers (for some reason) for arrays called input and output. Any tips on how to apply Sparsity to my example program here, I’d be very grateful.