\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 \cdot \cos \left(\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right) \cdot \left(\left(1 + 2 \cdot a\right) \cdot \left(\frac{b}{16} \cdot \left(\sqrt[3]{\sqrt[3]{t}} \cdot {\left(\sqrt[3]{\sqrt[3]{t}}\right)}^{2}\right)\right)\right)\right)double code(double x, double y, double z, double t, double a, double b) {
return ((double) (((double) (x * ((double) cos(((double) (((double) (((double) (((double) (((double) (y * 2.0)) + 1.0)) * z)) * t)) / 16.0)))))) * ((double) cos(((double) (((double) (((double) (((double) (((double) (a * 2.0)) + 1.0)) * b)) * t)) / 16.0))))));
}
double code(double x, double y, double z, double t, double a, double b) {
return ((double) (x * ((double) cos(((double) (((double) (((double) cbrt(t)) * ((double) cbrt(t)))) * ((double) (((double) (1.0 + ((double) (2.0 * a)))) * ((double) (((double) (b / 16.0)) * ((double) (((double) cbrt(((double) cbrt(t)))) * ((double) pow(((double) cbrt(((double) cbrt(t)))), 2.0))))))))))))));
}




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.7 |
|---|---|
| Target | 45.0 |
| Herbie | 45.9 |
Initial program 46.7
Simplified46.4
Taylor expanded around 0 46.0
rmApplied add-cube-cbrt46.0
Applied associate-*l*46.0
Simplified45.9
rmApplied add-cube-cbrt45.9
Simplified45.9
Final simplification45.9
herbie shell --seed 2020181
(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.0) (/ t (+ (- 1.0 (* a 2.0)) (pow (* a 2.0) 2.0))))))
(* (* x (cos (/ (* (* (+ (* y 2.0) 1.0) z) t) 16.0))) (cos (/ (* (* (+ (* a 2.0) 1.0) b) t) 16.0))))