 norm of the gradient.</p><pre class="codeinput">[M,M0,output] = cp_opt(X, R, <span class="string">'init'</span>, M_init);
</pre><pre class="codeoutput">Iter    10, f(x) = 1.239376e+03, ||grad||_infty = 9.52e+01
Iter    20, f(x) = 1.204254e+03, ||grad||_infty = 6.93e-01
Iter    30, f(x) = 1.033479e+03, ||grad||_infty = 9.20e+01
Iter    40, f(x) = 5.932041e+02, ||grad||_infty = 8.66e+00
Iter    50, f(x) = 5.891379e+02, ||grad||_infty = 2.58e-01
Iter    60, f(x) = 5.891340e+02, ||grad||_infty = 7.49e-03
Iter 