\frac{\left(-b\right) + \sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c}}{2 \cdot a}double f(double a, double b, double c) {
double r26885 = b;
double r26886 = -r26885;
double r26887 = r26885 * r26885;
double r26888 = 4.0;
double r26889 = a;
double r26890 = r26888 * r26889;
double r26891 = c;
double r26892 = r26890 * r26891;
double r26893 = r26887 - r26892;
double r26894 = sqrt(r26893);
double r26895 = r26886 + r26894;
double r26896 = 2.0;
double r26897 = r26896 * r26889;
double r26898 = r26895 / r26897;
return r26898;
}
Please include this information when filing a bug report:
herbie shell --seed 2019323 +o rules:numerics
(FPCore (a b c)
:name "Quadratic roots, narrow range"
:precision binary64
:pre (and (< 1.0536712127723509e-08 a 94906265.62425156) (< 1.0536712127723509e-08 b 94906265.62425156) (< 1.0536712127723509e-08 c 94906265.62425156))
(/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))
| dict-ref: no value for key: 'branch26884 in: '((a . binary64) (b . binary64) (c . binary64)) | L | C | |
|---|---|---|---|
| for-loop | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/points.rkt | 106 | 4 |
| make-sampler | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/points.rkt | 99 | 0 |
| prepare-points-intervals | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/points.rkt | 121 | 0 |
| for-loop | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/core/regimes.rkt | 233 | 9 |
| sindices->spoints | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/core/regimes.rkt | 194 | 0 |
| f33 | /opt/racket-7.2/collects/racket/match/compiler.rkt | 507 | 40 |
| get-final-combination | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/mainloop.rkt | 369 | 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 |