RR_IMAGE_SOURCE_SEPARATION - Comparison of different methods
Description
Reproducible research addendum for compressive source separation
THEORY AND METHODS FOR HYPERSPECTRAL IMAGING
Paper: Mohammad Golbabaee, Simon Arberet, and Pierre Vandergheynst
Demonstration matlab file: Perraudin Nathanael, Mohammad Golbabaee
EPFL -- August 2012
Dependencies
In order to use this matlab file you need the UNLocXbox toolbox. You can download it on https://lts2research.epfl.ch/unlocbox/
The problem
We would like to solve this problem:
\begin{equation*}
\begin{split} \operatorname{arg\,min}_S & \sum_j \|S_p(.,j)\|_{TV} \\\\ \text{ such that: } & \| \Phi ( S \cdot H^T ) - y \|_F < \epsilon \text{ Projection on a B2-Ball } \\\\ \text{and } & S(i,j) > 0 \text{ for all } i,j \text{ (positivity constraint) } \\\\ \text{and } & \sum_j S(i,j) = 1 \text{ for all } i \end{split}
\end{equation*}
- with
- \(S\) : Sources
- \(Y\) = \(\Phi(I_m) + Z\): the measurement (forward model)
- \(Z\) : noise
- \(I_m\): Original 3D image (64 x 64 x 64) in 'Data.mat'
- \(H\) : Mixing matrix
Uniform sampling
For UNIFORM sampling (Block_diag), a decorrelation step can be applied by projecting the measurements onto pinv(H), and solve a easier problem:
\begin{equation*}
\begin{split} \operatorname{arg\,min}_S & \sum_j \|S_p(.,j)\|_{TV} \\\\ \text{ such that: } & \| \Phi S_p - Y_p \|_F < \epsilon \text{ Projection on a B2-Ball } \\\\ \text{and } & S_p(i,j) > 0 \text{ for all } i,j \text{ (positivity constraint) } \\\\ \text{and } & \sum_j S_p(i,j) = 1 \text{ for all } i \\\\ \text{with } & Y_p = Y H^{+,T} (\text{ where } H^+ \text{ denotes the pseudo inverse of }H) \end{split}
\end{equation*}
- with
- \(Y\) = \(\Phi I_m + Z\): the measurement (forward model)
This step accelerates the computation and increase accuracy of source recovery. However, the use of the block diagonale (Block_diag) sampling is required!
Techniques for solving the problem
We will compare 4 different techniques to solve the problem
- Block_diag operator with TV regularization (Decorrelated measurements)
- Block_diag operator with Wavelets regularization (Decorrelated measurements)
- Dense operator with TV regularization (Correlated measurements)
- Dense operator with Wavelets regularization (Correlated measurements)
Results
Sprectral signature of the different sources
Orignial sources
Block_diag operator with TV regularization (Uncorrelated measurements)
Block_diag operator with Wavelets regularization (Uncorrelated measurements)
Dense operator with TV regularization (Correlated measurements)
Dense operator with Wavelets regularization (Correlated measurements)
This code produces the following output:
The time step is set manually to : 1
Algorithm selected: PPXA
Iteration 1:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 8.777459e-15, TOL_EPS, iter = 0
Prox_TV: obj = 0.000000e+00, rel_obj = NaN, MAX_IT, iter = 200
f(x^*) = 1.625159e+03, rel_eval = 1.000000e+00
Iteration 2:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.072550e-14, TOL_EPS, iter = 0
Prox_TV: obj = 1.969811e+02, rel_obj = 4.438895e-05, TOL_EPS, iter = 25
f(x^*) = 2.273401e+03, rel_eval = 2.851417e-01
Iteration 3:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.276980e-14, TOL_EPS, iter = 0
Prox_TV: obj = 5.878010e+02, rel_obj = 9.992657e-05, TOL_EPS, iter = 51
f(x^*) = 2.357874e+03, rel_eval = 3.582617e-02
Iteration 4:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.500319e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.956574e+02, rel_obj = 6.550449e-05, TOL_EPS, iter = 48
f(x^*) = 2.095839e+03, rel_eval = 1.250268e-01
Iteration 5:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.576912e-14, TOL_EPS, iter = 0
Prox_TV: obj = 1.271434e+03, rel_obj = 9.638162e-05, TOL_EPS, iter = 35
f(x^*) = 1.725926e+03, rel_eval = 2.143268e-01
Iteration 6:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.554978e-14, TOL_EPS, iter = 0
Prox_TV: obj = 1.389143e+03, rel_obj = 9.206145e-05, TOL_EPS, iter = 41
f(x^*) = 1.388456e+03, rel_eval = 2.430545e-01
Iteration 7:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.690903e-14, TOL_EPS, iter = 0
Prox_TV: obj = 1.390194e+03, rel_obj = 8.058345e-05, TOL_EPS, iter = 33
f(x^*) = 1.202921e+03, rel_eval = 1.542374e-01
Iteration 8:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.632558e-14, TOL_EPS, iter = 0
Prox_TV: obj = 1.312058e+03, rel_obj = 8.870730e-05, TOL_EPS, iter = 38
f(x^*) = 1.122423e+03, rel_eval = 7.171806e-02
Iteration 9:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.677427e-14, TOL_EPS, iter = 0
Prox_TV: obj = 1.204985e+03, rel_obj = 9.659891e-05, TOL_EPS, iter = 45
f(x^*) = 1.013007e+03, rel_eval = 1.080106e-01
Iteration 10:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.629780e-14, TOL_EPS, iter = 0
Prox_TV: obj = 1.107146e+03, rel_obj = 9.114429e-05, TOL_EPS, iter = 42
f(x^*) = 8.692454e+02, rel_eval = 1.653868e-01
Iteration 11:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.768578e-14, TOL_EPS, iter = 0
Prox_TV: obj = 1.032324e+03, rel_obj = 9.515797e-05, TOL_EPS, iter = 47
f(x^*) = 7.447572e+02, rel_eval = 1.671527e-01
Iteration 12:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.659520e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.856836e+02, rel_obj = 9.031932e-05, TOL_EPS, iter = 44
f(x^*) = 6.663420e+02, rel_eval = 1.176800e-01
Iteration 13:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.679560e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.615243e+02, rel_obj = 7.031575e-05, TOL_EPS, iter = 43
f(x^*) = 6.296317e+02, rel_eval = 5.830454e-02
Iteration 14:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.710756e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.525521e+02, rel_obj = 6.550090e-05, TOL_EPS, iter = 43
f(x^*) = 6.023504e+02, rel_eval = 4.529137e-02
Iteration 15:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.648123e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.525372e+02, rel_obj = 7.248738e-05, TOL_EPS, iter = 45
f(x^*) = 5.714816e+02, rel_eval = 5.401544e-02
Iteration 16:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.679333e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.570628e+02, rel_obj = 7.541219e-05, TOL_EPS, iter = 47
f(x^*) = 5.470100e+02, rel_eval = 4.473693e-02
Iteration 17:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.734636e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.629024e+02, rel_obj = 8.930399e-05, TOL_EPS, iter = 45
f(x^*) = 5.352792e+02, rel_eval = 2.191523e-02
Iteration 18:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.648757e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.670465e+02, rel_obj = 9.360061e-05, TOL_EPS, iter = 44
f(x^*) = 5.299098e+02, rel_eval = 1.013274e-02
Iteration 19:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.741253e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.684455e+02, rel_obj = 9.678243e-05, TOL_EPS, iter = 44
f(x^*) = 5.257404e+02, rel_eval = 7.930483e-03
Iteration 20:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.720721e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.675464e+02, rel_obj = 8.069933e-05, TOL_EPS, iter = 45
f(x^*) = 5.208072e+02, rel_eval = 9.472179e-03
Iteration 21:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.687121e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.657888e+02, rel_obj = 8.568492e-05, TOL_EPS, iter = 45
f(x^*) = 5.156501e+02, rel_eval = 1.000123e-02
Iteration 22:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.815002e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.643746e+02, rel_obj = 9.753904e-05, TOL_EPS, iter = 43
f(x^*) = 5.126174e+02, rel_eval = 5.916223e-03
Iteration 23:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.787413e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.641395e+02, rel_obj = 9.953912e-05, TOL_EPS, iter = 38
f(x^*) = 5.114311e+02, rel_eval = 2.319446e-03
Iteration 24:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.702688e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.633361e+02, rel_obj = 9.289655e-05, TOL_EPS, iter = 44
f(x^*) = 5.094405e+02, rel_eval = 3.907481e-03
Iteration 25:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.742122e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.639588e+02, rel_obj = 9.502711e-05, TOL_EPS, iter = 43
f(x^*) = 5.069093e+02, rel_eval = 4.993436e-03
Iteration 26:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.660721e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.655104e+02, rel_obj = 8.208210e-05, TOL_EPS, iter = 36
f(x^*) = 5.061889e+02, rel_eval = 1.423113e-03
Iteration 27:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.643814e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.658979e+02, rel_obj = 8.510157e-05, TOL_EPS, iter = 37
f(x^*) = 5.071753e+02, rel_eval = 1.944917e-03
Iteration 28:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.685600e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.661743e+02, rel_obj = 8.594692e-05, TOL_EPS, iter = 36
f(x^*) = 5.076564e+02, rel_eval = 9.477541e-04
Iteration 29:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.655590e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.660884e+02, rel_obj = 8.051686e-05, TOL_EPS, iter = 36
f(x^*) = 5.073344e+02, rel_eval = 6.348847e-04
Iteration 30:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.688142e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.658526e+02, rel_obj = 8.530702e-05, TOL_EPS, iter = 36
f(x^*) = 5.068857e+02, rel_eval = 8.850302e-04
Iteration 31:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.688433e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.641781e+02, rel_obj = 9.850714e-05, TOL_EPS, iter = 48
f(x^*) = 5.046407e+02, rel_eval = 4.448792e-03
Iteration 32:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.682173e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.639540e+02, rel_obj = 9.481559e-05, TOL_EPS, iter = 49
f(x^*) = 5.050078e+02, rel_eval = 7.268617e-04
Iteration 33:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.677029e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.654052e+02, rel_obj = 8.764792e-05, TOL_EPS, iter = 37
f(x^*) = 5.054268e+02, rel_eval = 8.289958e-04
Iteration 34:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.764682e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.656613e+02, rel_obj = 9.829042e-05, TOL_EPS, iter = 35
f(x^*) = 5.067527e+02, rel_eval = 2.616447e-03
Iteration 35:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.710003e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.656893e+02, rel_obj = 8.724792e-05, TOL_EPS, iter = 35
f(x^*) = 5.084887e+02, rel_eval = 3.414131e-03
Iteration 36:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.695152e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.657090e+02, rel_obj = 6.321279e-05, TOL_EPS, iter = 35
f(x^*) = 5.082035e+02, rel_eval = 5.611279e-04
Iteration 37:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.682356e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.656845e+02, rel_obj = 4.657418e-05, TOL_EPS, iter = 35
f(x^*) = 5.074231e+02, rel_eval = 1.538088e-03
Iteration 38:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.817315e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.656752e+02, rel_obj = 8.786710e-05, TOL_EPS, iter = 34
f(x^*) = 5.071711e+02, rel_eval = 4.967725e-04
Iteration 39:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.704567e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.655744e+02, rel_obj = 7.610909e-05, TOL_EPS, iter = 34
f(x^*) = 5.074016e+02, rel_eval = 4.543119e-04
Iteration 40:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.700117e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.654791e+02, rel_obj = 8.254583e-05, TOL_EPS, iter = 34
f(x^*) = 5.076152e+02, rel_eval = 4.207438e-04
Iteration 41:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.799048e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.653398e+02, rel_obj = 4.819717e-05, TOL_EPS, iter = 35
f(x^*) = 5.074052e+02, rel_eval = 4.139136e-04
Iteration 42:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.658532e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.653501e+02, rel_obj = 8.493400e-05, TOL_EPS, iter = 34
f(x^*) = 5.075034e+02, rel_eval = 1.935441e-04
Iteration 43:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.704545e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.653619e+02, rel_obj = 2.953633e-05, TOL_EPS, iter = 34
f(x^*) = 5.076786e+02, rel_eval = 3.450883e-04
Iteration 44:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.750627e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.653501e+02, rel_obj = 1.217152e-05, TOL_EPS, iter = 34
f(x^*) = 5.078605e+02, rel_eval = 3.581033e-04
Iteration 45:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.740216e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.653457e+02, rel_obj = 1.273960e-05, TOL_EPS, iter = 34
f(x^*) = 5.079543e+02, rel_eval = 1.846399e-04
Iteration 46:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.684522e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.654428e+02, rel_obj = 7.171570e-05, TOL_EPS, iter = 32
f(x^*) = 5.081546e+02, rel_eval = 3.941600e-04
Iteration 47:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.703563e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.653133e+02, rel_obj = 7.443116e-05, TOL_EPS, iter = 32
f(x^*) = 5.089314e+02, rel_eval = 1.526469e-03
Iteration 48:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.704651e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.650967e+02, rel_obj = 9.644243e-05, TOL_EPS, iter = 32
f(x^*) = 5.091420e+02, rel_eval = 4.136444e-04
Iteration 49:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.661744e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.649248e+02, rel_obj = 7.927922e-05, TOL_EPS, iter = 32
f(x^*) = 5.088920e+02, rel_eval = 4.912539e-04
Iteration 50:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.755470e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.648117e+02, rel_obj = 4.175398e-05, TOL_EPS, iter = 32
f(x^*) = 5.087945e+02, rel_eval = 1.916959e-04
Iteration 51:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.786111e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.646979e+02, rel_obj = 5.453823e-05, TOL_EPS, iter = 32
f(x^*) = 5.088493e+02, rel_eval = 1.076286e-04
Iteration 52:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.625270e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.646368e+02, rel_obj = 1.686300e-05, TOL_EPS, iter = 32
f(x^*) = 5.089939e+02, rel_eval = 2.840478e-04
Iteration 53:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.726339e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.645522e+02, rel_obj = 1.061202e-05, TOL_EPS, iter = 32
f(x^*) = 5.091439e+02, rel_eval = 2.946551e-04
Iteration 54:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.688572e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.644617e+02, rel_obj = 5.330352e-08, TOL_EPS, iter = 32
f(x^*) = 5.092535e+02, rel_eval = 2.152514e-04
Iteration 55:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.616389e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.643994e+02, rel_obj = 3.862874e-05, TOL_EPS, iter = 32
f(x^*) = 5.093549e+02, rel_eval = 1.990276e-04
Iteration 56:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.666791e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.643170e+02, rel_obj = 8.612587e-05, TOL_EPS, iter = 32
f(x^*) = 5.094653e+02, rel_eval = 2.167834e-04
Iteration 57:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.787552e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.642108e+02, rel_obj = 8.904680e-05, TOL_EPS, iter = 32
f(x^*) = 5.095209e+02, rel_eval = 1.090931e-04
Iteration 58:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.665885e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.640575e+02, rel_obj = 5.322693e-05, TOL_EPS, iter = 33
f(x^*) = 5.092216e+02, rel_eval = 5.876748e-04
Iteration 59:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.674536e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.640306e+02, rel_obj = 6.351518e-05, TOL_EPS, iter = 33
f(x^*) = 5.093938e+02, rel_eval = 3.379741e-04
Iteration 60:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.632609e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.640613e+02, rel_obj = 8.854250e-05, TOL_EPS, iter = 32
f(x^*) = 5.095625e+02, rel_eval = 3.310836e-04
Iteration 61:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.709584e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.639510e+02, rel_obj = 7.499677e-05, TOL_EPS, iter = 31
f(x^*) = 5.098260e+02, rel_eval = 5.168247e-04
Iteration 62:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.740104e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.638286e+02, rel_obj = 6.082638e-05, TOL_EPS, iter = 31
f(x^*) = 5.107226e+02, rel_eval = 1.755589e-03
Iteration 63:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.682239e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.636950e+02, rel_obj = 9.681383e-05, TOL_EPS, iter = 30
f(x^*) = 5.111641e+02, rel_eval = 8.636056e-04
Iteration 64:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.703079e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.631645e+02, rel_obj = 2.148582e-05, TOL_EPS, iter = 32
f(x^*) = 5.100070e+02, rel_eval = 2.268669e-03
Iteration 65:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.660419e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.629480e+02, rel_obj = 5.857607e-05, TOL_EPS, iter = 33
f(x^*) = 5.095512e+02, rel_eval = 8.946598e-04
Iteration 66:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.731754e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.628888e+02, rel_obj = 8.012750e-05, TOL_EPS, iter = 31
f(x^*) = 5.099085e+02, rel_eval = 7.008375e-04
Iteration 67:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.669935e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.628256e+02, rel_obj = 1.047561e-06, TOL_EPS, iter = 31
f(x^*) = 5.103936e+02, rel_eval = 9.504907e-04
Iteration 68:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.666067e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.627711e+02, rel_obj = 9.485841e-05, TOL_EPS, iter = 29
f(x^*) = 5.116249e+02, rel_eval = 2.406492e-03
Iteration 69:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.682854e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.622954e+02, rel_obj = 1.676968e-05, TOL_EPS, iter = 29
f(x^*) = 5.132937e+02, rel_eval = 3.251225e-03
Iteration 70:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.674792e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.616197e+02, rel_obj = 9.576158e-06, TOL_EPS, iter = 29
f(x^*) = 5.133986e+02, rel_eval = 2.043115e-04
Iteration 71:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.656947e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.608926e+02, rel_obj = 5.653960e-05, TOL_EPS, iter = 29
f(x^*) = 5.125206e+02, rel_eval = 1.713102e-03
Iteration 72:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.664374e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.602681e+02, rel_obj = 5.675119e-05, TOL_EPS, iter = 29
f(x^*) = 5.119761e+02, rel_eval = 1.063593e-03
Iteration 73:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.726224e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.597359e+02, rel_obj = 1.408013e-05, TOL_EPS, iter = 29
f(x^*) = 5.119166e+02, rel_eval = 1.162413e-04
Iteration 74:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.690016e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.592855e+02, rel_obj = 7.522378e-05, TOL_EPS, iter = 29
f(x^*) = 5.120614e+02, rel_eval = 2.828882e-04
Iteration 75:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.682644e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.589485e+02, rel_obj = 4.493032e-05, TOL_EPS, iter = 30
f(x^*) = 5.117427e+02, rel_eval = 6.226895e-04
Iteration 76:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.708431e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.586737e+02, rel_obj = 5.750271e-05, TOL_EPS, iter = 30
f(x^*) = 5.120414e+02, rel_eval = 5.833020e-04
Iteration 77:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.682970e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.582232e+02, rel_obj = 8.903049e-05, TOL_EPS, iter = 28
f(x^*) = 5.119438e+02, rel_eval = 1.906011e-04
Iteration 78:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.714459e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.578049e+02, rel_obj = 8.740493e-05, TOL_EPS, iter = 28
f(x^*) = 5.134614e+02, rel_eval = 2.955484e-03
Iteration 79:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.657831e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.571819e+02, rel_obj = 9.717739e-05, TOL_EPS, iter = 30
f(x^*) = 5.121690e+02, rel_eval = 2.523379e-03
Iteration 80:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.670293e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.568603e+02, rel_obj = 4.651045e-05, TOL_EPS, iter = 29
f(x^*) = 5.123649e+02, rel_eval = 3.823949e-04
Iteration 81:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.608134e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.565567e+02, rel_obj = 9.282292e-05, TOL_EPS, iter = 29
f(x^*) = 5.130037e+02, rel_eval = 1.245114e-03
Iteration 82:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.768741e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.561373e+02, rel_obj = 5.665340e-05, TOL_EPS, iter = 30
f(x^*) = 5.125562e+02, rel_eval = 8.730557e-04
Iteration 83:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.692308e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.559080e+02, rel_obj = 8.761303e-05, TOL_EPS, iter = 29
f(x^*) = 5.127862e+02, rel_eval = 4.486933e-04
Iteration 84:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.710148e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.554201e+02, rel_obj = 4.331473e-05, TOL_EPS, iter = 28
f(x^*) = 5.130666e+02, rel_eval = 5.463722e-04
Iteration 85:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.685404e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.486138e+02, rel_obj = 9.115555e-05, TOL_EPS, iter = 51
f(x^*) = 5.073338e+02, rel_eval = 1.129985e-02
Iteration 86:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.645391e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.545003e+02, rel_obj = 9.321397e-05, TOL_EPS, iter = 28
f(x^*) = 5.098026e+02, rel_eval = 4.842774e-03
Iteration 87:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.662888e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.541544e+02, rel_obj = 5.923399e-05, TOL_EPS, iter = 27
f(x^*) = 5.132396e+02, rel_eval = 6.696575e-03
Iteration 88:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.664247e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.532813e+02, rel_obj = 1.809506e-05, TOL_EPS, iter = 27
f(x^*) = 5.155573e+02, rel_eval = 4.495491e-03
Iteration 89:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.686404e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.522459e+02, rel_obj = 6.604846e-05, TOL_EPS, iter = 27
f(x^*) = 5.153542e+02, rel_eval = 3.939369e-04
Iteration 90:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.696841e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.514884e+02, rel_obj = 7.597329e-05, TOL_EPS, iter = 27
f(x^*) = 5.144504e+02, rel_eval = 1.756834e-03
Iteration 91:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.684266e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.443578e+02, rel_obj = 7.916613e-05, TOL_EPS, iter = 50
f(x^*) = 5.068072e+02, rel_eval = 1.508109e-02
Iteration 92:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.626499e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.439183e+02, rel_obj = 7.465175e-05, TOL_EPS, iter = 50
f(x^*) = 5.049207e+02, rel_eval = 3.736292e-03
Iteration 93:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.724381e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.505526e+02, rel_obj = 9.734100e-05, TOL_EPS, iter = 29
f(x^*) = 5.089589e+02, rel_eval = 7.934242e-03
Iteration 94:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.717038e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.506040e+02, rel_obj = 8.592485e-05, TOL_EPS, iter = 29
f(x^*) = 5.135867e+02, rel_eval = 9.010724e-03
Iteration 95:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.660466e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.504469e+02, rel_obj = 7.979441e-05, TOL_EPS, iter = 29
f(x^*) = 5.156362e+02, rel_eval = 3.974798e-03
Iteration 96:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.654211e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.501113e+02, rel_obj = 9.331138e-05, TOL_EPS, iter = 28
f(x^*) = 5.147059e+02, rel_eval = 1.807553e-03
Iteration 97:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.712684e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.429859e+02, rel_obj = 7.206344e-05, TOL_EPS, iter = 50
f(x^*) = 5.070321e+02, rel_eval = 1.513474e-02
Iteration 98:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.763534e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.427576e+02, rel_obj = 7.280402e-05, TOL_EPS, iter = 50
f(x^*) = 5.052304e+02, rel_eval = 3.566169e-03
Iteration 99:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.763595e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.428343e+02, rel_obj = 7.642108e-05, TOL_EPS, iter = 50
f(x^*) = 5.058642e+02, rel_eval = 1.253016e-03
Iteration 100:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.702635e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.428341e+02, rel_obj = 7.141435e-05, TOL_EPS, iter = 50
f(x^*) = 5.044328e+02, rel_eval = 2.837714e-03
Iteration 101:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.627313e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.495022e+02, rel_obj = 8.217591e-05, TOL_EPS, iter = 27
f(x^*) = 5.103841e+02, rel_eval = 1.166049e-02
Iteration 102:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.629177e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.423160e+02, rel_obj = 9.453390e-05, TOL_EPS, iter = 49
f(x^*) = 5.087814e+02, rel_eval = 3.150013e-03
Iteration 103:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.774158e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.417683e+02, rel_obj = 9.659827e-05, TOL_EPS, iter = 49
f(x^*) = 5.036407e+02, rel_eval = 1.020714e-02
Iteration 104:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.779286e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.486486e+02, rel_obj = 8.582642e-05, TOL_EPS, iter = 29
f(x^*) = 5.092642e+02, rel_eval = 1.104248e-02
Iteration 105:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.737590e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.417096e+02, rel_obj = 9.577651e-05, TOL_EPS, iter = 49
f(x^*) = 5.069787e+02, rel_eval = 4.508173e-03
Iteration 106:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.717261e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.414383e+02, rel_obj = 9.755700e-05, TOL_EPS, iter = 49
f(x^*) = 5.035135e+02, rel_eval = 6.882006e-03
Iteration 107:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.716658e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.481407e+02, rel_obj = 9.792177e-05, TOL_EPS, iter = 27
f(x^*) = 5.102316e+02, rel_eval = 1.316678e-02
Iteration 108:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.629473e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.483232e+02, rel_obj = 8.871759e-05, TOL_EPS, iter = 29
f(x^*) = 5.134716e+02, rel_eval = 6.309924e-03
Iteration 109:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.612028e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.477967e+02, rel_obj = 9.879167e-05, TOL_EPS, iter = 28
f(x^*) = 5.156777e+02, rel_eval = 4.278075e-03
Iteration 110:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.715577e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.402594e+02, rel_obj = 9.690510e-05, TOL_EPS, iter = 49
f(x^*) = 5.082462e+02, rel_eval = 1.462194e-02
Iteration 111:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.684734e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.399915e+02, rel_obj = 9.645877e-05, TOL_EPS, iter = 49
f(x^*) = 5.048037e+02, rel_eval = 6.819298e-03
Iteration 112:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.603687e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.400797e+02, rel_obj = 9.351505e-05, TOL_EPS, iter = 49
f(x^*) = 5.062025e+02, rel_eval = 2.763197e-03
Iteration 113:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.692092e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.400938e+02, rel_obj = 9.461755e-05, TOL_EPS, iter = 49
f(x^*) = 5.049139e+02, rel_eval = 2.552043e-03
Iteration 114:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.688792e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.467786e+02, rel_obj = 9.732117e-05, TOL_EPS, iter = 26
f(x^*) = 5.110860e+02, rel_eval = 1.207635e-02
Iteration 115:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.707829e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.461559e+02, rel_obj = 9.961725e-05, TOL_EPS, iter = 26
f(x^*) = 5.171340e+02, rel_eval = 1.169521e-02
Iteration 116:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.731259e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.450830e+02, rel_obj = 9.772022e-05, TOL_EPS, iter = 26
f(x^*) = 5.179732e+02, rel_eval = 1.620175e-03
Iteration 117:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.688260e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.442710e+02, rel_obj = 1.623547e-05, TOL_EPS, iter = 26
f(x^*) = 5.162678e+02, rel_eval = 3.303223e-03
Iteration 118:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.661862e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.365176e+02, rel_obj = 7.784416e-05, TOL_EPS, iter = 49
f(x^*) = 5.073716e+02, rel_eval = 1.753394e-02
Iteration 119:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.689149e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.361081e+02, rel_obj = 7.174079e-05, TOL_EPS, iter = 49
f(x^*) = 5.056150e+02, rel_eval = 3.474201e-03
Iteration 120:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.676421e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.361251e+02, rel_obj = 7.162849e-05, TOL_EPS, iter = 49
f(x^*) = 5.064115e+02, rel_eval = 1.572799e-03
Iteration 121:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.689845e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.360636e+02, rel_obj = 6.930210e-05, TOL_EPS, iter = 49
f(x^*) = 5.048064e+02, rel_eval = 3.179541e-03
Iteration 122:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.718640e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.358768e+02, rel_obj = 6.925172e-05, TOL_EPS, iter = 49
f(x^*) = 5.038804e+02, rel_eval = 1.837819e-03
Iteration 123:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.662115e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.427146e+02, rel_obj = 3.996131e-05, TOL_EPS, iter = 26
f(x^*) = 5.107747e+02, rel_eval = 1.349776e-02
Iteration 124:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.755962e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.353541e+02, rel_obj = 9.948048e-05, TOL_EPS, iter = 48
f(x^*) = 5.091521e+02, rel_eval = 3.186873e-03
Iteration 125:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.620866e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.419342e+02, rel_obj = 5.391066e-05, TOL_EPS, iter = 26
f(x^*) = 5.111336e+02, rel_eval = 3.876606e-03
Iteration 126:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.617377e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.346751e+02, rel_obj = 9.000759e-05, TOL_EPS, iter = 48
f(x^*) = 5.071134e+02, rel_eval = 7.927551e-03
Iteration 127:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.708255e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.343351e+02, rel_obj = 9.357795e-05, TOL_EPS, iter = 48
f(x^*) = 5.041455e+02, rel_eval = 5.887029e-03
Iteration 128:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.668925e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.342984e+02, rel_obj = 9.082321e-05, TOL_EPS, iter = 48
f(x^*) = 5.055145e+02, rel_eval = 2.708140e-03
Iteration 129:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.657285e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.342453e+02, rel_obj = 8.976582e-05, TOL_EPS, iter = 48
f(x^*) = 5.048738e+02, rel_eval = 1.269069e-03
Iteration 130:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.640060e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.412540e+02, rel_obj = 9.312129e-05, TOL_EPS, iter = 25
f(x^*) = 5.112004e+02, rel_eval = 1.237604e-02
Iteration 131:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.687132e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.407379e+02, rel_obj = 7.125908e-05, TOL_EPS, iter = 25
f(x^*) = 5.177799e+02, rel_eval = 1.270724e-02
Iteration 132:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.842083e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.398823e+02, rel_obj = 6.142927e-05, TOL_EPS, iter = 25
f(x^*) = 5.189342e+02, rel_eval = 2.224346e-03
Iteration 133:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.647545e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.392440e+02, rel_obj = 3.039531e-05, TOL_EPS, iter = 25
f(x^*) = 5.172352e+02, rel_eval = 3.284870e-03
Iteration 134:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.769539e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.310571e+02, rel_obj = 7.532765e-05, TOL_EPS, iter = 48
f(x^*) = 5.078518e+02, rel_eval = 1.847662e-02
Iteration 135:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.841949e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.306760e+02, rel_obj = 7.218525e-05, TOL_EPS, iter = 48
f(x^*) = 5.059853e+02, rel_eval = 3.688865e-03
Iteration 136:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.675863e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.387425e+02, rel_obj = 7.891800e-05, TOL_EPS, iter = 26
f(x^*) = 5.103269e+02, rel_eval = 8.507560e-03
Iteration 137:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.619767e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.306305e+02, rel_obj = 6.868328e-05, TOL_EPS, iter = 48
f(x^*) = 5.076302e+02, rel_eval = 5.312368e-03
Iteration 138:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.775519e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.301968e+02, rel_obj = 6.712617e-05, TOL_EPS, iter = 48
f(x^*) = 5.038791e+02, rel_eval = 7.444504e-03
Iteration 139:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.667778e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.300535e+02, rel_obj = 7.046140e-05, TOL_EPS, iter = 48
f(x^*) = 5.048002e+02, rel_eval = 1.824739e-03
Iteration 140:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.726234e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.300156e+02, rel_obj = 9.874628e-05, TOL_EPS, iter = 47
f(x^*) = 5.048961e+02, rel_eval = 1.900000e-04
Iteration 141:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.712220e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.298832e+02, rel_obj = 9.772865e-05, TOL_EPS, iter = 47
f(x^*) = 5.044205e+02, rel_eval = 9.429874e-04
Iteration 142:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.647839e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.374597e+02, rel_obj = 9.547916e-05, TOL_EPS, iter = 27
f(x^*) = 5.099375e+02, rel_eval = 1.081907e-02
Iteration 143:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.650270e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.377300e+02, rel_obj = 6.492428e-05, TOL_EPS, iter = 26
f(x^*) = 5.162121e+02, rel_eval = 1.215498e-02
Iteration 144:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.604842e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.375327e+02, rel_obj = 6.967735e-05, TOL_EPS, iter = 24
f(x^*) = 5.177130e+02, rel_eval = 2.899158e-03
Iteration 145:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.604422e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.371042e+02, rel_obj = 6.779801e-05, TOL_EPS, iter = 24
f(x^*) = 5.191798e+02, rel_eval = 2.825294e-03
Iteration 146:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.694268e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.283397e+02, rel_obj = 8.320834e-05, TOL_EPS, iter = 47
f(x^*) = 5.106836e+02, rel_eval = 1.663707e-02
Iteration 147:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.592520e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.279155e+02, rel_obj = 7.338307e-05, TOL_EPS, iter = 47
f(x^*) = 5.060198e+02, rel_eval = 9.216484e-03
Iteration 148:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.702442e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.279271e+02, rel_obj = 7.899135e-05, TOL_EPS, iter = 47
f(x^*) = 5.070866e+02, rel_eval = 2.103651e-03
Iteration 149:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.668272e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.278641e+02, rel_obj = 8.277234e-05, TOL_EPS, iter = 47
f(x^*) = 5.053579e+02, rel_eval = 3.420698e-03
Iteration 150:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.627429e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.276635e+02, rel_obj = 8.008706e-05, TOL_EPS, iter = 47
f(x^*) = 5.041055e+02, rel_eval = 2.484423e-03
Iteration 151:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.646271e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.356517e+02, rel_obj = 6.656945e-05, TOL_EPS, iter = 24
f(x^*) = 5.122493e+02, rel_eval = 1.589809e-02
Iteration 152:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.691188e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.270031e+02, rel_obj = 7.446383e-05, TOL_EPS, iter = 47
f(x^*) = 5.105725e+02, rel_eval = 3.284123e-03
Iteration 153:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.576616e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.347333e+02, rel_obj = 9.189528e-05, TOL_EPS, iter = 24
f(x^*) = 5.126542e+02, rel_eval = 4.060678e-03
Iteration 154:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.643746e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.262858e+02, rel_obj = 9.440255e-05, TOL_EPS, iter = 46
f(x^*) = 5.080317e+02, rel_eval = 9.098887e-03
Iteration 155:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.679116e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.259196e+02, rel_obj = 9.511982e-05, TOL_EPS, iter = 46
f(x^*) = 5.045127e+02, rel_eval = 6.975040e-03
Iteration 156:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.792177e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.258537e+02, rel_obj = 9.186389e-05, TOL_EPS, iter = 46
f(x^*) = 5.061116e+02, rel_eval = 3.159117e-03
Iteration 157:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.636041e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.257738e+02, rel_obj = 9.500786e-05, TOL_EPS, iter = 46
f(x^*) = 5.053354e+02, rel_eval = 1.535915e-03
Iteration 158:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.651028e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.256106e+02, rel_obj = 9.300290e-05, TOL_EPS, iter = 46
f(x^*) = 5.044901e+02, rel_eval = 1.675477e-03
Iteration 159:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.634811e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.254123e+02, rel_obj = 9.166895e-05, TOL_EPS, iter = 46
f(x^*) = 5.042984e+02, rel_eval = 3.801280e-04
Iteration 160:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.713735e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.334097e+02, rel_obj = 6.071582e-05, TOL_EPS, iter = 24
f(x^*) = 5.116104e+02, rel_eval = 1.429201e-02
Iteration 161:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.661434e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.249856e+02, rel_obj = 8.495566e-05, TOL_EPS, iter = 46
f(x^*) = 5.101983e+02, rel_eval = 2.767803e-03
Iteration 162:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.776660e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.329075e+02, rel_obj = 9.594540e-05, TOL_EPS, iter = 24
f(x^*) = 5.122506e+02, rel_eval = 4.006514e-03
Iteration 163:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.635454e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.244941e+02, rel_obj = 8.220083e-05, TOL_EPS, iter = 46
f(x^*) = 5.081424e+02, rel_eval = 8.084765e-03
Iteration 164:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.617189e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.242053e+02, rel_obj = 8.229282e-05, TOL_EPS, iter = 46
f(x^*) = 5.046703e+02, rel_eval = 6.879851e-03
Iteration 165:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.716248e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.241378e+02, rel_obj = 7.977098e-05, TOL_EPS, iter = 46
f(x^*) = 5.063139e+02, rel_eval = 3.246218e-03
Iteration 166:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.707778e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.240498e+02, rel_obj = 8.248997e-05, TOL_EPS, iter = 46
f(x^*) = 5.055233e+02, rel_eval = 1.563966e-03
Iteration 167:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.592021e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.238848e+02, rel_obj = 8.782100e-05, TOL_EPS, iter = 46
f(x^*) = 5.045907e+02, rel_eval = 1.848291e-03
Iteration 168:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.782662e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.236954e+02, rel_obj = 8.890357e-05, TOL_EPS, iter = 46
f(x^*) = 5.043368e+02, rel_eval = 5.032951e-04
Iteration 169:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.652495e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.318776e+02, rel_obj = 9.718835e-05, TOL_EPS, iter = 23
f(x^*) = 5.123837e+02, rel_eval = 1.570464e-02
Iteration 170:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.578833e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.315346e+02, rel_obj = 8.966175e-05, TOL_EPS, iter = 23
f(x^*) = 5.196565e+02, rel_eval = 1.399557e-02
Iteration 171:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.683384e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.310148e+02, rel_obj = 2.691544e-05, TOL_EPS, iter = 23
f(x^*) = 5.209687e+02, rel_eval = 2.518666e-03
Iteration 172:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.663225e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.218812e+02, rel_obj = 7.365451e-05, TOL_EPS, iter = 46
f(x^*) = 5.104812e+02, rel_eval = 2.054439e-02
Iteration 173:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.682264e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.302719e+02, rel_obj = 9.557333e-05, TOL_EPS, iter = 22
f(x^*) = 5.127663e+02, rel_eval = 4.456519e-03
Iteration 174:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.704557e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.304765e+02, rel_obj = 3.329160e-05, TOL_EPS, iter = 24
f(x^*) = 5.148706e+02, rel_eval = 4.086995e-03
Iteration 175:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.784753e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.299859e+02, rel_obj = 7.672830e-05, TOL_EPS, iter = 21
f(x^*) = 5.183968e+02, rel_eval = 6.802165e-03
Iteration 176:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.694454e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.290287e+02, rel_obj = 3.552863e-05, TOL_EPS, iter = 21
f(x^*) = 5.215400e+02, rel_eval = 6.026662e-03
Iteration 177:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.685766e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.283029e+02, rel_obj = 1.236127e-05, TOL_EPS, iter = 23
f(x^*) = 5.192210e+02, rel_eval = 4.466183e-03
Iteration 178:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.751613e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.280271e+02, rel_obj = 9.729774e-05, TOL_EPS, iter = 23
f(x^*) = 5.180346e+02, rel_eval = 2.290205e-03
Iteration 179:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.639082e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.181045e+02, rel_obj = 8.213758e-05, TOL_EPS, iter = 44
f(x^*) = 5.096391e+02, rel_eval = 1.647347e-02
Iteration 180:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.745451e-14, TOL_EPS, iter = 0
Prox_TV: obj = 9.279243e+02, rel_obj = 5.030567e-05, TOL_EPS, iter = 24
f(x^*) = 5.118180e+02, rel_eval = 4.257224e-03
PPXA:
f(x^*) = 5.118180e+02, rel_eval = 4.257224e-03
180 iterations
Stopping criterion: MAX_IT
estimation_time =
21.0500
Reconstruction_MSE =
0.0012
Sources_MSE =
0.0017
The time step is set manually to : 8.65781
Algorithm selected: PPXA
Iteration 1:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 8.507143e-15, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 0.000000e+00, REL_OBJ, iter = 1
f(x^*) = 8.324966e+02, rel_eval = 1.000000e+00
Iteration 2:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.046310e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 5.155743e+01, REL_OBJ, iter = 1
f(x^*) = 1.226197e+03, rel_eval = 3.210745e-01
Iteration 3:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.244523e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.360885e+02, REL_OBJ, iter = 1
f(x^*) = 1.503042e+03, rel_eval = 1.841897e-01
Iteration 4:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.326920e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.945618e+02, REL_OBJ, iter = 1
f(x^*) = 1.668201e+03, rel_eval = 9.900432e-02
Iteration 5:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.392998e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 4.621560e+02, REL_OBJ, iter = 1
f(x^*) = 1.696537e+03, rel_eval = 1.670214e-02
Iteration 6:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.157584e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 4.586981e+02, REL_OBJ, iter = 1
f(x^*) = 1.639505e+03, rel_eval = 3.478584e-02
Iteration 7:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.207639e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 4.243393e+02, REL_OBJ, iter = 1
f(x^*) = 1.573783e+03, rel_eval = 4.176064e-02
Iteration 8:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.171374e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.925352e+02, REL_OBJ, iter = 1
f(x^*) = 1.511921e+03, rel_eval = 4.091592e-02
Iteration 9:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.131759e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.731140e+02, REL_OBJ, iter = 1
f(x^*) = 1.446354e+03, rel_eval = 4.533306e-02
Iteration 10:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.114396e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.608897e+02, REL_OBJ, iter = 1
f(x^*) = 1.382333e+03, rel_eval = 4.631380e-02
Iteration 11:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.114491e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.503224e+02, REL_OBJ, iter = 1
f(x^*) = 1.329493e+03, rel_eval = 3.974434e-02
Iteration 12:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.153045e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.396193e+02, REL_OBJ, iter = 1
f(x^*) = 1.286633e+03, rel_eval = 3.331195e-02
Iteration 13:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.231835e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.297077e+02, REL_OBJ, iter = 1
f(x^*) = 1.249938e+03, rel_eval = 2.935718e-02
Iteration 14:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.287102e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.209717e+02, REL_OBJ, iter = 1
f(x^*) = 1.214782e+03, rel_eval = 2.894049e-02
Iteration 15:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.261321e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.132621e+02, REL_OBJ, iter = 1
f(x^*) = 1.182156e+03, rel_eval = 2.759834e-02
Iteration 16:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.259155e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.066014e+02, REL_OBJ, iter = 1
f(x^*) = 1.153681e+03, rel_eval = 2.468223e-02
Iteration 17:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.321697e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.007915e+02, REL_OBJ, iter = 1
f(x^*) = 1.130425e+03, rel_eval = 2.057232e-02
Iteration 18:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.263338e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.963839e+02, REL_OBJ, iter = 1
f(x^*) = 1.112075e+03, rel_eval = 1.650125e-02
Iteration 19:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.278689e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.928470e+02, REL_OBJ, iter = 1
f(x^*) = 1.095982e+03, rel_eval = 1.468364e-02
Iteration 20:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.278512e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.901999e+02, REL_OBJ, iter = 1
f(x^*) = 1.082835e+03, rel_eval = 1.214087e-02
Iteration 21:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.351759e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.883007e+02, REL_OBJ, iter = 1
f(x^*) = 1.072402e+03, rel_eval = 9.728230e-03
Iteration 22:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.359753e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.872513e+02, REL_OBJ, iter = 1
f(x^*) = 1.064918e+03, rel_eval = 7.028581e-03
Iteration 23:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.346758e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.865047e+02, REL_OBJ, iter = 1
f(x^*) = 1.059409e+03, rel_eval = 5.199939e-03
Iteration 24:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.368314e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.859210e+02, REL_OBJ, iter = 1
f(x^*) = 1.055089e+03, rel_eval = 4.094129e-03
Iteration 25:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.414030e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.855439e+02, REL_OBJ, iter = 1
f(x^*) = 1.051962e+03, rel_eval = 2.973016e-03
Iteration 26:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.435071e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.855441e+02, REL_OBJ, iter = 1
f(x^*) = 1.050442e+03, rel_eval = 1.446150e-03
Iteration 27:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.285195e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.859198e+02, REL_OBJ, iter = 1
f(x^*) = 1.049867e+03, rel_eval = 5.479442e-04
Iteration 28:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.383479e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.865002e+02, REL_OBJ, iter = 1
f(x^*) = 1.049804e+03, rel_eval = 6.018943e-05
Iteration 29:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.360231e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.871783e+02, REL_OBJ, iter = 1
f(x^*) = 1.050064e+03, rel_eval = 2.478197e-04
Iteration 30:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.432180e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.879159e+02, REL_OBJ, iter = 1
f(x^*) = 1.050585e+03, rel_eval = 4.958828e-04
Iteration 31:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.335577e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.887273e+02, REL_OBJ, iter = 1
f(x^*) = 1.051144e+03, rel_eval = 5.315579e-04
Iteration 32:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.406140e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.895480e+02, REL_OBJ, iter = 1
f(x^*) = 1.051793e+03, rel_eval = 6.167652e-04
Iteration 33:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.388073e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.903936e+02, REL_OBJ, iter = 1
f(x^*) = 1.052756e+03, rel_eval = 9.151709e-04
Iteration 34:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.351463e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.912854e+02, REL_OBJ, iter = 1
f(x^*) = 1.054279e+03, rel_eval = 1.444235e-03
Iteration 35:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.401191e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.920899e+02, REL_OBJ, iter = 1
f(x^*) = 1.055335e+03, rel_eval = 1.001009e-03
Iteration 36:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.402067e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.928526e+02, REL_OBJ, iter = 1
f(x^*) = 1.056054e+03, rel_eval = 6.807239e-04
Iteration 37:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.417990e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.936465e+02, REL_OBJ, iter = 1
f(x^*) = 1.057137e+03, rel_eval = 1.024760e-03
Iteration 38:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.338242e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.943948e+02, REL_OBJ, iter = 1
f(x^*) = 1.058199e+03, rel_eval = 1.002826e-03
Iteration 39:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.389635e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.950459e+02, REL_OBJ, iter = 1
f(x^*) = 1.059181e+03, rel_eval = 9.277650e-04
Iteration 40:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.428534e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.957137e+02, REL_OBJ, iter = 1
f(x^*) = 1.060181e+03, rel_eval = 9.434844e-04
Iteration 41:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.413689e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.962983e+02, REL_OBJ, iter = 1
f(x^*) = 1.060950e+03, rel_eval = 7.240300e-04
Iteration 42:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.413162e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.968388e+02, REL_OBJ, iter = 1
f(x^*) = 1.061510e+03, rel_eval = 5.280311e-04
Iteration 43:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.415205e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.973194e+02, REL_OBJ, iter = 1
f(x^*) = 1.062087e+03, rel_eval = 5.432180e-04
Iteration 44:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.437755e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.977769e+02, REL_OBJ, iter = 1
f(x^*) = 1.062633e+03, rel_eval = 5.134144e-04
Iteration 45:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.395852e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.982044e+02, REL_OBJ, iter = 1
f(x^*) = 1.063045e+03, rel_eval = 3.877610e-04
Iteration 46:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.367866e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.985725e+02, REL_OBJ, iter = 1
f(x^*) = 1.063198e+03, rel_eval = 1.442374e-04
Iteration 47:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.380284e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.989038e+02, REL_OBJ, iter = 1
f(x^*) = 1.063202e+03, rel_eval = 3.467588e-06
Iteration 48:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.375807e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.992382e+02, REL_OBJ, iter = 1
f(x^*) = 1.063146e+03, rel_eval = 5.235601e-05
Iteration 49:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.386149e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.995709e+02, REL_OBJ, iter = 1
f(x^*) = 1.063140e+03, rel_eval = 5.739794e-06
Iteration 50:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.448765e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 2.998905e+02, REL_OBJ, iter = 1
f(x^*) = 1.063345e+03, rel_eval = 1.930652e-04
Iteration 51:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.417419e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.002066e+02, REL_OBJ, iter = 1
f(x^*) = 1.063670e+03, rel_eval = 3.048238e-04
Iteration 52:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.435826e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.005183e+02, REL_OBJ, iter = 1
f(x^*) = 1.063994e+03, rel_eval = 3.048753e-04
Iteration 53:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.425987e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.007888e+02, REL_OBJ, iter = 1
f(x^*) = 1.064163e+03, rel_eval = 1.589238e-04
Iteration 54:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.395781e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.010237e+02, REL_OBJ, iter = 1
f(x^*) = 1.064252e+03, rel_eval = 8.325437e-05
Iteration 55:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.391959e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.012300e+02, REL_OBJ, iter = 1
f(x^*) = 1.064392e+03, rel_eval = 1.320897e-04
Iteration 56:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.427977e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.014269e+02, REL_OBJ, iter = 1
f(x^*) = 1.064528e+03, rel_eval = 1.276717e-04
Iteration 57:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.405086e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.016191e+02, REL_OBJ, iter = 1
f(x^*) = 1.064567e+03, rel_eval = 3.654754e-05
Iteration 58:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.439886e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.017944e+02, REL_OBJ, iter = 1
f(x^*) = 1.064552e+03, rel_eval = 1.409155e-05
Iteration 59:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.399174e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.019634e+02, REL_OBJ, iter = 1
f(x^*) = 1.064537e+03, rel_eval = 1.435018e-05
Iteration 60:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.440911e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.021224e+02, REL_OBJ, iter = 1
f(x^*) = 1.064552e+03, rel_eval = 1.413620e-05
Iteration 61:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.376744e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.022665e+02, REL_OBJ, iter = 1
f(x^*) = 1.064519e+03, rel_eval = 3.105970e-05
Iteration 62:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.447642e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.023856e+02, REL_OBJ, iter = 1
f(x^*) = 1.064486e+03, rel_eval = 3.120049e-05
Iteration 63:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.423614e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.024983e+02, REL_OBJ, iter = 1
f(x^*) = 1.064473e+03, rel_eval = 1.222721e-05
Iteration 64:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.427367e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.026103e+02, REL_OBJ, iter = 1
f(x^*) = 1.064492e+03, rel_eval = 1.823563e-05
Iteration 65:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.429346e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.027150e+02, REL_OBJ, iter = 1
f(x^*) = 1.064528e+03, rel_eval = 3.359791e-05
Iteration 66:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.392680e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.028205e+02, REL_OBJ, iter = 1
f(x^*) = 1.064522e+03, rel_eval = 5.760139e-06
Iteration 67:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.410144e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.029420e+02, REL_OBJ, iter = 1
f(x^*) = 1.064532e+03, rel_eval = 9.856490e-06
Iteration 68:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.521184e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.030642e+02, REL_OBJ, iter = 1
f(x^*) = 1.064552e+03, rel_eval = 1.838635e-05
Iteration 69:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.442261e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.031928e+02, REL_OBJ, iter = 1
f(x^*) = 1.064616e+03, rel_eval = 6.046057e-05
Iteration 70:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.401450e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.033379e+02, REL_OBJ, iter = 1
f(x^*) = 1.064696e+03, rel_eval = 7.467717e-05
Iteration 71:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.436518e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.034728e+02, REL_OBJ, iter = 1
f(x^*) = 1.064798e+03, rel_eval = 9.569798e-05
Iteration 72:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.438961e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.035930e+02, REL_OBJ, iter = 1
f(x^*) = 1.064882e+03, rel_eval = 7.905132e-05
Iteration 73:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.480084e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.037111e+02, REL_OBJ, iter = 1
f(x^*) = 1.064964e+03, rel_eval = 7.766965e-05
Iteration 74:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.458336e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.038242e+02, REL_OBJ, iter = 1
f(x^*) = 1.065052e+03, rel_eval = 8.251223e-05
Iteration 75:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.403431e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.039387e+02, REL_OBJ, iter = 1
f(x^*) = 1.065152e+03, rel_eval = 9.360627e-05
Iteration 76:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.431016e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.040507e+02, REL_OBJ, iter = 1
f(x^*) = 1.065249e+03, rel_eval = 9.118765e-05
Iteration 77:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.445941e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.041589e+02, REL_OBJ, iter = 1
f(x^*) = 1.065323e+03, rel_eval = 6.947463e-05
Iteration 78:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.459415e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.042516e+02, REL_OBJ, iter = 1
f(x^*) = 1.065383e+03, rel_eval = 5.609943e-05
Iteration 79:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.446165e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.043387e+02, REL_OBJ, iter = 1
f(x^*) = 1.065465e+03, rel_eval = 7.710514e-05
Iteration 80:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.469652e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.044111e+02, REL_OBJ, iter = 1
f(x^*) = 1.065518e+03, rel_eval = 4.955906e-05
Iteration 81:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.429385e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.044687e+02, REL_OBJ, iter = 1
f(x^*) = 1.065525e+03, rel_eval = 6.710305e-06
Iteration 82:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.438445e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.045276e+02, REL_OBJ, iter = 1
f(x^*) = 1.065551e+03, rel_eval = 2.424776e-05
Iteration 83:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.450750e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.045839e+02, REL_OBJ, iter = 1
f(x^*) = 1.065541e+03, rel_eval = 9.212922e-06
Iteration 84:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.408928e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.046285e+02, REL_OBJ, iter = 1
f(x^*) = 1.065501e+03, rel_eval = 3.750579e-05
Iteration 85:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.464229e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.046567e+02, REL_OBJ, iter = 1
f(x^*) = 1.065426e+03, rel_eval = 7.056023e-05
Iteration 86:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.432159e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.046700e+02, REL_OBJ, iter = 1
f(x^*) = 1.065297e+03, rel_eval = 1.214624e-04
Iteration 87:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.488635e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.046773e+02, REL_OBJ, iter = 1
f(x^*) = 1.065151e+03, rel_eval = 1.369636e-04
Iteration 88:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.406668e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.046858e+02, REL_OBJ, iter = 1
f(x^*) = 1.065008e+03, rel_eval = 1.337739e-04
Iteration 89:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.372665e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.046955e+02, REL_OBJ, iter = 1
f(x^*) = 1.064893e+03, rel_eval = 1.083353e-04
Iteration 90:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.441049e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.047125e+02, REL_OBJ, iter = 1
f(x^*) = 1.064814e+03, rel_eval = 7.354600e-05
Iteration 91:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.400841e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.047533e+02, REL_OBJ, iter = 1
f(x^*) = 1.064804e+03, rel_eval = 1.024571e-05
Iteration 92:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.474819e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.048123e+02, REL_OBJ, iter = 1
f(x^*) = 1.064879e+03, rel_eval = 7.085258e-05
Iteration 93:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.451247e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.048794e+02, REL_OBJ, iter = 1
f(x^*) = 1.064995e+03, rel_eval = 1.086992e-04
Iteration 94:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.430450e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.049537e+02, REL_OBJ, iter = 1
f(x^*) = 1.065143e+03, rel_eval = 1.388636e-04
Iteration 95:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.453730e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.050417e+02, REL_OBJ, iter = 1
f(x^*) = 1.065335e+03, rel_eval = 1.801305e-04
Iteration 96:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.459143e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.051237e+02, REL_OBJ, iter = 1
f(x^*) = 1.065510e+03, rel_eval = 1.642568e-04
Iteration 97:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.509093e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.052071e+02, REL_OBJ, iter = 1
f(x^*) = 1.065673e+03, rel_eval = 1.530028e-04
Iteration 98:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.429901e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.052858e+02, REL_OBJ, iter = 1
f(x^*) = 1.065813e+03, rel_eval = 1.312882e-04
Iteration 99:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.477817e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.053608e+02, REL_OBJ, iter = 1
f(x^*) = 1.065925e+03, rel_eval = 1.051880e-04
Iteration 100:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.454807e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.054399e+02, REL_OBJ, iter = 1
f(x^*) = 1.066035e+03, rel_eval = 1.033180e-04
Iteration 101:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.454427e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.055081e+02, REL_OBJ, iter = 1
f(x^*) = 1.066116e+03, rel_eval = 7.646394e-05
Iteration 102:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.417848e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.055629e+02, REL_OBJ, iter = 1
f(x^*) = 1.066183e+03, rel_eval = 6.211167e-05
Iteration 103:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.444682e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.056158e+02, REL_OBJ, iter = 1
f(x^*) = 1.066239e+03, rel_eval = 5.262533e-05
Iteration 104:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.442819e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.056639e+02, REL_OBJ, iter = 1
f(x^*) = 1.066276e+03, rel_eval = 3.462467e-05
Iteration 105:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.452390e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.057098e+02, REL_OBJ, iter = 1
f(x^*) = 1.066308e+03, rel_eval = 3.046528e-05
Iteration 106:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.433755e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.057534e+02, REL_OBJ, iter = 1
f(x^*) = 1.066354e+03, rel_eval = 4.323438e-05
Iteration 107:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.459357e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.057983e+02, REL_OBJ, iter = 1
f(x^*) = 1.066400e+03, rel_eval = 4.309871e-05
Iteration 108:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.418049e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.058440e+02, REL_OBJ, iter = 1
f(x^*) = 1.066445e+03, rel_eval = 4.230975e-05
Iteration 109:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.445951e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.058874e+02, REL_OBJ, iter = 1
f(x^*) = 1.066484e+03, rel_eval = 3.596375e-05
Iteration 110:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.423238e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.059271e+02, REL_OBJ, iter = 1
f(x^*) = 1.066535e+03, rel_eval = 4.784915e-05
Iteration 111:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.439503e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.059740e+02, REL_OBJ, iter = 1
f(x^*) = 1.066602e+03, rel_eval = 6.264916e-05
Iteration 112:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.444950e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.060247e+02, REL_OBJ, iter = 1
f(x^*) = 1.066725e+03, rel_eval = 1.152997e-04
Iteration 113:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.496133e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.060780e+02, REL_OBJ, iter = 1
f(x^*) = 1.066860e+03, rel_eval = 1.265420e-04
Iteration 114:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.414492e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.061331e+02, REL_OBJ, iter = 1
f(x^*) = 1.066983e+03, rel_eval = 1.156930e-04
Iteration 115:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.512818e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.061842e+02, REL_OBJ, iter = 1
f(x^*) = 1.067086e+03, rel_eval = 9.677803e-05
Iteration 116:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.396703e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.062272e+02, REL_OBJ, iter = 1
f(x^*) = 1.067156e+03, rel_eval = 6.552648e-05
Iteration 117:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.407691e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.062628e+02, REL_OBJ, iter = 1
f(x^*) = 1.067198e+03, rel_eval = 3.952609e-05
Iteration 118:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.429835e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.062941e+02, REL_OBJ, iter = 1
f(x^*) = 1.067229e+03, rel_eval = 2.872134e-05
Iteration 119:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.418982e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.063181e+02, REL_OBJ, iter = 1
f(x^*) = 1.067240e+03, rel_eval = 1.070987e-05
Iteration 120:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.473753e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.063406e+02, REL_OBJ, iter = 1
f(x^*) = 1.067224e+03, rel_eval = 1.502675e-05
Iteration 121:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.461782e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.063558e+02, REL_OBJ, iter = 1
f(x^*) = 1.067190e+03, rel_eval = 3.264267e-05
Iteration 122:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.434102e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.063689e+02, REL_OBJ, iter = 1
f(x^*) = 1.067149e+03, rel_eval = 3.766956e-05
Iteration 123:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.489292e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.063783e+02, REL_OBJ, iter = 1
f(x^*) = 1.067097e+03, rel_eval = 4.885456e-05
Iteration 124:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.485391e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.063842e+02, REL_OBJ, iter = 1
f(x^*) = 1.067027e+03, rel_eval = 6.584992e-05
Iteration 125:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.515321e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.063860e+02, REL_OBJ, iter = 1
f(x^*) = 1.066965e+03, rel_eval = 5.831218e-05
Iteration 126:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.457158e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.063892e+02, REL_OBJ, iter = 1
f(x^*) = 1.066916e+03, rel_eval = 4.532685e-05
Iteration 127:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.449073e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.063933e+02, REL_OBJ, iter = 1
f(x^*) = 1.066885e+03, rel_eval = 2.918032e-05
Iteration 128:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.420801e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.064016e+02, REL_OBJ, iter = 1
f(x^*) = 1.066866e+03, rel_eval = 1.807807e-05
Iteration 129:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.501085e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.064139e+02, REL_OBJ, iter = 1
f(x^*) = 1.066856e+03, rel_eval = 9.115549e-06
Iteration 130:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.482262e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.064311e+02, REL_OBJ, iter = 1
f(x^*) = 1.066857e+03, rel_eval = 1.185651e-06
Iteration 131:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.430672e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.064513e+02, REL_OBJ, iter = 1
f(x^*) = 1.066865e+03, rel_eval = 6.717695e-06
Iteration 132:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.457659e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.064716e+02, REL_OBJ, iter = 1
f(x^*) = 1.066877e+03, rel_eval = 1.144912e-05
Iteration 133:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.476788e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.064905e+02, REL_OBJ, iter = 1
f(x^*) = 1.066889e+03, rel_eval = 1.177225e-05
Iteration 134:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.482726e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.065090e+02, REL_OBJ, iter = 1
f(x^*) = 1.066899e+03, rel_eval = 9.399080e-06
Iteration 135:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.399861e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.065275e+02, REL_OBJ, iter = 1
f(x^*) = 1.066910e+03, rel_eval = 1.012861e-05
Iteration 136:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.525979e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.065410e+02, REL_OBJ, iter = 1
f(x^*) = 1.066912e+03, rel_eval = 1.692820e-06
Iteration 137:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.444053e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.065518e+02, REL_OBJ, iter = 1
f(x^*) = 1.066902e+03, rel_eval = 9.478910e-06
Iteration 138:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.481525e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.065588e+02, REL_OBJ, iter = 1
f(x^*) = 1.066879e+03, rel_eval = 2.188709e-05
Iteration 139:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.419255e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.065653e+02, REL_OBJ, iter = 1
f(x^*) = 1.066848e+03, rel_eval = 2.899126e-05
Iteration 140:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.421777e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.065715e+02, REL_OBJ, iter = 1
f(x^*) = 1.066806e+03, rel_eval = 3.870095e-05
Iteration 141:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.487848e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.065797e+02, REL_OBJ, iter = 1
f(x^*) = 1.066766e+03, rel_eval = 3.784288e-05
Iteration 142:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.454798e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.065870e+02, REL_OBJ, iter = 1
f(x^*) = 1.066728e+03, rel_eval = 3.593024e-05
Iteration 143:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.451519e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.065975e+02, REL_OBJ, iter = 1
f(x^*) = 1.066700e+03, rel_eval = 2.562918e-05
Iteration 144:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.541701e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.066138e+02, REL_OBJ, iter = 1
f(x^*) = 1.066695e+03, rel_eval = 4.686122e-06
Iteration 145:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.456043e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.066328e+02, REL_OBJ, iter = 1
f(x^*) = 1.066702e+03, rel_eval = 6.096987e-06
Iteration 146:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.429742e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.066547e+02, REL_OBJ, iter = 1
f(x^*) = 1.066720e+03, rel_eval = 1.701102e-05
Iteration 147:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.474613e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.066793e+02, REL_OBJ, iter = 1
f(x^*) = 1.066751e+03, rel_eval = 2.937593e-05
Iteration 148:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.460706e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.067037e+02, REL_OBJ, iter = 1
f(x^*) = 1.066790e+03, rel_eval = 3.665080e-05
Iteration 149:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.424009e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.067298e+02, REL_OBJ, iter = 1
f(x^*) = 1.066832e+03, rel_eval = 3.856028e-05
Iteration 150:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.476301e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.067544e+02, REL_OBJ, iter = 1
f(x^*) = 1.066881e+03, rel_eval = 4.644311e-05
Iteration 151:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.463090e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.067774e+02, REL_OBJ, iter = 1
f(x^*) = 1.066921e+03, rel_eval = 3.720220e-05
Iteration 152:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.490365e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.067966e+02, REL_OBJ, iter = 1
f(x^*) = 1.066952e+03, rel_eval = 2.911382e-05
Iteration 153:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.423637e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.068117e+02, REL_OBJ, iter = 1
f(x^*) = 1.066972e+03, rel_eval = 1.859371e-05
Iteration 154:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.434264e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.068230e+02, REL_OBJ, iter = 1
f(x^*) = 1.066981e+03, rel_eval = 8.598136e-06
Iteration 155:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.432444e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.068351e+02, REL_OBJ, iter = 1
f(x^*) = 1.066990e+03, rel_eval = 8.512206e-06
Iteration 156:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.405151e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.068457e+02, REL_OBJ, iter = 1
f(x^*) = 1.067000e+03, rel_eval = 9.666510e-06
Iteration 157:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.473686e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.068553e+02, REL_OBJ, iter = 1
f(x^*) = 1.067013e+03, rel_eval = 1.176776e-05
Iteration 158:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.502548e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.068663e+02, REL_OBJ, iter = 1
f(x^*) = 1.067032e+03, rel_eval = 1.773969e-05
Iteration 159:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.460885e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.068779e+02, REL_OBJ, iter = 1
f(x^*) = 1.067057e+03, rel_eval = 2.321841e-05
Iteration 160:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.458137e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.068908e+02, REL_OBJ, iter = 1
f(x^*) = 1.067085e+03, rel_eval = 2.681628e-05
Iteration 161:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.457584e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.069037e+02, REL_OBJ, iter = 1
f(x^*) = 1.067108e+03, rel_eval = 2.110177e-05
Iteration 162:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.423356e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.069195e+02, REL_OBJ, iter = 1
f(x^*) = 1.067126e+03, rel_eval = 1.734245e-05
Iteration 163:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.497593e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.069355e+02, REL_OBJ, iter = 1
f(x^*) = 1.067144e+03, rel_eval = 1.627408e-05
Iteration 164:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.461290e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.069521e+02, REL_OBJ, iter = 1
f(x^*) = 1.067159e+03, rel_eval = 1.436502e-05
Iteration 165:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.419311e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.069696e+02, REL_OBJ, iter = 1
f(x^*) = 1.067177e+03, rel_eval = 1.709565e-05
Iteration 166:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.470080e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.069869e+02, REL_OBJ, iter = 1
f(x^*) = 1.067204e+03, rel_eval = 2.548681e-05
Iteration 167:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.408482e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.070071e+02, REL_OBJ, iter = 1
f(x^*) = 1.067237e+03, rel_eval = 3.033719e-05
Iteration 168:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.483948e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.070306e+02, REL_OBJ, iter = 1
f(x^*) = 1.067279e+03, rel_eval = 3.999333e-05
Iteration 169:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.501275e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.070567e+02, REL_OBJ, iter = 1
f(x^*) = 1.067334e+03, rel_eval = 5.163102e-05
Iteration 170:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.464539e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.070832e+02, REL_OBJ, iter = 1
f(x^*) = 1.067387e+03, rel_eval = 4.897845e-05
Iteration 171:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.473450e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.071082e+02, REL_OBJ, iter = 1
f(x^*) = 1.067436e+03, rel_eval = 4.631763e-05
Iteration 172:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.472616e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.071356e+02, REL_OBJ, iter = 1
f(x^*) = 1.067489e+03, rel_eval = 4.947601e-05
Iteration 173:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.456063e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.071632e+02, REL_OBJ, iter = 1
f(x^*) = 1.067543e+03, rel_eval = 5.074503e-05
Iteration 174:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.498110e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.071918e+02, REL_OBJ, iter = 1
f(x^*) = 1.067601e+03, rel_eval = 5.445819e-05
Iteration 175:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.480379e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.072183e+02, REL_OBJ, iter = 1
f(x^*) = 1.067655e+03, rel_eval = 5.012765e-05
Iteration 176:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.509557e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.072416e+02, REL_OBJ, iter = 1
f(x^*) = 1.067703e+03, rel_eval = 4.464499e-05
Iteration 177:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.502884e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.072606e+02, REL_OBJ, iter = 1
f(x^*) = 1.067741e+03, rel_eval = 3.645046e-05
Iteration 178:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.498929e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.072786e+02, REL_OBJ, iter = 1
f(x^*) = 1.067775e+03, rel_eval = 3.136740e-05
Iteration 179:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.481977e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.072957e+02, REL_OBJ, iter = 1
f(x^*) = 1.067810e+03, rel_eval = 3.317855e-05
Iteration 180:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.467496e-14, TOL_EPS, iter = 0
prox_L1: ||A x-y||_1 = 3.073118e+02, REL_OBJ, iter = 1
f(x^*) = 1.067842e+03, rel_eval = 2.978034e-05
PPXA:
f(x^*) = 1.067842e+03, rel_eval = 2.978034e-05
180 iterations
Stopping criterion: MAX_IT
estimation_time =
4.4100
Reconstruction_MSE =
0.0528
Sources_MSE =
0.1347
The time step is set manually to : 1
Algorithm selected: PPXA
Iteration 1:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.386538e+00, MAX_IT, iter = 100
Prox_TV: obj = 0.000000e+00, rel_obj = NaN, MAX_IT, iter = 200
f(x^*) = 2.602217e+02, rel_eval = 1.000000e+00
Iteration 2:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.579161e+00, MAX_IT, iter = 100
Prox_TV: obj = 9.801191e+01, rel_obj = 6.933248e-05, TOL_EPS, iter = 14
f(x^*) = 4.823645e+02, rel_eval = 4.605290e-01
Iteration 3:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.493954e+00, MAX_IT, iter = 100
Prox_TV: obj = 1.712638e+02, rel_obj = 7.061559e-05, TOL_EPS, iter = 46
f(x^*) = 5.482431e+02, rel_eval = 1.201631e-01
Iteration 4:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.295418e+00, MAX_IT, iter = 100
Prox_TV: obj = 2.077407e+02, rel_obj = 2.293934e-06, TOL_EPS, iter = 51
f(x^*) = 5.371850e+02, rel_eval = 2.058528e-02
Iteration 5:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.109137e+00, MAX_IT, iter = 100
Prox_TV: obj = 1.975804e+02, rel_obj = 1.803965e-05, TOL_EPS, iter = 50
f(x^*) = 5.150468e+02, rel_eval = 4.298277e-02
Iteration 6:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 9.678108e-01, MAX_IT, iter = 100
Prox_TV: obj = 1.815432e+02, rel_obj = 8.644962e-05, TOL_EPS, iter = 49
f(x^*) = 5.111531e+02, rel_eval = 7.617629e-03
Iteration 7:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 8.591547e-01, MAX_IT, iter = 100
Prox_TV: obj = 1.765775e+02, rel_obj = 5.929779e-06, TOL_EPS, iter = 55
f(x^*) = 5.127840e+02, rel_eval = 3.180617e-03
Iteration 8:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 7.636999e-01, MAX_IT, iter = 100
Prox_TV: obj = 1.802669e+02, rel_obj = 2.149718e-05, TOL_EPS, iter = 71
f(x^*) = 5.127925e+02, rel_eval = 1.653152e-05
PPXA:
f(x^*) = 5.127925e+02, rel_eval = 1.653152e-05
8 iterations
Stopping criterion: REL_NORM_OBJ
estimation_time =
41.8900
Reconstruction_MSE =
0.0210
Sources_MSE =
0.2025
The time step is set manually to : 8.65781
Algorithm selected: PPXA
Iteration 1:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.394695e+00, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 0.000000e+00, REL_OBJ, iter = 1
f(x^*) = 3.619181e+02, rel_eval = 1.000000e+00
Iteration 2:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.591772e+00, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 8.460247e+01, REL_OBJ, iter = 1
f(x^*) = 6.668014e+02, rel_eval = 4.572324e-01
Iteration 3:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.516329e+00, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 1.940653e+02, REL_OBJ, iter = 1
f(x^*) = 8.339523e+02, rel_eval = 2.004322e-01
Iteration 4:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.325904e+00, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.506406e+02, REL_OBJ, iter = 1
f(x^*) = 8.949361e+02, rel_eval = 6.814316e-02
Iteration 5:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.144797e+00, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.508442e+02, REL_OBJ, iter = 1
f(x^*) = 9.008315e+02, rel_eval = 6.544486e-03
Iteration 6:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 1.018726e+00, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.392299e+02, REL_OBJ, iter = 1
f(x^*) = 8.992195e+02, rel_eval = 1.792731e-03
Iteration 7:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 9.241918e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.360853e+02, REL_OBJ, iter = 1
f(x^*) = 9.034586e+02, rel_eval = 4.692136e-03
Iteration 8:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 8.459033e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.420883e+02, REL_OBJ, iter = 1
f(x^*) = 9.153466e+02, rel_eval = 1.298737e-02
Iteration 9:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 7.775703e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.506803e+02, REL_OBJ, iter = 1
f(x^*) = 9.291197e+02, rel_eval = 1.482381e-02
Iteration 10:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 7.192597e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.578433e+02, REL_OBJ, iter = 1
f(x^*) = 9.414054e+02, rel_eval = 1.305038e-02
Iteration 11:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 6.730167e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.625463e+02, REL_OBJ, iter = 1
f(x^*) = 9.511286e+02, rel_eval = 1.022278e-02
Iteration 12:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 6.385809e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.659352e+02, REL_OBJ, iter = 1
f(x^*) = 9.594569e+02, rel_eval = 8.680238e-03
Iteration 13:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 6.122409e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.687991e+02, REL_OBJ, iter = 1
f(x^*) = 9.665447e+02, rel_eval = 7.333147e-03
Iteration 14:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 5.896034e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.712144e+02, REL_OBJ, iter = 1
f(x^*) = 9.726193e+02, rel_eval = 6.245633e-03
Iteration 15:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 5.685297e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.732245e+02, REL_OBJ, iter = 1
f(x^*) = 9.777085e+02, rel_eval = 5.205187e-03
Iteration 16:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 5.493263e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.750768e+02, REL_OBJ, iter = 1
f(x^*) = 9.822322e+02, rel_eval = 4.605554e-03
Iteration 17:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 5.325385e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.766984e+02, REL_OBJ, iter = 1
f(x^*) = 9.862913e+02, rel_eval = 4.115531e-03
Iteration 18:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 5.178787e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.782390e+02, REL_OBJ, iter = 1
f(x^*) = 9.902150e+02, rel_eval = 3.962498e-03
Iteration 19:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 5.045384e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.796931e+02, REL_OBJ, iter = 1
f(x^*) = 9.940809e+02, rel_eval = 3.888865e-03
Iteration 20:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 4.920804e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.811304e+02, REL_OBJ, iter = 1
f(x^*) = 9.979077e+02, rel_eval = 3.834850e-03
Iteration 21:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 4.803730e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.824750e+02, REL_OBJ, iter = 1
f(x^*) = 1.001461e+03, rel_eval = 3.547903e-03
Iteration 22:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 4.693881e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.836774e+02, REL_OBJ, iter = 1
f(x^*) = 1.004762e+03, rel_eval = 3.285374e-03
Iteration 23:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 4.590612e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.847719e+02, REL_OBJ, iter = 1
f(x^*) = 1.007863e+03, rel_eval = 3.077432e-03
Iteration 24:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 4.492646e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.858261e+02, REL_OBJ, iter = 1
f(x^*) = 1.010827e+03, rel_eval = 2.931410e-03
Iteration 25:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 4.397880e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.867612e+02, REL_OBJ, iter = 1
f(x^*) = 1.013557e+03, rel_eval = 2.693733e-03
Iteration 26:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 4.307040e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.876447e+02, REL_OBJ, iter = 1
f(x^*) = 1.016077e+03, rel_eval = 2.480700e-03
Iteration 27:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 4.219982e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.885083e+02, REL_OBJ, iter = 1
f(x^*) = 1.018414e+03, rel_eval = 2.294669e-03
Iteration 28:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 4.136396e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.893154e+02, REL_OBJ, iter = 1
f(x^*) = 1.020659e+03, rel_eval = 2.199417e-03
Iteration 29:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 4.056908e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.900779e+02, REL_OBJ, iter = 1
f(x^*) = 1.022849e+03, rel_eval = 2.140885e-03
Iteration 30:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.981500e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.908159e+02, REL_OBJ, iter = 1
f(x^*) = 1.024976e+03, rel_eval = 2.075439e-03
Iteration 31:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.910135e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.915272e+02, REL_OBJ, iter = 1
f(x^*) = 1.027081e+03, rel_eval = 2.049271e-03
Iteration 32:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.842740e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.922120e+02, REL_OBJ, iter = 1
f(x^*) = 1.029160e+03, rel_eval = 2.020166e-03
Iteration 33:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.778363e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.928682e+02, REL_OBJ, iter = 1
f(x^*) = 1.031187e+03, rel_eval = 1.965163e-03
Iteration 34:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.717164e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.935258e+02, REL_OBJ, iter = 1
f(x^*) = 1.033157e+03, rel_eval = 1.907262e-03
Iteration 35:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.659033e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.941670e+02, REL_OBJ, iter = 1
f(x^*) = 1.035086e+03, rel_eval = 1.863156e-03
Iteration 36:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.603964e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.947968e+02, REL_OBJ, iter = 1
f(x^*) = 1.036971e+03, rel_eval = 1.817902e-03
Iteration 37:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.551840e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.954012e+02, REL_OBJ, iter = 1
f(x^*) = 1.038770e+03, rel_eval = 1.732471e-03
Iteration 38:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.502506e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.959829e+02, REL_OBJ, iter = 1
f(x^*) = 1.040539e+03, rel_eval = 1.700229e-03
Iteration 39:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.455567e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.965501e+02, REL_OBJ, iter = 1
f(x^*) = 1.042326e+03, rel_eval = 1.713881e-03
Iteration 40:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.410743e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.971127e+02, REL_OBJ, iter = 1
f(x^*) = 1.044085e+03, rel_eval = 1.685108e-03
Iteration 41:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.367711e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.976808e+02, REL_OBJ, iter = 1
f(x^*) = 1.045848e+03, rel_eval = 1.685520e-03
Iteration 42:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.325985e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.982534e+02, REL_OBJ, iter = 1
f(x^*) = 1.047586e+03, rel_eval = 1.659143e-03
Iteration 43:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.285342e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.988104e+02, REL_OBJ, iter = 1
f(x^*) = 1.049303e+03, rel_eval = 1.635742e-03
Iteration 44:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.245976e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.993614e+02, REL_OBJ, iter = 1
f(x^*) = 1.051028e+03, rel_eval = 1.641922e-03
Iteration 45:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.208151e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 2.999024e+02, REL_OBJ, iter = 1
f(x^*) = 1.052754e+03, rel_eval = 1.639158e-03
Iteration 46:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.171809e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.004281e+02, REL_OBJ, iter = 1
f(x^*) = 1.054463e+03, rel_eval = 1.620983e-03
Iteration 47:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.136990e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.009372e+02, REL_OBJ, iter = 1
f(x^*) = 1.056134e+03, rel_eval = 1.582330e-03
Iteration 48:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.103774e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.014246e+02, REL_OBJ, iter = 1
f(x^*) = 1.057772e+03, rel_eval = 1.548300e-03
Iteration 49:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.072356e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.019059e+02, REL_OBJ, iter = 1
f(x^*) = 1.059352e+03, rel_eval = 1.491333e-03
Iteration 50:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.042816e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.023638e+02, REL_OBJ, iter = 1
f(x^*) = 1.060868e+03, rel_eval = 1.429155e-03
Iteration 51:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.014872e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.027972e+02, REL_OBJ, iter = 1
f(x^*) = 1.062317e+03, rel_eval = 1.363850e-03
Iteration 52:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.988586e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.032096e+02, REL_OBJ, iter = 1
f(x^*) = 1.063716e+03, rel_eval = 1.314939e-03
Iteration 53:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.963901e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.036096e+02, REL_OBJ, iter = 1
f(x^*) = 1.065047e+03, rel_eval = 1.250045e-03
Iteration 54:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.940849e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.039974e+02, REL_OBJ, iter = 1
f(x^*) = 1.066310e+03, rel_eval = 1.184245e-03
Iteration 55:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.919386e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.043746e+02, REL_OBJ, iter = 1
f(x^*) = 1.067515e+03, rel_eval = 1.129434e-03
Iteration 56:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.899406e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.047460e+02, REL_OBJ, iter = 1
f(x^*) = 1.068679e+03, rel_eval = 1.088366e-03
Iteration 57:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.880908e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.050964e+02, REL_OBJ, iter = 1
f(x^*) = 1.069804e+03, rel_eval = 1.052409e-03
Iteration 58:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.863758e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.054265e+02, REL_OBJ, iter = 1
f(x^*) = 1.070874e+03, rel_eval = 9.992030e-04
Iteration 59:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.847968e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.057411e+02, REL_OBJ, iter = 1
f(x^*) = 1.071914e+03, rel_eval = 9.694856e-04
Iteration 60:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.833435e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.060349e+02, REL_OBJ, iter = 1
f(x^*) = 1.072909e+03, rel_eval = 9.281277e-04
Iteration 61:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.820100e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.063121e+02, REL_OBJ, iter = 1
f(x^*) = 1.073848e+03, rel_eval = 8.735134e-04
Iteration 62:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.807780e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.065721e+02, REL_OBJ, iter = 1
f(x^*) = 1.074733e+03, rel_eval = 8.243744e-04
Iteration 63:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.796466e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.068188e+02, REL_OBJ, iter = 1
f(x^*) = 1.075550e+03, rel_eval = 7.595745e-04
Iteration 64:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.786153e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.070595e+02, REL_OBJ, iter = 1
f(x^*) = 1.076311e+03, rel_eval = 7.065569e-04
Iteration 65:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.776745e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.072926e+02, REL_OBJ, iter = 1
f(x^*) = 1.077060e+03, rel_eval = 6.956587e-04
Iteration 66:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.768222e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.075086e+02, REL_OBJ, iter = 1
f(x^*) = 1.077795e+03, rel_eval = 6.813985e-04
Iteration 67:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.760559e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.077197e+02, REL_OBJ, iter = 1
f(x^*) = 1.078490e+03, rel_eval = 6.443683e-04
Iteration 68:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.753540e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.079262e+02, REL_OBJ, iter = 1
f(x^*) = 1.079152e+03, rel_eval = 6.137376e-04
Iteration 69:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.747081e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.081209e+02, REL_OBJ, iter = 1
f(x^*) = 1.079772e+03, rel_eval = 5.740881e-04
Iteration 70:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.741125e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.083035e+02, REL_OBJ, iter = 1
f(x^*) = 1.080391e+03, rel_eval = 5.732328e-04
Iteration 71:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.735713e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.084830e+02, REL_OBJ, iter = 1
f(x^*) = 1.081002e+03, rel_eval = 5.649146e-04
Iteration 72:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.730907e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.086532e+02, REL_OBJ, iter = 1
f(x^*) = 1.081576e+03, rel_eval = 5.311010e-04
Iteration 73:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.726752e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.088138e+02, REL_OBJ, iter = 1
f(x^*) = 1.082110e+03, rel_eval = 4.929664e-04
Iteration 74:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.723194e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.089764e+02, REL_OBJ, iter = 1
f(x^*) = 1.082622e+03, rel_eval = 4.734194e-04
Iteration 75:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.720023e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.091419e+02, REL_OBJ, iter = 1
f(x^*) = 1.083129e+03, rel_eval = 4.680980e-04
Iteration 76:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.717340e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.093027e+02, REL_OBJ, iter = 1
f(x^*) = 1.083626e+03, rel_eval = 4.585864e-04
Iteration 77:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.715127e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.094542e+02, REL_OBJ, iter = 1
f(x^*) = 1.084114e+03, rel_eval = 4.502109e-04
Iteration 78:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.713413e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.095959e+02, REL_OBJ, iter = 1
f(x^*) = 1.084588e+03, rel_eval = 4.369448e-04
Iteration 79:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.712182e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.097379e+02, REL_OBJ, iter = 1
f(x^*) = 1.085068e+03, rel_eval = 4.423292e-04
Iteration 80:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.711360e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.098891e+02, REL_OBJ, iter = 1
f(x^*) = 1.085561e+03, rel_eval = 4.540313e-04
Iteration 81:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.710934e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.100483e+02, REL_OBJ, iter = 1
f(x^*) = 1.086071e+03, rel_eval = 4.693298e-04
Iteration 82:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.710941e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.102121e+02, REL_OBJ, iter = 1
f(x^*) = 1.086597e+03, rel_eval = 4.840994e-04
Iteration 83:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.711408e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.103714e+02, REL_OBJ, iter = 1
f(x^*) = 1.087119e+03, rel_eval = 4.805566e-04
Iteration 84:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.712343e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.105270e+02, REL_OBJ, iter = 1
f(x^*) = 1.087628e+03, rel_eval = 4.681228e-04
Iteration 85:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.713752e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.106789e+02, REL_OBJ, iter = 1
f(x^*) = 1.088136e+03, rel_eval = 4.669574e-04
Iteration 86:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.715627e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.108350e+02, REL_OBJ, iter = 1
f(x^*) = 1.088645e+03, rel_eval = 4.673185e-04
Iteration 87:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.717891e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.109937e+02, REL_OBJ, iter = 1
f(x^*) = 1.089152e+03, rel_eval = 4.652244e-04
Iteration 88:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.720559e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.111564e+02, REL_OBJ, iter = 1
f(x^*) = 1.089646e+03, rel_eval = 4.538486e-04
Iteration 89:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.723618e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.113154e+02, REL_OBJ, iter = 1
f(x^*) = 1.090131e+03, rel_eval = 4.446646e-04
Iteration 90:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.727167e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.114683e+02, REL_OBJ, iter = 1
f(x^*) = 1.090611e+03, rel_eval = 4.401886e-04
Iteration 91:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.731237e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.116284e+02, REL_OBJ, iter = 1
f(x^*) = 1.091095e+03, rel_eval = 4.430777e-04
Iteration 92:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.735794e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.117875e+02, REL_OBJ, iter = 1
f(x^*) = 1.091576e+03, rel_eval = 4.414288e-04
Iteration 93:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.740855e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.119463e+02, REL_OBJ, iter = 1
f(x^*) = 1.092059e+03, rel_eval = 4.418929e-04
Iteration 94:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.746383e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.120995e+02, REL_OBJ, iter = 1
f(x^*) = 1.092529e+03, rel_eval = 4.305712e-04
Iteration 95:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.752317e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.122486e+02, REL_OBJ, iter = 1
f(x^*) = 1.092982e+03, rel_eval = 4.141157e-04
Iteration 96:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.758641e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.123913e+02, REL_OBJ, iter = 1
f(x^*) = 1.093422e+03, rel_eval = 4.020110e-04
Iteration 97:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.765339e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.125265e+02, REL_OBJ, iter = 1
f(x^*) = 1.093855e+03, rel_eval = 3.962919e-04
Iteration 98:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.772371e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.126628e+02, REL_OBJ, iter = 1
f(x^*) = 1.094284e+03, rel_eval = 3.924009e-04
Iteration 99:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.779603e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.127990e+02, REL_OBJ, iter = 1
f(x^*) = 1.094722e+03, rel_eval = 3.993577e-04
Iteration 100:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.786989e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.129377e+02, REL_OBJ, iter = 1
f(x^*) = 1.095168e+03, rel_eval = 4.078362e-04
Iteration 101:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.794560e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.130830e+02, REL_OBJ, iter = 1
f(x^*) = 1.095621e+03, rel_eval = 4.129309e-04
Iteration 102:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.802363e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.132359e+02, REL_OBJ, iter = 1
f(x^*) = 1.096079e+03, rel_eval = 4.182355e-04
Iteration 103:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.810369e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.133814e+02, REL_OBJ, iter = 1
f(x^*) = 1.096544e+03, rel_eval = 4.235545e-04
Iteration 104:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.818629e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.135293e+02, REL_OBJ, iter = 1
f(x^*) = 1.097034e+03, rel_eval = 4.473043e-04
Iteration 105:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.827174e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.136841e+02, REL_OBJ, iter = 1
f(x^*) = 1.097542e+03, rel_eval = 4.624209e-04
Iteration 106:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.836024e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.138474e+02, REL_OBJ, iter = 1
f(x^*) = 1.098068e+03, rel_eval = 4.788310e-04
Iteration 107:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.845167e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.140106e+02, REL_OBJ, iter = 1
f(x^*) = 1.098595e+03, rel_eval = 4.803585e-04
Iteration 108:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.854591e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.141754e+02, REL_OBJ, iter = 1
f(x^*) = 1.099134e+03, rel_eval = 4.896878e-04
Iteration 109:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.864230e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.143385e+02, REL_OBJ, iter = 1
f(x^*) = 1.099655e+03, rel_eval = 4.742899e-04
Iteration 110:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.874114e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.145055e+02, REL_OBJ, iter = 1
f(x^*) = 1.100161e+03, rel_eval = 4.596147e-04
Iteration 111:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.884225e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.146709e+02, REL_OBJ, iter = 1
f(x^*) = 1.100663e+03, rel_eval = 4.559480e-04
Iteration 112:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.894553e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.148328e+02, REL_OBJ, iter = 1
f(x^*) = 1.101159e+03, rel_eval = 4.507368e-04
Iteration 113:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.905098e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.149907e+02, REL_OBJ, iter = 1
f(x^*) = 1.101657e+03, rel_eval = 4.522471e-04
Iteration 114:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.915856e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.151507e+02, REL_OBJ, iter = 1
f(x^*) = 1.102150e+03, rel_eval = 4.472335e-04
Iteration 115:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.926857e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.153178e+02, REL_OBJ, iter = 1
f(x^*) = 1.102647e+03, rel_eval = 4.510178e-04
Iteration 116:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.938080e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.154898e+02, REL_OBJ, iter = 1
f(x^*) = 1.103158e+03, rel_eval = 4.630679e-04
Iteration 117:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.949488e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.156583e+02, REL_OBJ, iter = 1
f(x^*) = 1.103673e+03, rel_eval = 4.667866e-04
Iteration 118:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.961120e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.158245e+02, REL_OBJ, iter = 1
f(x^*) = 1.104192e+03, rel_eval = 4.694114e-04
Iteration 119:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.973008e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.159965e+02, REL_OBJ, iter = 1
f(x^*) = 1.104719e+03, rel_eval = 4.768344e-04
Iteration 120:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.985185e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.161647e+02, REL_OBJ, iter = 1
f(x^*) = 1.105233e+03, rel_eval = 4.658224e-04
Iteration 121:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 2.997644e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.163211e+02, REL_OBJ, iter = 1
f(x^*) = 1.105721e+03, rel_eval = 4.406402e-04
Iteration 122:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.010383e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.164770e+02, REL_OBJ, iter = 1
f(x^*) = 1.106201e+03, rel_eval = 4.343933e-04
Iteration 123:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.023408e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.166224e+02, REL_OBJ, iter = 1
f(x^*) = 1.106665e+03, rel_eval = 4.191867e-04
Iteration 124:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.036697e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.167668e+02, REL_OBJ, iter = 1
f(x^*) = 1.107128e+03, rel_eval = 4.177709e-04
Iteration 125:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.050215e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.169183e+02, REL_OBJ, iter = 1
f(x^*) = 1.107603e+03, rel_eval = 4.294940e-04
Iteration 126:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.063962e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.170775e+02, REL_OBJ, iter = 1
f(x^*) = 1.108094e+03, rel_eval = 4.426459e-04
Iteration 127:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.077923e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.172384e+02, REL_OBJ, iter = 1
f(x^*) = 1.108590e+03, rel_eval = 4.474775e-04
Iteration 128:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.092099e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.173911e+02, REL_OBJ, iter = 1
f(x^*) = 1.109082e+03, rel_eval = 4.435239e-04
Iteration 129:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.106492e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.175507e+02, REL_OBJ, iter = 1
f(x^*) = 1.109591e+03, rel_eval = 4.591230e-04
Iteration 130:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.121090e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.177126e+02, REL_OBJ, iter = 1
f(x^*) = 1.110118e+03, rel_eval = 4.742078e-04
Iteration 131:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.135902e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.178837e+02, REL_OBJ, iter = 1
f(x^*) = 1.110651e+03, rel_eval = 4.798475e-04
Iteration 132:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.150909e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.180552e+02, REL_OBJ, iter = 1
f(x^*) = 1.111183e+03, rel_eval = 4.793884e-04
Iteration 133:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.166132e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.182152e+02, REL_OBJ, iter = 1
f(x^*) = 1.111702e+03, rel_eval = 4.665690e-04
Iteration 134:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.181595e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.183723e+02, REL_OBJ, iter = 1
f(x^*) = 1.112206e+03, rel_eval = 4.531063e-04
Iteration 135:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.197312e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.185286e+02, REL_OBJ, iter = 1
f(x^*) = 1.112694e+03, rel_eval = 4.387973e-04
Iteration 136:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.213265e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.186780e+02, REL_OBJ, iter = 1
f(x^*) = 1.113167e+03, rel_eval = 4.249764e-04
Iteration 137:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.229463e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.188346e+02, REL_OBJ, iter = 1
f(x^*) = 1.113663e+03, rel_eval = 4.448002e-04
Iteration 138:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.245867e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.190009e+02, REL_OBJ, iter = 1
f(x^*) = 1.114180e+03, rel_eval = 4.643957e-04
Iteration 139:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.262451e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.191759e+02, REL_OBJ, iter = 1
f(x^*) = 1.114709e+03, rel_eval = 4.750322e-04
Iteration 140:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.279187e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.193459e+02, REL_OBJ, iter = 1
f(x^*) = 1.115242e+03, rel_eval = 4.771907e-04
Iteration 141:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.296070e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.195137e+02, REL_OBJ, iter = 1
f(x^*) = 1.115771e+03, rel_eval = 4.744129e-04
Iteration 142:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.313120e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.196784e+02, REL_OBJ, iter = 1
f(x^*) = 1.116294e+03, rel_eval = 4.685864e-04
Iteration 143:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.330372e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.198379e+02, REL_OBJ, iter = 1
f(x^*) = 1.116812e+03, rel_eval = 4.641282e-04
Iteration 144:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.347852e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.200009e+02, REL_OBJ, iter = 1
f(x^*) = 1.117338e+03, rel_eval = 4.707477e-04
Iteration 145:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.365555e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.201684e+02, REL_OBJ, iter = 1
f(x^*) = 1.117871e+03, rel_eval = 4.765365e-04
Iteration 146:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.383473e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.203310e+02, REL_OBJ, iter = 1
f(x^*) = 1.118397e+03, rel_eval = 4.706861e-04
Iteration 147:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.401619e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.204885e+02, REL_OBJ, iter = 1
f(x^*) = 1.118903e+03, rel_eval = 4.515237e-04
Iteration 148:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.420001e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.206431e+02, REL_OBJ, iter = 1
f(x^*) = 1.119400e+03, rel_eval = 4.442816e-04
Iteration 149:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.438628e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.208033e+02, REL_OBJ, iter = 1
f(x^*) = 1.119895e+03, rel_eval = 4.421084e-04
Iteration 150:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.457464e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.209679e+02, REL_OBJ, iter = 1
f(x^*) = 1.120396e+03, rel_eval = 4.469231e-04
Iteration 151:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.476504e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.211269e+02, REL_OBJ, iter = 1
f(x^*) = 1.120898e+03, rel_eval = 4.478364e-04
Iteration 152:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.495758e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.212818e+02, REL_OBJ, iter = 1
f(x^*) = 1.121395e+03, rel_eval = 4.432494e-04
Iteration 153:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.515225e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.214296e+02, REL_OBJ, iter = 1
f(x^*) = 1.121880e+03, rel_eval = 4.320643e-04
Iteration 154:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.534896e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.215776e+02, REL_OBJ, iter = 1
f(x^*) = 1.122356e+03, rel_eval = 4.248707e-04
Iteration 155:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.554767e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.217221e+02, REL_OBJ, iter = 1
f(x^*) = 1.122825e+03, rel_eval = 4.175482e-04
Iteration 156:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.574869e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.218689e+02, REL_OBJ, iter = 1
f(x^*) = 1.123293e+03, rel_eval = 4.163989e-04
Iteration 157:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.595197e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.220154e+02, REL_OBJ, iter = 1
f(x^*) = 1.123764e+03, rel_eval = 4.194498e-04
Iteration 158:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.615760e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.221656e+02, REL_OBJ, iter = 1
f(x^*) = 1.124240e+03, rel_eval = 4.231867e-04
Iteration 159:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.636527e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.223157e+02, REL_OBJ, iter = 1
f(x^*) = 1.124720e+03, rel_eval = 4.266961e-04
Iteration 160:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.657476e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.224597e+02, REL_OBJ, iter = 1
f(x^*) = 1.125199e+03, rel_eval = 4.255534e-04
Iteration 161:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.678597e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.226007e+02, REL_OBJ, iter = 1
f(x^*) = 1.125672e+03, rel_eval = 4.206485e-04
Iteration 162:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.699921e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.227482e+02, REL_OBJ, iter = 1
f(x^*) = 1.126145e+03, rel_eval = 4.196783e-04
Iteration 163:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.721446e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.229007e+02, REL_OBJ, iter = 1
f(x^*) = 1.126622e+03, rel_eval = 4.230159e-04
Iteration 164:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.743137e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.230540e+02, REL_OBJ, iter = 1
f(x^*) = 1.127102e+03, rel_eval = 4.260999e-04
Iteration 165:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.765028e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.232024e+02, REL_OBJ, iter = 1
f(x^*) = 1.127572e+03, rel_eval = 4.170626e-04
Iteration 166:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.787123e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.233500e+02, REL_OBJ, iter = 1
f(x^*) = 1.128036e+03, rel_eval = 4.110909e-04
Iteration 167:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.809422e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.234999e+02, REL_OBJ, iter = 1
f(x^*) = 1.128504e+03, rel_eval = 4.144024e-04
Iteration 168:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.831915e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.236479e+02, REL_OBJ, iter = 1
f(x^*) = 1.128980e+03, rel_eval = 4.217401e-04
Iteration 169:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.854589e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.237944e+02, REL_OBJ, iter = 1
f(x^*) = 1.129454e+03, rel_eval = 4.200296e-04
Iteration 170:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.877439e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.239352e+02, REL_OBJ, iter = 1
f(x^*) = 1.129918e+03, rel_eval = 4.105259e-04
Iteration 171:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.900480e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.240743e+02, REL_OBJ, iter = 1
f(x^*) = 1.130373e+03, rel_eval = 4.023609e-04
Iteration 172:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.923722e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.242121e+02, REL_OBJ, iter = 1
f(x^*) = 1.130831e+03, rel_eval = 4.049984e-04
Iteration 173:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.947150e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.243519e+02, REL_OBJ, iter = 1
f(x^*) = 1.131295e+03, rel_eval = 4.105889e-04
Iteration 174:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.970754e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.244917e+02, REL_OBJ, iter = 1
f(x^*) = 1.131764e+03, rel_eval = 4.141261e-04
Iteration 175:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 3.994554e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.246326e+02, REL_OBJ, iter = 1
f(x^*) = 1.132241e+03, rel_eval = 4.214173e-04
Iteration 176:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 4.018546e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.247727e+02, REL_OBJ, iter = 1
f(x^*) = 1.132726e+03, rel_eval = 4.285072e-04
Iteration 177:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 4.042720e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.249198e+02, REL_OBJ, iter = 1
f(x^*) = 1.133225e+03, rel_eval = 4.395268e-04
Iteration 178:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 4.067081e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.250717e+02, REL_OBJ, iter = 1
f(x^*) = 1.133734e+03, rel_eval = 4.489149e-04
Iteration 179:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 4.091621e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.252313e+02, REL_OBJ, iter = 1
f(x^*) = 1.134257e+03, rel_eval = 4.619448e-04
Iteration 180:
Proj. B2: epsilon = 0.000000e+00, ||y-Ax||_2 = 4.116332e-01, MAX_IT, iter = 100
prox_L1: ||A x-y||_1 = 3.253928e+02, REL_OBJ, iter = 1
f(x^*) = 1.134794e+03, rel_eval = 4.727586e-04
PPXA:
f(x^*) = 1.134794e+03, rel_eval = 4.727586e-04
180 iterations
Stopping criterion: MAX_IT
estimation_time =
875.1300
Reconstruction_MSE =
0.0064
Sources_MSE =
0.0314
References:
M. Golbabaee, S. Arberet, P. Vandergheynst, et al. Multichannel compressed sensing via source separation for hyperspectral images. In Eusipco 2010, 2010.