Average Error: 0.0 → 0.0
Time: 1.1s
Precision: 64
\[re \cdot im + im \cdot re\]
\[im \cdot \left(re + re\right)\]
double f(double re, double im) {
        double r8320 = re;
        double r8321 = im;
        double r8322 = r8320 * r8321;
        double r8323 = r8321 * r8320;
        double r8324 = r8322 + r8323;
        return r8324;
}

double f(double re, double im) {
        double r8325 = im;
        double r8326 = re;
        double r8327 = r8326 + r8326;
        double r8328 = r8325 * r8327;
        return r8328;
}

re \cdot im + im \cdot re
im \cdot \left(re + re\right)

Error

Bits error versus re

Bits error versus im

Derivation

  1. Initial program 0.0

    \[re \cdot im + im \cdot re\]
  2. Simplified0.0

    \[\leadsto \color{blue}{im \cdot \left(re + re\right)}\]
  3. Final simplification0.0

    \[\leadsto im \cdot \left(re + re\right)\]

Reproduce

herbie shell --seed 2019101 
(FPCore (re im)
  :name "math.square on complex, imaginary part"
  (+ (* re im) (* im re)))