Mercurial > hg > absrec

comparison ABSmixed.py @ 21:d395461b92ae tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Lots and lots of modifications. Approximate recovery script working.
author Nic Cleju <nikcleju@gmail.com>
date Mon, 23 Apr 2012 10:54:57 +0300
parents 7fdf964f4edd
children
comparison
equal deleted inserted replaced
20:eccc7a5b9ee3 21:d395461b92ae
36 D = numpy.linalg.pinv(Omega) 36 D = numpy.linalg.pinv(Omega)
37 U,S,Vt = numpy.linalg.svd(D) 37 U,S,Vt = numpy.linalg.svd(D)
38 Aeps = numpy.dot(M,D) 38 Aeps = numpy.dot(M,D)
39 Aexact = Vt[-(N-n):,:] 39 Aexact = Vt[-(N-n):,:]
40 40
41 return numpy.dot(D, pyCSalgos.SL0.SL0_approx.SL0_approx_analysis(Aeps,Aexact,y,epsilon,sigma_min,sigma_decrease_factor,mu_0,L,Aeps_pinv,Aexact_pinv,true_s)) 41 #return numpy.dot(D, pyCSalgos.SL0.SL0_approx.SL0_approx_analysis(Aeps,Aexact,y,epsilon,sigma_min,sigma_decrease_factor,mu_0,L,Aeps_pinv,Aexact_pinv,true_s))
42 #return numpy.dot(D, pyCSalgos.SL0.SL0_approx.SL0_robust_analysis(Aeps,Aexact,y,epsilon,sigma_min,sigma_decrease_factor,mu_0,L,Aeps_pinv,Aexact_pinv,true_s))
43 #return numpy.dot(D, pyCSalgos.SL0.SL0_approx.SL0_approx_analysis_unconstrained(Aeps,Aexact,y,epsilon,sigma_min,sigma_decrease_factor,mu_0,L,Aeps_pinv,Aexact_pinv,true_s))
44 return numpy.dot(D, pyCSalgos.SL0.SL0_approx.SL0_approx_analysis_dai(Aeps,Aexact,y,epsilon,sigma_min,sigma_decrease_factor,mu_0,L,Aeps_pinv,Aexact_pinv,true_s))
42 45