disable non-IEEE 754 math

build.sh

@@ -413,7 +413,7 @@ all_args="$all_args -Wformat=2"
 all_args="$all_args -Wno-format-nonliteral"
 all_args="$all_args -Wno-parentheses"
 all_args="$all_args -fvisibility=hidden -fvisibility-inlines-hidden"
-all_args="$all_args -fno-math-errno -fno-trapping-math"
+all_args="$all_args -fno-math-errno -fno-trapping-math -fno-fast-math -ffp-contract=off"
 # TODO: once we have highly efficient tightly looped code, try no -fpic and see if that makes better code.  The compiler can save a register in this case. See https://akkadia.org/drepper/dsohowto.pdf
 # TODO: check no-plt compiler option
 all_args="$all_args -fpic"

shared/libebm/interpretable_numerics.cpp

@@ -1379,6 +1379,9 @@ static double Mean(const size_t cSamples,
    // https://stackoverflow.com/questions/895929/how-do-i-determine-the-standard-deviation-stddev-of-a-set-of-values
    // https://www.johndcook.com/blog/standard_deviation/
 
+   // do not put multiple floating point operations in the same statement since that can be optimized
+   // https://clang.llvm.org/docs/UsersManual.html
+
    double factor = 1.0;
    double mean;
    size_t cNaN;
@@ -1440,7 +1443,8 @@ static double Mean(const size_t cSamples,
                // if all the weights are zero, then weigh them all equally
                ratio = double{1} / static_cast<double>(cNormal);
             }
-            mean += numerator * ratio;
+            const double multiple = numerator * ratio;
+            mean += multiple;
          }
          if(nullptr != pWeight) {
             ++pWeight;