Merge pull request #133 from SciML/ap/infnorm

Add norm(_, Inf) to WorkPrecision

Author: ChrisRackauckas

Project.toml

@@ -1,7 +1,7 @@
 name = "DiffEqDevTools"
 uuid = "f3b72e0c-5b89-59e1-b016-84e28bfd966d"
 authors = ["Chris Rackauckas <[email protected]>"]
-version = "2.43.0"
+version = "2.44.0"
 
 [deps]
 DiffEqBase = "2b5f629d-d688-5b77-993f-72d75c75574e"

src/benchmark.jl

@@ -415,12 +415,11 @@ function WorkPrecision(prob::NonlinearProblem, alg, abstols, reltols, dts = noth
 
             stats[i] = sol.stats
 
+            err = appxsol === nothing ? sol.resid : (sol .- appxsol)
             if error_estimate == :l2
-                if isnothing(appxsol)
-                    errors[i] = Dict(error_estimate => sqrt(sum(abs2, sol.resid)))
-                else
-                    errors[i] = Dict(error_estimate => sqrt(sum(abs2, sol .- appxsol)))
-                end
+                errors[i] = Dict(error_estimate => norm(err, 2))
+            elseif error_estimate == :l∞ || error_estimate == :linf
+                errors[i] = Dict(error_estimate => norm(err, Inf))
             else
                 error("Unsupported norm used: $(error_estimate).")
             end