Fix ensemble indexing

Author: ChrisRackauckas

Project.toml

@@ -1,7 +1,7 @@
 name = "SciMLBase"
 uuid = "0bca4576-84f4-4d90-8ffe-ffa030f20462"
 authors = ["Chris Rackauckas <accounts@chrisrackauckas.com> and contributors"]
-version = "2.48.0"
+version = "2.48.1"
 
 [deps]
 ADTypes = "47edcb42-4c32-4615-8424-f2b9edc5f35b"

src/ensemble/ensemble_solutions.jl

@@ -116,7 +116,7 @@ function calculate_ensemble_errors(u; elapsedTime = 0.0, converged = false,
                            sol.W(densetimes[i])[1])
                        for i in eachindex(densetimes)] for sol in u]
         udense = [u[j](densetimes) for j in 1:length(u)]
-        dense_weak_errors = [mean([udense[j][i] - u_analytic[j][i] for j in 1:length(u)])
+        dense_weak_errors = [mean([udense[j].u[i] - u_analytic[j][i] for j in 1:length(u)])
                              for i in eachindex(densetimes)]
         dense_L2_errors = [sqrt.(sum(abs2, err) / length(err)) for err in dense_weak_errors]
         L2_tmp = sqrt(sum(abs2, dense_L2_errors) / length(dense_L2_errors))