Add in the images

SciML · Feb 27, 2023 · 6ab2cd8 · 6ab2cd8
1 parent 7133dd2
commit 6ab2cd8
Show file tree

Hide file tree

Showing 17 changed files with 380 additions and 369 deletions.
diff --git a/docs/src/mean_exit_time.md b/docs/src/mean_exit_time.md
@@ -100,6 +100,7 @@ using CairoMakie
 using OrdinaryDiffEq
 using SteadyStateDiffEq
 using LinearSolve
+using LinearAlgebra
 using FastLapackInterface
 
 R = 1.0
@@ -359,7 +360,7 @@ prob = FVMProblem(mesh, BCs;
     final_time=final_time,
     initial_condition=initial_guess,
     steady=true)
-sol = solve(prob, DynamicSS(TRBDF2(linsolve=KrylovJL_GMRES())))
+sol = solve(prob, DynamicSS(TRBDF2(linsolve=KrylovJL_GMRES())), parallel = true)
 fig = Figure(fontsize=38)
 ax = Axis(fig[1, 1], xlabel=L"x", ylabel=L"y", width=600, height=300)
 pt_mat = Matrix(get_points(tri)')

diff --git a/test/figures/annulus_mean_exit_time.png b/test/figures/annulus_mean_exit_time.png
diff --git a/test/figures/annulus_test.png b/test/figures/annulus_test.png
diff --git a/test/figures/ellipse_mean_exit_time.png b/test/figures/ellipse_mean_exit_time.png
diff --git a/test/figures/laplace_equation.png b/test/figures/laplace_equation.png
diff --git a/test/figures/perturbed_annulus_mean_exit_time.png b/test/figures/perturbed_annulus_mean_exit_time.png
diff --git a/test/figures/perturbed_circle_mean_exit_time.png b/test/figures/perturbed_circle_mean_exit_time.png
diff --git a/test/figures/perturbed_ellipse_mean_exit_time.png b/test/figures/perturbed_ellipse_mean_exit_time.png
diff --git a/test/mean_exit_time.jl b/test/mean_exit_time.jl
diff --git a/test/met_annulus.jl b/test/met_annulus.jl
@@ -0,0 +1,54 @@
+
+using ..FiniteVolumeMethod
+using DelaunayTriangulation
+using CairoMakie
+using OrdinaryDiffEq
+using SteadyStateDiffEq
+using Test
+using LinearAlgebra
+using StatsBase
+using LinearSolve
+using Krylov
+include("test_setup.jl")
+R₁ = 2.0
+R₂ = 3.0
+θ = LinRange(0, 2π, 250)
+x = [
+    [R₂ .* cos.(θ)],
+    [reverse(R₁ .* cos.(θ))] # inner boundaries are clockwise
+]
+y = [
+    [R₂ .* sin.(θ)],
+    [reverse(R₁ .* sin.(θ))] # inner boundaries are clockwise
+]
+tri = generate_mesh(x, y, 0.2; gmsh_path=GMSH_PATH)
+mesh = FVMGeometry(tri)
+outer_bc = (x, y, t, u, p) -> zero(u)
+inner_bc = (x, y, t, u, p) -> zero(u)
+type = [:D, :N]
+BCs = BoundaryConditions(mesh, [outer_bc, inner_bc], type)
+D = 6.25e-4
+initial_guess = zeros(num_points(tri))
+diffusion_function = (x, y, t, u, D) -> D
+diffusion_parameters = D
+reaction = (x, y, t, u, p) -> one(u)
+final_time = Inf
+prob = FVMProblem(mesh, BCs;
+    diffusion_function,
+    diffusion_parameters,
+    reaction_function=reaction,
+    final_time=final_time,
+    initial_condition=initial_guess,
+    steady=true)
+sol = solve(prob, DynamicSS(TRBDF2(linsolve=KrylovJL_GMRES())))
+fig = Figure(fontsize=38)
+ax = Axis(fig[1, 1], xlabel=L"x", ylabel=L"y", width=600, height=300)
+pt_mat = Matrix(get_points(tri)')
+T_mat = [T[j] for T in each_triangle(tri), j in 1:3]
+msh = mesh!(ax, pt_mat, T_mat, color=sol, colorrange=(0, 900))
+Colorbar(fig[1, 2], msh)
+resize_to_layout!(fig)
+SAVE_FIGURE && save("figures/annulus_mean_exit_time.png", fig)
+exact = [(R₂^2 - norm(p)^2) / (4D) + R₁^2 * log(norm(p) / R₂) / (2D) for p in each_point(tri)]
+error_fnc = (fvm_sol, exact_sol) -> 100abs(fvm_sol - exact_sol) / maximum(exact)
+@test median(error_fnc.(exact, sol)) < 0.1