
(FPCore im_sqr (re im) :precision binary64 (+ (* re im) (* im re)))
(FPCore im_sqr (re im) :precision binary64 (* re (+ im im)))
double im_sqr(double re, double im) {
return (re * im) + (im * re);
}
double im_sqr(double re, double im) {
return re * (im + im);
}
real(8) function im_sqr(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
im_sqr = (re * im) + (im * re)
end function
real(8) function im_sqr(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
im_sqr = re * (im + im)
end function
public static double im_sqr(double re, double im) {
return (re * im) + (im * re);
}
public static double im_sqr(double re, double im) {
return re * (im + im);
}
def im_sqr(re, im): return (re * im) + (im * re)
def im_sqr(re, im): return re * (im + im)
function im_sqr(re, im) return Float64(Float64(re * im) + Float64(im * re)) end
function im_sqr(re, im) return Float64(re * Float64(im + im)) end
function tmp = im_sqr(re, im) tmp = (re * im) + (im * re); end
function tmp = im_sqr(re, im) tmp = re * (im + im); end
im$95$sqr[re_, im_] := N[(N[(re * im), $MachinePrecision] + N[(im * re), $MachinePrecision]), $MachinePrecision]
im$95$sqr[re_, im_] := N[(re * N[(im + im), $MachinePrecision]), $MachinePrecision]
re \cdot im + im \cdot re
re \cdot \left(im + im\right)
Herbie found 1 alternatives:
| Alternative | Accuracy | Speedup |
|---|
Results
Initial program 100.0%
Simplified100.0%
[Start]100.0% | \[ re \cdot im + im \cdot re
\] |
|---|---|
*-commutative [=>]100.0% | \[ \color{blue}{im \cdot re} + im \cdot re
\] |
distribute-rgt-out [=>]100.0% | \[ \color{blue}{re \cdot \left(im + im\right)}
\] |
Final simplification100.0%
herbie shell --seed 2023178
(FPCore im_sqr (re im)
:name "math.square on complex, imaginary part"
:precision binary64
(+ (* re im) (* im re)))