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