diff --git a/rtengine/dual_demosaic_RT.cc b/rtengine/dual_demosaic_RT.cc
index b5839ee8b..0d9999e18 100644
--- a/rtengine/dual_demosaic_RT.cc
+++ b/rtengine/dual_demosaic_RT.cc
@@ -33,7 +33,7 @@
 
 #include "../rtgui/options.h"
 
-//#define BENCHMARK
+#define BENCHMARK
 #include "StopWatch.h"
 
 using namespace std;
@@ -66,14 +66,9 @@ void RawImageSource::dual_demosaic_RT(bool isBayer, const procparams::RAWParams
         return;
     }
 
-    array2D<float> redTmp(winw, winh);
-    array2D<float> greenTmp(winw, winh);
-    array2D<float> blueTmp(winw, winh);
     array2D<float> L(winw, winh);
 
     if (isBayer) {
-        vng4_demosaic(rawData, redTmp, greenTmp, blueTmp);
-
         if (raw.bayersensor.method == procparams::RAWParams::BayerSensor::getMethodString(procparams::RAWParams::BayerSensor::Method::AMAZEVNG4) || raw.bayersensor.method == procparams::RAWParams::BayerSensor::getMethodString(procparams::RAWParams::BayerSensor::Method::PIXELSHIFT)) {
             amaze_demosaic_RT(0, 0, winw, winh, rawData, red, green, blue, options.chunkSizeAMAZE, options.measure);
         } else if (raw.bayersensor.method == procparams::RAWParams::BayerSensor::getMethodString(procparams::RAWParams::BayerSensor::Method::DCBVNG4) ) {
@@ -87,7 +82,6 @@ void RawImageSource::dual_demosaic_RT(bool isBayer, const procparams::RAWParams
         } else {
             xtrans_interpolate (1, false, options.chunkSizeXT, options.measure);
         }
-        fast_xtrans_interpolate(rawData, redTmp, greenTmp, blueTmp);
     }
 
     const float xyz_rgb[3][3] = {          // XYZ from RGB
@@ -114,6 +108,17 @@ void RawImageSource::dual_demosaic_RT(bool isBayer, const procparams::RAWParams
     buildBlendMask(L, blend, winw, winh, contrastf, 1.f, autoContrast);
     contrast = contrastf * 100.f;
 
+    array2D<float>& redTmp = L; // L is not needed anymore => reuse it
+    array2D<float> greenTmp(winw, winh);
+    array2D<float> blueTmp(winw, winh);
+
+    if (isBayer) {
+        vng4_demosaic(rawData, redTmp, greenTmp, blueTmp);
+    } else {
+        fast_xtrans_interpolate(rawData, redTmp, greenTmp, blueTmp);
+    }
+
+
     // the following is split into 3 loops intentionally to avoid cache conflicts on CPUs with only 4-way cache
 #ifdef _OPENMP
     #pragma omp parallel for