docs: def point_in_prism(tri1, tri2, pt):
      '''
      point_in_prism(tri1, tri2, pt) yields True if the given point is inside the prism that stretches
        between triangle 1 and triangle 2. Will automatically thread over extended dimensions. If
        multiple triangles are given, then the vertices must be an earlier dimension than the
        coordinates; e.g., a 3 x 3 x n array will be assumed to organized such that element [0,1,k] is
        the y coordinate of the first vertex of the k'th triangle.
      '''
      bcs = prism_barycentric_coordinates(tri1, tri2, pt)
      return np.logical_not(np.isclose(np.sum(bcs, axis=0), 0))