\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2}0.5 \cdot \left(\left(y + \frac{x}{\frac{y}{x}}\right) - \frac{z}{\frac{y}{z}}\right)double f(double x, double y, double z) {
double r644646 = x;
double r644647 = r644646 * r644646;
double r644648 = y;
double r644649 = r644648 * r644648;
double r644650 = r644647 + r644649;
double r644651 = z;
double r644652 = r644651 * r644651;
double r644653 = r644650 - r644652;
double r644654 = 2.0;
double r644655 = r644648 * r644654;
double r644656 = r644653 / r644655;
return r644656;
}
double f(double x, double y, double z) {
double r644657 = 0.5;
double r644658 = y;
double r644659 = x;
double r644660 = r644658 / r644659;
double r644661 = r644659 / r644660;
double r644662 = r644658 + r644661;
double r644663 = z;
double r644664 = r644658 / r644663;
double r644665 = r644663 / r644664;
double r644666 = r644662 - r644665;
double r644667 = r644657 * r644666;
return r644667;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 29.3 |
|---|---|
| Target | 0.2 |
| Herbie | 0.1 |
Initial program 29.3
Taylor expanded around 0 12.8
Simplified12.8
rmApplied unpow212.8
Applied associate-/l*6.9
rmApplied unpow26.9
Applied associate-/l*0.1
Final simplification0.1
herbie shell --seed 2020036 +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)))