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) (((double) cos(phi1)) + ((double) (((double) cos(phi2)) * ((double) (((double) cos(lambda2)) * ((double) cos(lambda1)))))))) + ((double) (((double) cos(phi2)) * ((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.9
rmApplied sin-diff0.8
rmApplied cos-diff0.2
rmApplied distribute-lft-in0.2
Applied associate-+r+0.2
Final simplification0.2
herbie shell --seed 2020181
(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)))))))