\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 \frac{{\left(\frac{\sqrt[3]{\frac{1}{{t}^{1}}} \cdot \sqrt[3]{\frac{1}{{t}^{1}}}}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left(\left(\ell \cdot \cos k\right) \cdot {\left(\frac{\sqrt[3]{\frac{1}{{t}^{1}}}}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1}\right)}{\sin k \cdot \frac{\sin k}{\ell}}double f(double t, double l, double k) {
double r151131 = 2.0;
double r151132 = t;
double r151133 = 3.0;
double r151134 = pow(r151132, r151133);
double r151135 = l;
double r151136 = r151135 * r151135;
double r151137 = r151134 / r151136;
double r151138 = k;
double r151139 = sin(r151138);
double r151140 = r151137 * r151139;
double r151141 = tan(r151138);
double r151142 = r151140 * r151141;
double r151143 = 1.0;
double r151144 = r151138 / r151132;
double r151145 = pow(r151144, r151131);
double r151146 = r151143 + r151145;
double r151147 = r151146 - r151143;
double r151148 = r151142 * r151147;
double r151149 = r151131 / r151148;
return r151149;
}
double f(double t, double l, double k) {
double r151150 = 2.0;
double r151151 = 1.0;
double r151152 = t;
double r151153 = 1.0;
double r151154 = pow(r151152, r151153);
double r151155 = r151151 / r151154;
double r151156 = cbrt(r151155);
double r151157 = r151156 * r151156;
double r151158 = k;
double r151159 = 2.0;
double r151160 = r151150 / r151159;
double r151161 = pow(r151158, r151160);
double r151162 = r151157 / r151161;
double r151163 = pow(r151162, r151153);
double r151164 = l;
double r151165 = cos(r151158);
double r151166 = r151164 * r151165;
double r151167 = r151156 / r151161;
double r151168 = pow(r151167, r151153);
double r151169 = r151166 * r151168;
double r151170 = r151163 * r151169;
double r151171 = sin(r151158);
double r151172 = r151171 / r151164;
double r151173 = r151171 * r151172;
double r151174 = r151170 / r151173;
double r151175 = r151150 * r151174;
return r151175;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
Initial program 48.1
Simplified40.4
Taylor expanded around inf 22.1
rmApplied add-sqr-sqrt43.2
Applied unpow-prod-down43.2
Applied times-frac43.2
Simplified43.2
Simplified21.0
rmApplied frac-times20.1
Applied associate-*r/15.8
Simplified15.5
rmApplied sqr-pow15.5
Applied add-cube-cbrt15.8
Applied times-frac10.8
Applied unpow-prod-down10.8
Applied associate-*l*3.9
Simplified3.9
Final simplification3.9
herbie shell --seed 2019179
(FPCore (t l k)
:name "Toniolo and Linder, Equation (10-)"
(/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (- (+ 1.0 (pow (/ k t) 2.0)) 1.0))))