Average Error: 28.9 → 7.5
Time: 5.5m
Precision: 64
Internal Precision: 128
\[\sqrt{\left|\frac{a \cdot a - b \cdot b}{a \cdot a}\right|}\]
\[\sqrt{\left|\left(\frac{b}{a} + 1\right) \cdot \left(1 - \frac{b}{a}\right)\right|}\]

Error

Bits error versus a

Bits error versus b

Derivation

  1. Initial program 28.9

    \[\sqrt{\left|\frac{a \cdot a - b \cdot b}{a \cdot a}\right|}\]
  2. Simplified7.5

    \[\leadsto \color{blue}{\sqrt{\left|1 - \frac{b}{a} \cdot \frac{b}{a}\right|}}\]
  3. Using strategy rm
  4. Applied *-commutative7.5

    \[\leadsto \sqrt{\left|1 - \color{blue}{\frac{b}{a} \cdot \frac{b}{a}}\right|}\]
  5. Applied *-un-lft-identity7.5

    \[\leadsto \sqrt{\left|\color{blue}{1 \cdot 1} - \frac{b}{a} \cdot \frac{b}{a}\right|}\]
  6. Applied difference-of-squares7.5

    \[\leadsto \sqrt{\left|\color{blue}{\left(1 + \frac{b}{a}\right) \cdot \left(1 - \frac{b}{a}\right)}\right|}\]
  7. Final simplification7.5

    \[\leadsto \sqrt{\left|\left(\frac{b}{a} + 1\right) \cdot \left(1 - \frac{b}{a}\right)\right|}\]

Reproduce

herbie shell --seed 2019072 +o rules:numerics
(FPCore (a b)
  :name "Eccentricity of an ellipse"
  :pre (<= 0 b a 1)
  (sqrt (fabs (/ (- (* a a) (* b b)) (* a a)))))