Quadrature points for boundary

[MichaelLoibl]

Dear all,
I am working on a GridapEmbedded model. I compute certain boundary integrals which works fine. Now, I am interested to obtain the quadrature points along a boundary. I tried the code below which worked for domain integrals. CellQuadrature is apparently the wrong choice to obtain the quadrature rule for the boundary. What would be an appropriate function?

using Gridap
using Gridap.ReferenceFEs
using GridapEmbedded
using GridapEmbedded.CSG
using GridapEmbedded.LevelSetCutters

domain = (0,1,0,1)
n = 2
bgmodel = CartesianDiscreteModel(domain,(n,n))

reffe = ReferenceFE(lagrangian,Float64,1)
Ω_bg = Triangulation(bgmodel)
V_bg = FESpace(Ω_bg,reffe)
φh = interpolate(x->(x[1] - 1/2)^2 + (x[2] - 1/2)^2 - 1/25,V_bg)
geo = DiscreteGeometry(φh,bgmodel)

cutgeo = cut(bgmodel,geo)

Ω = Triangulation(cutgeo,PHYSICAL)
Γd = EmbeddedBoundary(cutgeo)

order = 1
degree = 2*order
dΩ = Measure(Ω,degree)
dΓd = Measure(Γd,degree)

Qₕ_domain = CellQuadrature(dΩ,degree)
Qₕ_bound = CellQuadrature(dΓd,degree)

[JordiManyer]

The measures have the CellQuadrature inside.

So both of these are unnecessary:

Qₕ_domain = CellQuadrature(dΩ,degree)
Qₕ_bound = CellQuadrature(dΓd,degree)

you can just retrieve the quadrature that the measures contain.