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) (((double) (((double) M_PI) / 2.0)) - ((double) asin(((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
rmApplied cos-diff3.9
rmApplied add-exp-log3.9
Simplified3.9
rmApplied acos-asin4.0
Final simplification4.0
herbie shell --seed 2020190
(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))