\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} \leq 0:\\
\;\;\;\;\frac{\cos \left(2 \cdot x\right)}{{c}^{\left(\frac{2}{2}\right)} \cdot \left(\left(x \cdot \left(x \cdot {s}^{\left(\frac{2}{2}\right)}\right)\right) \cdot \left({c}^{\left(\frac{2}{2}\right)} \cdot {s}^{\left(\frac{2}{2}\right)}\right)\right)}\\
\mathbf{elif}\;{s}^{2} \leq 1.30542547701306 \cdot 10^{+247}:\\
\;\;\;\;\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{\frac{\cos \left(2 \cdot x\right)}{{c}^{\left(\frac{2}{2}\right)}}}{x \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)}\\
\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)) <= 0.0)) {
VAR = (((double) cos(((double) (2.0 * x)))) / ((double) (((double) pow(c, (2.0 / 2.0))) * ((double) (((double) (x * ((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 {
double VAR_1;
if ((((double) pow(s, 2.0)) <= 1.30542547701306e+247)) {
VAR_1 = (((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_1 = ((((double) cos(((double) (2.0 * x)))) / ((double) pow(c, (2.0 / 2.0)))) / ((double) (x * ((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))))))))));
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus x



Bits error versus c



Bits error versus s
Results
if (pow s 2.0) < 0.0Initial program 64.0
rmApplied sqr-pow64.0
Applied associate-*l*64.0
Simplified64.0
rmApplied sqr-pow64.0
Applied associate-*r*36.8
rmApplied associate-*l*21.2
Simplified21.2
rmApplied associate-*r*18.9
if 0.0 < (pow s 2.0) < 1.30542547701306e247Initial program 21.6
rmApplied sqr-pow21.6
Applied associate-*l*15.7
Simplified8.3
rmApplied sqr-pow8.3
Applied associate-*r*8.1
rmApplied associate-*l*7.8
Simplified7.8
rmApplied associate-*r*3.9
Simplified3.9
if 1.30542547701306e247 < (pow s 2.0) Initial program 24.0
rmApplied sqr-pow24.0
Applied associate-*l*19.2
Simplified19.0
rmApplied sqr-pow19.0
Applied associate-*r*13.1
rmApplied associate-*l*8.5
Simplified8.5
rmApplied associate-/r*8.4
Final simplification7.6
herbie shell --seed 2020199
(FPCore (x c s)
:name "mixedcos"
:precision binary64
(/ (cos (* 2.0 x)) (* (pow c 2.0) (* (* x (pow s 2.0)) x))))