R \cdot \left(2 \cdot \tan^{-1}_* \frac{\sqrt{{\left(\sin \left(\frac{\phi_1 - \phi_2}{2}\right)\right)}^{2} + \left(\left(\cos \phi_1 \cdot \cos \phi_2\right) \cdot \sin \left(\frac{\lambda_1 - \lambda_2}{2}\right)\right) \cdot \sin \left(\frac{\lambda_1 - \lambda_2}{2}\right)}}{\sqrt{1 - \left({\left(\sin \left(\frac{\phi_1 - \phi_2}{2}\right)\right)}^{2} + \left(\left(\cos \phi_1 \cdot \cos \phi_2\right) \cdot \sin \left(\frac{\lambda_1 - \lambda_2}{2}\right)\right) \cdot \sin \left(\frac{\lambda_1 - \lambda_2}{2}\right)\right)}}\right)double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return ((double) (R * ((double) (2.0 * ((double) atan2(((double) sqrt(((double) (((double) pow(((double) sin(((double) (((double) (phi1 - phi2)) / 2.0)))), 2.0)) + ((double) (((double) (((double) (((double) cos(phi1)) * ((double) cos(phi2)))) * ((double) sin(((double) (((double) (lambda1 - lambda2)) / 2.0)))))) * ((double) sin(((double) (((double) (lambda1 - lambda2)) / 2.0)))))))))), ((double) sqrt(((double) (1.0 - ((double) (((double) pow(((double) sin(((double) (((double) (phi1 - phi2)) / 2.0)))), 2.0)) + ((double) (((double) (((double) (((double) cos(phi1)) * ((double) cos(phi2)))) * ((double) sin(((double) (((double) (lambda1 - lambda2)) / 2.0)))))) * ((double) sin(((double) (((double) (lambda1 - lambda2)) / 2.0))))))))))))))))));
}
Please include this information when filing a bug report:
herbie shell --seed 2020165
(FPCore (R lambda1 lambda2 phi1 phi2)
:name "Distance on a great circle"
:precision binary64
(* R (* 2.0 (atan2 (sqrt (+ (pow (sin (/ (- phi1 phi2) 2.0)) 2.0) (* (* (* (cos phi1) (cos phi2)) (sin (/ (- lambda1 lambda2) 2.0))) (sin (/ (- lambda1 lambda2) 2.0))))) (sqrt (- 1.0 (+ (pow (sin (/ (- phi1 phi2) 2.0)) 2.0) (* (* (* (cos phi1) (cos phi2)) (sin (/ (- lambda1 lambda2) 2.0))) (sin (/ (- lambda1 lambda2) 2.0))))))))))
| 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 |
herbie shell --seed 2020165
(FPCore (R lambda1 lambda2 phi1 phi2)
:name "Distance on a great circle"
:precision binary64
(* R (* 2.0 (atan2 (sqrt (+ (pow (sin (/ (- phi1 phi2) 2.0)) 2.0) (* (* (* (cos phi1) (cos phi2)) (sin (/ (- lambda1 lambda2) 2.0))) (sin (/ (- lambda1 lambda2) 2.0))))) (sqrt (- 1.0 (+ (pow (sin (/ (- phi1 phi2) 2.0)) 2.0) (* (* (* (cos phi1) (cos phi2)) (sin (/ (- lambda1 lambda2) 2.0))) (sin (/ (- lambda1 lambda2) 2.0))))))))))