\cos \left(\frac{K \cdot \left(m + n\right)}{2} - M\right) \cdot e^{\left(-{\left(\frac{m + n}{2} - M\right)}^{2}\right) - \left(\ell - \left|m - n\right|\right)}\sqrt{\frac{1}{e^{{\left(\frac{m + n}{2} - M\right)}^{2} + \left(\ell - \left|m - n\right|\right)}}} \cdot \sqrt{\frac{1}{e^{{\left(\frac{m + n}{2} - M\right)}^{2} + \left(\ell - \left|m - n\right|\right)}}}double code(double K, double m, double n, double M, double l) {
return (cos((((K * (m + n)) / 2.0) - M)) * exp((-pow((((m + n) / 2.0) - M), 2.0) - (l - fabs((m - n))))));
}
double code(double K, double m, double n, double M, double l) {
return (sqrt((1.0 / exp((pow((((m + n) / 2.0) - M), 2.0) + (l - fabs((m - n))))))) * sqrt((1.0 / exp((pow((((m + n) / 2.0) - M), 2.0) + (l - fabs((m - n))))))));
}



Bits error versus K



Bits error versus m



Bits error versus n



Bits error versus M



Bits error versus l
Results
Initial program 14.7
Simplified14.7
Taylor expanded around 0 1.4
rmApplied add-sqr-sqrt1.4
Final simplification1.4
herbie shell --seed 2020066
(FPCore (K m n M l)
:name "Maksimov and Kolovsky, Equation (32)"
:precision binary64
(* (cos (- (/ (* K (+ m n)) 2) M)) (exp (- (- (pow (- (/ (+ m n) 2) M) 2)) (- l (fabs (- m n)))))))