\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 r170423 = -2.0;
double r170424 = J;
double r170425 = r170423 * r170424;
double r170426 = K;
double r170427 = 2.0;
double r170428 = r170426 / r170427;
double r170429 = cos(r170428);
double r170430 = r170425 * r170429;
double r170431 = 1.0;
double r170432 = U;
double r170433 = r170427 * r170424;
double r170434 = r170433 * r170429;
double r170435 = r170432 / r170434;
double r170436 = pow(r170435, r170427);
double r170437 = r170431 + r170436;
double r170438 = sqrt(r170437);
double r170439 = r170430 * r170438;
return r170439;
}
double f(double J, double K, double U) {
double r170440 = -2.0;
double r170441 = J;
double r170442 = r170440 * r170441;
double r170443 = K;
double r170444 = 2.0;
double r170445 = r170443 / r170444;
double r170446 = cos(r170445);
double r170447 = 1.0;
double r170448 = sqrt(r170447);
double r170449 = U;
double r170450 = r170444 * r170441;
double r170451 = r170449 / r170450;
double r170452 = r170451 / r170446;
double r170453 = 2.0;
double r170454 = r170444 / r170453;
double r170455 = pow(r170452, r170454);
double r170456 = hypot(r170448, r170455);
double r170457 = r170446 * r170456;
double r170458 = r170442 * r170457;
return r170458;
}



Bits error versus J



Bits error versus K



Bits error versus U
Results
Initial program 18.0
rmApplied sqr-pow18.0
Applied add-sqr-sqrt18.0
Applied hypot-def7.9
rmApplied associate-*l*8.0
rmApplied associate-/r*7.9
Final simplification7.9
herbie shell --seed 2020056 +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)))))