Merge pull request #2505 from SciML/rewrap

Do not rewrap PDESolution with default algorithm

Author: ChrisRackauckas

lib/OrdinaryDiffEqDefault/src/default_alg.jl

@@ -40,10 +40,10 @@ function isdefaultalg(alg::CompositeAlgorithm{
 end
 
 function DiffEqBase.__init(prob::ODEProblem, ::Nothing, args...; kwargs...)
-    DiffEqBase.init(prob, DefaultODEAlgorithm(autodiff = false), args...; kwargs...)
+    DiffEqBase.init(prob, DefaultODEAlgorithm(autodiff = false), args...; wrap = Val(false), kwargs...)
 end
 function DiffEqBase.__solve(prob::ODEProblem, ::Nothing, args...; kwargs...)
-    DiffEqBase.solve(prob, DefaultODEAlgorithm(autodiff = false), args...; kwargs...)
+    DiffEqBase.solve(prob, DefaultODEAlgorithm(autodiff = false), args...; wrap = Val(false), kwargs...)
 end
 
 function is_stiff(integrator, alg, ntol, stol, is_stiffalg, current)