diff --git a/src/DehazingCE.cpp b/src/DehazingCE.cpp
index bcd6c1f..37e328e 100644
--- a/src/DehazingCE.cpp
+++ b/src/DehazingCE.cpp
@@ -5,7 +5,7 @@
 
 constexpr float SQRT_3 = 1.733f;
 
-dehazing::dehazing(int nW, int nH, int nBits, int nABlockSize, int nTBlockSize, float fTransInit, bool bPrevFlag, bool bPosFlag, float fL1, float fL2, int nGBlockSize)
+dehazing::dehazing(int nW, int nH, int nBits, int nABlockSize, int nTBlockSize, float fTransInit, bool bPrevFlag, bool bPosFlag, double dL1, float fL2, int nGBlockSize)
 {
     width = nW;
     height = nH;
@@ -17,7 +17,7 @@ dehazing::dehazing(int nW, int nH, int nBits, int nABlockSize, int nTBlockSize,
     m_PostFlag = bPosFlag;
 
     // Parameters for each cost (loss cost, temporal coherence cost)
-    Lambda1 = fL1;
+    Lambda1 = dL1;
     Lambda2 = fL2;  // only used in previous mode
 
     // Block size for transmission estimation
@@ -220,7 +220,7 @@ float dehazing::NFTrsEstimationColor(const T* pnImageR, const T* pnImageG, const
 {
     int nOutR, nOutG, nOutB;
     float fOptTrs;
-    float fCost, fMinCost, fMean;
+    double dCost, dMinCost, dMean;
 
     int nEndX = std::min(nStartX + TBlockSize, ref_width);
     int nEndY = std::min(nStartY + TBlockSize, ref_height);
@@ -232,10 +232,10 @@ float dehazing::NFTrsEstimationColor(const T* pnImageR, const T* pnImageG, const
 
     for (auto nCounter = 0; nCounter < bits - 1; nCounter++)
     {
-        int nSumofSLoss = 0;
+        long long int nSumofSLoss = 0;
         int nLossCount = 0;
-        int nSumofSquaredOuts = 0;
-        int nSumofOuts = 0;
+        long long int nSumofSquaredOuts = 0;
+        long long int nSumofOuts = 0;
 
         int half_peak = ((peak + 1) >> 1);
         for (auto y = nStartY; y < nEndY; y++)
@@ -282,13 +282,13 @@ float dehazing::NFTrsEstimationColor(const T* pnImageR, const T* pnImageG, const
             }
         }
 
-        fMean = (float)(nSumofOuts) / (float)(nNumberofPixels);
-        fCost = Lambda1 * (float)nSumofSLoss / (float)(nNumberofPixels)
-                -((float)nSumofSquaredOuts / (float)nNumberofPixels - fMean * fMean);
+        dMean = (double)nSumofOuts / nNumberofPixels;
+        dCost = Lambda1 * (double)nSumofSLoss / nNumberofPixels
+                -((double)nSumofSquaredOuts / nNumberofPixels - dMean * dMean);
 
-        if (nCounter == 0 || fMinCost > fCost)
+        if (nCounter == 0 || dMinCost > dCost)
         {
-            fMinCost = fCost;
+            dMinCost = dCost;
             fOptTrs = fTrans;
         }
 
diff --git a/src/DehazingCE.h b/src/DehazingCE.h
index 2abcb12..c88c98a 100644
--- a/src/DehazingCE.h
+++ b/src/DehazingCE.h
@@ -10,7 +10,7 @@
 class dehazing
 {
 public:
-    dehazing(int nW, int nH, int nBits, int nABlockSize, int nTBlockSize, float fTransInit, bool bPrevFlag, bool bPosFlag, float fL1, float fL2, int nGBlockSize);
+    dehazing(int nW, int nH, int nBits, int nABlockSize, int nTBlockSize, float fTransInit, bool bPrevFlag, bool bPosFlag, double dL1, float fL2, int nGBlockSize);
 
     ~dehazing();
 
@@ -68,7 +68,7 @@ class dehazing
     bool m_PreviousFlag;
     bool m_PostFlag;  // Flag for post processing(deblocking)
 
-    float Lambda1;
+    double Lambda1;
     float Lambda2;
 
     int TopLeftX;
diff --git a/src/main.cpp b/src/main.cpp
index 7e1490f..235f680 100644
--- a/src/main.cpp
+++ b/src/main.cpp
@@ -185,9 +185,9 @@ static void VS_CC filterCreate(const VSMap* in, VSMap* out, void* userData, VSCo
         if (err)
             TransInit = 0.3f;
 
-        float lamdaA = (float)(vsapi->propGetFloat(in, "lamda", 0, &err));
+        double lamdaA = vsapi->propGetFloat(in, "lamda", 0, &err);
         if (err)
-            lamdaA = 5.f;
+            lamdaA = 5.0;
 
 		float gamma = (float)(vsapi->propGetFloat(in, "gamma", 0, &err));
         if (err)