diff options
| author | Edoardo Pasca <edo.paskino@gmail.com> | 2017-10-24 14:41:20 +0100 |
|---|---|---|
| committer | Edoardo Pasca <edo.paskino@gmail.com> | 2017-10-24 14:41:20 +0100 |
| commit | 67280dbdd5dc30b6500ba663d5f35dd8d5aa1389 (patch) | |
| tree | 3cfa6e49088dfbde13452127e879fda6122a1687 /main_func | |
| parent | ac4408e8984be8ca23a46b2b75bb243a0a4720aa (diff) | |
| parent | c91436873e48d531b9313f9c10fa5f89bcb90ab6 (diff) | |
| download | regularization-67280dbdd5dc30b6500ba663d5f35dd8d5aa1389.tar.gz regularization-67280dbdd5dc30b6500ba663d5f35dd8d5aa1389.tar.bz2 regularization-67280dbdd5dc30b6500ba663d5f35dd8d5aa1389.tar.xz regularization-67280dbdd5dc30b6500ba663d5f35dd8d5aa1389.zip | |
Merge branch 'master' of https://github.com/vais-ral/CCPi-FISTA_Reconstruction into pythonize
Diffstat (limited to 'main_func')
| -rw-r--r-- | main_func/FISTA_REC.m | 4 | ||||
| -rw-r--r-- | main_func/regularizers_CPU/FGP_TV.c | 6 |
2 files changed, 7 insertions, 3 deletions
diff --git a/main_func/FISTA_REC.m b/main_func/FISTA_REC.m index 1e4228d..3d22b97 100644 --- a/main_func/FISTA_REC.m +++ b/main_func/FISTA_REC.m @@ -685,6 +685,10 @@ else counterInd = counterInd + numProjSub; end + if (i == 1) + r_old = r; + end + % working with a 'ring vector' if (lambdaR_L1 > 0) r = max(abs(r)-lambdaR_L1, 0).*sign(r); % soft-thresholding operator for ring vector diff --git a/main_func/regularizers_CPU/FGP_TV.c b/main_func/regularizers_CPU/FGP_TV.c index 66442c9..30cea1a 100644 --- a/main_func/regularizers_CPU/FGP_TV.c +++ b/main_func/regularizers_CPU/FGP_TV.c @@ -66,7 +66,7 @@ void mexFunction( A = (float *) mxGetData(prhs[0]); /*noisy image (2D/3D) */ lambda = (float) mxGetScalar(prhs[1]); /* regularization parameter */ iter = 50; /* default iterations number */ - epsil = 0.001; /* default tolerance constant */ + epsil = 0.0001; /* default tolerance constant */ methTV = 0; /* default isotropic TV penalty */ if ((nrhs == 3) || (nrhs == 4) || (nrhs == 5)) iter = (int) mxGetScalar(prhs[2]); /* iterations number */ @@ -89,7 +89,7 @@ void mexFunction( tk = 1.0f; tkp1=1.0f; - count = 1; + count = 0; re_old = 0.0f; if (number_of_dims == 2) { @@ -128,7 +128,7 @@ void mexFunction( } re = sqrt(re)/sqrt(re1); if (re < epsil) count++; - if (count > 3) { + if (count > 4) { Obj_func_CALC2D(A, D, funcvalA, lambda, dimX, dimY); break; } |