\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)}2 \cdot \left({\left(\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left({\left(\frac{\sqrt[3]{1}}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1} \cdot \left(\frac{\cos k}{\left|\sin k\right|} \cdot \frac{\ell}{\frac{\left|\sin k\right|}{\ell}}\right)\right)\right)double f(double t, double l, double k) {
double r94467 = 2.0;
double r94468 = t;
double r94469 = 3.0;
double r94470 = pow(r94468, r94469);
double r94471 = l;
double r94472 = r94471 * r94471;
double r94473 = r94470 / r94472;
double r94474 = k;
double r94475 = sin(r94474);
double r94476 = r94473 * r94475;
double r94477 = tan(r94474);
double r94478 = r94476 * r94477;
double r94479 = 1.0;
double r94480 = r94474 / r94468;
double r94481 = pow(r94480, r94467);
double r94482 = r94479 + r94481;
double r94483 = r94482 - r94479;
double r94484 = r94478 * r94483;
double r94485 = r94467 / r94484;
return r94485;
}
double f(double t, double l, double k) {
double r94486 = 2.0;
double r94487 = 1.0;
double r94488 = cbrt(r94487);
double r94489 = r94488 * r94488;
double r94490 = k;
double r94491 = 2.0;
double r94492 = r94486 / r94491;
double r94493 = pow(r94490, r94492);
double r94494 = r94489 / r94493;
double r94495 = 1.0;
double r94496 = pow(r94494, r94495);
double r94497 = t;
double r94498 = pow(r94497, r94495);
double r94499 = r94493 * r94498;
double r94500 = r94488 / r94499;
double r94501 = pow(r94500, r94495);
double r94502 = cos(r94490);
double r94503 = sin(r94490);
double r94504 = fabs(r94503);
double r94505 = r94502 / r94504;
double r94506 = l;
double r94507 = r94504 / r94506;
double r94508 = r94506 / r94507;
double r94509 = r94505 * r94508;
double r94510 = r94501 * r94509;
double r94511 = r94496 * r94510;
double r94512 = r94486 * r94511;
return r94512;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
Initial program 47.5
Simplified39.7
Taylor expanded around inf 21.8
rmApplied sqr-pow21.8
Applied associate-*l*19.7
rmApplied add-cube-cbrt19.7
Applied times-frac19.6
Applied unpow-prod-down19.6
Applied associate-*l*17.9
rmApplied add-sqr-sqrt17.9
Applied times-frac17.9
Simplified17.9
Simplified15.6
Final simplification15.6
herbie shell --seed 2020049
(FPCore (t l k)
:name "Toniolo and Linder, Equation (10-)"
:precision binary64
(/ 2 (* (* (* (/ (pow t 3) (* l l)) (sin k)) (tan k)) (- (+ 1 (pow (/ k t) 2)) 1))))