\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)}\frac{\frac{\sqrt{1}}{{\left(\sqrt[3]{t}\right)}^{3}} \cdot \left(\frac{\sqrt{2}}{{\left(\sqrt[3]{t}\right)}^{3}} \cdot \left(\frac{\sqrt{2}}{{\left(\sqrt[3]{t}\right)}^{3} \cdot \sin k} \cdot \ell\right)\right)}{\tan k} \cdot \frac{\ell}{\mathsf{fma}\left(2, 1, {\left(\frac{k}{t}\right)}^{2}\right)}double f(double t, double l, double k) {
double r137480 = 2.0;
double r137481 = t;
double r137482 = 3.0;
double r137483 = pow(r137481, r137482);
double r137484 = l;
double r137485 = r137484 * r137484;
double r137486 = r137483 / r137485;
double r137487 = k;
double r137488 = sin(r137487);
double r137489 = r137486 * r137488;
double r137490 = tan(r137487);
double r137491 = r137489 * r137490;
double r137492 = 1.0;
double r137493 = r137487 / r137481;
double r137494 = pow(r137493, r137480);
double r137495 = r137492 + r137494;
double r137496 = r137495 + r137492;
double r137497 = r137491 * r137496;
double r137498 = r137480 / r137497;
return r137498;
}
double f(double t, double l, double k) {
double r137499 = 1.0;
double r137500 = sqrt(r137499);
double r137501 = t;
double r137502 = cbrt(r137501);
double r137503 = 3.0;
double r137504 = pow(r137502, r137503);
double r137505 = r137500 / r137504;
double r137506 = 2.0;
double r137507 = sqrt(r137506);
double r137508 = r137507 / r137504;
double r137509 = k;
double r137510 = sin(r137509);
double r137511 = r137504 * r137510;
double r137512 = r137507 / r137511;
double r137513 = l;
double r137514 = r137512 * r137513;
double r137515 = r137508 * r137514;
double r137516 = r137505 * r137515;
double r137517 = tan(r137509);
double r137518 = r137516 / r137517;
double r137519 = 2.0;
double r137520 = 1.0;
double r137521 = r137509 / r137501;
double r137522 = pow(r137521, r137506);
double r137523 = fma(r137519, r137520, r137522);
double r137524 = r137513 / r137523;
double r137525 = r137518 * r137524;
return r137525;
}



Bits error versus t



Bits error versus l



Bits error versus k
Initial program 32.7
Simplified32.7
rmApplied *-un-lft-identity32.7
Applied times-frac31.9
Applied associate-*r*29.4
Simplified28.1
rmApplied add-cube-cbrt28.4
Applied unpow-prod-down28.4
Applied associate-*l*27.2
rmApplied add-sqr-sqrt27.2
Applied times-frac27.1
Applied associate-*l*24.1
rmApplied unpow-prod-down24.1
Applied *-un-lft-identity24.1
Applied sqrt-prod24.1
Applied times-frac23.9
Applied associate-*l*21.9
Final simplification21.9
herbie shell --seed 2020056 +o rules:numerics
(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))))