last try ...

SKopecz · Sep 8, 2024 · e18f22b · e18f22b
1 parent 9495976
commit e18f22b
Showing 1 changed file with 12 additions and 56 deletions.
diff --git a/docs/src/convergence.md b/docs/src/convergence.md
@@ -41,64 +41,20 @@ using PrettyTables
 # collect data and create headers
 data = Matrix{Float64}(undef, length(dts), 1 + 2*length(algs))
 data[:,1] = dts
-#=
-header = ["Δt"]
-subheader = [""]
-for i in eachindex(algs)
-    data = [data err[i] [NaN; eoc[i]]]
-    header = [header names[i] names[i]]
-    subheader = [subheader "Error" "EOC"]
-end
+
+#header = ["Δt"]
+#subheader = [""]
+#for i in eachindex(algs)
+#    data = [data err[i] [NaN; eoc[i]]]
+#    header = [header names[i] names[i]]
+#    subheader = [subheader "Error" "EOC"]
+#end
 
 # print table
-pretty_table(data, header = (header, subheader),
-             formatters = (ft_printf("%5.4e", [1, 2, 4, 6, 8]),
-                           ft_printf("%5.4f", [3, 5, 7, 9])))
-=#                           
+#pretty_table(data, header = (header, subheader),
+#             formatters = (ft_printf("%5.4e", [1, 2, 4, 6, 8]),
+#                           ft_printf("%5.4f", [3, 5, 7, 9])))
+                           
 ```
 
 ## Third order MPRK schemes
-
-```@example eoc
-#=
-#select schemes
-algs = [MPRK43I(1.0, 0.5)
-        MPRK43I(0.5, 0.75)
-        MPRK43II(0.5)
-        MPRK43II(2.0 / 3.0)
-        SSPMPRK43()]
-
-names = ["MPRK43I(1.0,0.5)"
-         "MPRK43I(0.5, 0.75)"
-         "MPRK43II(0.5)"
-         "MPRK43II(2.0/3.0)"
-         "SSPMPRK43()"]
-
-
-dts = 0.5 .^ (5:10)
-err = Vector{Vector{Float64}}(undef, length(algs))
-eoc = Vector{Vector{Float64}}(undef, length(algs))
-
-#compute errors and experimental order of convergence
-for i in eachindex(algs)
-    sim = test_convergence(dts, prob, algs[i])
-    err[i] = sim.errors[:l∞]
-    eoc[i] = -log2.(err[i][2:end] ./ err[i][1:(end - 1)])
-end
-
-# collect data and create headers
-data = Matrix{Float64}(undef, length(dts), 1 + 2*length(algs))
-data[:,1] = dts
-header = ["Δt"]
-subheader = [""]
-for i in eachindex(algs)
-    data = [data err[i] [NaN; eoc[i]]]
-    header = [header names[i] names[i]]
-    subheader = [subheader "Error" "EOC"]
-end
-
-pretty_table(data, header = (header, subheader),
-             formatters = (ft_printf("%5.4e", [1, 2, 4, 6, 8, 10]),
-                           ft_printf("%5.4f", [3, 5, 7, 9, 11])))
-=#                           
-```