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 2020128
(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))))))
| random-bits: contract violation expected: Positive-Integer given: 0 | L | C | |
|---|---|---|---|
| random-ranges | /data/pavpan/nightlies/herbie/separate-biginterval/src/common.rkt | 139 | 0 |
| sample-multi-bounded | /data/pavpan/nightlies/herbie/separate-biginterval/src/points.rkt | 14 | 0 |
| (unnamed) | /data/pavpan/nightlies/herbie/separate-biginterval/src/points.rkt | 127 | 6 |
| prepare-points-intervals | /data/pavpan/nightlies/herbie/separate-biginterval/src/points.rkt | 134 | 0 |
| pred | /data/pavpan/nightlies/herbie/separate-biginterval/src/core/regimes.rkt | 200 | 4 |
| for-loop | /data/pavpan/nightlies/herbie/separate-biginterval/src/core/regimes.rkt | 230 | 9 |
| sindices->spoints | /data/pavpan/nightlies/herbie/separate-biginterval/src/core/regimes.rkt | 190 | 0 |
| f37 | unknown | ||
| get-final-combination | /data/pavpan/nightlies/herbie/separate-biginterval/src/mainloop.rkt | 390 | 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 2020128
(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))))))