_divideComplex, real part

Average Error: 26.3 → 25.5

Time: 4.6s

Precision: binary64

Math

\[\]

↓

\[\]

Error

Bits error versus x.re

Bits error versus x.im

Bits error versus y.re

Bits error versus y.im

Derivation

1. Split input into 2 regimes

2. if (/ (+ (* x.re y.re) (* x.im y.im)) (+ (* y.re y.re) (* y.im y.im))) < 1.34211888946586264e301

   1. Initial program 14.3

      \[\]
   2. Using strategy rm

   3. Applied add-sqr-sqrt 14.3

      \[\leadsto \]

   4. Applied associate-/r* 14.2

      \[\leadsto \]

   if 1.34211888946586264e301 < (/ (+ (* x.re y.re) (* x.im y.im)) (+ (* y.re y.re) (* y.im y.im)))

   1. Initial program 63.5

      \[\]
   2. Using strategy rm

   3. Applied add-sqr-sqrt 63.5

      \[\leadsto \]

   4. Applied associate-/r* 63.5

      \[\leadsto \]

   5. Taylor expanded around -inf 60.5

      \[\leadsto \]

   6. Simplified 60.5

      \[\leadsto \]
3. Recombined 2 regimes into one program.

4. Final simplification 25.5

   \[\leadsto \]

Reproduce

herbie shell --seed 2020191 
(FPCore (x.re x.im y.re y.im)
  :name "_divideComplex, real part"
  :precision binary64
  (/ (+ (* x.re y.re) (* x.im y.im)) (+ (* y.re y.re) (* y.im y.im))))