\[1.00000000000000001 \cdot 10^{-150} \lt \left|x\right| \lt 9.99999999999999981 \cdot 10^{149}\]

\[\sqrt{0.5 \cdot \left(1 + \frac{x}{\sqrt{\left(4 \cdot p\right) \cdot p + x \cdot x}}\right)}\]

\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))))))))))))));
}

Reproduce

Please include this information when filing a bug report:

herbie shell --seed 2020165 
(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))))))))

Backtrace

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

Reproduce

herbie shell --seed 2020165 
(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))))))))