| Alternative 1 | |
|---|---|
| Accuracy | 95.9% |
| Cost | 13440 |
(FPCore (x c s) :precision binary64 (/ (cos (* 2.0 x)) (* (pow c 2.0) (* (* x (pow s 2.0)) x))))
(FPCore (x c s)
:precision binary64
(let* ((t_0 (cos (+ x x))) (t_1 (* x (* c s))))
(if (<= c -5.2e+212)
(* (/ t_0 t_1) (/ 1.0 t_1))
(if (<= c -8e-205)
(/ t_0 (pow (* c (* x s)) 2.0))
(* t_0 (pow (* s (* c x)) -2.0))))))double code(double x, double c, double s) {
return cos((2.0 * x)) / (pow(c, 2.0) * ((x * pow(s, 2.0)) * x));
}
double code(double x, double c, double s) {
double t_0 = cos((x + x));
double t_1 = x * (c * s);
double tmp;
if (c <= -5.2e+212) {
tmp = (t_0 / t_1) * (1.0 / t_1);
} else if (c <= -8e-205) {
tmp = t_0 / pow((c * (x * s)), 2.0);
} else {
tmp = t_0 * pow((s * (c * x)), -2.0);
}
return tmp;
}
real(8) function code(x, c, s)
real(8), intent (in) :: x
real(8), intent (in) :: c
real(8), intent (in) :: s
code = cos((2.0d0 * x)) / ((c ** 2.0d0) * ((x * (s ** 2.0d0)) * x))
end function
real(8) function code(x, c, s)
real(8), intent (in) :: x
real(8), intent (in) :: c
real(8), intent (in) :: s
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = cos((x + x))
t_1 = x * (c * s)
if (c <= (-5.2d+212)) then
tmp = (t_0 / t_1) * (1.0d0 / t_1)
else if (c <= (-8d-205)) then
tmp = t_0 / ((c * (x * s)) ** 2.0d0)
else
tmp = t_0 * ((s * (c * x)) ** (-2.0d0))
end if
code = tmp
end function
public static double code(double x, double c, double s) {
return Math.cos((2.0 * x)) / (Math.pow(c, 2.0) * ((x * Math.pow(s, 2.0)) * x));
}
public static double code(double x, double c, double s) {
double t_0 = Math.cos((x + x));
double t_1 = x * (c * s);
double tmp;
if (c <= -5.2e+212) {
tmp = (t_0 / t_1) * (1.0 / t_1);
} else if (c <= -8e-205) {
tmp = t_0 / Math.pow((c * (x * s)), 2.0);
} else {
tmp = t_0 * Math.pow((s * (c * x)), -2.0);
}
return tmp;
}
def code(x, c, s): return math.cos((2.0 * x)) / (math.pow(c, 2.0) * ((x * math.pow(s, 2.0)) * x))
def code(x, c, s): t_0 = math.cos((x + x)) t_1 = x * (c * s) tmp = 0 if c <= -5.2e+212: tmp = (t_0 / t_1) * (1.0 / t_1) elif c <= -8e-205: tmp = t_0 / math.pow((c * (x * s)), 2.0) else: tmp = t_0 * math.pow((s * (c * x)), -2.0) return tmp
function code(x, c, s) return Float64(cos(Float64(2.0 * x)) / Float64((c ^ 2.0) * Float64(Float64(x * (s ^ 2.0)) * x))) end
function code(x, c, s) t_0 = cos(Float64(x + x)) t_1 = Float64(x * Float64(c * s)) tmp = 0.0 if (c <= -5.2e+212) tmp = Float64(Float64(t_0 / t_1) * Float64(1.0 / t_1)); elseif (c <= -8e-205) tmp = Float64(t_0 / (Float64(c * Float64(x * s)) ^ 2.0)); else tmp = Float64(t_0 * (Float64(s * Float64(c * x)) ^ -2.0)); end return tmp end
function tmp = code(x, c, s) tmp = cos((2.0 * x)) / ((c ^ 2.0) * ((x * (s ^ 2.0)) * x)); end
function tmp_2 = code(x, c, s) t_0 = cos((x + x)); t_1 = x * (c * s); tmp = 0.0; if (c <= -5.2e+212) tmp = (t_0 / t_1) * (1.0 / t_1); elseif (c <= -8e-205) tmp = t_0 / ((c * (x * s)) ^ 2.0); else tmp = t_0 * ((s * (c * x)) ^ -2.0); end tmp_2 = tmp; end
code[x_, c_, s_] := N[(N[Cos[N[(2.0 * x), $MachinePrecision]], $MachinePrecision] / N[(N[Power[c, 2.0], $MachinePrecision] * N[(N[(x * N[Power[s, 2.0], $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[x_, c_, s_] := Block[{t$95$0 = N[Cos[N[(x + x), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(x * N[(c * s), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[c, -5.2e+212], N[(N[(t$95$0 / t$95$1), $MachinePrecision] * N[(1.0 / t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, -8e-205], N[(t$95$0 / N[Power[N[(c * N[(x * s), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(t$95$0 * N[Power[N[(s * N[(c * x), $MachinePrecision]), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision]]]]]
\frac{\cos \left(2 \cdot x\right)}{{c}^{2} \cdot \left(\left(x \cdot {s}^{2}\right) \cdot x\right)}
\begin{array}{l}
t_0 := \cos \left(x + x\right)\\
t_1 := x \cdot \left(c \cdot s\right)\\
\mathbf{if}\;c \leq -5.2 \cdot 10^{+212}:\\
\;\;\;\;\frac{t_0}{t_1} \cdot \frac{1}{t_1}\\
\mathbf{elif}\;c \leq -8 \cdot 10^{-205}:\\
\;\;\;\;\frac{t_0}{{\left(c \cdot \left(x \cdot s\right)\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;t_0 \cdot {\left(s \cdot \left(c \cdot x\right)\right)}^{-2}\\
\end{array}
Results
if c < -5.1999999999999997e212Initial program 58.1%
Simplified96.0%
[Start]58.1 | \[ \frac{\cos \left(2 \cdot x\right)}{{c}^{2} \cdot \left(\left(x \cdot {s}^{2}\right) \cdot x\right)}
\] |
|---|---|
*-commutative [=>]58.1 | \[ \frac{\cos \left(2 \cdot x\right)}{{c}^{2} \cdot \left(\color{blue}{\left({s}^{2} \cdot x\right)} \cdot x\right)}
\] |
associate-*l* [=>]52.5 | \[ \frac{\cos \left(2 \cdot x\right)}{{c}^{2} \cdot \color{blue}{\left({s}^{2} \cdot \left(x \cdot x\right)\right)}}
\] |
associate-*r* [=>]52.6 | \[ \frac{\cos \left(2 \cdot x\right)}{\color{blue}{\left({c}^{2} \cdot {s}^{2}\right) \cdot \left(x \cdot x\right)}}
\] |
*-commutative [=>]52.6 | \[ \frac{\cos \left(2 \cdot x\right)}{\color{blue}{\left(x \cdot x\right) \cdot \left({c}^{2} \cdot {s}^{2}\right)}}
\] |
unpow2 [=>]52.6 | \[ \frac{\cos \left(2 \cdot x\right)}{\left(x \cdot x\right) \cdot \left(\color{blue}{\left(c \cdot c\right)} \cdot {s}^{2}\right)}
\] |
unpow2 [=>]52.6 | \[ \frac{\cos \left(2 \cdot x\right)}{\left(x \cdot x\right) \cdot \left(\left(c \cdot c\right) \cdot \color{blue}{\left(s \cdot s\right)}\right)}
\] |
unswap-sqr [=>]76.0 | \[ \frac{\cos \left(2 \cdot x\right)}{\left(x \cdot x\right) \cdot \color{blue}{\left(\left(c \cdot s\right) \cdot \left(c \cdot s\right)\right)}}
\] |
unswap-sqr [=>]96.0 | \[ \frac{\cos \left(2 \cdot x\right)}{\color{blue}{\left(x \cdot \left(c \cdot s\right)\right) \cdot \left(x \cdot \left(c \cdot s\right)\right)}}
\] |
Applied egg-rr96.3%
if -5.1999999999999997e212 < c < -8e-205Initial program 60.2%
Simplified77.6%
[Start]60.2 | \[ \frac{\cos \left(2 \cdot x\right)}{{c}^{2} \cdot \left(\left(x \cdot {s}^{2}\right) \cdot x\right)}
\] |
|---|---|
associate-*r* [=>]64.2 | \[ \frac{\cos \left(2 \cdot x\right)}{\color{blue}{\left({c}^{2} \cdot \left(x \cdot {s}^{2}\right)\right) \cdot x}}
\] |
*-commutative [=>]64.2 | \[ \frac{\cos \left(2 \cdot x\right)}{\color{blue}{x \cdot \left({c}^{2} \cdot \left(x \cdot {s}^{2}\right)\right)}}
\] |
*-commutative [=>]64.2 | \[ \frac{\cos \left(2 \cdot x\right)}{x \cdot \left({c}^{2} \cdot \color{blue}{\left({s}^{2} \cdot x\right)}\right)}
\] |
associate-*r* [=>]60.3 | \[ \frac{\cos \left(2 \cdot x\right)}{x \cdot \color{blue}{\left(\left({c}^{2} \cdot {s}^{2}\right) \cdot x\right)}}
\] |
*-commutative [=>]60.3 | \[ \frac{\cos \left(2 \cdot x\right)}{x \cdot \color{blue}{\left(x \cdot \left({c}^{2} \cdot {s}^{2}\right)\right)}}
\] |
unpow2 [=>]60.3 | \[ \frac{\cos \left(2 \cdot x\right)}{x \cdot \left(x \cdot \left(\color{blue}{\left(c \cdot c\right)} \cdot {s}^{2}\right)\right)}
\] |
unpow2 [=>]60.3 | \[ \frac{\cos \left(2 \cdot x\right)}{x \cdot \left(x \cdot \left(\left(c \cdot c\right) \cdot \color{blue}{\left(s \cdot s\right)}\right)\right)}
\] |
unswap-sqr [=>]77.6 | \[ \frac{\cos \left(2 \cdot x\right)}{x \cdot \left(x \cdot \color{blue}{\left(\left(c \cdot s\right) \cdot \left(c \cdot s\right)\right)}\right)}
\] |
Taylor expanded in x around inf 54.7%
Simplified98.3%
[Start]54.7 | \[ \frac{\cos \left(2 \cdot x\right)}{{s}^{2} \cdot \left({c}^{2} \cdot {x}^{2}\right)}
\] |
|---|---|
count-2 [<=]54.7 | \[ \frac{\cos \color{blue}{\left(x + x\right)}}{{s}^{2} \cdot \left({c}^{2} \cdot {x}^{2}\right)}
\] |
associate-*r* [=>]54.1 | \[ \frac{\cos \left(x + x\right)}{\color{blue}{\left({s}^{2} \cdot {c}^{2}\right) \cdot {x}^{2}}}
\] |
associate-/r* [=>]54.0 | \[ \color{blue}{\frac{\frac{\cos \left(x + x\right)}{{s}^{2} \cdot {c}^{2}}}{{x}^{2}}}
\] |
*-commutative [<=]54.0 | \[ \frac{\frac{\cos \left(x + x\right)}{\color{blue}{{c}^{2} \cdot {s}^{2}}}}{{x}^{2}}
\] |
unpow2 [=>]54.0 | \[ \frac{\frac{\cos \left(x + x\right)}{\color{blue}{\left(c \cdot c\right)} \cdot {s}^{2}}}{{x}^{2}}
\] |
unpow2 [=>]54.0 | \[ \frac{\frac{\cos \left(x + x\right)}{\left(c \cdot c\right) \cdot \color{blue}{\left(s \cdot s\right)}}}{{x}^{2}}
\] |
swap-sqr [<=]67.4 | \[ \frac{\frac{\cos \left(x + x\right)}{\color{blue}{\left(c \cdot s\right) \cdot \left(c \cdot s\right)}}}{{x}^{2}}
\] |
unpow2 [<=]67.4 | \[ \frac{\frac{\cos \left(x + x\right)}{\color{blue}{{\left(c \cdot s\right)}^{2}}}}{{x}^{2}}
\] |
associate-/l/ [=>]67.6 | \[ \color{blue}{\frac{\cos \left(x + x\right)}{{x}^{2} \cdot {\left(c \cdot s\right)}^{2}}}
\] |
unpow2 [=>]67.6 | \[ \frac{\cos \left(x + x\right)}{\color{blue}{\left(x \cdot x\right)} \cdot {\left(c \cdot s\right)}^{2}}
\] |
unpow2 [=>]67.6 | \[ \frac{\cos \left(x + x\right)}{\left(x \cdot x\right) \cdot \color{blue}{\left(\left(c \cdot s\right) \cdot \left(c \cdot s\right)\right)}}
\] |
swap-sqr [<=]96.1 | \[ \frac{\cos \left(x + x\right)}{\color{blue}{\left(x \cdot \left(c \cdot s\right)\right) \cdot \left(x \cdot \left(c \cdot s\right)\right)}}
\] |
unpow2 [<=]96.1 | \[ \frac{\cos \left(x + x\right)}{\color{blue}{{\left(x \cdot \left(c \cdot s\right)\right)}^{2}}}
\] |
*-commutative [=>]96.1 | \[ \frac{\cos \left(x + x\right)}{{\color{blue}{\left(\left(c \cdot s\right) \cdot x\right)}}^{2}}
\] |
associate-*l* [=>]98.3 | \[ \frac{\cos \left(x + x\right)}{{\color{blue}{\left(c \cdot \left(s \cdot x\right)\right)}}^{2}}
\] |
if -8e-205 < c Initial program 46.0%
Simplified94.9%
[Start]46.0 | \[ \frac{\cos \left(2 \cdot x\right)}{{c}^{2} \cdot \left(\left(x \cdot {s}^{2}\right) \cdot x\right)}
\] |
|---|---|
*-commutative [=>]46.0 | \[ \frac{\cos \left(2 \cdot x\right)}{{c}^{2} \cdot \left(\color{blue}{\left({s}^{2} \cdot x\right)} \cdot x\right)}
\] |
associate-*l* [=>]41.5 | \[ \frac{\cos \left(2 \cdot x\right)}{{c}^{2} \cdot \color{blue}{\left({s}^{2} \cdot \left(x \cdot x\right)\right)}}
\] |
associate-*r* [=>]41.8 | \[ \frac{\cos \left(2 \cdot x\right)}{\color{blue}{\left({c}^{2} \cdot {s}^{2}\right) \cdot \left(x \cdot x\right)}}
\] |
*-commutative [=>]41.8 | \[ \frac{\cos \left(2 \cdot x\right)}{\color{blue}{\left(x \cdot x\right) \cdot \left({c}^{2} \cdot {s}^{2}\right)}}
\] |
unpow2 [=>]41.8 | \[ \frac{\cos \left(2 \cdot x\right)}{\left(x \cdot x\right) \cdot \left(\color{blue}{\left(c \cdot c\right)} \cdot {s}^{2}\right)}
\] |
unpow2 [=>]41.8 | \[ \frac{\cos \left(2 \cdot x\right)}{\left(x \cdot x\right) \cdot \left(\left(c \cdot c\right) \cdot \color{blue}{\left(s \cdot s\right)}\right)}
\] |
unswap-sqr [=>]66.5 | \[ \frac{\cos \left(2 \cdot x\right)}{\left(x \cdot x\right) \cdot \color{blue}{\left(\left(c \cdot s\right) \cdot \left(c \cdot s\right)\right)}}
\] |
unswap-sqr [=>]94.9 | \[ \frac{\cos \left(2 \cdot x\right)}{\color{blue}{\left(x \cdot \left(c \cdot s\right)\right) \cdot \left(x \cdot \left(c \cdot s\right)\right)}}
\] |
Taylor expanded in x around inf 41.5%
Simplified95.1%
[Start]41.5 | \[ \frac{\cos \left(2 \cdot x\right)}{{c}^{2} \cdot \left({s}^{2} \cdot {x}^{2}\right)}
\] |
|---|---|
count-2 [<=]41.5 | \[ \frac{\cos \color{blue}{\left(x + x\right)}}{{c}^{2} \cdot \left({s}^{2} \cdot {x}^{2}\right)}
\] |
associate-*r* [=>]41.8 | \[ \frac{\cos \left(x + x\right)}{\color{blue}{\left({c}^{2} \cdot {s}^{2}\right) \cdot {x}^{2}}}
\] |
associate-/r* [=>]42.0 | \[ \color{blue}{\frac{\frac{\cos \left(x + x\right)}{{c}^{2} \cdot {s}^{2}}}{{x}^{2}}}
\] |
unpow2 [=>]42.0 | \[ \frac{\frac{\cos \left(x + x\right)}{\color{blue}{\left(c \cdot c\right)} \cdot {s}^{2}}}{{x}^{2}}
\] |
unpow2 [=>]42.0 | \[ \frac{\frac{\cos \left(x + x\right)}{\left(c \cdot c\right) \cdot \color{blue}{\left(s \cdot s\right)}}}{{x}^{2}}
\] |
swap-sqr [<=]66.3 | \[ \frac{\frac{\cos \left(x + x\right)}{\color{blue}{\left(c \cdot s\right) \cdot \left(c \cdot s\right)}}}{{x}^{2}}
\] |
*-lft-identity [<=]66.3 | \[ \frac{\frac{\color{blue}{1 \cdot \cos \left(x + x\right)}}{\left(c \cdot s\right) \cdot \left(c \cdot s\right)}}{{x}^{2}}
\] |
times-frac [=>]66.4 | \[ \frac{\color{blue}{\frac{1}{c \cdot s} \cdot \frac{\cos \left(x + x\right)}{c \cdot s}}}{{x}^{2}}
\] |
associate-*l/ [<=]72.8 | \[ \color{blue}{\frac{\frac{1}{c \cdot s}}{{x}^{2}} \cdot \frac{\cos \left(x + x\right)}{c \cdot s}}
\] |
unpow2 [=>]72.8 | \[ \frac{\frac{1}{c \cdot s}}{\color{blue}{x \cdot x}} \cdot \frac{\cos \left(x + x\right)}{c \cdot s}
\] |
associate-/r* [<=]72.8 | \[ \color{blue}{\frac{1}{\left(c \cdot s\right) \cdot \left(x \cdot x\right)}} \cdot \frac{\cos \left(x + x\right)}{c \cdot s}
\] |
associate-*r* [=>]87.3 | \[ \frac{1}{\color{blue}{\left(\left(c \cdot s\right) \cdot x\right) \cdot x}} \cdot \frac{\cos \left(x + x\right)}{c \cdot s}
\] |
*-commutative [<=]87.3 | \[ \frac{1}{\color{blue}{\left(x \cdot \left(c \cdot s\right)\right)} \cdot x} \cdot \frac{\cos \left(x + x\right)}{c \cdot s}
\] |
associate-/r* [=>]87.7 | \[ \color{blue}{\frac{\frac{1}{x \cdot \left(c \cdot s\right)}}{x}} \cdot \frac{\cos \left(x + x\right)}{c \cdot s}
\] |
Final simplification97.0%
| Alternative 1 | |
|---|---|
| Accuracy | 95.9% |
| Cost | 13440 |
| Alternative 2 | |
|---|---|
| Accuracy | 95.3% |
| Cost | 7880 |
| Alternative 3 | |
|---|---|
| Accuracy | 95.3% |
| Cost | 7753 |
| Alternative 4 | |
|---|---|
| Accuracy | 84.9% |
| Cost | 7625 |
| Alternative 5 | |
|---|---|
| Accuracy | 89.6% |
| Cost | 7625 |
| Alternative 6 | |
|---|---|
| Accuracy | 93.8% |
| Cost | 7625 |
| Alternative 7 | |
|---|---|
| Accuracy | 95.0% |
| Cost | 7625 |
| Alternative 8 | |
|---|---|
| Accuracy | 88.0% |
| Cost | 7624 |
| Alternative 9 | |
|---|---|
| Accuracy | 95.7% |
| Cost | 7360 |
| Alternative 10 | |
|---|---|
| Accuracy | 73.2% |
| Cost | 6912 |
| Alternative 11 | |
|---|---|
| Accuracy | 73.3% |
| Cost | 6784 |
| Alternative 12 | |
|---|---|
| Accuracy | 72.1% |
| Cost | 1097 |
| Alternative 13 | |
|---|---|
| Accuracy | 71.3% |
| Cost | 964 |
| Alternative 14 | |
|---|---|
| Accuracy | 73.3% |
| Cost | 964 |
| Alternative 15 | |
|---|---|
| Accuracy | 73.9% |
| Cost | 964 |
| Alternative 16 | |
|---|---|
| Accuracy | 70.4% |
| Cost | 832 |
| Alternative 17 | |
|---|---|
| Accuracy | 73.2% |
| Cost | 832 |
herbie shell --seed 2023122
(FPCore (x c s)
:name "mixedcos"
:precision binary64
(/ (cos (* 2.0 x)) (* (pow c 2.0) (* (* x (pow s 2.0)) x))))