interpolation diffeqarray

SciML · Jul 1, 2017 · 04bace8 · 04bace8
1 parent 911ec37
commit 04bace8
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 4 deletions.
diff --git a/src/interpolation.jl b/src/interpolation.jl
@@ -63,7 +63,7 @@ interp_summary(sol::DESolution) = interp_summary(sol.interp)
       end
     end
   end
-  vals
+  DiffEqArray(vals,tvals)
 end
 
 """

diff --git a/src/solutions/ode_solutions.jl b/src/solutions/ode_solutions.jl
@@ -84,9 +84,9 @@ function calculate_solution_errors!(sol::AbstractODESolution;fill_uanalytic=true
       if sol.dense && dense_errors
         densetimes = collect(linspace(sol.t[1],sol.t[end],100))
         interp_u = sol(densetimes)
-        interp_analytic = [f(Val{:analytic},t,sol.u[1]) for t in densetimes]
-        sol.errors[:L∞] = maximum(vecvecapply((x)->abs.(x),interp_u-interp_analytic))
-        sol.errors[:L2] = sqrt(recursive_mean(vecvecapply((x)->float.(x).^2,interp_u-interp_analytic)))
+        interp_analytic = VectorOfArray([f(Val{:analytic},t,sol.prob.u0) for t in densetimes])
+        sol.errors[:L∞] = maximum(abs.(interp_u.-interp_analytic))
+        sol.errors[:L2] = sqrt(mean((interp_u.-interp_analytic).^2))
       end
     end
   end