1 - \left(\frac{1}{1 + 0.3275911000000000239396058532292954623699 \cdot \left|x\right|} \cdot \left(0.2548295919999999936678136691625695675611 + \frac{1}{1 + 0.3275911000000000239396058532292954623699 \cdot \left|x\right|} \cdot \left(-0.2844967359999999723108032867457950487733 + \frac{1}{1 + 0.3275911000000000239396058532292954623699 \cdot \left|x\right|} \cdot \left(1.421413741000000063863240029604639858007 + \frac{1}{1 + 0.3275911000000000239396058532292954623699 \cdot \left|x\right|} \cdot \left(-1.453152027000000012790792425221297889948 + \frac{1}{1 + 0.3275911000000000239396058532292954623699 \cdot \left|x\right|} \cdot 1.061405428999999900341322245367337018251\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|}double f(double x) {
double r8833 = 1.0;
double r8834 = 0.3275911;
double r8835 = x;
double r8836 = fabs(r8835);
double r8837 = r8834 * r8836;
double r8838 = r8833 + r8837;
double r8839 = r8833 / r8838;
double r8840 = 0.254829592;
double r8841 = -0.284496736;
double r8842 = 1.421413741;
double r8843 = -1.453152027;
double r8844 = 1.061405429;
double r8845 = r8839 * r8844;
double r8846 = r8843 + r8845;
double r8847 = r8839 * r8846;
double r8848 = r8842 + r8847;
double r8849 = r8839 * r8848;
double r8850 = r8841 + r8849;
double r8851 = r8839 * r8850;
double r8852 = r8840 + r8851;
double r8853 = r8839 * r8852;
double r8854 = r8836 * r8836;
double r8855 = -r8854;
double r8856 = exp(r8855);
double r8857 = r8853 * r8856;
double r8858 = r8833 - r8857;
return r8858;
}
Please include this information when filing a bug report:
herbie shell --seed 2019173 +o rules:numerics
(FPCore (x)
:name "Jmat.Real.erf"
(- 1.0 (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))))
| atab-pick-alt: contract violation expected: alt? given: '(#<alt (λ (x) (- 1.0 (/ (/ (fma (fma (/ 1.0 (fma 0.3275911 (fabs x) 1.0)) (fma (/ 1.0 (fma 0.3275911 (fabs x) 1.0)) (fma (/ 1.0 (fma 0.3275911 (fabs x) 1.0)) 1.061405429 -1.453152027) 1.421413741) -0.284496736) (/ 1.0 (fma 0.3275911 (fabs x) 1.0)) 0.2... 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 |