x + \left(\tan \left(y + z\right) - \tan a\right)
x + \left(\frac{\tan y + \tan z}{1 - {\left(\tan y \cdot \tan z\right)}^{3}} \cdot \left(1 + \left(\left(\tan y \cdot \tan z\right) \cdot \frac{\sin y \cdot \sin z}{\cos y \cdot \cos z} + 1 \cdot \left(\tan y \cdot \tan z\right)\right)\right) - \tan a\right)double f(double x, double y, double z, double a) {
double r218153 = x;
double r218154 = y;
double r218155 = z;
double r218156 = r218154 + r218155;
double r218157 = tan(r218156);
double r218158 = a;
double r218159 = tan(r218158);
double r218160 = r218157 - r218159;
double r218161 = r218153 + r218160;
return r218161;
}
double f(double x, double y, double z, double a) {
double r218162 = x;
double r218163 = y;
double r218164 = tan(r218163);
double r218165 = z;
double r218166 = tan(r218165);
double r218167 = r218164 + r218166;
double r218168 = 1.0;
double r218169 = r218164 * r218166;
double r218170 = 3.0;
double r218171 = pow(r218169, r218170);
double r218172 = r218168 - r218171;
double r218173 = r218167 / r218172;
double r218174 = sin(r218163);
double r218175 = sin(r218165);
double r218176 = r218174 * r218175;
double r218177 = cos(r218163);
double r218178 = cos(r218165);
double r218179 = r218177 * r218178;
double r218180 = r218176 / r218179;
double r218181 = r218169 * r218180;
double r218182 = r218168 * r218169;
double r218183 = r218181 + r218182;
double r218184 = r218168 + r218183;
double r218185 = r218173 * r218184;
double r218186 = a;
double r218187 = tan(r218186);
double r218188 = r218185 - r218187;
double r218189 = r218162 + r218188;
return r218189;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus a
Results
Initial program 13.8
rmApplied tan-sum0.2
rmApplied flip3--0.2
Applied associate-/r/0.2
Simplified0.2
rmApplied tan-quot0.2
Applied tan-quot0.2
Applied frac-times0.2
Final simplification0.2
herbie shell --seed 2019344
(FPCore (x y z a)
:name "(+ x (- (tan (+ y z)) (tan a)))"
:precision binary64
:pre (and (or (== x 0.0) (<= 0.5884142 x 505.5909)) (or (<= -1.796658e+308 y -9.425585e-310) (<= 1.284938e-309 y 1.7512240000000001e+308)) (or (<= -1.7767070000000002e+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.7512240000000001e+308)))
(+ x (- (tan (+ y z)) (tan a))))