def point_in_prism(tri1, tri2, pt):
bcs = prism_barycentric_coordinates(tri1, tri2, pt)<EOL>return np.logical_not(np.isclose(np.sum(bcs, axis=<NUM_LIT:0>), <NUM_LIT:0>))<EOL>
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.