\frac{\tan \left(\frac{x}{y \cdot 2}\right)}{\sin \left(\frac{x}{y \cdot 2}\right)}\begin{array}{l}
\mathbf{if}\;\frac{\tan \left(\frac{x}{y \cdot 2}\right)}{\sin \left(\frac{x}{y \cdot 2}\right)} \leq 2.542199372979687:\\
\;\;\;\;\sqrt[3]{{\left(\frac{\tan \left(\frac{x}{y \cdot 2}\right)}{\sin \left(\frac{x}{y \cdot 2}\right)}\right)}^{3}}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}(FPCore (x y) :precision binary64 (/ (tan (/ x (* y 2.0))) (sin (/ x (* y 2.0)))))
(FPCore (x y) :precision binary64 (if (<= (/ (tan (/ x (* y 2.0))) (sin (/ x (* y 2.0)))) 2.542199372979687) (cbrt (pow (/ (tan (/ x (* y 2.0))) (sin (/ x (* y 2.0)))) 3.0)) 1.0))
double code(double x, double y) {
return tan(x / (y * 2.0)) / sin(x / (y * 2.0));
}
double code(double x, double y) {
double tmp;
if ((tan(x / (y * 2.0)) / sin(x / (y * 2.0))) <= 2.542199372979687) {
tmp = cbrt(pow((tan(x / (y * 2.0)) / sin(x / (y * 2.0))), 3.0));
} else {
tmp = 1.0;
}
return tmp;
}




Bits error versus x




Bits error versus y
Results
| Original | 36.3 |
|---|---|
| Target | 29.3 |
| Herbie | 28.1 |
if (/.f64 (tan.f64 (/.f64 x (*.f64 y 2))) (sin.f64 (/.f64 x (*.f64 y 2)))) < 2.54219937297968679Initial program 25.6
rmApplied add-cbrt-cube_binary64_1989225.6
Simplified25.6
if 2.54219937297968679 < (/.f64 (tan.f64 (/.f64 x (*.f64 y 2))) (sin.f64 (/.f64 x (*.f64 y 2)))) Initial program 62.5
Taylor expanded around 0 34.1
Final simplification28.1
herbie shell --seed 2020344
(FPCore (x y)
:name "Diagrams.TwoD.Layout.CirclePacking:approxRadius from diagrams-contrib-1.3.0.5"
:precision binary64
:herbie-target
(if (< y -1.2303690911306994e+114) 1.0 (if (< y -9.102852406811914e-222) (/ (sin (/ x (* y 2.0))) (* (sin (/ x (* y 2.0))) (log (exp (cos (/ x (* y 2.0))))))) 1.0))
(/ (tan (/ x (* y 2.0))) (sin (/ x (* y 2.0)))))