\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2}0.5 \cdot \left(\left(y + \frac{{x}^{\left(\frac{2}{2}\right)}}{\frac{y}{x}}\right) - \left|z\right| \cdot \frac{\left|z\right|}{y}\right)double f(double x, double y, double z) {
double r624848 = x;
double r624849 = r624848 * r624848;
double r624850 = y;
double r624851 = r624850 * r624850;
double r624852 = r624849 + r624851;
double r624853 = z;
double r624854 = r624853 * r624853;
double r624855 = r624852 - r624854;
double r624856 = 2.0;
double r624857 = r624850 * r624856;
double r624858 = r624855 / r624857;
return r624858;
}
double f(double x, double y, double z) {
double r624859 = 0.5;
double r624860 = y;
double r624861 = x;
double r624862 = 2.0;
double r624863 = r624862 / r624862;
double r624864 = pow(r624861, r624863);
double r624865 = r624860 / r624861;
double r624866 = r624864 / r624865;
double r624867 = r624860 + r624866;
double r624868 = z;
double r624869 = fabs(r624868);
double r624870 = r624869 / r624860;
double r624871 = r624869 * r624870;
double r624872 = r624867 - r624871;
double r624873 = r624859 * r624872;
return r624873;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 28.5 |
|---|---|
| Target | 0.2 |
| Herbie | 0.1 |
Initial program 28.5
Taylor expanded around 0 12.7
Simplified12.7
rmApplied sqr-pow12.7
Applied associate-/l*6.9
Simplified6.9
rmApplied *-un-lft-identity6.9
Applied add-sqr-sqrt6.9
Applied times-frac6.9
Simplified6.9
Simplified0.1
Final simplification0.1
herbie shell --seed 2020020
(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)))