x + \left(\tan \left(y + z\right) - \tan a\right)
\left(\frac{\frac{\cos z \cdot \left(\tan y \cdot \sin y\right) - \left(\sin z \cdot \tan z\right) \cdot \cos y}{\left(\cos z \cdot \cos y\right) \cdot \left(\tan y - \tan z\right)}}{1 - \frac{\tan y \cdot \sin z}{\cos z}} - \tan a\right) + xdouble f(double x, double y, double z, double a) {
double r31326926 = x;
double r31326927 = y;
double r31326928 = z;
double r31326929 = r31326927 + r31326928;
double r31326930 = tan(r31326929);
double r31326931 = a;
double r31326932 = tan(r31326931);
double r31326933 = r31326930 - r31326932;
double r31326934 = r31326926 + r31326933;
return r31326934;
}
double f(double x, double y, double z, double a) {
double r31326935 = z;
double r31326936 = cos(r31326935);
double r31326937 = y;
double r31326938 = tan(r31326937);
double r31326939 = sin(r31326937);
double r31326940 = r31326938 * r31326939;
double r31326941 = r31326936 * r31326940;
double r31326942 = sin(r31326935);
double r31326943 = tan(r31326935);
double r31326944 = r31326942 * r31326943;
double r31326945 = cos(r31326937);
double r31326946 = r31326944 * r31326945;
double r31326947 = r31326941 - r31326946;
double r31326948 = r31326936 * r31326945;
double r31326949 = r31326938 - r31326943;
double r31326950 = r31326948 * r31326949;
double r31326951 = r31326947 / r31326950;
double r31326952 = 1.0;
double r31326953 = r31326938 * r31326942;
double r31326954 = r31326953 / r31326936;
double r31326955 = r31326952 - r31326954;
double r31326956 = r31326951 / r31326955;
double r31326957 = a;
double r31326958 = tan(r31326957);
double r31326959 = r31326956 - r31326958;
double r31326960 = x;
double r31326961 = r31326959 + r31326960;
return r31326961;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus a
Results
Initial program 13.3
rmApplied tan-sum0.2
rmApplied tan-quot0.2
Applied associate-*r/0.2
rmApplied flip-+0.2
rmApplied tan-quot0.2
Applied associate-*l/0.2
Applied tan-quot0.2
Applied associate-*r/0.2
Applied frac-sub0.2
Applied associate-/l/0.2
Final simplification0.2
herbie shell --seed 2019125
(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))))