summaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
authorDaniil Kazantsev <dkazanc@hotmail.com>2018-07-23 09:52:43 +0100
committerDaniil Kazantsev <dkazanc@hotmail.com>2018-07-23 09:52:43 +0100
commit06695c2e1afceef81395dc56c7bf01ff98c486bb (patch)
tree00c227de88936e55713ed586b0678b7422cb508b
parent6d01bb2d3236698efaa824d09d10064e9cbf8ca1 (diff)
downloadregularization-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.c4
-rw-r--r--Core/regularisers_GPU/LLT_ROF_GPU_core.cu4
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]));
}
}