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



Bits error versus lambda1



Bits error versus lambda2



Bits error versus phi1



Bits error versus phi2
Results
Initial program 13.1
rmApplied sin-diff6.8
Simplified6.8
rmApplied cos-diff0.2
Applied distribute-lft-in0.2
Simplified0.2
Simplified0.2
rmApplied add-log-exp0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2020180
(FPCore (lambda1 lambda2 phi1 phi2)
:name "Bearing on a great circle"
:precision binary64
(atan2 (* (sin (- lambda1 lambda2)) (cos phi2)) (- (* (cos phi1) (sin phi2)) (* (* (sin phi1) (cos phi2)) (cos (- lambda1 lambda2))))))