Merge branch 'moment-based-reffes' of github.com:Antoinemarteau/Tutorials into v0.20

JordiManyer · JordiManyer · commit 9692b46bdd9b · 2026-03-24T18:57:28.000+11:00
diff --git a/Project.toml b/Project.toml
@@ -40,6 +40,7 @@ BlockArrays = "1.6.3"
 DataStructures = "0.18.22"
 Gridap = "0.20"
 GridapDistributed = "0.4"
+GridapEmbedded = "0.9.9"
 GridapGmsh = "0.7"
 GridapP4est = "0.3"
 GridapPETSc = "0.5"
diff --git a/src/darcy.jl b/src/darcy.jl
@@ -63,14 +63,18 @@ model = CartesianDiscreteModel(domain,partition)
 # Next, we build the FE spaces. We consider the first order RT space for the flux and the discontinuous pressure space as described above.  This mixed FE pair satisfies the inf-sup condition and, thus, it is stable.
 
 order = 1
+h = 1/100
 
 V = FESpace(model, ReferenceFE(raviart_thomas,Float64,order),
-      conformity=:HDiv, dirichlet_tags=[5,6])
+      conformity=:HDiv, dirichlet_tags=[5,6], scale_dof=true, global_meshsize=h)
 
 Q = FESpace(model, ReferenceFE(lagrangian,Float64,order),
-      conformity=:L2)
+      conformity=:L2, scale_dof=true, global_meshsize=h)
 
-# Note that the Dirichlet boundary for the flux are the bottom and top sides of the squared domain (identified with the boundary tags 5, and 6 respectively), whereas no Dirichlet data can be imposed on the pressure space. We select `conformity=:HDiv` for the flux (i.e., shape functions with $H^1(\mathrm{div};\Omega)$ regularity) and `conformity=:L2` for the pressure (i.e. discontinuous shape functions).
+# We select `conformity=:HDiv` for the flux (i.e., shape functions with $H^1(\mathrm{div};\Omega)$ regularity) and `conformity=:L2` for the pressure (i.e. discontinuous shape functions).
+# Note that the Dirichlet boundary for the flux are the bottom and top sides of the squared domain (identified with the boundary tags 5, and 6 respectively), whereas no Dirichlet data can be imposed on the pressure space.
+#
+# The `scale_dof=true` argument is given to cancel the scaling of the DoFs with the meshsize `h` that is introduced by the piola map, e.g. `h`$^{2-1}$ (2 is the dimension) for the div-conforming contra-variant Piola map. This is usefull for mixed elements of different conformity that use small or non-uniform meshsize. On a uniform Cartesian mesh like `model` used here, we also passed `global_meshsize=1/100`. This argument fixes the local meshsize estimate needed to re-scale the DoFs, avoiding the computation of the volume of each face of the mesh that own a DoF. But this should only be used on quasi-uniform and shape-regular meshes.
 #
 # From these objects, we construct the trial spaces. Note that we impose homogeneous boundary conditions for the flux.
 
diff --git a/src/inc_navier_stokes.jl b/src/inc_navier_stokes.jl
@@ -105,6 +105,8 @@ P = TrialFESpace(Q)
 Y = MultiFieldFESpace([V, Q])
 X = MultiFieldFESpace([U, P])
 
+# `scale_dof` is not used here, because currently $H^1$ and $L^2$ conforming elements both use the same Piola map in Gridap (the trivial one).
+
 # ## Triangulation and integration quadrature
 #
 # From the discrete model we can define the triangulation and integration measure

Original file line number	Diff line number	Diff line change
`@@ -105,6 +105,8 @@ P = TrialFESpace(Q)`
`105`	`105`	`Y = MultiFieldFESpace([V, Q])`
`106`	`106`	`X = MultiFieldFESpace([U, P])`
`107`	`107`
	`108`	+# `scale_dof` is not used here, because currently $H^1$ and $L^2$ conforming elements both use the same Piola map in Gridap (the trivial one).
	`109`	`+`
`108`	`110`	`# ## Triangulation and integration quadrature`
`109`	`111`	`#`
`110`	`112`	`# From the discrete model we can define the triangulation and integration measure`