def point_in_tetrahedron(tetra, pt):
bcs = tetrahedral_barycentric_coordinates(tetra, pt)<EOL>return np.logical_not(np.all(np.isclose(bcs, <NUM_LIT:0>), axis=<NUM_LIT:0>))<EOL>
point_in_tetrahedron(tetrahedron, point) yields True if the given point is in the given 
  tetrahedron. If either tetrahedron or point (or both) are lists of shapes/points, then this
  calculation is automatically threaded over all the given arguments.