From 2d5b96a88069960dd250f7ccc968f55dd07ea855 Mon Sep 17 00:00:00 2001
From: Oliver Duis <oduis@hotmail.com>
Date: Sun, 31 Jul 2011 22:53:40 +0200
Subject: [PATCH] Fixed RL deconvolution problem with negative lightness
 channel see issue 868

---
 rtengine/ipsharpen.cc | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/rtengine/ipsharpen.cc b/rtengine/ipsharpen.cc
index f6431a634..f877f31a0 100644
--- a/rtengine/ipsharpen.cc
+++ b/rtengine/ipsharpen.cc
@@ -39,19 +39,19 @@ namespace rtengine {
 extern Settings* settings;
 
 void ImProcFunctions::dcdamping (float** aI, float** aO, float damping, int W, int H) {
-
+    const float dampingFac=2.0/(damping*damping);
 #ifdef _OPENMP
 #pragma omp for
 #endif
 	for (int i=0; i<H; i++)
         for (int j=0; j<W; j++) {
             float I = aI[i][j];
-            float O = (float)aO[i][j];
-            if (O==0.0 || I==0.0) {
+            float O = aO[i][j];
+            if (O<=0.0 || I<=0.0) {
                 aI[i][j] = 0.0;
                 continue;
             }
-            float U = -(O * log(I/O) - I + O) * 2.0 / (damping*damping);
+            float U = -(O * log(I/O) - I + O) * dampingFac;
             U = MIN(U,1.0);
             U = U*U*U*U*(5.0-U*4.0);
             aI[i][j] = (O - I) / I * U + 1.0;