\frac{\cos \left(2 \cdot x\right)}{{c}^{2} \cdot \left(\left(x \cdot {s}^{2}\right) \cdot x\right)}\begin{array}{l}
\mathbf{if}\;c \leq 7.440582287819774 \cdot 10^{-222}:\\
\;\;\;\;\frac{\cos \left(2 \cdot x\right)}{\left(x \cdot {c}^{\left(\frac{2}{2}\right)}\right) \cdot \left(\left(x \cdot {s}^{\left(\frac{2}{2}\right)}\right) \cdot \left({c}^{\left(\frac{2}{2}\right)} \cdot {s}^{\left(\frac{2}{2}\right)}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\cos \left(2 \cdot x\right)}{{c}^{\left(\frac{2}{2}\right)} \cdot \left(\left|\left(x \cdot {s}^{\left(\frac{2}{2}\right)}\right) \cdot {c}^{\left(\frac{2}{4}\right)}\right| \cdot \left|\left(x \cdot {s}^{\left(\frac{2}{2}\right)}\right) \cdot {c}^{\left(\frac{2}{4}\right)}\right|\right)}\\
\end{array}(FPCore (x c s) :precision binary64 (/ (cos (* 2.0 x)) (* (pow c 2.0) (* (* x (pow s 2.0)) x))))
(FPCore (x c s)
:precision binary64
(if (<= c 7.440582287819774e-222)
(/
(cos (* 2.0 x))
(*
(* x (pow c (/ 2.0 2.0)))
(*
(* x (pow s (/ 2.0 2.0)))
(* (pow c (/ 2.0 2.0)) (pow s (/ 2.0 2.0))))))
(/
(cos (* 2.0 x))
(*
(pow c (/ 2.0 2.0))
(*
(fabs (* (* x (pow s (/ 2.0 2.0))) (pow c (/ 2.0 4.0))))
(fabs (* (* x (pow s (/ 2.0 2.0))) (pow c (/ 2.0 4.0)))))))))double code(double x, double c, double s) {
return (((double) cos(((double) (2.0 * x)))) / ((double) (((double) pow(c, 2.0)) * ((double) (((double) (x * ((double) pow(s, 2.0)))) * x)))));
}
double code(double x, double c, double s) {
double VAR;
if ((c <= 7.440582287819774e-222)) {
VAR = (((double) cos(((double) (2.0 * x)))) / ((double) (((double) (x * ((double) pow(c, (2.0 / 2.0))))) * ((double) (((double) (x * ((double) pow(s, (2.0 / 2.0))))) * ((double) (((double) pow(c, (2.0 / 2.0))) * ((double) pow(s, (2.0 / 2.0))))))))));
} else {
VAR = (((double) cos(((double) (2.0 * x)))) / ((double) (((double) pow(c, (2.0 / 2.0))) * ((double) (((double) fabs(((double) (((double) (x * ((double) pow(s, (2.0 / 2.0))))) * ((double) pow(c, (2.0 / 4.0))))))) * ((double) fabs(((double) (((double) (x * ((double) pow(s, (2.0 / 2.0))))) * ((double) pow(c, (2.0 / 4.0))))))))))));
}
return VAR;
}



Bits error versus x



Bits error versus c



Bits error versus s
Results
if c < 7.4405822878197738e-222Initial program 29.8
rmApplied sqr-pow29.8
Applied associate-*l*24.2
Simplified19.8
rmApplied sqr-pow19.8
Applied associate-*r*14.6
rmApplied associate-*l*10.3
Simplified10.3
rmApplied associate-*r*10.0
Simplified10.0
if 7.4405822878197738e-222 < c Initial program 25.1
rmApplied sqr-pow25.1
Applied associate-*l*21.6
Simplified18.6
rmApplied sqr-pow18.6
Applied associate-*r*12.1
rmApplied add-sqr-sqrt12.2
Simplified12.2
Simplified5.4
Final simplification7.9
herbie shell --seed 2020198
(FPCore (x c s)
:name "mixedcos"
:precision binary64
(/ (cos (* 2.0 x)) (* (pow c 2.0) (* (* x (pow s 2.0)) x))))