Mesh connectivity in Fenicsx

dolfinx
5

I’ve made an issue regarding this at:

A workaround would be to use:

Ve = fem.FunctionSpace(domain, ("DG", 0))
q_degree = 3
vol2 = fem.Function(Ve)
fem.petsc.assemble_vector(vol2.vector, fem.form(
    1*ufl.TestFunction(Ve)*ufl.dx(metadata={"quadrature_degree": q_degree})))
vol2.x.scatter_forward()
print(vol2.x.array - elmVol.x.array)
assert np.allclose(vol2.x.array, elmVol.x.array)
vol2.name = "ElementVolume (assembly)"
1 Like