def tetra_vol(pts):
    '''Volume of a tetrahedron'''
    tetra = pts[1:] - pts[0]
    return np.abs(np.dot(tetra[0], np.cross(tetra[1], tetra[2]))) / 6
dataset/ETHPy150Open gallantlab/pycortex/cortex/polyutils.py/tetra_vol/original