Quotient of sum of exps

Average Error: 0.6 → 0.6

Time: 9.6s

Precision: 64

• Falling back on racket egraph because egg-herbie package not installed

• could not determine a ground truth for program body

  1. a = -2.3686796720332175e+206
  2. b = -5.202796025710632e+132

Math

\[\frac{e^{a}}{e^{a} + e^{b}}\]

↓

\[\frac{e^{a}}{\mathsf{fma}\left(\sqrt{e^{a}}, \sqrt{e^{a}}, e^{b}\right)}\]

C

double f(double a, double b) {
        double r123637 = a;
        double r123638 = exp(r123637);
        double r123639 = b;
        double r123640 = exp(r123639);
        double r123641 = r123638 + r123640;
        double r123642 = r123638 / r123641;
        return r123642;
}

↓

double f(double a, double b) {
        double r123643 = a;
        double r123644 = exp(r123643);
        double r123645 = sqrt(r123644);
        double r123646 = b;
        double r123647 = exp(r123646);
        double r123648 = fma(r123645, r123645, r123647);
        double r123649 = r123644 / r123648;
        return r123649;
}

Error

Bits error versus a

Bits error versus b

Target

Original	0.6
Target	0.0
Herbie	0.6

\[\frac{1}{1 + e^{b - a}}\]

Derivation

1. Initial program 0.6

   \[\frac{e^{a}}{e^{a} + e^{b}}\]
2. Using strategy rm

3. Applied add-sqr-sqrt 0.6

   \[\leadsto \frac{e^{a}}{\color{blue}{\sqrt{e^{a}} \cdot \sqrt{e^{a}}} + e^{b}}\]

4. Applied fma-def 0.6

   \[\leadsto \frac{e^{a}}{\color{blue}{\mathsf{fma}\left(\sqrt{e^{a}}, \sqrt{e^{a}}, e^{b}\right)}}\]

5. Final simplification 0.6

   \[\leadsto \frac{e^{a}}{\mathsf{fma}\left(\sqrt{e^{a}}, \sqrt{e^{a}}, e^{b}\right)}\]

Reproduce

herbie shell --seed 2020047 +o rules:numerics
(FPCore (a b)
  :name "Quotient of sum of exps"
  :precision binary64

  :herbie-target
  (/ 1 (+ 1 (exp (- b a))))

  (/ (exp a) (+ (exp a) (exp b))))