\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2}0.5 \cdot \left(\left(y + \frac{\frac{{x}^{1}}{y}}{\frac{1}{x}}\right) - \frac{z}{\frac{y}{z}}\right)double f(double x, double y, double z) {
double r638078 = x;
double r638079 = r638078 * r638078;
double r638080 = y;
double r638081 = r638080 * r638080;
double r638082 = r638079 + r638081;
double r638083 = z;
double r638084 = r638083 * r638083;
double r638085 = r638082 - r638084;
double r638086 = 2.0;
double r638087 = r638080 * r638086;
double r638088 = r638085 / r638087;
return r638088;
}
double f(double x, double y, double z) {
double r638089 = 0.5;
double r638090 = y;
double r638091 = x;
double r638092 = 1.0;
double r638093 = pow(r638091, r638092);
double r638094 = r638093 / r638090;
double r638095 = r638092 / r638091;
double r638096 = r638094 / r638095;
double r638097 = r638090 + r638096;
double r638098 = z;
double r638099 = r638090 / r638098;
double r638100 = r638098 / r638099;
double r638101 = r638097 - r638100;
double r638102 = r638089 * r638101;
return r638102;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 28.7 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 28.7
Taylor expanded around 0 12.3
Simplified12.3
rmApplied unpow212.3
Applied associate-/l*6.7
rmApplied sqr-pow6.7
Applied associate-/l*0.2
Simplified0.2
rmApplied div-inv0.2
Applied associate-/r*0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2020089 +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)))