\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)double f(double z) {
double r8586 = atan2(1.0, 0.0);
double r8587 = z;
double r8588 = r8586 * r8587;
double r8589 = sin(r8588);
double r8590 = r8586 / r8589;
double r8591 = 2.0;
double r8592 = r8586 * r8591;
double r8593 = sqrt(r8592);
double r8594 = 1.0;
double r8595 = r8594 - r8587;
double r8596 = r8595 - r8594;
double r8597 = 7.0;
double r8598 = r8596 + r8597;
double r8599 = 0.5;
double r8600 = r8598 + r8599;
double r8601 = r8596 + r8599;
double r8602 = pow(r8600, r8601);
double r8603 = r8593 * r8602;
double r8604 = -r8600;
double r8605 = exp(r8604);
double r8606 = r8603 * r8605;
double r8607 = 0.9999999999998099;
double r8608 = 676.5203681218851;
double r8609 = r8596 + r8594;
double r8610 = r8608 / r8609;
double r8611 = r8607 + r8610;
double r8612 = -1259.1392167224028;
double r8613 = r8596 + r8591;
double r8614 = r8612 / r8613;
double r8615 = r8611 + r8614;
double r8616 = 771.3234287776531;
double r8617 = 3.0;
double r8618 = r8596 + r8617;
double r8619 = r8616 / r8618;
double r8620 = r8615 + r8619;
double r8621 = -176.6150291621406;
double r8622 = 4.0;
double r8623 = r8596 + r8622;
double r8624 = r8621 / r8623;
double r8625 = r8620 + r8624;
double r8626 = 12.507343278686905;
double r8627 = 5.0;
double r8628 = r8596 + r8627;
double r8629 = r8626 / r8628;
double r8630 = r8625 + r8629;
double r8631 = -0.13857109526572012;
double r8632 = 6.0;
double r8633 = r8596 + r8632;
double r8634 = r8631 / r8633;
double r8635 = r8630 + r8634;
double r8636 = 9.984369578019572e-06;
double r8637 = r8636 / r8598;
double r8638 = r8635 + r8637;
double r8639 = 1.5056327351493116e-07;
double r8640 = 8.0;
double r8641 = r8596 + r8640;
double r8642 = r8639 / r8641;
double r8643 = r8638 + r8642;
double r8644 = r8606 * r8643;
double r8645 = r8590 * r8644;
return r8645;
}
Please include this information when filing a bug report:
herbie shell --seed 2019173 +o rules:numerics
(FPCore (z)
:name "Jmat.Real.gamma, branch z less than 0.5"
(* (/ PI (sin (* PI z))) (* (* (* (sqrt (* PI 2.0)) (pow (+ (+ (- (- 1.0 z) 1.0) 7.0) 0.5) (+ (- (- 1.0 z) 1.0) 0.5))) (exp (- (+ (+ (- (- 1.0 z) 1.0) 7.0) 0.5)))) (+ (+ (+ (+ (+ (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 (+ (- (- 1.0 z) 1.0) 1.0))) (/ -1259.1392167224028 (+ (- (- 1.0 z) 1.0) 2.0))) (/ 771.3234287776531 (+ (- (- 1.0 z) 1.0) 3.0))) (/ -176.6150291621406 (+ (- (- 1.0 z) 1.0) 4.0))) (/ 12.507343278686905 (+ (- (- 1.0 z) 1.0) 5.0))) (/ -0.13857109526572012 (+ (- (- 1.0 z) 1.0) 6.0))) (/ 9.984369578019572e-06 (+ (- (- 1.0 z) 1.0) 7.0))) (/ 1.5056327351493116e-07 (+ (- (- 1.0 z) 1.0) 8.0))))))
| atab-pick-alt: contract violation expected: alt? given: '(#<alt (λ (z) (* (/ PI (sin (* PI z))) (* (* (* (sqrt (* PI 2.0)) (pow (+ (+ (- (- 1.0 z) 1.0) 7.0) 0.5) (+ (- (- 1.0 z) 1.0) 0.5))) (exp (- (+ (+ (- (- 1.0 z) 1.0) 7.0) 0.5)))) (+ (+ (+ (+ (+ (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 (+ (- (- ... in: the range of the #:picking-func argument of (-> alt-table? #:only-fresh boolean? #:picking-func (-> (listof alt?) alt?) (values alt? alt-table?)) contract from: /data/pavpan/nightlies/herbie/no-costs-2/src/core/alt-table.rkt blaming: /data/pavpan/nightlies/herbie/no-costs-2/src/mainloop.rkt (assuming the contract is correct) at: /data/pavpan/nightlies/herbie/no-costs-2/src/core/alt-table.rkt:13.3 | L | C | |
|---|---|---|---|
| raise-blame-error16 | /opt/racket-7.0/collects/racket/contract/private/blame.rkt | 347 | 0 |
| (unnamed) | /data/pavpan/nightlies/herbie/no-costs-2/src/core/alt-table.rkt | 59 | 0 |
| (unnamed) | /opt/racket-7.0/collects/racket/contract/private/arrow-val-first.rkt | 388 | 18 |
| choose-best-alt! | /data/pavpan/nightlies/herbie/no-costs-2/src/mainloop.rkt | 110 | 0 |
| run-iter! | /data/pavpan/nightlies/herbie/no-costs-2/src/mainloop.rkt | 319 | 0 |
| for-loop | /data/pavpan/nightlies/herbie/no-costs-2/src/mainloop.rkt | 353 | 4 |
| run-improve43 | /data/pavpan/nightlies/herbie/no-costs-2/src/mainloop.rkt | 338 | 0 |
| (unnamed) | /opt/racket-7.0/collects/racket/private/more-scheme.rkt | 261 | 28 |
| run | /opt/racket-7.0/share/pkgs/profile-lib/main.rkt | 39 | 2 |
| profile-thunk16 | /opt/racket-7.0/share/pkgs/profile-lib/main.rkt | 9 | 0 |
| (unnamed) | /opt/racket-7.0/collects/racket/private/more-scheme.rkt | 261 | 28 |