\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2}\frac{y + \frac{x + z}{\frac{y}{x - z}}}{2}double f(double x, double y, double z) {
double r392974 = x;
double r392975 = r392974 * r392974;
double r392976 = y;
double r392977 = r392976 * r392976;
double r392978 = r392975 + r392977;
double r392979 = z;
double r392980 = r392979 * r392979;
double r392981 = r392978 - r392980;
double r392982 = 2.0;
double r392983 = r392976 * r392982;
double r392984 = r392981 / r392983;
return r392984;
}
double f(double x, double y, double z) {
double r392985 = y;
double r392986 = x;
double r392987 = z;
double r392988 = r392986 + r392987;
double r392989 = r392986 - r392987;
double r392990 = r392985 / r392989;
double r392991 = r392988 / r392990;
double r392992 = r392985 + r392991;
double r392993 = 2.0;
double r392994 = r392992 / r392993;
return r392994;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 27.8 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 27.8
Simplified12.6
rmApplied difference-of-squares12.6
Applied associate-/l*0.2
Final simplification0.2
herbie shell --seed 2019208 +o rules:numerics
(FPCore (x y z)
:name "Diagrams.TwoD.Apollonian:initialConfig from diagrams-contrib-1.3.0.5, A"
:precision binary64
:herbie-target
(- (* y 0.5) (* (* (/ 0.5 y) (+ z x)) (- z x)))
(/ (- (+ (* x x) (* y y)) (* z z)) (* y 2)))