| Alternative 1 | |
|---|---|
| Error | 0.2 |
| Cost | 1352 |
(FPCore (x.re x.im) :precision binary64 (- (* (- (* x.re x.re) (* x.im x.im)) x.re) (* (+ (* x.re x.im) (* x.im x.re)) x.im)))
(FPCore (x.re x.im)
:precision binary64
(if (<= x.im -5.5e+124)
(* (* x.im (* x.re -3.0)) x.im)
(if (<= x.im 5.2e+116)
(-
(* (+ x.re x.im) (* (- x.re x.im) x.re))
(* x.re (* (+ x.im x.im) x.im)))
(* (* x.re x.im) (* -3.0 x.im)))))double code(double x_46_re, double x_46_im) {
return (((x_46_re * x_46_re) - (x_46_im * x_46_im)) * x_46_re) - (((x_46_re * x_46_im) + (x_46_im * x_46_re)) * x_46_im);
}
double code(double x_46_re, double x_46_im) {
double tmp;
if (x_46_im <= -5.5e+124) {
tmp = (x_46_im * (x_46_re * -3.0)) * x_46_im;
} else if (x_46_im <= 5.2e+116) {
tmp = ((x_46_re + x_46_im) * ((x_46_re - x_46_im) * x_46_re)) - (x_46_re * ((x_46_im + x_46_im) * x_46_im));
} else {
tmp = (x_46_re * x_46_im) * (-3.0 * x_46_im);
}
return tmp;
}
real(8) function code(x_46re, x_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
code = (((x_46re * x_46re) - (x_46im * x_46im)) * x_46re) - (((x_46re * x_46im) + (x_46im * x_46re)) * x_46im)
end function
real(8) function code(x_46re, x_46im)
real(8), intent (in) :: x_46re
real(8), intent (in) :: x_46im
real(8) :: tmp
if (x_46im <= (-5.5d+124)) then
tmp = (x_46im * (x_46re * (-3.0d0))) * x_46im
else if (x_46im <= 5.2d+116) then
tmp = ((x_46re + x_46im) * ((x_46re - x_46im) * x_46re)) - (x_46re * ((x_46im + x_46im) * x_46im))
else
tmp = (x_46re * x_46im) * ((-3.0d0) * x_46im)
end if
code = tmp
end function
public static double code(double x_46_re, double x_46_im) {
return (((x_46_re * x_46_re) - (x_46_im * x_46_im)) * x_46_re) - (((x_46_re * x_46_im) + (x_46_im * x_46_re)) * x_46_im);
}
public static double code(double x_46_re, double x_46_im) {
double tmp;
if (x_46_im <= -5.5e+124) {
tmp = (x_46_im * (x_46_re * -3.0)) * x_46_im;
} else if (x_46_im <= 5.2e+116) {
tmp = ((x_46_re + x_46_im) * ((x_46_re - x_46_im) * x_46_re)) - (x_46_re * ((x_46_im + x_46_im) * x_46_im));
} else {
tmp = (x_46_re * x_46_im) * (-3.0 * x_46_im);
}
return tmp;
}
def code(x_46_re, x_46_im): return (((x_46_re * x_46_re) - (x_46_im * x_46_im)) * x_46_re) - (((x_46_re * x_46_im) + (x_46_im * x_46_re)) * x_46_im)
def code(x_46_re, x_46_im): tmp = 0 if x_46_im <= -5.5e+124: tmp = (x_46_im * (x_46_re * -3.0)) * x_46_im elif x_46_im <= 5.2e+116: tmp = ((x_46_re + x_46_im) * ((x_46_re - x_46_im) * x_46_re)) - (x_46_re * ((x_46_im + x_46_im) * x_46_im)) else: tmp = (x_46_re * x_46_im) * (-3.0 * x_46_im) return tmp
function code(x_46_re, x_46_im) return Float64(Float64(Float64(Float64(x_46_re * x_46_re) - Float64(x_46_im * x_46_im)) * x_46_re) - Float64(Float64(Float64(x_46_re * x_46_im) + Float64(x_46_im * x_46_re)) * x_46_im)) end
function code(x_46_re, x_46_im) tmp = 0.0 if (x_46_im <= -5.5e+124) tmp = Float64(Float64(x_46_im * Float64(x_46_re * -3.0)) * x_46_im); elseif (x_46_im <= 5.2e+116) tmp = Float64(Float64(Float64(x_46_re + x_46_im) * Float64(Float64(x_46_re - x_46_im) * x_46_re)) - Float64(x_46_re * Float64(Float64(x_46_im + x_46_im) * x_46_im))); else tmp = Float64(Float64(x_46_re * x_46_im) * Float64(-3.0 * x_46_im)); end return tmp end
function tmp = code(x_46_re, x_46_im) tmp = (((x_46_re * x_46_re) - (x_46_im * x_46_im)) * x_46_re) - (((x_46_re * x_46_im) + (x_46_im * x_46_re)) * x_46_im); end
function tmp_2 = code(x_46_re, x_46_im) tmp = 0.0; if (x_46_im <= -5.5e+124) tmp = (x_46_im * (x_46_re * -3.0)) * x_46_im; elseif (x_46_im <= 5.2e+116) tmp = ((x_46_re + x_46_im) * ((x_46_re - x_46_im) * x_46_re)) - (x_46_re * ((x_46_im + x_46_im) * x_46_im)); else tmp = (x_46_re * x_46_im) * (-3.0 * x_46_im); end tmp_2 = tmp; end
code[x$46$re_, x$46$im_] := N[(N[(N[(N[(x$46$re * x$46$re), $MachinePrecision] - N[(x$46$im * x$46$im), $MachinePrecision]), $MachinePrecision] * x$46$re), $MachinePrecision] - N[(N[(N[(x$46$re * x$46$im), $MachinePrecision] + N[(x$46$im * x$46$re), $MachinePrecision]), $MachinePrecision] * x$46$im), $MachinePrecision]), $MachinePrecision]
code[x$46$re_, x$46$im_] := If[LessEqual[x$46$im, -5.5e+124], N[(N[(x$46$im * N[(x$46$re * -3.0), $MachinePrecision]), $MachinePrecision] * x$46$im), $MachinePrecision], If[LessEqual[x$46$im, 5.2e+116], N[(N[(N[(x$46$re + x$46$im), $MachinePrecision] * N[(N[(x$46$re - x$46$im), $MachinePrecision] * x$46$re), $MachinePrecision]), $MachinePrecision] - N[(x$46$re * N[(N[(x$46$im + x$46$im), $MachinePrecision] * x$46$im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x$46$re * x$46$im), $MachinePrecision] * N[(-3.0 * x$46$im), $MachinePrecision]), $MachinePrecision]]]
\left(x.re \cdot x.re - x.im \cdot x.im\right) \cdot x.re - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
\begin{array}{l}
\mathbf{if}\;x.im \leq -5.5 \cdot 10^{+124}:\\
\;\;\;\;\left(x.im \cdot \left(x.re \cdot -3\right)\right) \cdot x.im\\
\mathbf{elif}\;x.im \leq 5.2 \cdot 10^{+116}:\\
\;\;\;\;\left(x.re + x.im\right) \cdot \left(\left(x.re - x.im\right) \cdot x.re\right) - x.re \cdot \left(\left(x.im + x.im\right) \cdot x.im\right)\\
\mathbf{else}:\\
\;\;\;\;\left(x.re \cdot x.im\right) \cdot \left(-3 \cdot x.im\right)\\
\end{array}
Results
| Original | 7.2 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
if x.im < -5.49999999999999977e124Initial program 44.8
Simplified44.9
[Start]44.8 | \[ \left(x.re \cdot x.re - x.im \cdot x.im\right) \cdot x.re - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
\] |
|---|---|
rational.json-simplify-2 [=>]44.8 | \[ \color{blue}{x.re \cdot \left(x.re \cdot x.re - x.im \cdot x.im\right)} - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
\] |
rational.json-simplify-34 [=>]44.8 | \[ x.re \cdot \color{blue}{\left(\left(x.im + x.re\right) \cdot \left(x.re - x.im\right)\right)} - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
\] |
rational.json-simplify-43 [=>]0.4 | \[ \color{blue}{\left(x.im + x.re\right) \cdot \left(\left(x.re - x.im\right) \cdot x.re\right)} - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
\] |
rational.json-simplify-1 [=>]0.4 | \[ \color{blue}{\left(x.re + x.im\right)} \cdot \left(\left(x.re - x.im\right) \cdot x.re\right) - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
\] |
rational.json-simplify-2 [=>]0.4 | \[ \left(x.re + x.im\right) \cdot \left(\left(x.re - x.im\right) \cdot x.re\right) - \color{blue}{x.im \cdot \left(x.re \cdot x.im + x.im \cdot x.re\right)}
\] |
rational.json-simplify-51 [=>]0.4 | \[ \left(x.re + x.im\right) \cdot \left(\left(x.re - x.im\right) \cdot x.re\right) - x.im \cdot \color{blue}{\left(x.re \cdot \left(x.im + x.im\right)\right)}
\] |
rational.json-simplify-43 [=>]44.9 | \[ \left(x.re + x.im\right) \cdot \left(\left(x.re - x.im\right) \cdot x.re\right) - \color{blue}{x.re \cdot \left(\left(x.im + x.im\right) \cdot x.im\right)}
\] |
Taylor expanded in x.re around 0 44.9
Simplified44.9
[Start]44.9 | \[ \left(x.re + x.im\right) \cdot \left(-1 \cdot \left(x.re \cdot x.im\right)\right) - x.re \cdot \left(\left(x.im + x.im\right) \cdot x.im\right)
\] |
|---|---|
rational.json-simplify-43 [=>]44.9 | \[ \left(x.re + x.im\right) \cdot \color{blue}{\left(x.re \cdot \left(x.im \cdot -1\right)\right)} - x.re \cdot \left(\left(x.im + x.im\right) \cdot x.im\right)
\] |
rational.json-simplify-9 [=>]44.9 | \[ \left(x.re + x.im\right) \cdot \left(x.re \cdot \color{blue}{\left(-x.im\right)}\right) - x.re \cdot \left(\left(x.im + x.im\right) \cdot x.im\right)
\] |
Applied egg-rr0.4
Taylor expanded in x.re around 0 0.4
Simplified0.4
[Start]0.4 | \[ \left(\left(-1 \cdot x.im - 2 \cdot x.im\right) \cdot x.re\right) \cdot x.im
\] |
|---|---|
rational.json-simplify-2 [=>]0.4 | \[ \color{blue}{\left(x.re \cdot \left(-1 \cdot x.im - 2 \cdot x.im\right)\right)} \cdot x.im
\] |
rational.json-simplify-2 [=>]0.4 | \[ \left(x.re \cdot \left(\color{blue}{x.im \cdot -1} - 2 \cdot x.im\right)\right) \cdot x.im
\] |
rational.json-simplify-52 [=>]0.4 | \[ \left(x.re \cdot \color{blue}{\left(x.im \cdot \left(-1 - 2\right)\right)}\right) \cdot x.im
\] |
rational.json-simplify-2 [=>]0.4 | \[ \left(x.re \cdot \color{blue}{\left(\left(-1 - 2\right) \cdot x.im\right)}\right) \cdot x.im
\] |
rational.json-simplify-43 [<=]0.4 | \[ \color{blue}{\left(x.im \cdot \left(x.re \cdot \left(-1 - 2\right)\right)\right)} \cdot x.im
\] |
metadata-eval [=>]0.4 | \[ \left(x.im \cdot \left(x.re \cdot \color{blue}{-3}\right)\right) \cdot x.im
\] |
if -5.49999999999999977e124 < x.im < 5.19999999999999973e116Initial program 0.2
Simplified0.2
[Start]0.2 | \[ \left(x.re \cdot x.re - x.im \cdot x.im\right) \cdot x.re - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
\] |
|---|---|
rational.json-simplify-2 [=>]0.2 | \[ \color{blue}{x.re \cdot \left(x.re \cdot x.re - x.im \cdot x.im\right)} - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
\] |
rational.json-simplify-34 [=>]0.2 | \[ x.re \cdot \color{blue}{\left(\left(x.im + x.re\right) \cdot \left(x.re - x.im\right)\right)} - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
\] |
rational.json-simplify-43 [=>]0.2 | \[ \color{blue}{\left(x.im + x.re\right) \cdot \left(\left(x.re - x.im\right) \cdot x.re\right)} - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
\] |
rational.json-simplify-1 [=>]0.2 | \[ \color{blue}{\left(x.re + x.im\right)} \cdot \left(\left(x.re - x.im\right) \cdot x.re\right) - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
\] |
rational.json-simplify-2 [=>]0.2 | \[ \left(x.re + x.im\right) \cdot \left(\left(x.re - x.im\right) \cdot x.re\right) - \color{blue}{x.im \cdot \left(x.re \cdot x.im + x.im \cdot x.re\right)}
\] |
rational.json-simplify-51 [=>]0.2 | \[ \left(x.re + x.im\right) \cdot \left(\left(x.re - x.im\right) \cdot x.re\right) - x.im \cdot \color{blue}{\left(x.re \cdot \left(x.im + x.im\right)\right)}
\] |
rational.json-simplify-43 [=>]0.2 | \[ \left(x.re + x.im\right) \cdot \left(\left(x.re - x.im\right) \cdot x.re\right) - \color{blue}{x.re \cdot \left(\left(x.im + x.im\right) \cdot x.im\right)}
\] |
if 5.19999999999999973e116 < x.im Initial program 41.0
Simplified41.1
[Start]41.0 | \[ \left(x.re \cdot x.re - x.im \cdot x.im\right) \cdot x.re - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
\] |
|---|---|
rational.json-simplify-2 [=>]41.0 | \[ \color{blue}{x.re \cdot \left(x.re \cdot x.re - x.im \cdot x.im\right)} - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
\] |
rational.json-simplify-34 [=>]41.0 | \[ x.re \cdot \color{blue}{\left(\left(x.im + x.re\right) \cdot \left(x.re - x.im\right)\right)} - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
\] |
rational.json-simplify-43 [=>]0.4 | \[ \color{blue}{\left(x.im + x.re\right) \cdot \left(\left(x.re - x.im\right) \cdot x.re\right)} - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
\] |
rational.json-simplify-1 [=>]0.4 | \[ \color{blue}{\left(x.re + x.im\right)} \cdot \left(\left(x.re - x.im\right) \cdot x.re\right) - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im
\] |
rational.json-simplify-2 [=>]0.4 | \[ \left(x.re + x.im\right) \cdot \left(\left(x.re - x.im\right) \cdot x.re\right) - \color{blue}{x.im \cdot \left(x.re \cdot x.im + x.im \cdot x.re\right)}
\] |
rational.json-simplify-51 [=>]0.4 | \[ \left(x.re + x.im\right) \cdot \left(\left(x.re - x.im\right) \cdot x.re\right) - x.im \cdot \color{blue}{\left(x.re \cdot \left(x.im + x.im\right)\right)}
\] |
rational.json-simplify-43 [=>]41.1 | \[ \left(x.re + x.im\right) \cdot \left(\left(x.re - x.im\right) \cdot x.re\right) - \color{blue}{x.re \cdot \left(\left(x.im + x.im\right) \cdot x.im\right)}
\] |
Taylor expanded in x.re around 0 41.1
Simplified41.1
[Start]41.1 | \[ \left(x.re + x.im\right) \cdot \left(-1 \cdot \left(x.re \cdot x.im\right)\right) - x.re \cdot \left(\left(x.im + x.im\right) \cdot x.im\right)
\] |
|---|---|
rational.json-simplify-43 [=>]41.1 | \[ \left(x.re + x.im\right) \cdot \color{blue}{\left(x.re \cdot \left(x.im \cdot -1\right)\right)} - x.re \cdot \left(\left(x.im + x.im\right) \cdot x.im\right)
\] |
rational.json-simplify-9 [=>]41.1 | \[ \left(x.re + x.im\right) \cdot \left(x.re \cdot \color{blue}{\left(-x.im\right)}\right) - x.re \cdot \left(\left(x.im + x.im\right) \cdot x.im\right)
\] |
Applied egg-rr0.4
Taylor expanded in x.re around 0 0.4
Simplified0.4
[Start]0.4 | \[ \left(\left(-2 \cdot x.im - x.im\right) \cdot x.re\right) \cdot x.im
\] |
|---|---|
rational.json-simplify-2 [<=]0.4 | \[ \left(\left(\color{blue}{x.im \cdot -2} - x.im\right) \cdot x.re\right) \cdot x.im
\] |
rational.json-simplify-6 [<=]0.4 | \[ \left(\left(x.im \cdot -2 - \color{blue}{1 \cdot x.im}\right) \cdot x.re\right) \cdot x.im
\] |
rational.json-simplify-52 [=>]0.4 | \[ \left(\color{blue}{\left(x.im \cdot \left(-2 - 1\right)\right)} \cdot x.re\right) \cdot x.im
\] |
metadata-eval [=>]0.4 | \[ \left(\left(x.im \cdot \color{blue}{-3}\right) \cdot x.re\right) \cdot x.im
\] |
Applied egg-rr0.4
Simplified0.4
[Start]0.4 | \[ x.im \cdot \left(x.im \cdot \left(-3 \cdot x.re\right)\right) + 0
\] |
|---|---|
rational.json-simplify-4 [=>]0.4 | \[ \color{blue}{x.im \cdot \left(x.im \cdot \left(-3 \cdot x.re\right)\right)}
\] |
rational.json-simplify-43 [=>]0.4 | \[ x.im \cdot \color{blue}{\left(-3 \cdot \left(x.re \cdot x.im\right)\right)}
\] |
rational.json-simplify-2 [=>]0.4 | \[ x.im \cdot \color{blue}{\left(\left(x.re \cdot x.im\right) \cdot -3\right)}
\] |
rational.json-simplify-43 [=>]0.4 | \[ \color{blue}{\left(x.re \cdot x.im\right) \cdot \left(-3 \cdot x.im\right)}
\] |
Final simplification0.2
| Alternative 1 | |
|---|---|
| Error | 0.2 |
| Cost | 1352 |
| Alternative 2 | |
|---|---|
| Error | 0.2 |
| Cost | 1352 |
| Alternative 3 | |
|---|---|
| Error | 0.2 |
| Cost | 1216 |
| Alternative 4 | |
|---|---|
| Error | 0.2 |
| Cost | 1088 |
| Alternative 5 | |
|---|---|
| Error | 0.2 |
| Cost | 1088 |
| Alternative 6 | |
|---|---|
| Error | 0.3 |
| Cost | 968 |
| Alternative 7 | |
|---|---|
| Error | 18.9 |
| Cost | 448 |
| Alternative 8 | |
|---|---|
| Error | 18.9 |
| Cost | 448 |
| Alternative 9 | |
|---|---|
| Error | 18.9 |
| Cost | 448 |
| Alternative 10 | |
|---|---|
| Error | 46.0 |
| Cost | 384 |
herbie shell --seed 2023069
(FPCore (x.re x.im)
:name "math.cube on complex, real part"
:precision binary64
:herbie-target
(+ (* (* x.re x.re) (- x.re x.im)) (* (* x.re x.im) (- x.re (* 3.0 x.im))))
(- (* (- (* x.re x.re) (* x.im x.im)) x.re) (* (+ (* x.re x.im) (* x.im x.re)) x.im)))