x + \left(\tan \left(y + z\right) - \tan a\right)
\sqrt[3]{\left(\left(\frac{\tan y + \tan z}{1 - \tan z \cdot \tan y} - \tan a\right) \cdot \left(\frac{\tan y + \tan z}{1 - \tan z \cdot \tan y} - \tan a\right)\right) \cdot \left(\frac{\tan y + \tan z}{1 - \frac{\sin y \cdot \tan z}{\cos y}} - \tan a\right)} + xdouble f(double x, double y, double z, double a) {
double r2585849 = x;
double r2585850 = y;
double r2585851 = z;
double r2585852 = r2585850 + r2585851;
double r2585853 = tan(r2585852);
double r2585854 = a;
double r2585855 = tan(r2585854);
double r2585856 = r2585853 - r2585855;
double r2585857 = r2585849 + r2585856;
return r2585857;
}
double f(double x, double y, double z, double a) {
double r2585858 = y;
double r2585859 = tan(r2585858);
double r2585860 = z;
double r2585861 = tan(r2585860);
double r2585862 = r2585859 + r2585861;
double r2585863 = 1.0;
double r2585864 = r2585861 * r2585859;
double r2585865 = r2585863 - r2585864;
double r2585866 = r2585862 / r2585865;
double r2585867 = a;
double r2585868 = tan(r2585867);
double r2585869 = r2585866 - r2585868;
double r2585870 = r2585869 * r2585869;
double r2585871 = sin(r2585858);
double r2585872 = r2585871 * r2585861;
double r2585873 = cos(r2585858);
double r2585874 = r2585872 / r2585873;
double r2585875 = r2585863 - r2585874;
double r2585876 = r2585862 / r2585875;
double r2585877 = r2585876 - r2585868;
double r2585878 = r2585870 * r2585877;
double r2585879 = cbrt(r2585878);
double r2585880 = x;
double r2585881 = r2585879 + r2585880;
return r2585881;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus a
Results
Initial program 13.1
rmApplied tan-sum0.2
rmApplied add-log-exp0.2
rmApplied add-cbrt-cube0.3
Simplified0.3
rmApplied tan-quot0.3
Applied associate-*r/0.3
Final simplification0.3
herbie shell --seed 2019128 +o rules:numerics
(FPCore (x y z a)
:name "(+ x (- (tan (+ y z)) (tan a)))"
:pre (and (or (== x 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))))