\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 10.123476350845582:\\
\;\;\;\;\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 r614258 = x;
double r614259 = y;
double r614260 = 2.0;
double r614261 = r614259 * r614260;
double r614262 = r614258 / r614261;
double r614263 = tan(r614262);
double r614264 = sin(r614262);
double r614265 = r614263 / r614264;
return r614265;
}
double f(double x, double y) {
double r614266 = x;
double r614267 = y;
double r614268 = 2.0;
double r614269 = r614267 * r614268;
double r614270 = r614266 / r614269;
double r614271 = tan(r614270);
double r614272 = sin(r614270);
double r614273 = r614271 / r614272;
double r614274 = 10.123476350845582;
bool r614275 = r614273 <= r614274;
double r614276 = cbrt(r614273);
double r614277 = r614276 * r614276;
double r614278 = r614277 * r614276;
double r614279 = 1.0;
double r614280 = r614275 ? r614278 : r614279;
return r614280;
}




Bits error versus x




Bits error versus y
Results
| Original | 35.4 |
|---|---|
| Target | 28.8 |
| Herbie | 27.4 |
if (/ (tan (/ x (* y 2.0))) (sin (/ x (* y 2.0)))) < 10.123476350845582Initial program 26.1
rmApplied add-cube-cbrt26.1
if 10.123476350845582 < (/ (tan (/ x (* y 2.0))) (sin (/ x (* y 2.0)))) Initial program 63.5
Taylor expanded around 0 31.4
Final simplification27.4
herbie shell --seed 2020064 +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)))))