fenics

Please describe the problem you are trying to solve.
The current implementation of the barycentric coordinates are overly complicated and explicitly require the normals of triangles or edges. There are simpler ways to implement this.

Describe the solution you propose.

• Simplify math/barycenter.[ch]pp
• Remove Edge::computeNormal() as it is unnecessary.

**Describe alternativ

The actual ghostUpdate is truly needed only in very few places, for most local operations one knows the value for ghosts without communication.

This

u.vector.copy(result=u0.vector)
u0.vector.ghostUpdate(addv=PETSc.InsertMode.INSERT, mode=PETSc.ScatterMode.FORWARD)

could become

with u0.vector.localForm() as loc0, u.vector.localForm() as loc:
    loc.copy(loc0)
``