x - \frac{y}{1 + \frac{x \cdot y}{2}}x - \frac{y}{1 + \frac{x \cdot y}{2}}(FPCore (x y) :precision binary64 (- x (/ y (+ 1.0 (/ (* x y) 2.0)))))
(FPCore (x y) :precision binary64 (- x (/ y (+ 1.0 (/ (* x y) 2.0)))))
double code(double x, double y) {
return ((double) (x - (y / ((double) (1.0 + (((double) (x * y)) / 2.0))))));
}
double code(double x, double y) {
return ((double) (x - (y / ((double) (1.0 + (((double) (x * y)) / 2.0))))));
}



Bits error versus x



Bits error versus y
Results
Initial program 0.0
Final simplification0.0
herbie shell --seed 2020198
(FPCore (x y)
:name "Data.Number.Erf:$cinvnormcdf from erf-2.0.0.0, B"
:precision binary64
(- x (/ y (+ 1.0 (/ (* x y) 2.0)))))