@verbose
def _points_outside_surface(rr, surf, n_jobs=1, verbose=None):
    rr = np.atleast_2d(rr)
    assert rr.shape[1] == 3
    parallel, p_fun, _ = parallel_func(_get_solids, n_jobs)
    tot_angles = parallel(p_fun(surf['rr'][tris], rr)
                          for tris in np.array_split(surf['tris'], n_jobs))
    return np.abs(np.sum(tot_angles, axis=0) / (2 * np.pi) - 1.0) > 1e-5
Check whether points are outside a surface
dataset/ETHPy150Open mne-tools/mne-python/mne/source_space.py/_points_outside_surface