(log (/ (+ (* (exp a) an) (* (exp c) cn)) (+ an cn)))

Average Error: 52.0 → 52.0

Time: 5.2s

Precision: binary64

• could not determine a ground truth for program body

  1. a = -9.426156606175285e+27
  2. an = -2.6305831255475505e-188
  3. c = -4.6986329538044373e+55
  4. cn = 7.960518986233747e+131

Math

\[\log \left(\frac{e^{a} \cdot an + e^{c} \cdot cn}{an + cn}\right)\]

↓

\[\log \left(\frac{e^{a} \cdot an + e^{c} \cdot cn}{an + cn}\right)\]

C

double code(double a, double an, double c, double cn) {
	return ((double) log(((double) (((double) (((double) (((double) exp(a)) * an)) + ((double) (((double) exp(c)) * cn)))) / ((double) (an + cn))))));
}

↓

double code(double a, double an, double c, double cn) {
	return ((double) log(((double) (((double) (((double) (((double) exp(a)) * an)) + ((double) (((double) exp(c)) * cn)))) / ((double) (an + cn))))));
}

Error

Bits error versus a

Bits error versus an

Bits error versus c

Bits error versus cn

Derivation

1. Initial program 52.0

   \[\log \left(\frac{e^{a} \cdot an + e^{c} \cdot cn}{an + cn}\right)\]

2. Final simplification 52.0

   \[\leadsto \log \left(\frac{e^{a} \cdot an + e^{c} \cdot cn}{an + cn}\right)\]

Reproduce

herbie shell --seed 2020153 
(FPCore (a an c cn)
  :name "(log (/ (+ (* (exp a) an) (* (exp c) cn)) (+ an cn)))"
  :precision binary64
  (log (/ (+ (* (exp a) an) (* (exp c) cn)) (+ an cn))))