In [1]:
from __future__ import division
import numpy as np
import scipy.linalg as lin
import pylab as pl
from locore.identifiability import crit_l1_analysis
from locore.operators import finite_diff_1d
Dimension of the problem
In [2]:
n = 100
Denoising setting
In [3]:
A = np.eye(n)
First order dictionary (TV)
In [4]:
D = finite_diff_1d(n)
Generate a test vector
In [5]:
x = np.ones((n,1))
x[n//2:n//2+10] = 2
Compute IC
In [6]:
crit_l1_analysis(D, A, x)
Out[6]: