\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)} \le 2.094986369407834114753086396376602351665:\\
\;\;\;\;\left(\sqrt[3]{\frac{\tan \left(\frac{x}{y \cdot 2}\right)}{\sin \left(\frac{x}{y \cdot 2}\right)}} \cdot \sqrt[3]{\frac{\tan \left(\frac{x}{y \cdot 2}\right)}{\sin \left(\frac{x}{y \cdot 2}\right)}}\right) \cdot \sqrt[3]{\frac{\tan \left(\frac{x}{y \cdot 2}\right)}{\sin \left(\frac{x}{y \cdot 2}\right)}}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}double f(double x, double y) {
double r476887 = x;
double r476888 = y;
double r476889 = 2.0;
double r476890 = r476888 * r476889;
double r476891 = r476887 / r476890;
double r476892 = tan(r476891);
double r476893 = sin(r476891);
double r476894 = r476892 / r476893;
return r476894;
}
double f(double x, double y) {
double r476895 = x;
double r476896 = y;
double r476897 = 2.0;
double r476898 = r476896 * r476897;
double r476899 = r476895 / r476898;
double r476900 = tan(r476899);
double r476901 = sin(r476899);
double r476902 = r476900 / r476901;
double r476903 = 2.094986369407834;
bool r476904 = r476902 <= r476903;
double r476905 = cbrt(r476902);
double r476906 = r476905 * r476905;
double r476907 = r476906 * r476905;
double r476908 = 1.0;
double r476909 = r476904 ? r476907 : r476908;
return r476909;
}




Bits error versus x




Bits error versus y
Results
| Original | 36.1 |
|---|---|
| Target | 29.3 |
| Herbie | 28.0 |
if (/ (tan (/ x (* y 2.0))) (sin (/ x (* y 2.0)))) < 2.094986369407834Initial program 25.3
rmApplied add-cube-cbrt25.3
if 2.094986369407834 < (/ (tan (/ x (* y 2.0))) (sin (/ x (* y 2.0)))) Initial program 62.1
Taylor expanded around 0 34.4
Final simplification28.0
herbie shell --seed 2019303
(FPCore (x y)
:name "Diagrams.TwoD.Layout.CirclePacking:approxRadius from diagrams-contrib-1.3.0.5"
:precision binary64
:herbie-target
(if (< y -1.23036909113069936e114) 1 (if (< y -9.1028524068119138e-222) (/ (sin (/ x (* y 2))) (* (sin (/ x (* y 2))) (log (exp (cos (/ x (* y 2))))))) 1))
(/ (tan (/ x (* y 2))) (sin (/ x (* y 2)))))