\left(x \cdot \cos \left(\frac{\left(\left(y \cdot 2 + 1\right) \cdot z\right) \cdot t}{16}\right)\right) \cdot \cos \left(\frac{\left(\left(a \cdot 2 + 1\right) \cdot b\right) \cdot t}{16}\right)x
double f(double x, double y, double z, double t, double a, double b) {
double r717929 = x;
double r717930 = y;
double r717931 = 2.0;
double r717932 = r717930 * r717931;
double r717933 = 1.0;
double r717934 = r717932 + r717933;
double r717935 = z;
double r717936 = r717934 * r717935;
double r717937 = t;
double r717938 = r717936 * r717937;
double r717939 = 16.0;
double r717940 = r717938 / r717939;
double r717941 = cos(r717940);
double r717942 = r717929 * r717941;
double r717943 = a;
double r717944 = r717943 * r717931;
double r717945 = r717944 + r717933;
double r717946 = b;
double r717947 = r717945 * r717946;
double r717948 = r717947 * r717937;
double r717949 = r717948 / r717939;
double r717950 = cos(r717949);
double r717951 = r717942 * r717950;
return r717951;
}
double f(double x, double __attribute__((unused)) y, double __attribute__((unused)) z, double __attribute__((unused)) t, double __attribute__((unused)) a, double __attribute__((unused)) b) {
double r717952 = x;
return r717952;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a




Bits error versus b
Results
| Original | 46.5 |
|---|---|
| Target | 44.6 |
| Herbie | 44.4 |
Initial program 46.5
Taylor expanded around 0 45.7
Taylor expanded around 0 44.4
Final simplification44.4
herbie shell --seed 2019304
(FPCore (x y z t a b)
:name "Codec.Picture.Jpg.FastDct:referenceDct from JuicyPixels-3.2.6.1"
:precision binary64
:herbie-target
(* x (cos (* (/ b 16) (/ t (+ (- 1 (* a 2)) (pow (* a 2) 2))))))
(* (* x (cos (/ (* (* (+ (* y 2) 1) z) t) 16))) (cos (/ (* (* (+ (* a 2) 1) b) t) 16))))