\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(\left(\frac{\sqrt[3]{\sqrt[3]{\frac{2}{t}} \cdot \sqrt[3]{\frac{2}{t}}}}{\frac{1}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}} \cdot \frac{\frac{\sqrt[3]{\sqrt[3]{\frac{2}{t}}}}{\sqrt[3]{\tan k}}}{\frac{k}{\sqrt[3]{\ell}}}\right) \cdot \frac{\frac{\sqrt[3]{\frac{2}{t}}}{\sqrt[3]{\tan k}}}{\sin k}\right) \cdot \frac{\frac{\sqrt[3]{\frac{2}{t}}}{\sqrt[3]{\tan k}}}{\frac{k}{\ell}}double f(double t, double l, double k) {
double r11587841 = 2.0;
double r11587842 = t;
double r11587843 = 3.0;
double r11587844 = pow(r11587842, r11587843);
double r11587845 = l;
double r11587846 = r11587845 * r11587845;
double r11587847 = r11587844 / r11587846;
double r11587848 = k;
double r11587849 = sin(r11587848);
double r11587850 = r11587847 * r11587849;
double r11587851 = tan(r11587848);
double r11587852 = r11587850 * r11587851;
double r11587853 = 1.0;
double r11587854 = r11587848 / r11587842;
double r11587855 = pow(r11587854, r11587841);
double r11587856 = r11587853 + r11587855;
double r11587857 = r11587856 - r11587853;
double r11587858 = r11587852 * r11587857;
double r11587859 = r11587841 / r11587858;
return r11587859;
}
double f(double t, double l, double k) {
double r11587860 = 2.0;
double r11587861 = t;
double r11587862 = r11587860 / r11587861;
double r11587863 = cbrt(r11587862);
double r11587864 = r11587863 * r11587863;
double r11587865 = cbrt(r11587864);
double r11587866 = 1.0;
double r11587867 = l;
double r11587868 = cbrt(r11587867);
double r11587869 = r11587868 * r11587868;
double r11587870 = r11587866 / r11587869;
double r11587871 = r11587865 / r11587870;
double r11587872 = cbrt(r11587863);
double r11587873 = k;
double r11587874 = tan(r11587873);
double r11587875 = cbrt(r11587874);
double r11587876 = r11587872 / r11587875;
double r11587877 = r11587873 / r11587868;
double r11587878 = r11587876 / r11587877;
double r11587879 = r11587871 * r11587878;
double r11587880 = r11587863 / r11587875;
double r11587881 = sin(r11587873);
double r11587882 = r11587880 / r11587881;
double r11587883 = r11587879 * r11587882;
double r11587884 = r11587873 / r11587867;
double r11587885 = r11587880 / r11587884;
double r11587886 = r11587883 * r11587885;
return r11587886;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
Initial program 47.2
Simplified30.8
rmApplied associate-/r/30.7
Applied times-frac30.6
Applied associate-/l*28.0
Simplified15.0
rmApplied add-cube-cbrt15.2
Applied add-cube-cbrt15.3
Applied times-frac15.3
Applied times-frac14.7
Simplified2.7
rmApplied associate-*l*1.0
rmApplied add-cube-cbrt1.1
Applied *-un-lft-identity1.1
Applied times-frac1.1
Applied *-un-lft-identity1.1
Applied add-cube-cbrt1.1
Applied cbrt-prod1.1
Applied times-frac1.1
Applied times-frac1.1
Final simplification1.1
herbie shell --seed 2019164
(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))))