\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[3]{\frac{2}{t}}}{\sqrt[3]{\tan k}} \cdot \left(\frac{\sqrt[3]{\frac{2}{t}}}{\sqrt[3]{\tan k}} \cdot \left(\sqrt[3]{\cos k} \cdot \frac{\frac{\sqrt[3]{\frac{2}{t}}}{\sqrt[3]{\sin k}}}{\frac{k}{\ell}}\right)\right)}{\frac{k}{\ell} \cdot \sin k}double f(double t, double l, double k) {
double r10130105 = 2.0;
double r10130106 = t;
double r10130107 = 3.0;
double r10130108 = pow(r10130106, r10130107);
double r10130109 = l;
double r10130110 = r10130109 * r10130109;
double r10130111 = r10130108 / r10130110;
double r10130112 = k;
double r10130113 = sin(r10130112);
double r10130114 = r10130111 * r10130113;
double r10130115 = tan(r10130112);
double r10130116 = r10130114 * r10130115;
double r10130117 = 1.0;
double r10130118 = r10130112 / r10130106;
double r10130119 = pow(r10130118, r10130105);
double r10130120 = r10130117 + r10130119;
double r10130121 = r10130120 - r10130117;
double r10130122 = r10130116 * r10130121;
double r10130123 = r10130105 / r10130122;
return r10130123;
}
double f(double t, double l, double k) {
double r10130124 = 2.0;
double r10130125 = t;
double r10130126 = r10130124 / r10130125;
double r10130127 = cbrt(r10130126);
double r10130128 = k;
double r10130129 = tan(r10130128);
double r10130130 = cbrt(r10130129);
double r10130131 = r10130127 / r10130130;
double r10130132 = cos(r10130128);
double r10130133 = cbrt(r10130132);
double r10130134 = sin(r10130128);
double r10130135 = cbrt(r10130134);
double r10130136 = r10130127 / r10130135;
double r10130137 = l;
double r10130138 = r10130128 / r10130137;
double r10130139 = r10130136 / r10130138;
double r10130140 = r10130133 * r10130139;
double r10130141 = r10130131 * r10130140;
double r10130142 = r10130131 * r10130141;
double r10130143 = r10130138 * r10130134;
double r10130144 = r10130142 / r10130143;
return r10130144;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
Initial program 47.3
Simplified30.3
rmApplied associate-/l*28.0
Simplified15.3
rmApplied add-cube-cbrt15.5
Applied add-cube-cbrt15.5
Applied times-frac15.5
Applied times-frac14.9
Simplified2.7
rmApplied associate-*r/2.7
Applied frac-times1.0
rmApplied associate-/r/1.0
Applied tan-quot0.9
Applied cbrt-div0.9
Applied associate-/r/0.9
Applied times-frac0.9
Final simplification0.9
herbie shell --seed 2019162 +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))))