wj - \frac{wj \cdot e^{wj} - x}{e^{wj} + wj \cdot e^{wj}}double code(double wj, double x) {
return ((double) (wj - ((double) (((double) (((double) (wj * ((double) exp(wj)))) - x)) / ((double) (((double) exp(wj)) + ((double) (wj * ((double) exp(wj))))))))));
}
Please include this information when filing a bug report:
herbie shell --seed 2020169
(FPCore (wj x)
:name "Jmat.Real.lambertw, newton loop step"
:precision binary64
:herbie-target
(- wj (- (/ wj (+ wj 1.0)) (/ x (+ (exp wj) (* wj (exp wj))))))
(- wj (/ (- (* wj (exp wj)) x) (+ (exp wj) (* wj (exp wj))))))
| match-define: no matching clause for '(== branch457 -1.7792498876631063e-12) | L | C | |
|---|---|---|---|
| match:error | /opt/racket-7.5/collects/racket/match/runtime.rkt | 24 | 0 |
| program-body | /data/pavpan/nightlies/herbie/ival-find-range/src/programs.rkt | 25 | 0 |
| (unnamed) | /opt/racket-7.5/collects/racket/contract/private/arrow-higher-order.rkt | 379 | 33 |
| prepare-points | /data/pavpan/nightlies/herbie/ival-find-range/src/points.rkt | 130 | 0 |
| pred | /data/pavpan/nightlies/herbie/ival-find-range/src/core/regimes.rkt | 204 | 4 |
| .../more-scheme.rkt:261:28 | unknown | ||
| for-loop | /data/pavpan/nightlies/herbie/ival-find-range/src/core/regimes.rkt | 240 | 3 |
| sindices->spoints | /data/pavpan/nightlies/herbie/ival-find-range/src/core/regimes.rkt | 194 | 0 |
| f37 | unknown | ||
| get-final-combination | /data/pavpan/nightlies/herbie/ival-find-range/src/mainloop.rkt | 389 | 0 |
| .../more-scheme.rkt:261:28 | unknown | ||
| run | /opt/racket-7.5/share/pkgs/profile-lib/main.rkt | 39 | 2 |
| profile-thunk16 | /opt/racket-7.5/share/pkgs/profile-lib/main.rkt | 9 | 0 |
| (unnamed) | /opt/racket-7.5/collects/racket/private/more-scheme.rkt | 261 | 28 |
herbie shell --seed 2020169
(FPCore (wj x)
:name "Jmat.Real.lambertw, newton loop step"
:precision binary64
:herbie-target
(- wj (- (/ wj (+ wj 1.0)) (/ x (+ (exp wj) (* wj (exp wj))))))
(- wj (/ (- (* wj (exp wj)) x) (+ (exp wj) (* wj (exp wj))))))