x + \left(\tan \left(y + z\right) - \tan a\right)
\frac{\left(x + \left(\frac{\tan y + \tan z}{1 - \tan y \cdot \tan z} - \tan a\right)\right) \cdot \left(x - \left(\frac{\tan y + \tan z}{1 - \tan y \cdot \tan z} - \tan a\right)\right)}{x - \left(\frac{\tan y + \tan z}{1 - \tan y \cdot \tan z} - \tan a\right)}double f(double x, double y, double z, double a) {
double r128610 = x;
double r128611 = y;
double r128612 = z;
double r128613 = r128611 + r128612;
double r128614 = tan(r128613);
double r128615 = a;
double r128616 = tan(r128615);
double r128617 = r128614 - r128616;
double r128618 = r128610 + r128617;
return r128618;
}
double f(double x, double y, double z, double a) {
double r128619 = x;
double r128620 = y;
double r128621 = tan(r128620);
double r128622 = z;
double r128623 = tan(r128622);
double r128624 = r128621 + r128623;
double r128625 = 1.0;
double r128626 = r128621 * r128623;
double r128627 = r128625 - r128626;
double r128628 = r128624 / r128627;
double r128629 = a;
double r128630 = tan(r128629);
double r128631 = r128628 - r128630;
double r128632 = r128619 + r128631;
double r128633 = r128619 - r128631;
double r128634 = r128632 * r128633;
double r128635 = r128634 / r128633;
return r128635;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus a
Results
Initial program 13.5
rmApplied tan-sum0.2
rmApplied flip-+0.4
rmApplied difference-of-squares0.3
Final simplification0.3
herbie shell --seed 2019235
(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))))