\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 r196986 = -2.0;
double r196987 = J;
double r196988 = r196986 * r196987;
double r196989 = K;
double r196990 = 2.0;
double r196991 = r196989 / r196990;
double r196992 = cos(r196991);
double r196993 = r196988 * r196992;
double r196994 = 1.0;
double r196995 = U;
double r196996 = r196990 * r196987;
double r196997 = r196996 * r196992;
double r196998 = r196995 / r196997;
double r196999 = pow(r196998, r196990);
double r197000 = r196994 + r196999;
double r197001 = sqrt(r197000);
double r197002 = r196993 * r197001;
return r197002;
}
double f(double J, double K, double U) {
double r197003 = -2.0;
double r197004 = J;
double r197005 = r197003 * r197004;
double r197006 = K;
double r197007 = 2.0;
double r197008 = r197006 / r197007;
double r197009 = cos(r197008);
double r197010 = 1.0;
double r197011 = sqrt(r197010);
double r197012 = U;
double r197013 = r197007 * r197004;
double r197014 = r197012 / r197013;
double r197015 = r197014 / r197009;
double r197016 = 2.0;
double r197017 = r197007 / r197016;
double r197018 = pow(r197015, r197017);
double r197019 = hypot(r197011, r197018);
double r197020 = r197009 * r197019;
double r197021 = r197005 * r197020;
return r197021;
}



Bits error versus J



Bits error versus K



Bits error versus U
Results
Initial program 17.8
rmApplied sqr-pow17.8
Applied add-sqr-sqrt17.8
Applied hypot-def8.2
rmApplied associate-*l*8.3
rmApplied associate-/r*8.2
Final simplification8.2
herbie shell --seed 2020083 +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)))))