\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)}{\left(\sqrt[3]{\cos M}\right)}^{2} \cdot \left(\sqrt[3]{\cos M} \cdot e^{\left|m - n\right| - \left(\ell + {\left(\frac{m + n}{2} - M\right)}^{2}\right)}\right)(FPCore (K m n M l) :precision binary64 (* (cos (- (/ (* K (+ m n)) 2.0) M)) (exp (- (- (pow (- (/ (+ m n) 2.0) M) 2.0)) (- l (fabs (- m n)))))))
(FPCore (K m n M l) :precision binary64 (* (pow (cbrt (cos M)) 2.0) (* (cbrt (cos M)) (exp (- (fabs (- m n)) (+ l (pow (- (/ (+ m n) 2.0) M) 2.0)))))))
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 pow(cbrt(cos(M)), 2.0) * (cbrt(cos(M)) * exp(fabs(m - n) - (l + pow((((m + n) / 2.0) - M), 2.0))));
}



Bits error versus K



Bits error versus m



Bits error versus n



Bits error versus M



Bits error versus l
Results
Initial program 14.9
Simplified14.9
Taylor expanded around 0 1.3
Simplified1.3
rmApplied add-cube-cbrt_binary64_7951.3
Applied associate-*l*_binary64_7011.3
Simplified1.3
rmApplied pow2_binary64_8411.3
Final simplification1.3
herbie shell --seed 2021024
(FPCore (K m n M l)
:name "Maksimov and Kolovsky, Equation (32)"
:precision binary64
(* (cos (- (/ (* K (+ m n)) 2.0) M)) (exp (- (- (pow (- (/ (+ m n) 2.0) M) 2.0)) (- l (fabs (- m n)))))))