1 - x \cdot \left(0.253 + x \cdot 0.12\right)
1 - \frac{1}{\frac{1}{x \cdot \mathsf{fma}\left(0.12, x, 0.253\right)}}double code(double x) {
return (1.0 - (x * (0.253 + (x * 0.12))));
}
double code(double x) {
return (1.0 - (1.0 / (1.0 / (x * fma(0.12, x, 0.253)))));
}



Bits error versus x
Results
Initial program 0.1
rmApplied flip-+0.1
Applied associate-*r/7.4
Simplified7.3
rmApplied clear-num7.4
Simplified0.1
Final simplification0.1
herbie shell --seed 2020071 +o rules:numerics
(FPCore (x)
:name "Numeric.SpecFunctions:invIncompleteGamma from math-functions-0.1.5.2, A"
:precision binary64
(- 1 (* x (+ 0.253 (* x 0.12)))))