#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 r22324 = lambda1; float r22325 = theta; float r22326 = sin(r22325); float r22327 = delta; float r22328 = sin(r22327); float r22329 = r22326 * r22328; float r22330 = phi1; float r22331 = cos(r22330); float r22332 = r22329 * r22331; float r22333 = cos(r22327); float r22334 = sin(r22330); float r22335 = r22334 * r22333; float r22336 = r22331 * r22328; float r22337 = cos(r22325); float r22338 = r22336 * r22337; float r22339 = r22335 + r22338; float r22340 = asin(r22339); float r22341 = sin(r22340); float r22342 = r22334 * r22341; float r22343 = r22333 - r22342; float r22344 = atan2(r22332, r22343); float r22345 = r22324 + r22344; return r22345; } double f_id(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { double r22346 = lambda1; double r22347 = theta; double r22348 = sin(r22347); double r22349 = delta; double r22350 = sin(r22349); double r22351 = r22348 * r22350; double r22352 = phi1; double r22353 = cos(r22352); double r22354 = r22351 * r22353; double r22355 = cos(r22349); double r22356 = sin(r22352); double r22357 = r22356 * r22355; double r22358 = r22353 * r22350; double r22359 = cos(r22347); double r22360 = r22358 * r22359; double r22361 = r22357 + r22360; double r22362 = asin(r22361); double r22363 = sin(r22362); double r22364 = r22356 * r22363; double r22365 = r22355 - r22364; double r22366 = atan2(r22354, r22365); double r22367 = r22346 + r22366; return r22367; } double f_of(float lambda1, float phi1, float __attribute__((unused)) phi2, float delta, float theta) { float r22368 = lambda1; float r22369 = theta; float r22370 = sin(r22369); float r22371 = delta; float r22372 = sin(r22371); float r22373 = r22370 * r22372; float r22374 = phi1; float r22375 = cos(r22374); float r22376 = r22373 * r22375; float r22377 = cos(r22371); float r22378 = sin(r22374); float r22379 = atan2(1.0, 0.0); float r22380 = 2; float r22381 = r22379 / r22380; float r22382 = cos(r22369); float r22383 = r22372 * r22382; float r22384 = r22377 * r22378; float r22385 = fma(r22375, r22383, r22384); float r22386 = acos(r22385); float r22387 = r22381 - r22386; float r22388 = sin(r22387); float r22389 = r22378 * r22388; float r22390 = r22377 - r22389; float r22391 = expm1(r22390); float r22392 = log1p(r22391); float r22393 = atan2(r22376, r22392); float r22394 = r22368 + r22393; return r22394; } double f_od(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { double r22395 = lambda1; double r22396 = theta; double r22397 = sin(r22396); double r22398 = delta; double r22399 = sin(r22398); double r22400 = r22397 * r22399; double r22401 = phi1; double r22402 = cos(r22401); double r22403 = r22400 * r22402; double r22404 = cos(r22398); double r22405 = sin(r22401); double r22406 = atan2(1.0, 0.0); double r22407 = 2; double r22408 = r22406 / r22407; double r22409 = cos(r22396); double r22410 = r22399 * r22409; double r22411 = r22404 * r22405; double r22412 = fma(r22402, r22410, r22411); double r22413 = acos(r22412); double r22414 = r22408 - r22413; double r22415 = sin(r22414); double r22416 = r22405 * r22415; double r22417 = r22404 - r22416; double r22418 = expm1(r22417); double r22419 = log1p(r22418); double r22420 = atan2(r22403, r22419); double r22421 = r22395 + r22420; return r22421; } 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 r22422, r22423, r22424, r22425, r22426, r22427, r22428, r22429, r22430, r22431, r22432, r22433, r22434, r22435, r22436, r22437, r22438, r22439, r22440, r22441, r22442, r22443; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22422); mpfr_init(r22423); mpfr_init(r22424); mpfr_init(r22425); mpfr_init(r22426); mpfr_init(r22427); mpfr_init(r22428); mpfr_init(r22429); mpfr_init(r22430); mpfr_init(r22431); mpfr_init(r22432); mpfr_init(r22433); mpfr_init(r22434); mpfr_init(r22435); mpfr_init(r22436); mpfr_init(r22437); mpfr_init(r22438); mpfr_init(r22439); mpfr_init(r22440); mpfr_init(r22441); mpfr_init(r22442); mpfr_init(r22443); } double f_im(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { mpfr_set_d(r22422, lambda1, MPFR_RNDN); mpfr_set_d(r22423, theta, MPFR_RNDN); mpfr_sin(r22424, r22423, MPFR_RNDN); mpfr_set_d(r22425, delta, MPFR_RNDN); mpfr_sin(r22426, r22425, MPFR_RNDN); mpfr_mul(r22427, r22424, r22426, MPFR_RNDN); mpfr_set_d(r22428, phi1, MPFR_RNDN); mpfr_cos(r22429, r22428, MPFR_RNDN); mpfr_mul(r22430, r22427, r22429, MPFR_RNDN); mpfr_cos(r22431, r22425, MPFR_RNDN); mpfr_sin(r22432, r22428, MPFR_RNDN); mpfr_mul(r22433, r22432, r22431, MPFR_RNDN); mpfr_mul(r22434, r22429, r22426, MPFR_RNDN); mpfr_cos(r22435, r22423, MPFR_RNDN); mpfr_mul(r22436, r22434, r22435, MPFR_RNDN); mpfr_add(r22437, r22433, r22436, MPFR_RNDN); mpfr_asin(r22438, r22437, MPFR_RNDN); mpfr_sin(r22439, r22438, MPFR_RNDN); mpfr_mul(r22440, r22432, r22439, MPFR_RNDN); mpfr_sub(r22441, r22431, r22440, MPFR_RNDN); mpfr_atan2(r22442, r22430, r22441, MPFR_RNDN); mpfr_add(r22443, r22422, r22442, MPFR_RNDN); return mpfr_get_d(r22443, MPFR_RNDN); } static mpfr_t r22444, r22445, r22446, r22447, r22448, r22449, r22450, r22451, r22452, r22453, r22454, r22455, r22456, r22457, r22458, r22459, r22460, r22461, r22462, r22463, r22464, r22465, r22466, r22467, r22468, r22469, r22470; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22444); mpfr_init(r22445); mpfr_init(r22446); mpfr_init(r22447); mpfr_init(r22448); mpfr_init(r22449); mpfr_init(r22450); mpfr_init(r22451); mpfr_init(r22452); mpfr_init(r22453); mpfr_init(r22454); mpfr_init(r22455); mpfr_init_set_str(r22456, "2", 10, MPFR_RNDN); mpfr_init(r22457); mpfr_init(r22458); mpfr_init(r22459); mpfr_init(r22460); mpfr_init(r22461); mpfr_init(r22462); mpfr_init(r22463); mpfr_init(r22464); mpfr_init(r22465); mpfr_init(r22466); mpfr_init(r22467); mpfr_init(r22468); mpfr_init(r22469); mpfr_init(r22470); } double f_fm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { mpfr_set_d(r22444, lambda1, MPFR_RNDN); mpfr_set_d(r22445, theta, MPFR_RNDN); mpfr_sin(r22446, r22445, MPFR_RNDN); mpfr_set_d(r22447, delta, MPFR_RNDN); mpfr_sin(r22448, r22447, MPFR_RNDN); mpfr_mul(r22449, r22446, r22448, MPFR_RNDN); mpfr_set_d(r22450, phi1, MPFR_RNDN); mpfr_cos(r22451, r22450, MPFR_RNDN); mpfr_mul(r22452, r22449, r22451, MPFR_RNDN); mpfr_cos(r22453, r22447, MPFR_RNDN); mpfr_sin(r22454, r22450, MPFR_RNDN); mpfr_const_pi(r22455, MPFR_RNDN); ; mpfr_div(r22457, r22455, r22456, MPFR_RNDN); mpfr_cos(r22458, r22445, MPFR_RNDN); mpfr_mul(r22459, r22448, r22458, MPFR_RNDN); mpfr_mul(r22460, r22453, r22454, MPFR_RNDN); mpfr_fma(r22461, r22451, r22459, r22460, MPFR_RNDN); mpfr_acos(r22462, r22461, MPFR_RNDN); mpfr_sub(r22463, r22457, r22462, MPFR_RNDN); mpfr_sin(r22464, r22463, MPFR_RNDN); mpfr_mul(r22465, r22454, r22464, MPFR_RNDN); mpfr_sub(r22466, r22453, r22465, MPFR_RNDN); mpfr_expm1(r22467, r22466, MPFR_RNDN); mpfr_log1p(r22468, r22467, MPFR_RNDN); mpfr_atan2(r22469, r22452, r22468, MPFR_RNDN); mpfr_add(r22470, r22444, r22469, MPFR_RNDN); return mpfr_get_d(r22470, MPFR_RNDN); } static mpfr_t r22471, r22472, r22473, r22474, r22475, r22476, r22477, r22478, r22479, r22480, r22481, r22482, r22483, r22484, r22485, r22486, r22487, r22488, r22489, r22490, r22491, r22492, r22493, r22494, r22495, r22496, r22497; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22471); mpfr_init(r22472); mpfr_init(r22473); mpfr_init(r22474); mpfr_init(r22475); mpfr_init(r22476); mpfr_init(r22477); mpfr_init(r22478); mpfr_init(r22479); mpfr_init(r22480); mpfr_init(r22481); mpfr_init(r22482); mpfr_init_set_str(r22483, "2", 10, MPFR_RNDN); 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); } double f_dm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { mpfr_set_d(r22471, lambda1, MPFR_RNDN); mpfr_set_d(r22472, theta, MPFR_RNDN); mpfr_sin(r22473, r22472, MPFR_RNDN); mpfr_set_d(r22474, delta, MPFR_RNDN); mpfr_sin(r22475, r22474, MPFR_RNDN); mpfr_mul(r22476, r22473, r22475, MPFR_RNDN); mpfr_set_d(r22477, phi1, MPFR_RNDN); mpfr_cos(r22478, r22477, MPFR_RNDN); mpfr_mul(r22479, r22476, r22478, MPFR_RNDN); mpfr_cos(r22480, r22474, MPFR_RNDN); mpfr_sin(r22481, r22477, MPFR_RNDN); mpfr_const_pi(r22482, MPFR_RNDN); ; mpfr_div(r22484, r22482, r22483, MPFR_RNDN); mpfr_cos(r22485, r22472, MPFR_RNDN); mpfr_mul(r22486, r22475, r22485, MPFR_RNDN); mpfr_mul(r22487, r22480, r22481, MPFR_RNDN); mpfr_fma(r22488, r22478, r22486, r22487, MPFR_RNDN); mpfr_acos(r22489, r22488, MPFR_RNDN); mpfr_sub(r22490, r22484, r22489, MPFR_RNDN); mpfr_sin(r22491, r22490, MPFR_RNDN); mpfr_mul(r22492, r22481, r22491, MPFR_RNDN); mpfr_sub(r22493, r22480, r22492, MPFR_RNDN); mpfr_expm1(r22494, r22493, MPFR_RNDN); mpfr_log1p(r22495, r22494, MPFR_RNDN); mpfr_atan2(r22496, r22479, r22495, MPFR_RNDN); mpfr_add(r22497, r22471, r22496, MPFR_RNDN); return mpfr_get_d(r22497, MPFR_RNDN); }