\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 \frac{\frac{\cos k}{\frac{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}{\ell}}}{{\left(\sqrt[3]{\sin k}\right)}^{2}}\right)\right)double f(double t, double l, double k) {
double r87700 = 2.0;
double r87701 = t;
double r87702 = 3.0;
double r87703 = pow(r87701, r87702);
double r87704 = l;
double r87705 = r87704 * r87704;
double r87706 = r87703 / r87705;
double r87707 = k;
double r87708 = sin(r87707);
double r87709 = r87706 * r87708;
double r87710 = tan(r87707);
double r87711 = r87709 * r87710;
double r87712 = 1.0;
double r87713 = r87707 / r87701;
double r87714 = pow(r87713, r87700);
double r87715 = r87712 + r87714;
double r87716 = r87715 - r87712;
double r87717 = r87711 * r87716;
double r87718 = r87700 / r87717;
return r87718;
}
double f(double t, double l, double k) {
double r87719 = 2.0;
double r87720 = 1.0;
double r87721 = cbrt(r87720);
double r87722 = r87721 * r87721;
double r87723 = k;
double r87724 = 2.0;
double r87725 = r87719 / r87724;
double r87726 = pow(r87723, r87725);
double r87727 = r87722 / r87726;
double r87728 = 1.0;
double r87729 = pow(r87727, r87728);
double r87730 = t;
double r87731 = pow(r87730, r87728);
double r87732 = r87726 * r87731;
double r87733 = r87721 / r87732;
double r87734 = pow(r87733, r87728);
double r87735 = cos(r87723);
double r87736 = sin(r87723);
double r87737 = cbrt(r87736);
double r87738 = 4.0;
double r87739 = pow(r87737, r87738);
double r87740 = l;
double r87741 = r87739 / r87740;
double r87742 = r87741 / r87740;
double r87743 = r87735 / r87742;
double r87744 = pow(r87737, r87724);
double r87745 = r87743 / r87744;
double r87746 = r87734 * r87745;
double r87747 = r87729 * r87746;
double r87748 = r87719 * r87747;
return r87748;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
Initial program 48.6
Simplified40.4
Taylor expanded around inf 22.2
rmApplied sqr-pow22.2
Applied associate-*l*19.8
rmApplied add-cube-cbrt19.8
Applied times-frac19.6
Applied unpow-prod-down19.6
Applied associate-*l*18.0
rmApplied add-cube-cbrt18.4
Applied unpow-prod-down18.4
Applied associate-/r*18.1
Simplified15.7
Final simplification15.7
herbie shell --seed 2020001
(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))))