def run_exact_tst(y, M, Omega):
"""
Wrapper for TST algorithm (with default optimized params) for exact analysis recovery
"""
nsweep = 300
tol = 1e-5
return ABSexact.tst_recom(y, M, Omega, nsweep, tol)
def run_exact_tst(y, M, Omega):
"""
Wrapper for TST algorithm (with default optimized params) for exact analysis recovery
"""
nsweep = 300
tol = 1e-5
return ABSexact.tst_recom(y, M, Omega, nsweep, tol)
def run_exact_sl0(y, M, Omega):
"""
Wrapper for SL0 algorithm for exact analysis recovery
"""
sigma_min = 1e-12
sigma_decrease_factor = 0.5
mu_0 = 2
L = 20
return ABSexact.sl0(y, M, Omega, sigma_min, sigma_decrease_factor, mu_0, L)
def run_exact_sl0(y, M, Omega):
"""
Wrapper for SL0 algorithm for exact analysis recovery
"""
sigma_min = 1e-12
sigma_decrease_factor = 0.5
mu_0 = 2
L = 20
return ABSexact.sl0(y, M, Omega, sigma_min, sigma_decrease_factor, mu_0, L)
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)
def run_exact_ompeps(y, M, Omega):
"""
Wrapper for OMP algorithm for exact analysis recovery, with stopping criterion = epsilon
"""
return ABSexact.ompeps(y, M, Omega, 1e-9)
def run_exact_bp_cvxopt(y, M, Omega):
"""
Wrapper for BP algorithm for exact analysis recovery
Algorithm implementation is using cvxopt linear programming
"""
return ABSexact.bp_cvxopt(y, M, Omega)
def run_exact_ompeps(y, M, Omega):
"""
Wrapper for OMP algorithm for exact analysis recovery, with stopping criterion = epsilon
"""
return ABSexact.ompeps(y, M, Omega, 1e-9)
def run_exact_bp_cvxopt(y, M, Omega):
"""
Wrapper for BP algorithm for exact analysis recovery
Algorithm implementation is using cvxopt linear programming
"""
return ABSexact.bp_cvxopt(y, M, Omega)
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)
