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))))))))))));
}
double code(double lambda1, double lambda2, double phi1, double phi2) {
return ((double) (lambda1 + ((double) atan2(((double) (((double) cos(phi2)) * ((double) (((double) (((double) sin(lambda1)) * ((double) cos(lambda2)))) - ((double) (((double) cos(lambda1)) * ((double) sin(lambda2)))))))), ((double) (((double) cos(phi1)) + ((double) log(((double) pow(((double) exp(((double) cos(phi2)))), ((double) (((double) (((double) cos(lambda2)) * ((double) cos(lambda1)))) + ((double) (((double) sin(lambda1)) * ((double) sin(lambda2))))))))))))))));
}



Bits error versus lambda1



Bits error versus lambda2



Bits error versus phi1



Bits error versus phi2
Results
Initial program 0.7
rmApplied sin-diff0.7
rmApplied cos-diff0.2
rmApplied add-log-exp0.2
Simplified0.3
Final simplification0.3
herbie shell --seed 2020191
(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)))))))