1 + \frac{4 \cdot \left(\left(x + y \cdot 0.25\right) - z\right)}{y}1 + \frac{4 \cdot \left(\left(x + y \cdot 0.25\right) - z\right)}{y}double code(double x, double y, double z) {
return (1.0 + ((4.0 * ((x + (y * 0.25)) - z)) / y));
}
double code(double x, double y, double z) {
return (1.0 + ((4.0 * ((x + (y * 0.25)) - z)) / y));
}



Bits error versus x



Bits error versus y



Bits error versus z
Results
Initial program 0.1
Final simplification0.1
herbie shell --seed 2020066 +o rules:numerics
(FPCore (x y z)
:name "Data.Array.Repa.Algorithms.ColorRamp:rampColorHotToCold from repa-algorithms-3.4.0.1, C"
:precision binary64
(+ 1 (/ (* 4 (- (+ x (* y 0.25)) z)) y)))