Skip to content

Commit

Permalink
More tests
Browse files Browse the repository at this point in the history
  • Loading branch information
JordiManyer committed Oct 12, 2024
1 parent eae50d2 commit 05eb506
Show file tree
Hide file tree
Showing 3 changed files with 24 additions and 7 deletions.
14 changes: 7 additions & 7 deletions test/AdaptivityTests/FineToCoarseFieldsTests.jl
Original file line number Diff line number Diff line change
Expand Up @@ -91,13 +91,13 @@ eh2 = u_c - u_fc
eh3 = u_c - u_fc2
@test sum((eh3eh3)*dΩ_c) < 1.e-12

modelH=CartesianDiscreteModel((0,1,0,1),(1,1))
modelh=refine(modelH,2)
reffe=LagrangianRefFE(Float64,QUAD,1)
XH = TestFESpace(modelH,reffe)
xH = get_fe_basis(XH)
modelH = CartesianDiscreteModel((0,1,0,1),(1,1))
modelh = refine(modelH,2)
reffe = LagrangianRefFE(Float64,QUAD,1)
XH = TestFESpace(modelH,reffe)
xH = get_fe_basis(XH)
xHh = change_domain(xH,get_triangulation(modelh),ReferenceDomain())
evaluate(Gridap.CellData.get_data(xHh)[1],
[Point(0.0,0.0),Point(0.5,0.5)])
evaluate(Gridap.CellData.get_data(xHh)[1],[Point(0.0,0.0),Point(0.5,0.5)])
evaluate(Gridap.CellData.get_data(xHh)[1],Point(0.5,0.5))

end
5 changes: 5 additions & 0 deletions test/AdaptivityTests/MacroFETests.jl
Original file line number Diff line number Diff line change
Expand Up @@ -13,7 +13,12 @@ function test_macro_reffe(model,fmodel,rrule,order)
reffe = LagrangianRefFE(Float64,poly,order)
sub_reffes = Fill(reffe,num_subcells(rrule))
macro_reffe = Adaptivity.MacroReferenceFE(rrule,sub_reffes)
ReferenceFEs.test_reference_fe(macro_reffe)

macro_quad = Quadrature(poly,Adaptivity.CompositeQuadrature(),rrule,2*order)
macro_quad_bis = Adaptivity.CompositeQuadrature(Quadrature(poly,4*order),rrule)
ReferenceFEs.test_quadrature(macro_quad)
ReferenceFEs.test_quadrature(macro_quad_bis)

Ω = Triangulation(model)
Ωf = Triangulation(fmodel)
Expand Down
12 changes: 12 additions & 0 deletions test/AdaptivityTests/RefinementRulesTests.jl
Original file line number Diff line number Diff line change
Expand Up @@ -24,4 +24,16 @@ for poly in polys
end
end

rr_bc2 = Adaptivity.BarycentricRefinementRule(TRI)
Adaptivity.test_refinement_rule(rr_bc2)

rr_bc3 = Adaptivity.BarycentricRefinementRule(TET)
Adaptivity.test_refinement_rule(rr_bc3)

rr_ps2 = Adaptivity.PowellSabinRefinementRule(TRI)
Adaptivity.test_refinement_rule(rr_ps2)

rr_ps3 = Adaptivity.PowellSabinRefinementRule(TET)
Adaptivity.test_refinement_rule(rr_ps3)

end

0 comments on commit 05eb506

Please sign in to comment.