Piecewise constants on the bounday of 2-D domain

Main Forum Old-Forum - do not post here
1

So I can create a piecewise linear FE space on a boundary of a 2-D square. But there does not seem to be a way to create piecewise constants. Is this just not doable?

2

for that we have the space of element-wise polynomials on surfaces:

fes = SurfaceL2(mesh, order=0)
3

Yes, but does this work on a sub-domain of co-dimension 1? Do I need to add an argument, like “definedon=”?

4

For this, you can use facetfespace, and use definedon=“”

5

Thanks. Actually, it needs to be:

bdy_L2 = FacetFESpace(mesh, order=0, definedon=“”, definedonbound=“outer”)