x + \left(\tan \left(y + z\right) - \tan a\right)
\left(\mathsf{fma}\left(\mathsf{fma}\left(\sin y, \cos z, \cos y \cdot \sin z\right), \frac{\frac{1}{\cos z \cdot \cos y}}{1 - \tan y \cdot \tan z}, -\tan a\right) + x\right) + \mathsf{fma}\left(-\sqrt[3]{\tan a}, \sqrt[3]{\tan a} \cdot \sqrt[3]{\tan a}, \sqrt[3]{\tan a} \cdot \left(\sqrt[3]{\tan a} \cdot \sqrt[3]{\tan a}\right)\right)double f(double x, double y, double z, double a) {
double r106739 = x;
double r106740 = y;
double r106741 = z;
double r106742 = r106740 + r106741;
double r106743 = tan(r106742);
double r106744 = a;
double r106745 = tan(r106744);
double r106746 = r106743 - r106745;
double r106747 = r106739 + r106746;
return r106747;
}
double f(double x, double y, double z, double a) {
double r106748 = y;
double r106749 = sin(r106748);
double r106750 = z;
double r106751 = cos(r106750);
double r106752 = cos(r106748);
double r106753 = sin(r106750);
double r106754 = r106752 * r106753;
double r106755 = fma(r106749, r106751, r106754);
double r106756 = 1.0;
double r106757 = r106751 * r106752;
double r106758 = r106756 / r106757;
double r106759 = tan(r106748);
double r106760 = tan(r106750);
double r106761 = r106759 * r106760;
double r106762 = r106756 - r106761;
double r106763 = r106758 / r106762;
double r106764 = a;
double r106765 = tan(r106764);
double r106766 = -r106765;
double r106767 = fma(r106755, r106763, r106766);
double r106768 = x;
double r106769 = r106767 + r106768;
double r106770 = cbrt(r106765);
double r106771 = -r106770;
double r106772 = r106770 * r106770;
double r106773 = r106770 * r106772;
double r106774 = fma(r106771, r106772, r106773);
double r106775 = r106769 + r106774;
return r106775;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus a
Initial program 13.2
rmApplied tan-sum0.2
rmApplied tan-quot0.2
Applied tan-quot0.2
Applied frac-add0.2
Simplified0.2
Simplified0.2
rmApplied add-cube-cbrt0.3
Applied *-un-lft-identity0.3
Applied div-inv0.3
Applied times-frac0.3
Applied prod-diff0.3
Applied associate-+r+0.3
Simplified0.2
Final simplification0.2
herbie shell --seed 2019209 +o rules:numerics
(FPCore (x y z a)
:name "(+ x (- (tan (+ y z)) (tan a)))"
:precision binary64
:pre (and (or (== x 0.0) (<= 0.588414199999999998 x 505.590899999999976)) (or (<= -1.79665800000000009e308 y -9.425585000000013e-310) (<= 1.284938e-309 y 1.75122399999999993e308)) (or (<= -1.776707e308 z -8.59979600000002e-310) (<= 3.29314499999998e-311 z 1.72515400000000009e308)) (or (<= -1.79665800000000009e308 a -9.425585000000013e-310) (<= 1.284938e-309 a 1.75122399999999993e308)))
(+ x (- (tan (+ y z)) (tan a))))