Fix initial guess tests

Author: ErikQQY

lib/BoundaryValueDiffEqCore/src/utils.jl

@@ -625,7 +625,7 @@ end
 end
 @inline function __initial_guess_on_mesh(u₀::SciMLBase.ODESolution, mesh, p; tune_parameters = false)
     tune_parameters && return VectorOfArray([vcat(vec(u), __tunable_part(p)) for u in u₀.u])
-    return VectorOfArray([deepcopy(u) for u in u₀.u])
+    return VectorOfArray(deepcopy(u₀.u))
 end
 @inline function __initial_guess_on_mesh(u₀::AbstractArray, mesh, p; tune_parameters = false)
     tune_parameters && return VectorOfArray([vcat(vec(u₀), __tunable_part(p)) for _ in mesh])

lib/BoundaryValueDiffEqMIRK/test/mirk_basic_tests.jl

@@ -714,5 +714,5 @@ end
 
     # Ensure the initial guess is not modified in-place
     sol2 = solve(bvp2, MIRK4(), dt = 0.1, adaptive = false)
-    @test bvp2.u0 == u_guess
+    @test bvp2.u0.u == u_guess
 end