Exception for Midpoint on a great circle

Math

\[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\]

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	return ((double) (lambda1 + ((double) atan2(((double) (((double) cos(phi2)) * ((double) sin(((double) (lambda1 - lambda2)))))), ((double) (((double) cos(phi1)) + ((double) (((double) cos(phi2)) * ((double) cos(((double) (lambda1 - lambda2))))))))))));
}

Reproduce

Please include this information when filing a bug report:

herbie shell --seed 2020165 
(FPCore (lambda1 lambda2 phi1 phi2)
  :name "Midpoint on a great circle"
  :precision binary64
  (+ lambda1 (atan2 (* (cos phi2) (sin (- lambda1 lambda2))) (+ (cos phi1) (* (cos phi2) (cos (- lambda1 lambda2)))))))

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 (lambda1 lambda2 phi1 phi2)
  :name "Midpoint on a great circle"
  :precision binary64
  (+ lambda1 (atan2 (* (cos phi2) (sin (- lambda1 lambda2))) (+ (cos phi1) (* (cos phi2) (cos (- lambda1 lambda2)))))))