\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]{\frac{2}{t}}}{\sqrt[3]{\tan k}}}{\frac{k}{\ell}} \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 r11213757 = 2.0;
double r11213758 = t;
double r11213759 = 3.0;
double r11213760 = pow(r11213758, r11213759);
double r11213761 = l;
double r11213762 = r11213761 * r11213761;
double r11213763 = r11213760 / r11213762;
double r11213764 = k;
double r11213765 = sin(r11213764);
double r11213766 = r11213763 * r11213765;
double r11213767 = tan(r11213764);
double r11213768 = r11213766 * r11213767;
double r11213769 = 1.0;
double r11213770 = r11213764 / r11213758;
double r11213771 = pow(r11213770, r11213757);
double r11213772 = r11213769 + r11213771;
double r11213773 = r11213772 - r11213769;
double r11213774 = r11213768 * r11213773;
double r11213775 = r11213757 / r11213774;
return r11213775;
}
double f(double t, double l, double k) {
double r11213776 = 2.0;
double r11213777 = t;
double r11213778 = r11213776 / r11213777;
double r11213779 = cbrt(r11213778);
double r11213780 = k;
double r11213781 = tan(r11213780);
double r11213782 = cbrt(r11213781);
double r11213783 = r11213779 / r11213782;
double r11213784 = l;
double r11213785 = r11213780 / r11213784;
double r11213786 = r11213783 / r11213785;
double r11213787 = sin(r11213780);
double r11213788 = r11213783 / r11213787;
double r11213789 = r11213786 * r11213788;
double r11213790 = r11213789 * r11213786;
return r11213790;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
Initial program 47.1
Simplified30.5
rmApplied associate-/l*28.3
Simplified15.4
rmApplied add-cube-cbrt15.6
Applied add-cube-cbrt15.7
Applied times-frac15.7
Applied times-frac15.1
Simplified2.8
rmApplied associate-*l*1.0
Final simplification1.0
herbie shell --seed 2019168 +o rules:numerics
(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))))