\frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\left(\frac{\frac{\sqrt[3]{2}}{\sqrt{\frac{k}{t} \cdot \frac{k}{t} + 2}}}{\frac{\sqrt[3]{\tan k}}{\frac{1}{t}}} \cdot \frac{\frac{\ell}{t} \cdot \frac{\sqrt[3]{2} \cdot \sqrt[3]{2}}{\sqrt{\frac{k}{t} \cdot \frac{k}{t} + 2}}}{\sqrt[3]{\tan k} \cdot \sqrt[3]{\tan k}}\right) \cdot \frac{\frac{\ell}{t}}{\sin k}double f(double t, double l, double k) {
double r11793327 = 2.0;
double r11793328 = t;
double r11793329 = 3.0;
double r11793330 = pow(r11793328, r11793329);
double r11793331 = l;
double r11793332 = r11793331 * r11793331;
double r11793333 = r11793330 / r11793332;
double r11793334 = k;
double r11793335 = sin(r11793334);
double r11793336 = r11793333 * r11793335;
double r11793337 = tan(r11793334);
double r11793338 = r11793336 * r11793337;
double r11793339 = 1.0;
double r11793340 = r11793334 / r11793328;
double r11793341 = pow(r11793340, r11793327);
double r11793342 = r11793339 + r11793341;
double r11793343 = r11793342 + r11793339;
double r11793344 = r11793338 * r11793343;
double r11793345 = r11793327 / r11793344;
return r11793345;
}
double f(double t, double l, double k) {
double r11793346 = 2.0;
double r11793347 = cbrt(r11793346);
double r11793348 = k;
double r11793349 = t;
double r11793350 = r11793348 / r11793349;
double r11793351 = r11793350 * r11793350;
double r11793352 = r11793351 + r11793346;
double r11793353 = sqrt(r11793352);
double r11793354 = r11793347 / r11793353;
double r11793355 = tan(r11793348);
double r11793356 = cbrt(r11793355);
double r11793357 = 1.0;
double r11793358 = r11793357 / r11793349;
double r11793359 = r11793356 / r11793358;
double r11793360 = r11793354 / r11793359;
double r11793361 = l;
double r11793362 = r11793361 / r11793349;
double r11793363 = r11793347 * r11793347;
double r11793364 = r11793363 / r11793353;
double r11793365 = r11793362 * r11793364;
double r11793366 = r11793356 * r11793356;
double r11793367 = r11793365 / r11793366;
double r11793368 = r11793360 * r11793367;
double r11793369 = sin(r11793348);
double r11793370 = r11793362 / r11793369;
double r11793371 = r11793368 * r11793370;
return r11793371;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
Initial program 31.9
Simplified18.3
rmApplied times-frac16.8
Applied *-un-lft-identity16.8
Applied times-frac16.1
Applied *-un-lft-identity16.1
Applied times-frac14.7
Simplified14.7
rmApplied div-inv14.7
Applied *-un-lft-identity14.7
Applied times-frac12.5
Applied *-un-lft-identity12.5
Applied times-frac12.0
Simplified11.9
rmApplied *-un-lft-identity11.9
Applied add-cube-cbrt11.9
Applied times-frac11.9
Applied add-cube-cbrt12.1
Applied times-frac12.1
Applied add-sqr-sqrt12.1
Applied add-cube-cbrt12.2
Applied times-frac12.2
Applied times-frac12.1
Applied associate-*r*11.4
Simplified11.4
Final simplification11.4
herbie shell --seed 2019130
(FPCore (t l k)
:name "Toniolo and Linder, Equation (10+)"
(/ 2 (* (* (* (/ (pow t 3) (* l l)) (sin k)) (tan k)) (+ (+ 1 (pow (/ k t) 2)) 1))))