\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 21.37040012626098928194551263004541397095:\\
\;\;\;\;\frac{\frac{1}{\sqrt[3]{\frac{\sin \left(\frac{x}{y \cdot 2}\right)}{\tan \left(\frac{x}{y \cdot 2}\right)}} \cdot \sqrt[3]{\frac{\sin \left(\frac{x}{y \cdot 2}\right)}{\tan \left(\frac{x}{y \cdot 2}\right)}}}}{\sqrt[3]{\frac{\sin \left(\frac{x}{y \cdot 2}\right)}{\tan \left(\frac{x}{y \cdot 2}\right)}}}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}double f(double x, double y) {
double r638269 = x;
double r638270 = y;
double r638271 = 2.0;
double r638272 = r638270 * r638271;
double r638273 = r638269 / r638272;
double r638274 = tan(r638273);
double r638275 = sin(r638273);
double r638276 = r638274 / r638275;
return r638276;
}
double f(double x, double y) {
double r638277 = x;
double r638278 = y;
double r638279 = 2.0;
double r638280 = r638278 * r638279;
double r638281 = r638277 / r638280;
double r638282 = tan(r638281);
double r638283 = sin(r638281);
double r638284 = r638282 / r638283;
double r638285 = 21.37040012626099;
bool r638286 = r638284 <= r638285;
double r638287 = 1.0;
double r638288 = r638283 / r638282;
double r638289 = cbrt(r638288);
double r638290 = r638289 * r638289;
double r638291 = r638287 / r638290;
double r638292 = r638291 / r638289;
double r638293 = 1.0;
double r638294 = r638286 ? r638292 : r638293;
return r638294;
}




Bits error versus x




Bits error versus y
Results
| Original | 35.5 |
|---|---|
| Target | 28.7 |
| Herbie | 27.3 |
if (/ (tan (/ x (* y 2.0))) (sin (/ x (* y 2.0)))) < 21.37040012626099Initial program 26.4
rmApplied clear-num26.4
rmApplied add-cube-cbrt26.5
Applied associate-/r*26.5
if 21.37040012626099 < (/ (tan (/ x (* y 2.0))) (sin (/ x (* y 2.0)))) Initial program 63.8
Taylor expanded around 0 30.0
Final simplification27.3
herbie shell --seed 2019354 +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)))))