Average Error: 0.6 → 0.5
Time: 5.5s
Precision: 64
\[\sqrt{\left(\frac{\left(re \cdot re\right)}{\left(im \cdot im\right)}\right)}\]
↓
\[\sqrt{\left(\mathsf{qma}\left(\left(\left(re \cdot re\right)\right), im, im\right)\right)}\]
\sqrt{\left(\frac{\left(re \cdot re\right)}{\left(im \cdot im\right)}\right)}
↓
\sqrt{\left(\mathsf{qma}\left(\left(\left(re \cdot re\right)\right), im, im\right)\right)}
double f(double re, double im) {
        double r3220495 = re;
        double r3220496 = r3220495 * r3220495;
        double r3220497 = im;
        double r3220498 = r3220497 * r3220497;
        double r3220499 = r3220496 + r3220498;
        double r3220500 = sqrt(r3220499);
        return r3220500;
}
↓
double f(double re, double im) {
        double r3220501 = re;
        double r3220502 = r3220501 * r3220501;
        double r3220503 = /*Error: no posit support in C */;
        double r3220504 = im;
        double r3220505 = /*Error: no posit support in C */;
        double r3220506 = /*Error: no posit support in C */;
        double r3220507 = sqrt(r3220506);
        return r3220507;
}

Error

Bits error versus re

Bits error versus im

Derivation

  1. Initial program 0.6

    \[\sqrt{\left(\frac{\left(re \cdot re\right)}{\left(im \cdot im\right)}\right)}\]
  2. Using strategy rm

  3. Applied introduce-quire0.6

    \[\leadsto \sqrt{\left(\frac{\color{blue}{\left(\left(\left(re \cdot re\right)\right)\right)}}{\left(im \cdot im\right)}\right)}\]

  4. Applied insert-quire-fdp-add0.5

    \[\leadsto \sqrt{\color{blue}{\left(\left(\mathsf{qma}\left(\left(\left(re \cdot re\right)\right), im, im\right)\right)\right)}}\]

  5. Final simplification0.5

    \[\leadsto \sqrt{\left(\mathsf{qma}\left(\left(\left(re \cdot re\right)\right), im, im\right)\right)}\]

Reproduce

herbie shell --seed 2019151 
(FPCore (re im)
  :name "math.abs on complex"
  (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))))