\left(\left(-2 \cdot J\right) \cdot \cos \left(\frac{K}{2}\right)\right) \cdot \sqrt{1 + {\left(\frac{U}{\left(2 \cdot J\right) \cdot \cos \left(\frac{K}{2}\right)}\right)}^{2}}\left(-2 \cdot J\right) \cdot \left(\cos \left(\frac{K}{2}\right) \cdot \mathsf{hypot}\left(\sqrt{1}, {\left(\frac{\frac{U}{2 \cdot J}}{\cos \left(\frac{K}{2}\right)}\right)}^{\left(\frac{2}{2}\right)}\right)\right)double f(double J, double K, double U) {
double r178170 = -2.0;
double r178171 = J;
double r178172 = r178170 * r178171;
double r178173 = K;
double r178174 = 2.0;
double r178175 = r178173 / r178174;
double r178176 = cos(r178175);
double r178177 = r178172 * r178176;
double r178178 = 1.0;
double r178179 = U;
double r178180 = r178174 * r178171;
double r178181 = r178180 * r178176;
double r178182 = r178179 / r178181;
double r178183 = pow(r178182, r178174);
double r178184 = r178178 + r178183;
double r178185 = sqrt(r178184);
double r178186 = r178177 * r178185;
return r178186;
}
double f(double J, double K, double U) {
double r178187 = -2.0;
double r178188 = J;
double r178189 = r178187 * r178188;
double r178190 = K;
double r178191 = 2.0;
double r178192 = r178190 / r178191;
double r178193 = cos(r178192);
double r178194 = 1.0;
double r178195 = sqrt(r178194);
double r178196 = U;
double r178197 = r178191 * r178188;
double r178198 = r178196 / r178197;
double r178199 = r178198 / r178193;
double r178200 = 2.0;
double r178201 = r178191 / r178200;
double r178202 = pow(r178199, r178201);
double r178203 = hypot(r178195, r178202);
double r178204 = r178193 * r178203;
double r178205 = r178189 * r178204;
return r178205;
}



Bits error versus J



Bits error versus K



Bits error versus U
Results
Initial program 18.3
rmApplied sqr-pow18.3
Applied add-sqr-sqrt18.3
Applied hypot-def8.5
rmApplied associate-*l*8.5
rmApplied associate-/r*8.5
Final simplification8.5
herbie shell --seed 2019351 +o rules:numerics
(FPCore (J K U)
:name "Maksimov and Kolovsky, Equation (3)"
:precision binary64
(* (* (* -2 J) (cos (/ K 2))) (sqrt (+ 1 (pow (/ U (* (* 2 J) (cos (/ K 2)))) 2)))))