\frac{\tan \left(\frac{x}{y \cdot 2}\right)}{\sin \left(\frac{x}{y \cdot 2}\right)}\mathsf{log1p}\left(\mathsf{expm1}\left(\sqrt[3]{\frac{1}{\cos \left(\frac{x}{2 \cdot y}\right)} \cdot \left(\frac{1}{\cos \left(\frac{x}{2 \cdot y}\right)} \cdot \frac{1}{\cos \left(\frac{x}{2 \cdot y}\right)}\right)}\right)\right)double f(double x, double y) {
double r23371703 = x;
double r23371704 = y;
double r23371705 = 2.0;
double r23371706 = r23371704 * r23371705;
double r23371707 = r23371703 / r23371706;
double r23371708 = tan(r23371707);
double r23371709 = sin(r23371707);
double r23371710 = r23371708 / r23371709;
return r23371710;
}
double f(double x, double y) {
double r23371711 = 1.0;
double r23371712 = x;
double r23371713 = 2.0;
double r23371714 = y;
double r23371715 = r23371713 * r23371714;
double r23371716 = r23371712 / r23371715;
double r23371717 = cos(r23371716);
double r23371718 = r23371711 / r23371717;
double r23371719 = r23371718 * r23371718;
double r23371720 = r23371718 * r23371719;
double r23371721 = cbrt(r23371720);
double r23371722 = expm1(r23371721);
double r23371723 = log1p(r23371722);
return r23371723;
}




Bits error versus x




Bits error versus y
Results
| Original | 35.2 |
|---|---|
| Target | 28.5 |
| Herbie | 28.2 |
Initial program 35.2
rmApplied tan-quot35.2
rmApplied log1p-expm1-u35.3
Simplified28.2
rmApplied add-cbrt-cube28.2
Final simplification28.2
herbie shell --seed 2019172 +o rules:numerics
(FPCore (x y)
:name "Diagrams.TwoD.Layout.CirclePacking:approxRadius from diagrams-contrib-1.3.0.5"
: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)))))