\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{\frac{x}{y}}{2}\right)} \cdot \left(\frac{1}{\cos \left(\frac{\frac{x}{y}}{2}\right)} \cdot \frac{1}{\cos \left(\frac{\frac{x}{y}}{2}\right)}\right)}\right)\right)double f(double x, double y) {
double r28300038 = x;
double r28300039 = y;
double r28300040 = 2.0;
double r28300041 = r28300039 * r28300040;
double r28300042 = r28300038 / r28300041;
double r28300043 = tan(r28300042);
double r28300044 = sin(r28300042);
double r28300045 = r28300043 / r28300044;
return r28300045;
}
double f(double x, double y) {
double r28300046 = 1.0;
double r28300047 = x;
double r28300048 = y;
double r28300049 = r28300047 / r28300048;
double r28300050 = 2.0;
double r28300051 = r28300049 / r28300050;
double r28300052 = cos(r28300051);
double r28300053 = r28300046 / r28300052;
double r28300054 = r28300053 * r28300053;
double r28300055 = r28300053 * r28300054;
double r28300056 = cbrt(r28300055);
double r28300057 = expm1(r28300056);
double r28300058 = log1p(r28300057);
return r28300058;
}




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)))))