Issue with vector-valued FacetSpace

For the latest nightly build (self-compiled)
gf.Set((x*y,y), dual=True)
works, whereas
gf.Set(x*y, dual=True)
does not work as expected.

You can try using
Bdy_Space = VectorFacetFESpace(Om_mesh, order=0, definedon="", definedonbound="Gamma")
instead of
Bdy_Space = FacetFESpace(Om_mesh, order=0, dim=2, definedon="", definedonbound="Gamma")

VectorFacetFESpace has automatically the dimension of the mesh, i.e. 2 in your example.

Best,
Michael

Thank you for this.  I think the reason it did not work is because I didn't have the latest build.  The other forum thread told me what to do.

Note: I need to use the "dim" argument, because in 3-D, I need to set dim=5 (not the spatial dimension).