From f8b40b659152e266b1face04974e803257cbfc71 Mon Sep 17 00:00:00 2001 From: "Alberto F. Martin" Date: Fri, 8 Nov 2024 18:39:27 +1100 Subject: [PATCH] Minor --- test/mpi/RefinedPatchDarcyCubedSphereTests.jl | 12 +++--------- 1 file changed, 3 insertions(+), 9 deletions(-) diff --git a/test/mpi/RefinedPatchDarcyCubedSphereTests.jl b/test/mpi/RefinedPatchDarcyCubedSphereTests.jl index e42e37d..7a4cc85 100644 --- a/test/mpi/RefinedPatchDarcyCubedSphereTests.jl +++ b/test/mpi/RefinedPatchDarcyCubedSphereTests.jl @@ -74,7 +74,7 @@ module RefinedPatchDarcyCubedSphereTestsMPI uh,ph=xh eph = ph-p_ex euh = uh-u_ex - degree=10 #2*order + degree=2*order Ω = Triangulation(model) dΩ = Measure(Ω, degree) err_p = sum(∫(eph*eph)dΩ) @@ -102,13 +102,6 @@ module RefinedPatchDarcyCubedSphereTestsMPI coarse_model, cell_panels, coarse_cell_wise_vertex_coordinates= parse_cubed_sphere_coarse_model("connectivity-gridapgeo.txt","geometry-gridapgeo.txt") - - #println(coarse_cell_wise_vertex_coordinates_data) - #println(coarse_cell_wise_vertex_coordinates_ptrs) - # model,_=adapt_model(ranks,model) - # writevtk(model,"model_adapted_1") - # model,_=adapt_model(ranks,model) - # writevtk(model,"model_adapted_2") GridapPETSc.with(args=split(petsc_mumps_options())) do order=0 num_uniform_refinements=0 @@ -121,8 +114,9 @@ module RefinedPatchDarcyCubedSphereTestsMPI radius=1.0, adaptive=true, order=1) - op=assemble_darcy_problem(model, order) + op=assemble_darcy_problem(model, order) xh=solve_darcy_problem(op) + println(compute_darcy_errors(model, order, xh)) for num_uniform_refinements=1:3 model,_=adapt_model(ranks,model) op=assemble_darcy_problem(model, order)