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 2020165
(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))))))
| ival-hi: contract violation expected: ival? given: #<procedure:interval> | L | C | |
|---|---|---|---|
| for-loop | /data/pavpan/nightlies/herbie/ival-find-range/src/sampling.rkt | 37 | 11 |
| hyperrect-weight | /data/pavpan/nightlies/herbie/ival-find-range/src/sampling.rkt | 36 | 0 |
| total-weight | /data/pavpan/nightlies/herbie/ival-find-range/src/sampling.rkt | 113 | 2 |
| log-space-improvement | /data/pavpan/nightlies/herbie/ival-find-range/src/sampling.rkt | 110 | 0 |
| (unnamed) | /opt/racket-7.5/collects/racket/contract/private/arrow-higher-order.rkt | 375 | 33 |
| get-hyperrects | /data/pavpan/nightlies/herbie/ival-find-range/src/sampling.rkt | 123 | 0 |
| make-sampler | /data/pavpan/nightlies/herbie/ival-find-range/src/sampling.rkt | 147 | 0 |
| prepare-points-intervals | /data/pavpan/nightlies/herbie/ival-find-range/src/points.rkt | 90 | 0 |
| run-improve47 | /data/pavpan/nightlies/herbie/ival-find-range/src/mainloop.rkt | 365 | 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 2020165
(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))))))