diff options
author | Eugeniy E. Mikhailov <evgmik@gmail.com> | 2021-01-29 16:29:19 -0500 |
---|---|---|
committer | Eugeniy E. Mikhailov <evgmik@gmail.com> | 2021-01-29 16:43:25 -0500 |
commit | 08d3232bdcbada58ffe9d20c6924b7e55aeeafda (patch) | |
tree | e3c94b2ed8bb95f79dd0f06ad09644f8a912ab43 | |
parent | 3983eb46023c1edd00617729ba929057fda8d0bd (diff) | |
download | l1magic-08d3232bdcbada58ffe9d20c6924b7e55aeeafda.tar.gz l1magic-08d3232bdcbada58ffe9d20c6924b7e55aeeafda.zip |
Small fix to make code Matlab 2020 compatible
variables and function names should be different
-rw-r--r-- | Optimization/tvdantzig_newton.m | 4 | ||||
-rw-r--r-- | Optimization/tvqc_newton.m | 4 |
2 files changed, 4 insertions, 4 deletions
diff --git a/Optimization/tvdantzig_newton.m b/Optimization/tvdantzig_newton.m index 68d7148..b60b173 100644 --- a/Optimization/tvdantzig_newton.m +++ b/Optimization/tvdantzig_newton.m @@ -88,7 +88,7 @@ while (~done) w11 = ntgx - Dh'*(Dhx.*(sig12./sig22).*ntgt) - Dv'*(Dvx.*(sig12./sig22).*ntgt); if (largescale) - h11pfun = @(w) H11p(w, A, At, Dh, Dv, Dhx, Dvx, sigb, ft, siga); + h11pfun = @(w) H11pFun(w, A, At, Dh, Dv, Dhx, Dvx, sigb, ft, siga); [dx, cgres, cgiter] = cgsolve(h11pfun, w11, cgtol, cgmaxiter, 0); if (cgres > 1/2) disp('Cannot solve system. Returning previous iterate. (See Section 4 of notes for more information.)'); @@ -173,7 +173,7 @@ end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % H11p auxiliary function -function y = H11p(v, A, At, Dh, Dv, Dhx, Dvx, sigb, ft, siga) +function y = H11pFun(v, A, At, Dh, Dv, Dhx, Dvx, sigb, ft, siga) Dhv = Dh*v; Dvv = Dv*v; diff --git a/Optimization/tvqc_newton.m b/Optimization/tvqc_newton.m index febe8ff..f7ca330 100644 --- a/Optimization/tvqc_newton.m +++ b/Optimization/tvqc_newton.m @@ -79,7 +79,7 @@ while (~done) w1p = ntgx - Dh'*(Dhx.*(sig12./sig22).*ntgt) - Dv'*(Dvx.*(sig12./sig22).*ntgt); if (largescale) - h11pfun = @(z) H11p(z, A, At, Dh, Dv, Dhx, Dvx, sigb, ft, fe, Atr); + h11pfun = @(z) H11pFun(z, A, At, Dh, Dv, Dhx, Dvx, sigb, ft, fe, Atr); [dx, cgres, cgiter] = cgsolve(h11pfun, w1p, cgtol, cgmaxiter, 0); if (cgres > 1/2) disp('Cannot solve system. Returning previous iterate. (See Section 4 of notes for more information.)'); @@ -164,7 +164,7 @@ end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % H11p auxiliary function -function y = H11p(v, A, At, Dh, Dv, Dhx, Dvx, sigb, ft, fe, atr) +function y = H11pFun(v, A, At, Dh, Dv, Dhx, Dvx, sigb, ft, fe, atr) Dhv = Dh*v; Dvv = Dv*v; |