x - \frac{y}{1 + \frac{x \cdot y}{2}}x - \frac{y}{1 + \frac{x \cdot y}{2}}double f(double x, double y) {
double r256274 = x;
double r256275 = y;
double r256276 = 1.0;
double r256277 = r256274 * r256275;
double r256278 = 2.0;
double r256279 = r256277 / r256278;
double r256280 = r256276 + r256279;
double r256281 = r256275 / r256280;
double r256282 = r256274 - r256281;
return r256282;
}
double f(double x, double y) {
double r256283 = x;
double r256284 = y;
double r256285 = 1.0;
double r256286 = r256283 * r256284;
double r256287 = 2.0;
double r256288 = r256286 / r256287;
double r256289 = r256285 + r256288;
double r256290 = r256284 / r256289;
double r256291 = r256283 - r256290;
return r256291;
}



Bits error versus x



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