\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 \le 1.16359825181547585 \cdot 10^{-212}:\\
\;\;\;\;\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)}{\sqrt{\left|{c}^{\left(\frac{2}{2}\right)} \cdot \left(x \cdot {s}^{\left(\frac{2}{2}\right)}\right)\right|}} \cdot {\left(\frac{1}{\sqrt{\left|{c}^{\left(\frac{2}{2}\right)} \cdot \left(x \cdot {s}^{\left(\frac{2}{2}\right)}\right)\right|}}\right)}^{3}\\
\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 ((c <= 1.1635982518154759e-212)) {
VAR = ((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 = ((double) (((double) (((double) cos(((double) (2.0 * x)))) / ((double) sqrt(((double) fabs(((double) (((double) pow(c, ((double) (2.0 / 2.0)))) * ((double) (x * ((double) pow(s, ((double) (2.0 / 2.0)))))))))))))) * ((double) pow(((double) (1.0 / ((double) sqrt(((double) fabs(((double) (((double) pow(c, ((double) (2.0 / 2.0)))) * ((double) (x * ((double) pow(s, ((double) (2.0 / 2.0)))))))))))))), 3.0))));
}
return VAR;
}



Bits error versus x



Bits error versus c



Bits error versus s
Results
if c < 1.16359825181547585e-212Initial program 30.8
rmApplied sqr-pow30.8
Applied associate-*r*25.4
rmApplied add-sqr-sqrt25.4
Simplified25.4
Simplified3.6
Taylor expanded around inf 3.3
if 1.16359825181547585e-212 < c Initial program 24.8
rmApplied sqr-pow24.8
Applied associate-*r*17.9
rmApplied add-sqr-sqrt17.9
Simplified17.9
Simplified2.1
rmApplied associate-/r*1.8
rmApplied add-sqr-sqrt1.9
Applied div-inv1.9
Applied times-frac1.9
Simplified2.1
Final simplification2.7
herbie shell --seed 2020173
(FPCore (x c s)
:name "mixedcos"
:precision binary64
(/ (cos (* 2.0 x)) (* (pow c 2.0) (* (* x (pow s 2.0)) x))))