\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 r593236 = x;
double r593237 = y;
double r593238 = 2.0;
double r593239 = r593237 * r593238;
double r593240 = 1.0;
double r593241 = r593239 + r593240;
double r593242 = z;
double r593243 = r593241 * r593242;
double r593244 = t;
double r593245 = r593243 * r593244;
double r593246 = 16.0;
double r593247 = r593245 / r593246;
double r593248 = cos(r593247);
double r593249 = r593236 * r593248;
double r593250 = a;
double r593251 = r593250 * r593238;
double r593252 = r593251 + r593240;
double r593253 = b;
double r593254 = r593252 * r593253;
double r593255 = r593254 * r593244;
double r593256 = r593255 / r593246;
double r593257 = cos(r593256);
double r593258 = r593249 * r593257;
return r593258;
}
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 r593259 = x;
return r593259;
}




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))))