diff options
author | Daniil Kazantsev <dkazanc@hotmail.com> | 2018-07-23 09:52:43 +0100 |
---|---|---|
committer | Daniil Kazantsev <dkazanc@hotmail.com> | 2018-07-23 09:52:43 +0100 |
commit | 06695c2e1afceef81395dc56c7bf01ff98c486bb (patch) | |
tree | 00c227de88936e55713ed586b0678b7422cb508b | |
parent | 6d01bb2d3236698efaa824d09d10064e9cbf8ca1 (diff) | |
download | regularization-06695c2e1afceef81395dc56c7bf01ff98c486bb.tar.gz regularization-06695c2e1afceef81395dc56c7bf01ff98c486bb.tar.bz2 regularization-06695c2e1afceef81395dc56c7bf01ff98c486bb.tar.xz regularization-06695c2e1afceef81395dc56c7bf01ff98c486bb.zip |
minor fix of regul parameter for LLT
-rw-r--r-- | Core/regularisers_CPU/LLT_ROF_core.c | 4 | ||||
-rw-r--r-- | Core/regularisers_GPU/LLT_ROF_GPU_core.cu | 4 |
2 files changed, 4 insertions, 4 deletions
diff --git a/Core/regularisers_CPU/LLT_ROF_core.c b/Core/regularisers_CPU/LLT_ROF_core.c index 6dcf018..650c0af 100644 --- a/Core/regularisers_CPU/LLT_ROF_core.c +++ b/Core/regularisers_CPU/LLT_ROF_core.c @@ -363,7 +363,7 @@ float Update2D_LLT_ROF(float *U0, float *U, float *D1_LLT, float *D2_LLT, float div = dv1 + dv2; /*build Divirgent*/ /*combine all into one cost function to minimise */ - U[index] += tau*(lambdaROF*(div) - lambdaLLT*(laplc) - (U[index] - U0[index])); + U[index] += tau*(2.0f*lambdaROF*(div) - lambdaLLT*(laplc) - (U[index] - U0[index])); } } return *U; @@ -400,7 +400,7 @@ float Update3D_LLT_ROF(float *U0, float *U, float *D1_LLT, float *D2_LLT, float div = dv1 + dv2 + dv3; /*build Divirgent*/ /*combine all into one cost function to minimise */ - U[index] += tau*(lambdaROF*(div) - lambdaLLT*(laplc) - (U[index] - U0[index])); + U[index] += tau*(2.0f*lambdaROF*(div) - lambdaLLT*(laplc) - (U[index] - U0[index])); } } } diff --git a/Core/regularisers_GPU/LLT_ROF_GPU_core.cu b/Core/regularisers_GPU/LLT_ROF_GPU_core.cu index 70c9295..e71ca12 100644 --- a/Core/regularisers_GPU/LLT_ROF_GPU_core.cu +++ b/Core/regularisers_GPU/LLT_ROF_GPU_core.cu @@ -357,7 +357,7 @@ __global__ void Update2D_LLT_ROF_kernel(float *U0, float *U, float *D1_LLT, floa div = dv1 + dv2; /*build Divirgent*/ /*combine all into one cost function to minimise */ - U[index] += tau*(lambdaROF*(div) - lambdaLLT*(laplc) - (U[index] - U0[index])); + U[index] += tau*(2.0f*lambdaROF*(div) - lambdaLLT*(laplc) - (U[index] - U0[index])); } } @@ -395,7 +395,7 @@ __global__ void Update3D_LLT_ROF_kernel(float *U0, float *U, float *D1_LLT, floa div = dv1 + dv2 + dv3; /*build Divirgent*/ /*combine all into one cost function to minimise */ - U[index] += tau*(lambdaROF*(div) - lambdaLLT*(laplc) - (U[index] - U0[index])); + U[index] += tau*(2.0f*lambdaROF*(div) - lambdaLLT*(laplc) - (U[index] - U0[index])); } } |