Exception for math.sqrt on complex, real part

Math

\[0.5 \cdot \sqrt{2 \cdot \left(\sqrt{re \cdot re + im \cdot im} + re\right)}\]

C

double code(double re, double im) {
	return ((double) (0.5 * ((double) sqrt(((double) (2.0 * ((double) (((double) sqrt(((double) (((double) (re * re)) + ((double) (im * im)))))) + re))))))));
}

Reproduce

Please include this information when filing a bug report:

herbie shell --seed 2020128 
(FPCore (re im)
  :name "math.sqrt on complex, real part"
  :precision binary64

  :herbie-target
  (if (< re 0.0) (* 0.5 (* (sqrt 2.0) (sqrt (/ (* im im) (- (sqrt (+ (* re re) (* im im))) re))))) (* 0.5 (sqrt (* 2.0 (+ (sqrt (+ (* re re) (* im im))) re)))))

  (* 0.5 (sqrt (* 2.0 (+ (sqrt (+ (* re re) (* im im))) re)))))

Backtrace

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

Reproduce

herbie shell --seed 2020128 
(FPCore (re im)
  :name "math.sqrt on complex, real part"
  :precision binary64

  :herbie-target
  (if (< re 0.0) (* 0.5 (* (sqrt 2.0) (sqrt (/ (* im im) (- (sqrt (+ (* re re) (* im im))) re))))) (* 0.5 (sqrt (* 2.0 (+ (sqrt (+ (* re re) (* im im))) re)))))

  (* 0.5 (sqrt (* 2.0 (+ (sqrt (+ (* re re) (* im im))) re)))))