\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2}0.5 \cdot \left(\left(y + \left|x\right| \cdot \frac{\left|x\right|}{y}\right) - \left|z\right| \cdot \frac{\left|z\right|}{y}\right)double f(double x, double y, double z) {
double r652523 = x;
double r652524 = r652523 * r652523;
double r652525 = y;
double r652526 = r652525 * r652525;
double r652527 = r652524 + r652526;
double r652528 = z;
double r652529 = r652528 * r652528;
double r652530 = r652527 - r652529;
double r652531 = 2.0;
double r652532 = r652525 * r652531;
double r652533 = r652530 / r652532;
return r652533;
}
double f(double x, double y, double z) {
double r652534 = 0.5;
double r652535 = y;
double r652536 = x;
double r652537 = fabs(r652536);
double r652538 = r652537 / r652535;
double r652539 = r652537 * r652538;
double r652540 = r652535 + r652539;
double r652541 = z;
double r652542 = fabs(r652541);
double r652543 = r652542 / r652535;
double r652544 = r652542 * r652543;
double r652545 = r652540 - r652544;
double r652546 = r652534 * r652545;
return r652546;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 28.9 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 28.9
Taylor expanded around 0 12.9
Simplified12.9
rmApplied *-un-lft-identity12.9
Applied add-sqr-sqrt12.9
Applied times-frac12.9
Simplified12.9
Simplified7.1
rmApplied *-un-lft-identity7.1
Applied add-sqr-sqrt7.1
Applied times-frac7.1
Simplified7.1
Simplified0.2
Final simplification0.2
herbie shell --seed 2020046 +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)))