\lambda_1 + \tan^{-1}_* \frac{\left(\sin theta \cdot \sin delta\right) \cdot \cos \phi_1}{\cos delta - \sin \phi_1 \cdot \sin \sin^{-1} \left(\sin \phi_1 \cdot \cos delta + \left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta\right)}\lambda_1 + \tan^{-1}_* \frac{\left(\sin theta \cdot \sin delta\right) \cdot \cos \phi_1}{\frac{{\cos delta}^{3} - \log \left(e^{{\left(\sin \phi_1 \cdot \sin \sin^{-1} \left(\cos delta \cdot \sin \phi_1 + \left(\sin delta \cdot \cos \phi_1\right) \cdot \cos theta\right)\right)}^{3}}\right)}{\cos delta \cdot \cos delta + \sin \phi_1 \cdot \left(\sin \sin^{-1} \left(\cos delta \cdot \sin \phi_1 + \left(\sin delta \cdot \cos \phi_1\right) \cdot \cos theta\right) \cdot \left(\cos delta + \sin \phi_1 \cdot \sin \sin^{-1} \left(\cos delta \cdot \sin \phi_1 + \left(\sin delta \cdot \cos \phi_1\right) \cdot \cos theta\right)\right)\right)}}(FPCore (lambda1 phi1 phi2 delta theta)
:precision binary64
(+
lambda1
(atan2
(* (* (sin theta) (sin delta)) (cos phi1))
(-
(cos delta)
(*
(sin phi1)
(sin
(asin
(+
(* (sin phi1) (cos delta))
(* (* (cos phi1) (sin delta)) (cos theta))))))))))(FPCore (lambda1 phi1 phi2 delta theta)
:precision binary64
(+
lambda1
(atan2
(* (* (sin theta) (sin delta)) (cos phi1))
(/
(-
(pow (cos delta) 3.0)
(log
(exp
(pow
(*
(sin phi1)
(sin
(asin
(+
(* (cos delta) (sin phi1))
(* (* (sin delta) (cos phi1)) (cos theta))))))
3.0))))
(+
(* (cos delta) (cos delta))
(*
(sin phi1)
(*
(sin
(asin
(+
(* (cos delta) (sin phi1))
(* (* (sin delta) (cos phi1)) (cos theta)))))
(+
(cos delta)
(*
(sin phi1)
(sin
(asin
(+
(* (cos delta) (sin phi1))
(* (* (sin delta) (cos phi1)) (cos theta))))))))))))))double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1 + atan2(((sin(theta) * sin(delta)) * cos(phi1)), (cos(delta) - (sin(phi1) * sin(asin((sin(phi1) * cos(delta)) + ((cos(phi1) * sin(delta)) * cos(theta)))))));
}
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1 + atan2(((sin(theta) * sin(delta)) * cos(phi1)), ((pow(cos(delta), 3.0) - log(exp(pow((sin(phi1) * sin(asin((cos(delta) * sin(phi1)) + ((sin(delta) * cos(phi1)) * cos(theta))))), 3.0)))) / ((cos(delta) * cos(delta)) + (sin(phi1) * (sin(asin((cos(delta) * sin(phi1)) + ((sin(delta) * cos(phi1)) * cos(theta)))) * (cos(delta) + (sin(phi1) * sin(asin((cos(delta) * sin(phi1)) + ((sin(delta) * cos(phi1)) * cos(theta)))))))))));
}



Bits error versus lambda1



Bits error versus phi1



Bits error versus phi2



Bits error versus delta



Bits error versus theta
Results
Initial program 0.1
rmApplied flip3--_binary640.2
Simplified0.2
Simplified0.2
rmApplied add-log-exp_binary640.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2020232
(FPCore (lambda1 phi1 phi2 delta theta)
:name "Destination given bearing on a great circle"
:precision binary64
(+ lambda1 (atan2 (* (* (sin theta) (sin delta)) (cos phi1)) (- (cos delta) (* (sin phi1) (sin (asin (+ (* (sin phi1) (cos delta)) (* (* (cos phi1) (sin delta)) (cos theta))))))))))