In a UFL file, what's the best way to normalise a vector, i.e.,
h = VectorConstant("tetrahedron") h0, h1, h2 = h.split() normh = sqrt(h0**2 + h1**2 + h2**2) unit_h = h/normh
Is this the best way?