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



Bits error versus R



Bits error versus lambda1



Bits error versus lambda2



Bits error versus phi1



Bits error versus phi2
Results
Initial program 16.7
Simplified16.7
rmApplied cos-diff4.0
rmApplied add-exp-log4.0
rmApplied add-log-exp4.0
Final simplification4.0
herbie shell --seed 2020191
(FPCore (R lambda1 lambda2 phi1 phi2)
:name "Spherical law of cosines"
:precision binary64
(* (acos (+ (* (sin phi1) (sin phi2)) (* (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2))))) R))