\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}\;{s}^{2} \le 3.49050698809960419 \cdot 10^{211}:\\
\;\;\;\;\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{\sqrt[3]{\cos \left(2 \cdot x\right)} \cdot \sqrt[3]{\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 \frac{\sqrt[3]{\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|}\\
\end{array}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 ((((double) pow(s, 2.0)) <= 3.490506988099604e+211)) {
VAR = (((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 = ((double) ((((double) (((double) cbrt(((double) cos(((double) (2.0 * x)))))) * ((double) cbrt(((double) cos(((double) (2.0 * x)))))))) / ((double) fabs(((double) (((double) pow(c, (2.0 / 2.0))) * ((double) (x * ((double) pow(s, (2.0 / 2.0)))))))))) * (((double) cbrt(((double) cos(((double) (2.0 * x)))))) / ((double) fabs(((double) (((double) pow(c, (2.0 / 2.0))) * ((double) (x * ((double) pow(s, (2.0 / 2.0))))))))))));
}
return VAR;
}



Bits error versus x



Bits error versus c



Bits error versus s
Results
if (pow s 2.0) < 3.49050698809960419e211Initial program 31.9
rmApplied sqr-pow31.9
Applied associate-*r*25.3
rmApplied add-sqr-sqrt25.3
Simplified25.3
Simplified3.5
Taylor expanded around inf 3.2
if 3.49050698809960419e211 < (pow s 2.0) Initial program 23.2
rmApplied sqr-pow23.2
Applied associate-*r*17.6
rmApplied add-sqr-sqrt17.6
Simplified17.6
Simplified2.7
rmApplied add-cube-cbrt2.7
Applied times-frac2.5
Final simplification2.9
herbie shell --seed 2020182
(FPCore (x c s)
:name "mixedcos"
:precision binary64
(/ (cos (* 2.0 x)) (* (pow c 2.0) (* (* x (pow s 2.0)) x))))