\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}\;x \le -2.84998239742904699 \cdot 10^{218}:\\
\;\;\;\;\frac{\frac{\cos \left(2 \cdot x\right)}{\left|{c}^{\left(\frac{2}{2}\right)} \cdot \left(x \cdot {s}^{\left(\frac{2}{2}\right)}\right)\right|}}{\left|\left({c}^{\left(\frac{2}{2}\right)} \cdot x\right) \cdot {s}^{\left(\frac{2}{2}\right)}\right|}\\
\mathbf{elif}\;x \le -2.42422156316495082 \cdot 10^{-275}:\\
\;\;\;\;\frac{\cos \left(2 \cdot x\right)}{{\left(\left|{\left({s}^{1} \cdot {c}^{1}\right)}^{1} \cdot x\right|\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\cos \left(2 \cdot x\right)}{\left|{c}^{\left(\frac{2}{2}\right)} \cdot \left(x \cdot {s}^{\left(\frac{2}{2}\right)}\right)\right| \cdot \left|{c}^{\left(\frac{2}{2}\right)} \cdot \left(x \cdot {s}^{\left(\frac{2}{2}\right)}\right)\right|}\\
\end{array}double code(double x, double c, double s) {
return ((double) (((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 ((x <= -2.849982397429047e+218)) {
VAR = ((double) (((double) (((double) cos(((double) (2.0 * x)))) / ((double) fabs(((double) (((double) pow(c, ((double) (2.0 / 2.0)))) * ((double) (x * ((double) pow(s, ((double) (2.0 / 2.0)))))))))))) / ((double) fabs(((double) (((double) (((double) pow(c, ((double) (2.0 / 2.0)))) * x)) * ((double) pow(s, ((double) (2.0 / 2.0))))))))));
} else {
double VAR_1;
if ((x <= -2.4242215631649508e-275)) {
VAR_1 = ((double) (((double) cos(((double) (2.0 * x)))) / ((double) pow(((double) fabs(((double) (((double) pow(((double) (((double) pow(s, 1.0)) * ((double) pow(c, 1.0)))), 1.0)) * x)))), 2.0))));
} else {
VAR_1 = ((double) (((double) cos(((double) (2.0 * x)))) / ((double) (((double) fabs(((double) (((double) pow(c, ((double) (2.0 / 2.0)))) * ((double) (x * ((double) pow(s, ((double) (2.0 / 2.0)))))))))) * ((double) fabs(((double) (((double) pow(c, ((double) (2.0 / 2.0)))) * ((double) (x * ((double) pow(s, ((double) (2.0 / 2.0))))))))))))));
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus x



Bits error versus c



Bits error versus s
Results
if x < -2.84998239742904699e218Initial program 27.1
rmApplied sqr-pow27.1
Applied associate-*r*17.2
rmApplied add-sqr-sqrt17.2
Simplified17.2
Simplified2.6
rmApplied associate-/r*2.5
rmApplied associate-*r*3.9
if -2.84998239742904699e218 < x < -2.42422156316495082e-275Initial program 27.2
rmApplied sqr-pow27.2
Applied associate-*r*22.3
rmApplied add-sqr-sqrt22.3
Simplified22.3
Simplified2.9
Taylor expanded around 0 1.8
if -2.42422156316495082e-275 < x Initial program 27.7
rmApplied sqr-pow27.7
Applied associate-*r*21.5
rmApplied add-sqr-sqrt21.5
Simplified21.5
Simplified2.9
Final simplification2.6
herbie shell --seed 2020155
(FPCore (x c s)
:name "mixedcos"
:precision binary64
(/ (cos (* 2.0 x)) (* (pow c 2.0) (* (* x (pow s 2.0)) x))))