grid = np.vstack((grid, np.ones((1, grid.shape[1]))))
barycoords = np.dot(barytransform, grid)
barycoords = barycoords[:, (np.all(barycoords >= 0, axis=0))]
Повышение эффективности расчета барицентрических координат в python