Exception for Harley's example

• could not determine a ground truth for program body

  1. c_p = 1524059487.3410294
  2. c_n = 7.908113942053237e-177
  3. t = 2.816439611861168e-154
  4. s = -1.0385456294046283e-09

\[0.0 \lt c_p \land 0.0 \lt c_n\]

Math

\[\frac{{\left(\frac{1}{1 + e^{-s}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-s}}\right)}^{c_n}}{{\left(\frac{1}{1 + e^{-t}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-t}}\right)}^{c_n}}\]

C

double code(double c_p, double c_n, double t, double s) {
	return ((pow((1.0 / (1.0 + exp(-s))), c_p) * pow((1.0 - (1.0 / (1.0 + exp(-s)))), c_n)) / (pow((1.0 / (1.0 + exp(-t))), c_p) * pow((1.0 - (1.0 / (1.0 + exp(-t)))), c_n)));
}

Reproduce

herbie shell --seed 2020091 +o rules:numerics
(FPCore (c_p c_n t s)
  :name "Harley's example"
  :precision binary64
  :pre (and (< 0.0 c_p) (< 0.0 c_n))

  :herbie-target
  (* (pow (/ (+ 1 (exp (- t))) (+ 1 (exp (- s)))) c_p) (pow (/ (+ 1 (exp t)) (+ 1 (exp s))) c_n))

  (/ (* (pow (/ 1 (+ 1 (exp (- s)))) c_p) (pow (- 1 (/ 1 (+ 1 (exp (- s))))) c_n)) (* (pow (/ 1 (+ 1 (exp (- t)))) c_p) (pow (- 1 (/ 1 (+ 1 (exp (- t))))) c_n))))