\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 r91824 = 2.0;
double r91825 = t;
double r91826 = 3.0;
double r91827 = pow(r91825, r91826);
double r91828 = l;
double r91829 = r91828 * r91828;
double r91830 = r91827 / r91829;
double r91831 = k;
double r91832 = sin(r91831);
double r91833 = r91830 * r91832;
double r91834 = tan(r91831);
double r91835 = r91833 * r91834;
double r91836 = 1.0;
double r91837 = r91831 / r91825;
double r91838 = pow(r91837, r91824);
double r91839 = r91836 + r91838;
double r91840 = r91839 - r91836;
double r91841 = r91835 * r91840;
double r91842 = r91824 / r91841;
return r91842;
}
double f(double t, double l, double k) {
double r91843 = 2.0;
double r91844 = 1.0;
double r91845 = cbrt(r91844);
double r91846 = r91845 * r91845;
double r91847 = k;
double r91848 = 2.0;
double r91849 = r91843 / r91848;
double r91850 = pow(r91847, r91849);
double r91851 = r91846 / r91850;
double r91852 = 1.0;
double r91853 = pow(r91851, r91852);
double r91854 = t;
double r91855 = pow(r91854, r91852);
double r91856 = r91850 * r91855;
double r91857 = r91845 / r91856;
double r91858 = pow(r91857, r91852);
double r91859 = cos(r91847);
double r91860 = sin(r91847);
double r91861 = cbrt(r91860);
double r91862 = 4.0;
double r91863 = pow(r91861, r91862);
double r91864 = l;
double r91865 = r91863 / r91864;
double r91866 = r91865 / r91864;
double r91867 = r91859 / r91866;
double r91868 = pow(r91861, r91848);
double r91869 = r91867 / r91868;
double r91870 = r91858 * r91869;
double r91871 = r91853 * r91870;
double r91872 = r91843 * r91871;
return r91872;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
Initial program 48.5
Simplified41.2
Taylor expanded around inf 22.4
rmApplied sqr-pow22.4
Applied associate-*l*20.0
rmApplied add-cube-cbrt20.0
Applied times-frac19.9
Applied unpow-prod-down19.9
Applied associate-*l*18.0
rmApplied add-cube-cbrt18.3
Applied unpow-prod-down18.3
Applied associate-/r*18.0
Simplified15.6
Final simplification15.6
herbie shell --seed 2020024
(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))))