x + \left(\tan \left(y + z\right) - \tan a\right)
\left(0 \cdot \tan a + \left(\mathsf{fma}\left(\mathsf{fma}\left(\tan z \cdot \tan y, \tan z \cdot \tan y, \tan z \cdot \tan y\right), \frac{\tan y + \tan z}{1 - \left(\left(\tan z \cdot \tan y\right) \cdot \left(\tan z \cdot \tan y\right)\right) \cdot \left(\tan z \cdot \tan y\right)}, \frac{\tan y + \tan z}{1 - \sqrt[3]{\left(\left(\left(\tan z \cdot \tan y\right) \cdot \left(\tan z \cdot \tan y\right)\right) \cdot \left(\tan z \cdot \tan y\right)\right) \cdot \left(\left(\left(\tan z \cdot \tan y\right) \cdot \left(\tan z \cdot \tan y\right)\right) \cdot \left(\tan z \cdot \tan y\right)\right)} \cdot \left(\tan z \cdot \tan y\right)}\right) - \tan a\right)\right) + xdouble f(double x, double y, double z, double a) {
double r4729652 = x;
double r4729653 = y;
double r4729654 = z;
double r4729655 = r4729653 + r4729654;
double r4729656 = tan(r4729655);
double r4729657 = a;
double r4729658 = tan(r4729657);
double r4729659 = r4729656 - r4729658;
double r4729660 = r4729652 + r4729659;
return r4729660;
}
double f(double x, double y, double z, double a) {
double r4729661 = 0.0;
double r4729662 = a;
double r4729663 = tan(r4729662);
double r4729664 = r4729661 * r4729663;
double r4729665 = z;
double r4729666 = tan(r4729665);
double r4729667 = y;
double r4729668 = tan(r4729667);
double r4729669 = r4729666 * r4729668;
double r4729670 = fma(r4729669, r4729669, r4729669);
double r4729671 = r4729668 + r4729666;
double r4729672 = 1.0;
double r4729673 = r4729669 * r4729669;
double r4729674 = r4729673 * r4729669;
double r4729675 = r4729672 - r4729674;
double r4729676 = r4729671 / r4729675;
double r4729677 = r4729674 * r4729674;
double r4729678 = cbrt(r4729677);
double r4729679 = r4729678 * r4729669;
double r4729680 = r4729672 - r4729679;
double r4729681 = r4729671 / r4729680;
double r4729682 = fma(r4729670, r4729676, r4729681);
double r4729683 = r4729682 - r4729663;
double r4729684 = r4729664 + r4729683;
double r4729685 = x;
double r4729686 = r4729684 + r4729685;
return r4729686;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus a
Initial program 13.1
rmApplied tan-sum0.2
rmApplied add-cube-cbrt0.3
Applied flip3--0.3
Applied associate-/r/0.3
Applied prod-diff0.3
Simplified0.2
Simplified0.2
rmApplied add-cbrt-cube0.2
Applied add-cbrt-cube0.2
Applied cbrt-unprod0.2
Final simplification0.2
herbie shell --seed 2019142 +o rules:numerics
(FPCore (x y z a)
:name "(+ x (- (tan (+ y z)) (tan a)))"
:pre (and (or (== x 0) (<= 0.5884142 x 505.5909)) (or (<= -1.796658e+308 y -9.425585e-310) (<= 1.284938e-309 y 1.751224e+308)) (or (<= -1.776707e+308 z -8.599796e-310) (<= 3.293145e-311 z 1.725154e+308)) (or (<= -1.796658e+308 a -9.425585e-310) (<= 1.284938e-309 a 1.751224e+308)))
(+ x (- (tan (+ y z)) (tan a))))