Plot solutions with failures

CliMA · Nov 18, 2024 · e25f26f · e25f26f
1 parent 5801f31
commit e25f26f
Showing 1 changed file with 8 additions and 1 deletion.
diff --git a/docs/src/dev/compute_convergence.jl b/docs/src/dev/compute_convergence.jl
@@ -274,7 +274,14 @@ function compute_convergence!(
         @show default_dt
         @show count(isnan, plot2_values.u[end])
         @show length(plot2_values.u[end])
-        error("NaN found in plot2_values in problem $(test_name)")
+        @show count(isnan, plot2_values.u[1])
+        @show length(plot2_values.u[1])
+        out_path = joinpath("output", "convergence_failure")
+        mkpath(out_path)
+        fname(i) = "$(key1)_$(key2)_step_$(i).png"
+        for (i, u) in enumerate(plot2_values.u)
+            Plots.png(Plots.plot(u), joinpath(out_path, fname(i)))
+        end
     end
     out_dict[key1][key2]["plot2_values"] = plot2_values