\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(\sqrt{\mathsf{hypot}\left({\left(\frac{U}{\left(\cos \left(\frac{K}{2}\right) \cdot J\right) \cdot 2}\right)}^{\left(\frac{2}{2}\right)}, \sqrt{1}\right)} \cdot \left(\left(\cos \left(\frac{K}{2}\right) \cdot J\right) \cdot -2\right)\right) \cdot \sqrt{\mathsf{hypot}\left({\left(\frac{U}{\left(\cos \left(\frac{K}{2}\right) \cdot J\right) \cdot 2}\right)}^{\left(\frac{2}{2}\right)}, \sqrt{1}\right)}double f(double J, double K, double U) {
double r7725132 = -2.0;
double r7725133 = J;
double r7725134 = r7725132 * r7725133;
double r7725135 = K;
double r7725136 = 2.0;
double r7725137 = r7725135 / r7725136;
double r7725138 = cos(r7725137);
double r7725139 = r7725134 * r7725138;
double r7725140 = 1.0;
double r7725141 = U;
double r7725142 = r7725136 * r7725133;
double r7725143 = r7725142 * r7725138;
double r7725144 = r7725141 / r7725143;
double r7725145 = pow(r7725144, r7725136);
double r7725146 = r7725140 + r7725145;
double r7725147 = sqrt(r7725146);
double r7725148 = r7725139 * r7725147;
return r7725148;
}
double f(double J, double K, double U) {
double r7725149 = U;
double r7725150 = K;
double r7725151 = 2.0;
double r7725152 = r7725150 / r7725151;
double r7725153 = cos(r7725152);
double r7725154 = J;
double r7725155 = r7725153 * r7725154;
double r7725156 = r7725155 * r7725151;
double r7725157 = r7725149 / r7725156;
double r7725158 = 2.0;
double r7725159 = r7725151 / r7725158;
double r7725160 = pow(r7725157, r7725159);
double r7725161 = 1.0;
double r7725162 = sqrt(r7725161);
double r7725163 = hypot(r7725160, r7725162);
double r7725164 = sqrt(r7725163);
double r7725165 = -2.0;
double r7725166 = r7725155 * r7725165;
double r7725167 = r7725164 * r7725166;
double r7725168 = r7725167 * r7725164;
return r7725168;
}



Bits error versus J



Bits error versus K



Bits error versus U
Results
Initial program 17.9
Simplified17.9
rmApplied add-sqr-sqrt17.9
Applied sqr-pow17.9
Applied hypot-def8.4
rmApplied add-sqr-sqrt8.5
Applied associate-*l*8.5
Final simplification8.5
herbie shell --seed 2019168 +o rules:numerics
(FPCore (J K U)
:name "Maksimov and Kolovsky, Equation (3)"
(* (* (* -2.0 J) (cos (/ K 2.0))) (sqrt (+ 1.0 (pow (/ U (* (* 2.0 J) (cos (/ K 2.0)))) 2.0)))))