Skip to content

Commit

Permalink
Adding a mock test in FineToCoarseFieldsTests
Browse files Browse the repository at this point in the history
  • Loading branch information
amartinhuertas committed Jan 10, 2025
1 parent b4c3619 commit c7fe919
Showing 1 changed file with 87 additions and 0 deletions.
87 changes: 87 additions & 0 deletions test/AdaptivityTests/FineToCoarseFieldsTests.jl
Original file line number Diff line number Diff line change
Expand Up @@ -100,4 +100,91 @@ 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.5,0.5))

function test_PR_1074()

function setup_coarse_discrete_model()
ptr = [ 1, 5 ]
data = [ 1,2,3,4 ]
cell_vertex_lids = Gridap.Arrays.Table(data,ptr)
node_coordinates = Vector{Point{2,Float64}}(undef,4)
node_coordinates[1]=Point{2,Float64}(0.0,0.0)
node_coordinates[2]=Point{2,Float64}(1.0,0.0)
node_coordinates[3]=Point{2,Float64}(0.0,1.0)
node_coordinates[4]=Point{2,Float64}(1.0,1.0)

polytope=QUAD
scalar_reffe=Gridap.ReferenceFEs.ReferenceFE(polytope,Gridap.ReferenceFEs.lagrangian,Float64,1)
cell_types=collect(Fill(1,length(cell_vertex_lids)))
cell_reffes=[scalar_reffe]
grid = Gridap.Geometry.UnstructuredGrid(node_coordinates,
cell_vertex_lids,
cell_reffes,
cell_types,
Gridap.Geometry.NonOriented())
Gridap.Geometry.UnstructuredDiscreteModel(grid)
end

function setup_adapted_discrete_model(parent)
ptr = [ 1, 5, 9 ]
data = [ 1,2,3,4, 2,5,4,6 ]
cell_vertex_lids = Gridap.Arrays.Table(data,ptr)
node_coordinates = Vector{Point{2,Float64}}(undef,6)
node_coordinates[1]=Point{2,Float64}(0.0,0.5)
node_coordinates[2]=Point{2,Float64}(0.5,0.5)
node_coordinates[3]=Point{2,Float64}(0.0,1.0)
node_coordinates[4]=Point{2,Float64}(0.5,1.0)
node_coordinates[5]=Point{2,Float64}(1.0,0.5)
node_coordinates[6]=Point{2,Float64}(1.0,1.0)

polytope=QUAD
scalar_reffe=Gridap.ReferenceFEs.ReferenceFE(polytope,Gridap.ReferenceFEs.lagrangian,Float64,1)
cell_types=collect(Fill(1,length(cell_vertex_lids)))
cell_reffes=[scalar_reffe]
grid = Gridap.Geometry.UnstructuredGrid(node_coordinates,
cell_vertex_lids,
cell_reffes,
cell_types,
Gridap.Geometry.NonOriented())
model=Gridap.Geometry.UnstructuredDiscreteModel(grid)

n2o_faces_map=Vector{Vector{Int64}}(undef,3)
n2o_faces_map[3]=[1,1]
n2o_cell_to_child_id=[2,3]

ref_rules = [Gridap.Adaptivity.RefinementRule(Gridap.ReferenceFEs.LagrangianRefFE(Float64,QUAD,1),2)]

glue=Gridap.Adaptivity.AdaptivityGlue(Gridap.Adaptivity.RefinementGlue(),
n2o_faces_map,
n2o_cell_to_child_id,
ref_rules)
Gridap.Adaptivity.AdaptedDiscreteModel(model,parent,glue)
end

coarse_model=setup_coarse_discrete_model()
fine_model=setup_adapted_discrete_model(coarse_model)

order=0

VH = FESpace(coarse_model,
ReferenceFE(raviart_thomas,Float64,order),
conformity=:Hdiv)

UH = TrialFESpace(VH)

Vh = FESpace(fine_model,
ReferenceFE(raviart_thomas,Float64,order),
conformity=:Hdiv)

Uh = TrialFESpace(Vh)

u(x) = VectorValue(x[1],x[2])
uh = interpolate(u,Uh)
uH = interpolate(u,UH)
uhH = interpolate(uh,UH)

@test get_free_dof_values(uhH)[2] 1.0
end

test_PR_1074()

end

0 comments on commit c7fe919

Please sign in to comment.