WIP/HACK: Save output using maximum range for our worst precision (float16)

Author: Entropy512

src/DngFloatWriter.cpp

@@ -263,11 +263,11 @@ void DngFloatWriter::createRawIFD() {
     uint16_t cblack[cfaRows * cfaCols];
     for (int row = 0; row < cfaRows; ++row) {
         for (int col = 0; col < cfaCols; ++col) {
-            cblack[row*cfaCols + col] = params->blackAt(col, row);
+            cblack[row*cfaCols + col] = 0;
         }
     }
     rawIFD.addEntry(BLACKLEVEL, IFD::SHORT, cfaRows * cfaCols, cblack);
-    rawIFD.addEntry(WHITELEVEL, IFD::SHORT, params->max);
+    rawIFD.addEntry(WHITELEVEL, IFD::SHORT, 65504); //Largest possible non-inf float16 value
     rawIFD.addEntry(SAMPLESPERPIXEL, IFD::SHORT, 1);
     rawIFD.addEntry(BITSPERSAMPLE, IFD::SHORT, bps);
     if (bps == 24) {

src/ImageStack.cpp

@@ -458,13 +458,13 @@ Array2D<float> ImageStack::compose(const RawParameters & params, int featherRadi
     }
 
     dst.displace(params.leftMargin, params.topMargin);
-    // Scale to params.max and recover the black levels
-    float mult = (params.max - params.maxBlack) / max;
+    // Scale to float16 limit
+    float mult = 65504.0 / max; //65504 is the largest non-inf float16
     #pragma omp parallel for
     for (size_t y = 0; y < params.rawHeight; ++y) {
         for (size_t x = 0; x < params.rawWidth; ++x) {
             dst(x, y) *= mult;
-            dst(x, y) += params.blackAt(x - params.leftMargin, y - params.topMargin);
+            //dst(x, y) += params.blackAt(x - params.leftMargin, y - params.topMargin);
         }
     }