def run_exact_bp(y, M, Omega): """ Wrapper for BP algorithm for exact analysis recovery Algorithm implementation is l1eq_pd() from l1-magic toolbox """ return ABSexact.bp(y, M, Omega, numpy.zeros(Omega.shape[0]), pdtol=1e-5, pdmaxiter=100)