Exception for Jmat.Real.gamma, branch z less than 0.5

• could not determine a ground truth for program body

  1. z = -8.464915614153485e+164

Math

\[\frac{\pi}{\sin \left(\pi \cdot z\right)} \cdot \left(\left(\left(\sqrt{\pi \cdot 2} \cdot {\left(\left(\left(\left(1 - z\right) - 1\right) + 7\right) + 0.5\right)}^{\left(\left(\left(1 - z\right) - 1\right) + 0.5\right)}\right) \cdot e^{-\left(\left(\left(\left(1 - z\right) - 1\right) + 7\right) + 0.5\right)}\right) \cdot \left(\left(\left(\left(\left(\left(\left(\left(0.9999999999998099298181841732002794742584 + \frac{676.5203681218850988443591631948947906494}{\left(\left(1 - z\right) - 1\right) + 1}\right) + \frac{-1259.139216722402807135949842631816864014}{\left(\left(1 - z\right) - 1\right) + 2}\right) + \frac{771.3234287776531346025876700878143310547}{\left(\left(1 - z\right) - 1\right) + 3}\right) + \frac{-176.6150291621405870046146446838974952698}{\left(\left(1 - z\right) - 1\right) + 4}\right) + \frac{12.50734327868690520801919774385169148445}{\left(\left(1 - z\right) - 1\right) + 5}\right) + \frac{-0.1385710952657201178173096423051902092993}{\left(\left(1 - z\right) - 1\right) + 6}\right) + \frac{9.984369578019571583242346146658263705831 \cdot 10^{-6}}{\left(\left(1 - z\right) - 1\right) + 7}\right) + \frac{1.505632735149311617592788074479481785772 \cdot 10^{-7}}{\left(\left(1 - z\right) - 1\right) + 8}\right)\right)\]

C

double f(double z) {
        double r8500 = atan2(1.0, 0.0);
        double r8501 = z;
        double r8502 = r8500 * r8501;
        double r8503 = sin(r8502);
        double r8504 = r8500 / r8503;
        double r8505 = 2.0;
        double r8506 = r8500 * r8505;
        double r8507 = sqrt(r8506);
        double r8508 = 1.0;
        double r8509 = r8508 - r8501;
        double r8510 = r8509 - r8508;
        double r8511 = 7.0;
        double r8512 = r8510 + r8511;
        double r8513 = 0.5;
        double r8514 = r8512 + r8513;
        double r8515 = r8510 + r8513;
        double r8516 = pow(r8514, r8515);
        double r8517 = r8507 * r8516;
        double r8518 = -r8514;
        double r8519 = exp(r8518);
        double r8520 = r8517 * r8519;
        double r8521 = 0.9999999999998099;
        double r8522 = 676.5203681218851;
        double r8523 = r8510 + r8508;
        double r8524 = r8522 / r8523;
        double r8525 = r8521 + r8524;
        double r8526 = -1259.1392167224028;
        double r8527 = r8510 + r8505;
        double r8528 = r8526 / r8527;
        double r8529 = r8525 + r8528;
        double r8530 = 771.3234287776531;
        double r8531 = 3.0;
        double r8532 = r8510 + r8531;
        double r8533 = r8530 / r8532;
        double r8534 = r8529 + r8533;
        double r8535 = -176.6150291621406;
        double r8536 = 4.0;
        double r8537 = r8510 + r8536;
        double r8538 = r8535 / r8537;
        double r8539 = r8534 + r8538;
        double r8540 = 12.507343278686905;
        double r8541 = 5.0;
        double r8542 = r8510 + r8541;
        double r8543 = r8540 / r8542;
        double r8544 = r8539 + r8543;
        double r8545 = -0.13857109526572012;
        double r8546 = 6.0;
        double r8547 = r8510 + r8546;
        double r8548 = r8545 / r8547;
        double r8549 = r8544 + r8548;
        double r8550 = 9.984369578019572e-06;
        double r8551 = r8550 / r8512;
        double r8552 = r8549 + r8551;
        double r8553 = 1.5056327351493116e-07;
        double r8554 = 8.0;
        double r8555 = r8510 + r8554;
        double r8556 = r8553 / r8555;
        double r8557 = r8552 + r8556;
        double r8558 = r8520 * r8557;
        double r8559 = r8504 * r8558;
        return r8559;
}

Reproduce

Please include this information when filing a bug report:

herbie shell --seed 2019305 
(FPCore (z)
  :name "Jmat.Real.gamma, branch z less than 0.5"
  :precision binary64
  (* (/ PI (sin (* PI z))) (* (* (* (sqrt (* PI 2)) (pow (+ (+ (- (- 1 z) 1) 7) 0.5) (+ (- (- 1 z) 1) 0.5))) (exp (- (+ (+ (- (- 1 z) 1) 7) 0.5)))) (+ (+ (+ (+ (+ (+ (+ (+ 0.99999999999980993 (/ 676.520368121885099 (+ (- (- 1 z) 1) 1))) (/ -1259.13921672240281 (+ (- (- 1 z) 1) 2))) (/ 771.32342877765313 (+ (- (- 1 z) 1) 3))) (/ -176.615029162140587 (+ (- (- 1 z) 1) 4))) (/ 12.5073432786869052 (+ (- (- 1 z) 1) 5))) (/ -0.138571095265720118 (+ (- (- 1 z) 1) 6))) (/ 9.98436957801957158e-6 (+ (- (- 1 z) 1) 7))) (/ 1.50563273514931162e-7 (+ (- (- 1 z) 1) 8))))))

Backtrace

exact?: contract violation expected: number? given: 'h1		L	C
egg-parsed->expr	/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/eggmath.rkt	65	0
for-loop	/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/eggmath.rkt	70	6
egg-parsed->expr	/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/eggmath.rkt	65	0
for-loop	/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/eggmath.rkt	70	6
egg-parsed->expr	/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/eggmath.rkt	65	0
for-loop	/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/eggmath.rkt	70	6
egg-parsed->expr	/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/eggmath.rkt	65	0
for-loop	/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/eggmath.rkt	70	6
egg-parsed->expr	/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/eggmath.rkt	65	0
(unnamed)	/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/simplify.rkt	49	8
egraph-add-exprs	/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/eggmath.rkt	122	0
egraph-run	/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/eggmath.rkt	55	0
(unnamed)	/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/simplify.rkt	38	0
simplify!	/data/pavpan/nightlies/herbie/eggmath-testversion/src/mainloop.rkt	223	0
run-improve47	/data/pavpan/nightlies/herbie/eggmath-testversion/src/mainloop.rkt	345	0
(unnamed)	/opt/racket-7.2/collects/racket/private/more-scheme.rkt	261	28
run	/opt/racket-7.2/share/pkgs/profile-lib/main.rkt	39	2
profile-thunk16	/opt/racket-7.2/share/pkgs/profile-lib/main.rkt	9	0
(unnamed)	/opt/racket-7.2/collects/racket/private/more-scheme.rkt	261	28