x + \left(\tan \left(y + z\right) - \tan a\right)
\frac{\cos a \cdot \left(\tan y + \tan z\right) - \sin a \cdot \left(1 - \frac{\sin z \cdot \sin y}{\cos z \cdot \cos y}\right)}{\left(1 - \frac{\sin z \cdot \sin y}{\cos z \cdot \cos y}\right) \cdot \cos a} + xdouble f(double x, double y, double z, double a) {
double r5412034 = x;
double r5412035 = y;
double r5412036 = z;
double r5412037 = r5412035 + r5412036;
double r5412038 = tan(r5412037);
double r5412039 = a;
double r5412040 = tan(r5412039);
double r5412041 = r5412038 - r5412040;
double r5412042 = r5412034 + r5412041;
return r5412042;
}
double f(double x, double y, double z, double a) {
double r5412043 = a;
double r5412044 = cos(r5412043);
double r5412045 = y;
double r5412046 = tan(r5412045);
double r5412047 = z;
double r5412048 = tan(r5412047);
double r5412049 = r5412046 + r5412048;
double r5412050 = r5412044 * r5412049;
double r5412051 = sin(r5412043);
double r5412052 = 1.0;
double r5412053 = sin(r5412047);
double r5412054 = sin(r5412045);
double r5412055 = r5412053 * r5412054;
double r5412056 = cos(r5412047);
double r5412057 = cos(r5412045);
double r5412058 = r5412056 * r5412057;
double r5412059 = r5412055 / r5412058;
double r5412060 = r5412052 - r5412059;
double r5412061 = r5412051 * r5412060;
double r5412062 = r5412050 - r5412061;
double r5412063 = r5412060 * r5412044;
double r5412064 = r5412062 / r5412063;
double r5412065 = x;
double r5412066 = r5412064 + r5412065;
return r5412066;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus a
Results
Initial program 13.0
rmApplied tan-sum0.2
Taylor expanded around inf 0.2
rmApplied tan-quot0.2
Applied frac-sub0.2
Final simplification0.2
herbie shell --seed 2019169 +o rules:numerics
(FPCore (x y z a)
:name "(+ x (- (tan (+ y z)) (tan a)))"
:pre (and (or (== x 0.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))))