\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 r85270 = 2.0;
double r85271 = t;
double r85272 = 3.0;
double r85273 = pow(r85271, r85272);
double r85274 = l;
double r85275 = r85274 * r85274;
double r85276 = r85273 / r85275;
double r85277 = k;
double r85278 = sin(r85277);
double r85279 = r85276 * r85278;
double r85280 = tan(r85277);
double r85281 = r85279 * r85280;
double r85282 = 1.0;
double r85283 = r85277 / r85271;
double r85284 = pow(r85283, r85270);
double r85285 = r85282 + r85284;
double r85286 = r85285 - r85282;
double r85287 = r85281 * r85286;
double r85288 = r85270 / r85287;
return r85288;
}
double f(double t, double l, double k) {
double r85289 = 2.0;
double r85290 = 1.0;
double r85291 = cbrt(r85290);
double r85292 = r85291 * r85291;
double r85293 = k;
double r85294 = 2.0;
double r85295 = r85289 / r85294;
double r85296 = pow(r85293, r85295);
double r85297 = r85292 / r85296;
double r85298 = 1.0;
double r85299 = pow(r85297, r85298);
double r85300 = t;
double r85301 = pow(r85300, r85298);
double r85302 = r85296 * r85301;
double r85303 = r85291 / r85302;
double r85304 = pow(r85303, r85298);
double r85305 = cos(r85293);
double r85306 = sin(r85293);
double r85307 = cbrt(r85306);
double r85308 = 4.0;
double r85309 = pow(r85307, r85308);
double r85310 = l;
double r85311 = r85309 / r85310;
double r85312 = r85311 / r85310;
double r85313 = r85305 / r85312;
double r85314 = pow(r85307, r85294);
double r85315 = r85313 / r85314;
double r85316 = r85304 * r85315;
double r85317 = r85299 * r85316;
double r85318 = r85289 * r85317;
return r85318;
}



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))))