\lambda_1 + \tan^{-1}_* \frac{\left(\sin theta \cdot \sin delta\right) \cdot \cos \phi_1}{\cos delta - \sin \phi_1 \cdot \sin \left(\sin^{-1} \left(\sin \phi_1 \cdot \cos delta + \left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta\right)\right)}\lambda_1 + \tan^{-1}_* \frac{\left(\sin theta \cdot \sin delta\right) \cdot \cos \phi_1}{\frac{{\left(\cos delta\right)}^{3} - \left(\left(\sin \phi_1 \cdot \sin \phi_1\right) \cdot \sin \phi_1\right) \cdot \left(\left(\sin \left(\sin^{-1} \left(\sin \phi_1 \cdot \cos delta + \left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta\right)\right) \cdot \sin \left(\sin^{-1} \left(\sin \phi_1 \cdot \cos delta + \left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta\right)\right)\right) \cdot \sin \left(\sin^{-1} \left(\sin \phi_1 \cdot \cos delta + \left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta\right)\right)\right)}{\left(\sin \phi_1 \cdot \sin \left(\sin^{-1} \left(\sin \phi_1 \cdot \cos delta + \left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta\right)\right)\right) \cdot \left(\sin \phi_1 \cdot \sin \left(\sin^{-1} \left(\sin \phi_1 \cdot \cos delta + \left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta\right)\right) + \cos delta\right) + \cos delta \cdot \cos delta}}double f(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
double r90272 = lambda1;
double r90273 = theta;
double r90274 = sin(r90273);
double r90275 = delta;
double r90276 = sin(r90275);
double r90277 = r90274 * r90276;
double r90278 = phi1;
double r90279 = cos(r90278);
double r90280 = r90277 * r90279;
double r90281 = cos(r90275);
double r90282 = sin(r90278);
double r90283 = r90282 * r90281;
double r90284 = r90279 * r90276;
double r90285 = cos(r90273);
double r90286 = r90284 * r90285;
double r90287 = r90283 + r90286;
double r90288 = asin(r90287);
double r90289 = sin(r90288);
double r90290 = r90282 * r90289;
double r90291 = r90281 - r90290;
double r90292 = atan2(r90280, r90291);
double r90293 = r90272 + r90292;
return r90293;
}
double f(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
double r90294 = lambda1;
double r90295 = theta;
double r90296 = sin(r90295);
double r90297 = delta;
double r90298 = sin(r90297);
double r90299 = r90296 * r90298;
double r90300 = phi1;
double r90301 = cos(r90300);
double r90302 = r90299 * r90301;
double r90303 = cos(r90297);
double r90304 = 3.0;
double r90305 = pow(r90303, r90304);
double r90306 = sin(r90300);
double r90307 = r90306 * r90306;
double r90308 = r90307 * r90306;
double r90309 = r90306 * r90303;
double r90310 = r90301 * r90298;
double r90311 = cos(r90295);
double r90312 = r90310 * r90311;
double r90313 = r90309 + r90312;
double r90314 = asin(r90313);
double r90315 = sin(r90314);
double r90316 = r90315 * r90315;
double r90317 = r90316 * r90315;
double r90318 = r90308 * r90317;
double r90319 = r90305 - r90318;
double r90320 = r90306 * r90315;
double r90321 = r90320 + r90303;
double r90322 = r90320 * r90321;
double r90323 = r90303 * r90303;
double r90324 = r90322 + r90323;
double r90325 = r90319 / r90324;
double r90326 = atan2(r90302, r90325);
double r90327 = r90294 + r90326;
return r90327;
}



Bits error versus lambda1



Bits error versus phi1



Bits error versus phi2



Bits error versus delta



Bits error versus theta
Results
Initial program 0.2
rmApplied flip3--0.2
Simplified0.2
rmApplied add-cbrt-cube0.2
Applied add-cbrt-cube0.2
Applied cbrt-unprod0.2
Applied rem-cube-cbrt0.2
Final simplification0.2
herbie shell --seed 2020002
(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))))))))))