x + \frac{y - x}{z}
x - \frac{x - y}{z}
(FPCore (x y z) :precision binary64 (+ x (/ (- y x) z)))
(FPCore (x y z) :precision binary64 (- x (/ (- x y) z)))
double code(double x, double y, double z) {
return x + ((y - x) / z);
}
double code(double x, double y, double z) {
return x - ((x - y) / z);
}



Bits error versus x



Bits error versus y



Bits error versus z
Results
Initial program 0.0
Taylor expanded in y around 0 0.0
Taylor expanded in z around -inf 0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2022077
(FPCore (x y z)
:name "Statistics.Sample:$swelfordMean from math-functions-0.1.5.2"
:precision binary64
(+ x (/ (- y x) z)))