mymesh.quality.tet_circumradius

mymesh.quality.tet_circumradius(NodeCoords, NodeConn, V=None)

Circumradii for elements in a tetrahedral mesh.

\[V = \frac{\sqrt{(aA + bB + cC)(aA + bB - cC)(aA - bB + cC)(-aA + bB + cC)}}{24 V}\]

where \(a\), \(b\), and \(c\), are the lengths of the three edges meeting at a vertex and \(A\), \(B\), and \(C\) are the lengths of their opposite edges.

Parameters:

• NodeCoords (np.ndarray) – Node coordinates (shape=(n,3))

• NodeConn (np.ndarray) – Node connectivity (shape=(m,4), dtype=int)

• V (array_like, optional) – Volume of tetrahedra. If not provided, the volumes will be calculated.

Returns:

R – Circumradii of tetrahedra

Return type:

np.ndarray