def test1(): import numpy as np import pylab from scipy import sparse from regreg.algorithms import FISTA from regreg.atoms import l1norm from regreg.container import container from regreg.smooth import quadratic Y = np.random.standard_normal(500); Y[100:150] += 7; Y[250:300] += 14 sparsity = l1norm(500, lagrange=1.0) #Create D D = (np.identity(500) + np.diag([-1]*499,k=1))[:-1] D = sparse.csr_matrix(D) fused = l1norm.linear(D, lagrange=19.5) loss = quadratic.shift(-Y, lagrange=0.5) p = container(loss, sparsity, fused) soln1 = blockwise([sparsity, fused], Y) solver = FISTA(p) solver.fit(max_its=800,tol=1e-10) soln2 = solver.composite.coefs #plot solution pylab.figure(num=1) pylab.clf() pylab.scatter(np.arange(Y.shape[0]), Y, c='r') pylab.plot(soln1, c='y', linewidth=6) pylab.plot(soln2, c='b', linewidth=2)
import numpy as np import pylab from scipy import sparse from regreg.algorithms import FISTA from regreg.atoms import l1norm from regreg.container import container from regreg.smooth import l2normsq Y = np.random.standard_normal(500); Y[100:150] += 7; Y[250:300] += 14 loss = l2normsq.shift(-Y, coef=0.5) sparsity = l1norm(len(Y), 1.4) # TODO should make a module to compute typical Ds D = sparse.csr_matrix((np.identity(500) + np.diag([-1]*499,k=1))[:-1]) fused = l1norm.linear(D, 25.5) problem = container(loss, sparsity, fused) solver = FISTA(problem.composite()) solver.fit(max_its=100, tol=1e-10) solution = solver.composite.coefs delta1 = np.fabs(D * solution).sum() delta2 = np.fabs(solution).sum() fused_constraint = l1norm.linear(D, bound=delta1) sparsity_constraint = l1norm(500, bound=delta2) constrained_problem = container(loss, fused_constraint, sparsity_constraint) constrained_solver = FISTA(constrained_problem.composite()) constrained_solver.composite.lipshitz = 1.01
from scipy import sparse from regreg.algorithms import FISTA from regreg.atoms import l1norm from regreg.container import container from regreg.smooth import l2normsq Y = np.random.standard_normal(500) Y[100:150] += 7 Y[250:300] += 14 loss = l2normsq.shift(-Y, coef=0.5) sparsity = l1norm(len(Y), 1.4) # TODO should make a module to compute typical Ds D = sparse.csr_matrix((np.identity(500) + np.diag([-1] * 499, k=1))[:-1]) fused = l1norm.linear(D, 25.5) problem = container(loss, sparsity, fused) solver = FISTA(problem.composite()) solver.fit(max_its=100, tol=1e-10) solution = solver.composite.coefs delta1 = np.fabs(D * solution).sum() delta2 = np.fabs(solution).sum() fused_constraint = l1norm.linear(D, bound=delta1) sparsity_constraint = l1norm(500, bound=delta2) constrained_problem = container(loss, fused_constraint, sparsity_constraint) constrained_solver = FISTA(constrained_problem.composite()) constrained_solver.composite.lipshitz = 1.01
import numpy as np import pylab from scipy import sparse from regreg.algorithms import FISTA from regreg.atoms import l1norm from regreg.container import container from regreg.smooth import signal_approximator Y = np.random.standard_normal(500); Y[100:150] += 7; Y[250:300] += 14 sparsity = l1norm(500, l=1.3) #Create D D = (np.identity(500) + np.diag([-1]*499,k=1))[:-1] D = sparse.csr_matrix(D) fused = l1norm.linear(D, l=25.5) loss = signal_approximator(Y) problem = container(loss, sparsity, fused) solver = FISTA(problem.problem()) solver.fit(max_its=800,tol=1e-10) soln = solver.problem.coefs #plot solution pylab.figure(num=1) pylab.clf() pylab.plot(soln, c='g') pylab.scatter(np.arange(Y.shape[0]), Y)
import pylab from scipy import sparse from regreg.algorithms import FISTA from regreg.atoms import l1norm from regreg.container import container from regreg.smooth import signal_approximator Y = np.random.standard_normal(500) Y[100:150] += 7 Y[250:300] += 14 sparsity = l1norm(500, l=1.3) #Create D D = (np.identity(500) + np.diag([-1] * 499, k=1))[:-1] D = sparse.csr_matrix(D) fused = l1norm.linear(D, l=25.5) loss = signal_approximator(Y) problem = container(loss, sparsity, fused) solver = FISTA(problem.problem()) solver.fit(max_its=800, tol=1e-10) soln = solver.problem.coefs #plot solution pylab.figure(num=1) pylab.clf() pylab.plot(soln, c='g') pylab.scatter(np.arange(Y.shape[0]), Y)