def METHOD_NAME(p1, p2, p3, p4):  # pragma: no cover
    """doesn't compute the elemental for a tet element"""
    ## R vector joins midpoints of edges G1-G2 and G3-G4.
    p12 = (p1 + p2) / 2.
    p34 = (p3 + p4) / 2.
    r = p34 - p12
    ## S vector joins midpoints of edges G1-G3 and G2-G4.
    p13 = (p1 + p3) / 2.
    p24 = (p2 + p4) / 2.
    s = p24 - p13
    ## T vector joins midpoints of edges G1-G4 and G2-G3.
    p14 = (p1 + p4) / 2.
    p23 = (p2 + p3) / 2.
    t = p23 - p14
    origin = p1
    # The element coordinate system is chosen as close as possible to the
    # R, S, and T vectors and points in the same general direction.
    #
    # Mathematically speaking, the coordinate system is computed in such
    # a way that, if the R, S, and T vectors are described in the element
    # coordinate system, a 3x3 positive definite symmetric matrix would
    # be produced.
    #
    # ???
    return r, s, t, origin