\sqrt{0.5 \cdot \left(1 + \frac{x}{\sqrt{\left(4 \cdot p\right) \cdot p + x \cdot x}}\right)}double code(double p, double x) {
return ((double) sqrt(((double) (0.5 * ((double) (1.0 + ((double) (x / ((double) sqrt(((double) (((double) (((double) (4.0 * p)) * p)) + ((double) (x * x))))))))))))));
}
Please include this information when filing a bug report:
herbie shell --seed 2020155
(FPCore (p x)
:name "Given's Rotation SVD example"
:precision binary64
:pre (< 1e-150 (fabs x) 1e+150)
:herbie-target
(sqrt (+ 0.5 (/ (copysign 0.5 x) (hypot 1.0 (/ (* 2.0 p) x)))))
(sqrt (* 0.5 (+ 1.0 (/ x (sqrt (+ (* (* 4.0 p) p) (* x x))))))))
| hash-ref: no value found for key key: 'p | L | C | |
|---|---|---|---|
| range-table->hyperrects | /data/pavpan/nightlies/herbie/ival-find-range/src/findroot.rkt | 58 | 0 |
| make-sampler | /data/pavpan/nightlies/herbie/ival-find-range/src/points.rkt | 124 | 0 |
| prepare-points-intervals | /data/pavpan/nightlies/herbie/ival-find-range/src/points.rkt | 155 | 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 2020155
(FPCore (p x)
:name "Given's Rotation SVD example"
:precision binary64
:pre (< 1e-150 (fabs x) 1e+150)
:herbie-target
(sqrt (+ 0.5 (/ (copysign 0.5 x) (hypot 1.0 (/ (* 2.0 p) x)))))
(sqrt (* 0.5 (+ 1.0 (/ x (sqrt (+ (* (* 4.0 p) p) (* x x))))))))