x + \left(\tan \left(y + z\right) - \tan a\right)
\frac{\mathsf{fma}\left(\mathsf{fma}\left(\tan z, \tan y, -1\right), \sin a, \left(\tan y + \tan z\right) \cdot \cos a\right)}{\cos a \cdot \left(1 - {\left(\frac{\sin z \cdot \tan y}{\cos z}\right)}^{3}\right)} \cdot \left(\left(\frac{\sin z \cdot \tan y}{\cos z} \cdot \frac{\sin z \cdot \tan y}{\cos z} + \frac{\sin z \cdot \tan y}{\cos z}\right) + 1\right) + xdouble f(double x, double y, double z, double a) {
double r124265 = x;
double r124266 = y;
double r124267 = z;
double r124268 = r124266 + r124267;
double r124269 = tan(r124268);
double r124270 = a;
double r124271 = tan(r124270);
double r124272 = r124269 - r124271;
double r124273 = r124265 + r124272;
return r124273;
}
double f(double x, double y, double z, double a) {
double r124274 = z;
double r124275 = tan(r124274);
double r124276 = y;
double r124277 = tan(r124276);
double r124278 = -1.0;
double r124279 = fma(r124275, r124277, r124278);
double r124280 = a;
double r124281 = sin(r124280);
double r124282 = r124277 + r124275;
double r124283 = cos(r124280);
double r124284 = r124282 * r124283;
double r124285 = fma(r124279, r124281, r124284);
double r124286 = 1.0;
double r124287 = sin(r124274);
double r124288 = r124287 * r124277;
double r124289 = cos(r124274);
double r124290 = r124288 / r124289;
double r124291 = 3.0;
double r124292 = pow(r124290, r124291);
double r124293 = r124286 - r124292;
double r124294 = r124283 * r124293;
double r124295 = r124285 / r124294;
double r124296 = r124290 * r124290;
double r124297 = r124296 + r124290;
double r124298 = r124297 + r124286;
double r124299 = r124295 * r124298;
double r124300 = x;
double r124301 = r124299 + r124300;
return r124301;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus a
Initial program 13.4
rmApplied tan-quot13.4
Applied tan-sum0.2
Applied frac-sub0.2
Simplified0.2
rmApplied tan-quot0.2
Applied associate-*r/0.2
Simplified0.2
rmApplied flip3--0.2
Applied associate-*l/0.2
Applied associate-/r/0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2019304 +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))))