gridap / gridap/GridapEmbedded.jl
Periodic background mesh
- Dominant language
- Julia
- Stars
- 50
- Forks
- 16
- PR merge metrics
- No merged PRs in 30d
Description
Hi,
Is there support for periodic background meshes? The MWE below does not work.
```
using Gridap
using GridapEmbedded,GridapEmbedded.LevelSetCutters
import Gridap.Geometry: get_node_coordinates,collect1d
model = CartesianDiscreteModel((0,1,0,1),(100,100),isperiodic=(true,true));
V_φ = TestFESpace(model,ReferenceFE(lagrangian,Float64,1))
φh = interpolate(x->-cos(4π*x[1])*cos(4*pi*x[2])/4-0.2/4,V_φ)
point_to_coords = model.grid_topology.vertex_coordinates
geo = DiscreteGeometry(get_free_dof_values(φh),point_to_coords,name="")
cutgeo = cut(model,geo)
Ω = Triangulation(cutgeo,PHYSICAL)
writevtk(Ω,"discrete_geo_periodic",cellfields=["φh"=>φh])
```
Based on some testing I've conducted, it seems that analytic geometries do respect periodicity of the background mesh though.
Contributor guide
No contributing guide indexed for this repository
Assessment
This issue has not been assessed yet.