Exception for normal distribution

\[0.0 \le u1 \le 1 \land 0.0 \le u2 \le 1\]

Math

\[\left(\frac{1}{6} \cdot {\left(-2 \cdot \log u1\right)}^{0.5}\right) \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right) + 0.5\]

C

double code(double u1, double u2) {
	return ((double) (((double) (((double) (((double) (1.0 / 6.0)) * ((double) pow(((double) (-2.0 * ((double) log(u1)))), 0.5)))) * ((double) cos(((double) (((double) (2.0 * ((double) M_PI))) * u2)))))) + 0.5));
}

Reproduce

Please include this information when filing a bug report:

herbie shell --seed 2020165 
(FPCore (u1 u2)
  :name "normal distribution"
  :precision binary64
  :pre (and (<= 0.0 u1 1.0) (<= 0.0 u2 1.0))
  (+ (* (* (/ 1.0 6.0) (pow (* -2.0 (log u1)) 0.5)) (cos (* (* 2.0 PI) u2))) 0.5))

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 (u1 u2)
  :name "normal distribution"
  :precision binary64
  :pre (and (<= 0.0 u1 1.0) (<= 0.0 u2 1.0))
  (+ (* (* (/ 1.0 6.0) (pow (* -2.0 (log u1)) 0.5)) (cos (* (* 2.0 PI) u2))) 0.5))