\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{\ell}{k} \cdot \frac{\frac{\sqrt[3]{\frac{2}{t}}}{\sqrt[3]{\tan k}} \cdot \frac{\sqrt[3]{\frac{2}{t}}}{\sqrt[3]{\tan k}}}{\frac{\sin k}{\frac{\ell}{k}}}\right) \cdot \frac{\sqrt[3]{\frac{2}{t}}}{\sqrt[3]{\tan k}}double f(double t, double l, double k) {
double r3022242 = 2.0;
double r3022243 = t;
double r3022244 = 3.0;
double r3022245 = pow(r3022243, r3022244);
double r3022246 = l;
double r3022247 = r3022246 * r3022246;
double r3022248 = r3022245 / r3022247;
double r3022249 = k;
double r3022250 = sin(r3022249);
double r3022251 = r3022248 * r3022250;
double r3022252 = tan(r3022249);
double r3022253 = r3022251 * r3022252;
double r3022254 = 1.0;
double r3022255 = r3022249 / r3022243;
double r3022256 = pow(r3022255, r3022242);
double r3022257 = r3022254 + r3022256;
double r3022258 = r3022257 - r3022254;
double r3022259 = r3022253 * r3022258;
double r3022260 = r3022242 / r3022259;
return r3022260;
}
double f(double t, double l, double k) {
double r3022261 = l;
double r3022262 = k;
double r3022263 = r3022261 / r3022262;
double r3022264 = 2.0;
double r3022265 = t;
double r3022266 = r3022264 / r3022265;
double r3022267 = cbrt(r3022266);
double r3022268 = tan(r3022262);
double r3022269 = cbrt(r3022268);
double r3022270 = r3022267 / r3022269;
double r3022271 = r3022270 * r3022270;
double r3022272 = sin(r3022262);
double r3022273 = r3022272 / r3022263;
double r3022274 = r3022271 / r3022273;
double r3022275 = r3022263 * r3022274;
double r3022276 = r3022275 * r3022270;
return r3022276;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
Initial program 47.4
Simplified22.1
rmApplied associate-*r*18.2
rmApplied *-un-lft-identity18.2
Applied times-frac18.2
Applied times-frac18.0
Simplified14.4
Taylor expanded around inf 23.0
Simplified7.7
rmApplied add-cube-cbrt8.0
Applied add-cube-cbrt8.2
Applied times-frac8.2
Applied associate-*r*7.7
Simplified2.5
Final simplification2.5
herbie shell --seed 2019151 +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))))