x + \left(\tan \left(y + z\right) - \tan a\right)
x + \frac{\left(\tan y + \tan z\right) \cdot \cos a - \left(1 - \tan y \cdot \tan z\right) \cdot \sin a}{\left(1 - \left(\sin y \cdot \tan z\right) \cdot \frac{1}{\cos y}\right) \cdot \cos a}double f(double x, double y, double z, double a) {
double r182148 = x;
double r182149 = y;
double r182150 = z;
double r182151 = r182149 + r182150;
double r182152 = tan(r182151);
double r182153 = a;
double r182154 = tan(r182153);
double r182155 = r182152 - r182154;
double r182156 = r182148 + r182155;
return r182156;
}
double f(double x, double y, double z, double a) {
double r182157 = x;
double r182158 = y;
double r182159 = tan(r182158);
double r182160 = z;
double r182161 = tan(r182160);
double r182162 = r182159 + r182161;
double r182163 = a;
double r182164 = cos(r182163);
double r182165 = r182162 * r182164;
double r182166 = 1.0;
double r182167 = r182159 * r182161;
double r182168 = r182166 - r182167;
double r182169 = sin(r182163);
double r182170 = r182168 * r182169;
double r182171 = r182165 - r182170;
double r182172 = sin(r182158);
double r182173 = r182172 * r182161;
double r182174 = cos(r182158);
double r182175 = r182166 / r182174;
double r182176 = r182173 * r182175;
double r182177 = r182166 - r182176;
double r182178 = r182177 * r182164;
double r182179 = r182171 / r182178;
double r182180 = r182157 + r182179;
return r182180;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus a
Results
Initial program 13.2
rmApplied tan-quot13.2
Applied tan-sum0.2
Applied frac-sub0.2
rmApplied tan-quot0.2
Applied associate-*l/0.2
rmApplied div-inv0.2
Final simplification0.2
herbie shell --seed 2020001 +o rules:numerics
(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))))