\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.388236677338849:\\
\;\;\;\;\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 r659303 = x;
double r659304 = y;
double r659305 = 2.0;
double r659306 = r659304 * r659305;
double r659307 = r659303 / r659306;
double r659308 = tan(r659307);
double r659309 = sin(r659307);
double r659310 = r659308 / r659309;
return r659310;
}
double f(double x, double y) {
double r659311 = x;
double r659312 = y;
double r659313 = 2.0;
double r659314 = r659312 * r659313;
double r659315 = r659311 / r659314;
double r659316 = tan(r659315);
double r659317 = sin(r659315);
double r659318 = r659316 / r659317;
double r659319 = 2.388236677338849;
bool r659320 = r659318 <= r659319;
double r659321 = cbrt(r659318);
double r659322 = r659321 * r659321;
double r659323 = r659322 * r659321;
double r659324 = 1.0;
double r659325 = r659320 ? r659323 : r659324;
return r659325;
}




Bits error versus x




Bits error versus y
Results
| Original | 35.8 |
|---|---|
| Target | 29.0 |
| Herbie | 27.9 |
if (/ (tan (/ x (* y 2.0))) (sin (/ x (* y 2.0)))) < 2.388236677338849Initial program 24.9
rmApplied add-cube-cbrt25.0
if 2.388236677338849 < (/ (tan (/ x (* y 2.0))) (sin (/ x (* y 2.0)))) Initial program 62.4
Taylor expanded around 0 35.0
Final simplification27.9
herbie shell --seed 2020047 +o rules:numerics
(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 (if (< y -9.102852406811914e-222) (/ (sin (/ x (* y 2))) (* (sin (/ x (* y 2))) (log (exp (cos (/ x (* y 2))))))) 1))
(/ (tan (/ x (* y 2))) (sin (/ x (* y 2)))))