#include #include #include #include #include char *name = "Destination given bearing on a great circle"; double f_if(float lambda1, float phi1, float __attribute__((unused)) phi2, float delta, float theta) { float r22396 = lambda1; float r22397 = theta; float r22398 = sin(r22397); float r22399 = delta; float r22400 = sin(r22399); float r22401 = r22398 * r22400; float r22402 = phi1; float r22403 = cos(r22402); float r22404 = r22401 * r22403; float r22405 = cos(r22399); float r22406 = sin(r22402); float r22407 = r22406 * r22405; float r22408 = r22403 * r22400; float r22409 = cos(r22397); float r22410 = r22408 * r22409; float r22411 = r22407 + r22410; float r22412 = asin(r22411); float r22413 = sin(r22412); float r22414 = r22406 * r22413; float r22415 = r22405 - r22414; float r22416 = atan2(r22404, r22415); float r22417 = r22396 + r22416; return r22417; } double f_id(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { double r22418 = lambda1; double r22419 = theta; double r22420 = sin(r22419); double r22421 = delta; double r22422 = sin(r22421); double r22423 = r22420 * r22422; double r22424 = phi1; double r22425 = cos(r22424); double r22426 = r22423 * r22425; double r22427 = cos(r22421); double r22428 = sin(r22424); double r22429 = r22428 * r22427; double r22430 = r22425 * r22422; double r22431 = cos(r22419); double r22432 = r22430 * r22431; double r22433 = r22429 + r22432; double r22434 = asin(r22433); double r22435 = sin(r22434); double r22436 = r22428 * r22435; double r22437 = r22427 - r22436; double r22438 = atan2(r22426, r22437); double r22439 = r22418 + r22438; return r22439; } double f_of(float lambda1, float phi1, float __attribute__((unused)) phi2, float delta, float theta) { float r22440 = delta; float r22441 = sin(r22440); float r22442 = phi1; float r22443 = cos(r22442); float r22444 = theta; float r22445 = sin(r22444); float r22446 = r22443 * r22445; float r22447 = r22441 * r22446; float r22448 = cos(r22440); float r22449 = sin(r22442); float r22450 = cos(r22444); float r22451 = r22443 * r22450; float r22452 = r22449 * r22448; float r22453 = fma(r22451, r22441, r22452); float r22454 = r22449 * r22453; float r22455 = r22448 - r22454; float r22456 = atan2(r22447, r22455); float r22457 = lambda1; float r22458 = r22456 + r22457; return r22458; } double f_od(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { double r22459 = delta; double r22460 = sin(r22459); double r22461 = phi1; double r22462 = cos(r22461); double r22463 = theta; double r22464 = sin(r22463); double r22465 = r22462 * r22464; double r22466 = r22460 * r22465; double r22467 = cos(r22459); double r22468 = sin(r22461); double r22469 = cos(r22463); double r22470 = r22462 * r22469; double r22471 = r22468 * r22467; double r22472 = fma(r22470, r22460, r22471); double r22473 = r22468 * r22472; double r22474 = r22467 - r22473; double r22475 = atan2(r22466, r22474); double r22476 = lambda1; double r22477 = r22475 + r22476; return r22477; } void mpfr_fmod2(mpfr_t r, mpfr_t n, mpfr_t d, mpfr_rnd_t rmd) { mpfr_fmod(r, n, d, rmd); if (mpfr_cmp_ui(r, 0) < 0) mpfr_add(r, r, d, rmd); } static mpfr_t r22478, r22479, r22480, r22481, r22482, r22483, r22484, r22485, r22486, r22487, r22488, r22489, r22490, r22491, r22492, r22493, r22494, r22495, r22496, r22497, r22498, r22499; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22478); mpfr_init(r22479); mpfr_init(r22480); mpfr_init(r22481); mpfr_init(r22482); mpfr_init(r22483); mpfr_init(r22484); mpfr_init(r22485); mpfr_init(r22486); mpfr_init(r22487); mpfr_init(r22488); mpfr_init(r22489); mpfr_init(r22490); mpfr_init(r22491); mpfr_init(r22492); mpfr_init(r22493); mpfr_init(r22494); mpfr_init(r22495); mpfr_init(r22496); mpfr_init(r22497); mpfr_init(r22498); mpfr_init(r22499); } double f_im(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { mpfr_set_d(r22478, lambda1, MPFR_RNDN); mpfr_set_d(r22479, theta, MPFR_RNDN); mpfr_sin(r22480, r22479, MPFR_RNDN); mpfr_set_d(r22481, delta, MPFR_RNDN); mpfr_sin(r22482, r22481, MPFR_RNDN); mpfr_mul(r22483, r22480, r22482, MPFR_RNDN); mpfr_set_d(r22484, phi1, MPFR_RNDN); mpfr_cos(r22485, r22484, MPFR_RNDN); mpfr_mul(r22486, r22483, r22485, MPFR_RNDN); mpfr_cos(r22487, r22481, MPFR_RNDN); mpfr_sin(r22488, r22484, MPFR_RNDN); mpfr_mul(r22489, r22488, r22487, MPFR_RNDN); mpfr_mul(r22490, r22485, r22482, MPFR_RNDN); mpfr_cos(r22491, r22479, MPFR_RNDN); mpfr_mul(r22492, r22490, r22491, MPFR_RNDN); mpfr_add(r22493, r22489, r22492, MPFR_RNDN); mpfr_asin(r22494, r22493, MPFR_RNDN); mpfr_sin(r22495, r22494, MPFR_RNDN); mpfr_mul(r22496, r22488, r22495, MPFR_RNDN); mpfr_sub(r22497, r22487, r22496, MPFR_RNDN); mpfr_atan2(r22498, r22486, r22497, MPFR_RNDN); mpfr_add(r22499, r22478, r22498, MPFR_RNDN); return mpfr_get_d(r22499, MPFR_RNDN); } static mpfr_t r22500, r22501, r22502, r22503, r22504, r22505, r22506, r22507, r22508, r22509, r22510, r22511, r22512, r22513, r22514, r22515, r22516, r22517, r22518; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22500); mpfr_init(r22501); mpfr_init(r22502); mpfr_init(r22503); mpfr_init(r22504); mpfr_init(r22505); mpfr_init(r22506); mpfr_init(r22507); mpfr_init(r22508); mpfr_init(r22509); mpfr_init(r22510); mpfr_init(r22511); mpfr_init(r22512); mpfr_init(r22513); mpfr_init(r22514); mpfr_init(r22515); mpfr_init(r22516); mpfr_init(r22517); mpfr_init(r22518); } double f_fm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { mpfr_set_d(r22500, delta, MPFR_RNDN); mpfr_sin(r22501, r22500, MPFR_RNDN); mpfr_set_d(r22502, phi1, MPFR_RNDN); mpfr_cos(r22503, r22502, MPFR_RNDN); mpfr_set_d(r22504, theta, MPFR_RNDN); mpfr_sin(r22505, r22504, MPFR_RNDN); mpfr_mul(r22506, r22503, r22505, MPFR_RNDN); mpfr_mul(r22507, r22501, r22506, MPFR_RNDN); mpfr_cos(r22508, r22500, MPFR_RNDN); mpfr_sin(r22509, r22502, MPFR_RNDN); mpfr_cos(r22510, r22504, MPFR_RNDN); mpfr_mul(r22511, r22503, r22510, MPFR_RNDN); mpfr_mul(r22512, r22509, r22508, MPFR_RNDN); mpfr_fma(r22513, r22511, r22501, r22512, MPFR_RNDN); mpfr_mul(r22514, r22509, r22513, MPFR_RNDN); mpfr_sub(r22515, r22508, r22514, MPFR_RNDN); mpfr_atan2(r22516, r22507, r22515, MPFR_RNDN); mpfr_set_d(r22517, lambda1, MPFR_RNDN); mpfr_add(r22518, r22516, r22517, MPFR_RNDN); return mpfr_get_d(r22518, MPFR_RNDN); } static mpfr_t r22519, r22520, r22521, r22522, r22523, r22524, r22525, r22526, r22527, r22528, r22529, r22530, r22531, r22532, r22533, r22534, r22535, r22536, r22537; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22519); mpfr_init(r22520); mpfr_init(r22521); mpfr_init(r22522); mpfr_init(r22523); mpfr_init(r22524); mpfr_init(r22525); mpfr_init(r22526); mpfr_init(r22527); mpfr_init(r22528); mpfr_init(r22529); mpfr_init(r22530); mpfr_init(r22531); mpfr_init(r22532); mpfr_init(r22533); mpfr_init(r22534); mpfr_init(r22535); mpfr_init(r22536); mpfr_init(r22537); } double f_dm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { mpfr_set_d(r22519, delta, MPFR_RNDN); mpfr_sin(r22520, r22519, MPFR_RNDN); mpfr_set_d(r22521, phi1, MPFR_RNDN); mpfr_cos(r22522, r22521, MPFR_RNDN); mpfr_set_d(r22523, theta, MPFR_RNDN); mpfr_sin(r22524, r22523, MPFR_RNDN); mpfr_mul(r22525, r22522, r22524, MPFR_RNDN); mpfr_mul(r22526, r22520, r22525, MPFR_RNDN); mpfr_cos(r22527, r22519, MPFR_RNDN); mpfr_sin(r22528, r22521, MPFR_RNDN); mpfr_cos(r22529, r22523, MPFR_RNDN); mpfr_mul(r22530, r22522, r22529, MPFR_RNDN); mpfr_mul(r22531, r22528, r22527, MPFR_RNDN); mpfr_fma(r22532, r22530, r22520, r22531, MPFR_RNDN); mpfr_mul(r22533, r22528, r22532, MPFR_RNDN); mpfr_sub(r22534, r22527, r22533, MPFR_RNDN); mpfr_atan2(r22535, r22526, r22534, MPFR_RNDN); mpfr_set_d(r22536, lambda1, MPFR_RNDN); mpfr_add(r22537, r22535, r22536, MPFR_RNDN); return mpfr_get_d(r22537, MPFR_RNDN); }