\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}^{2} \le 1.2069967661290195 \cdot 10^{-118}:\\
\;\;\;\;\frac{\frac{\cos \left(2 \cdot x\right)}{x \cdot {c}^{\left(\frac{2}{2}\right)}}}{{c}^{\left(\frac{2}{2}\right)} \cdot \left(\left(x \cdot {\left(\sqrt[3]{s} \cdot \sqrt[3]{s}\right)}^{2}\right) \cdot {\left(\sqrt[3]{s}\right)}^{2}\right)}\\
\mathbf{elif}\;{c}^{2} \le 5.4405840560867589 \cdot 10^{247}:\\
\;\;\;\;\frac{\cos \left(2 \cdot x\right)}{{c}^{2} \cdot \left(x \cdot \left({s}^{\left(\frac{2}{2}\right)} \cdot \left(x \cdot {s}^{\left(\frac{2}{2}\right)}\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{x \cdot {c}^{\left(\frac{2}{2}\right)}} \cdot \frac{\cos \left(2 \cdot x\right)}{x \cdot \left({c}^{\left(\frac{2}{2}\right)} \cdot {s}^{2}\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 ((((double) pow(c, 2.0)) <= 1.2069967661290195e-118)) {
VAR = ((double) (((double) (((double) cos(((double) (2.0 * x)))) / ((double) (x * ((double) pow(c, ((double) (2.0 / 2.0)))))))) / ((double) (((double) pow(c, ((double) (2.0 / 2.0)))) * ((double) (((double) (x * ((double) pow(((double) (((double) cbrt(s)) * ((double) cbrt(s)))), 2.0)))) * ((double) pow(((double) cbrt(s)), 2.0))))))));
} else {
double VAR_1;
if ((((double) pow(c, 2.0)) <= 5.440584056086759e+247)) {
VAR_1 = ((double) (((double) cos(((double) (2.0 * x)))) / ((double) (((double) pow(c, 2.0)) * ((double) (x * ((double) (((double) pow(s, ((double) (2.0 / 2.0)))) * ((double) (x * ((double) pow(s, ((double) (2.0 / 2.0))))))))))))));
} else {
VAR_1 = ((double) (((double) (1.0 / ((double) (x * ((double) pow(c, ((double) (2.0 / 2.0)))))))) * ((double) (((double) cos(((double) (2.0 * x)))) / ((double) (x * ((double) (((double) pow(c, ((double) (2.0 / 2.0)))) * ((double) pow(s, 2.0))))))))));
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus x



Bits error versus c



Bits error versus s
Results
if (pow c 2.0) < 1.2069967661290195e-118Initial program 44.0
rmApplied sqr-pow44.0
Applied associate-*l*33.2
Simplified26.1
rmApplied associate-*r*25.0
Simplified25.0
rmApplied associate-/r*25.0
Simplified25.0
rmApplied add-cube-cbrt25.5
Applied unpow-prod-down25.5
Applied associate-*r*21.4
if 1.2069967661290195e-118 < (pow c 2.0) < 5.4405840560867589e247Initial program 18.4
rmApplied sqr-pow18.4
Applied associate-*r*8.3
if 5.4405840560867589e247 < (pow c 2.0) Initial program 24.6
rmApplied sqr-pow24.6
Applied associate-*l*20.6
Simplified17.2
rmApplied associate-*r*16.7
Simplified16.7
rmApplied *-un-lft-identity16.7
Applied times-frac16.6
Simplified16.2
Final simplification15.0
herbie shell --seed 2020185
(FPCore (x c s)
:name "mixedcos"
:precision binary64
(/ (cos (* 2.0 x)) (* (pow c 2.0) (* (* x (pow s 2.0)) x))))