\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 r646283 = x;
double r646284 = r646283 * r646283;
double r646285 = y;
double r646286 = r646285 * r646285;
double r646287 = r646284 + r646286;
double r646288 = z;
double r646289 = r646288 * r646288;
double r646290 = r646287 - r646289;
double r646291 = 2.0;
double r646292 = r646285 * r646291;
double r646293 = r646290 / r646292;
return r646293;
}
double f(double x, double y, double z) {
double r646294 = 0.5;
double r646295 = y;
double r646296 = x;
double r646297 = fabs(r646296);
double r646298 = r646297 / r646295;
double r646299 = r646297 * r646298;
double r646300 = r646295 + r646299;
double r646301 = z;
double r646302 = fabs(r646301);
double r646303 = r646302 / r646295;
double r646304 = r646302 * r646303;
double r646305 = r646300 - r646304;
double r646306 = r646294 * r646305;
return r646306;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 28.8 |
|---|---|
| Target | 0.2 |
| Herbie | 0.1 |
Initial program 28.8
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.0
rmApplied *-un-lft-identity7.0
Applied add-sqr-sqrt7.0
Applied times-frac7.0
Simplified7.0
Simplified0.1
Final simplification0.1
herbie shell --seed 2020081
(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)))