#include #include #include #include #include char *name = "Bearing on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r26431 = lambda1; float r26432 = lambda2; float r26433 = r26431 - r26432; float r26434 = sin(r26433); float r26435 = phi2; float r26436 = cos(r26435); float r26437 = r26434 * r26436; float r26438 = phi1; float r26439 = cos(r26438); float r26440 = sin(r26435); float r26441 = r26439 * r26440; float r26442 = sin(r26438); float r26443 = r26442 * r26436; float r26444 = cos(r26433); float r26445 = r26443 * r26444; float r26446 = r26441 - r26445; float r26447 = atan2(r26437, r26446); return r26447; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r26448 = lambda1; double r26449 = lambda2; double r26450 = r26448 - r26449; double r26451 = sin(r26450); double r26452 = phi2; double r26453 = cos(r26452); double r26454 = r26451 * r26453; double r26455 = phi1; double r26456 = cos(r26455); double r26457 = sin(r26452); double r26458 = r26456 * r26457; double r26459 = sin(r26455); double r26460 = r26459 * r26453; double r26461 = cos(r26450); double r26462 = r26460 * r26461; double r26463 = r26458 - r26462; double r26464 = atan2(r26454, r26463); return r26464; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r26465 = lambda1; float r26466 = sin(r26465); float r26467 = lambda2; float r26468 = cos(r26467); float r26469 = r26466 * r26468; float r26470 = cos(r26465); float r26471 = -r26467; float r26472 = sin(r26471); float r26473 = r26470 * r26472; float r26474 = cbrt(r26473); float r26475 = r26474 * r26474; float r26476 = r26475 * r26474; float r26477 = r26469 + r26476; float r26478 = phi2; float r26479 = cos(r26478); float r26480 = r26477 * r26479; float r26481 = phi1; float r26482 = cos(r26481); float r26483 = sin(r26478); float r26484 = r26482 * r26483; float r26485 = sin(r26481); float r26486 = r26485 * r26479; float r26487 = r26470 * r26468; float r26488 = r26466 * r26472; float r26489 = r26487 - r26488; float r26490 = r26486 * r26489; float r26491 = r26484 - r26490; float r26492 = atan2(r26480, r26491); return r26492; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r26493 = lambda1; double r26494 = sin(r26493); double r26495 = lambda2; double r26496 = cos(r26495); double r26497 = r26494 * r26496; double r26498 = cos(r26493); double r26499 = -r26495; double r26500 = sin(r26499); double r26501 = r26498 * r26500; double r26502 = cbrt(r26501); double r26503 = r26502 * r26502; double r26504 = r26503 * r26502; double r26505 = r26497 + r26504; double r26506 = phi2; double r26507 = cos(r26506); double r26508 = r26505 * r26507; double r26509 = phi1; double r26510 = cos(r26509); double r26511 = sin(r26506); double r26512 = r26510 * r26511; double r26513 = sin(r26509); double r26514 = r26513 * r26507; double r26515 = r26498 * r26496; double r26516 = r26494 * r26500; double r26517 = r26515 - r26516; double r26518 = r26514 * r26517; double r26519 = r26512 - r26518; double r26520 = atan2(r26508, r26519); return r26520; } 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 r26521, r26522, r26523, r26524, r26525, r26526, r26527, r26528, r26529, r26530, r26531, r26532, r26533, r26534, r26535, r26536, r26537; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r26521); mpfr_init(r26522); mpfr_init(r26523); mpfr_init(r26524); mpfr_init(r26525); mpfr_init(r26526); mpfr_init(r26527); mpfr_init(r26528); mpfr_init(r26529); mpfr_init(r26530); mpfr_init(r26531); mpfr_init(r26532); mpfr_init(r26533); mpfr_init(r26534); mpfr_init(r26535); mpfr_init(r26536); mpfr_init(r26537); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r26521, lambda1, MPFR_RNDN); mpfr_set_d(r26522, lambda2, MPFR_RNDN); mpfr_sub(r26523, r26521, r26522, MPFR_RNDN); mpfr_sin(r26524, r26523, MPFR_RNDN); mpfr_set_d(r26525, phi2, MPFR_RNDN); mpfr_cos(r26526, r26525, MPFR_RNDN); mpfr_mul(r26527, r26524, r26526, MPFR_RNDN); mpfr_set_d(r26528, phi1, MPFR_RNDN); mpfr_cos(r26529, r26528, MPFR_RNDN); mpfr_sin(r26530, r26525, MPFR_RNDN); mpfr_mul(r26531, r26529, r26530, MPFR_RNDN); mpfr_sin(r26532, r26528, MPFR_RNDN); mpfr_mul(r26533, r26532, r26526, MPFR_RNDN); mpfr_cos(r26534, r26523, MPFR_RNDN); mpfr_mul(r26535, r26533, r26534, MPFR_RNDN); mpfr_sub(r26536, r26531, r26535, MPFR_RNDN); mpfr_atan2(r26537, r26527, r26536, MPFR_RNDN); return mpfr_get_d(r26537, MPFR_RNDN); } static mpfr_t r26538, r26539, r26540, r26541, r26542, r26543, r26544, r26545, r26546, r26547, r26548, r26549, r26550, r26551, r26552, r26553, r26554, r26555, r26556, r26557, r26558, r26559, r26560, r26561, r26562, r26563, r26564, r26565; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r26538); mpfr_init(r26539); mpfr_init(r26540); mpfr_init(r26541); mpfr_init(r26542); mpfr_init(r26543); mpfr_init(r26544); mpfr_init(r26545); mpfr_init(r26546); mpfr_init(r26547); mpfr_init(r26548); mpfr_init(r26549); mpfr_init(r26550); mpfr_init(r26551); mpfr_init(r26552); mpfr_init(r26553); mpfr_init(r26554); mpfr_init(r26555); mpfr_init(r26556); mpfr_init(r26557); mpfr_init(r26558); mpfr_init(r26559); mpfr_init(r26560); mpfr_init(r26561); mpfr_init(r26562); mpfr_init(r26563); mpfr_init(r26564); mpfr_init(r26565); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r26538, lambda1, MPFR_RNDN); mpfr_sin(r26539, r26538, MPFR_RNDN); mpfr_set_d(r26540, lambda2, MPFR_RNDN); mpfr_cos(r26541, r26540, MPFR_RNDN); mpfr_mul(r26542, r26539, r26541, MPFR_RNDN); mpfr_cos(r26543, r26538, MPFR_RNDN); mpfr_neg(r26544, r26540, MPFR_RNDN); mpfr_sin(r26545, r26544, MPFR_RNDN); mpfr_mul(r26546, r26543, r26545, MPFR_RNDN); mpfr_cbrt(r26547, r26546, MPFR_RNDN); mpfr_mul(r26548, r26547, r26547, MPFR_RNDN); mpfr_mul(r26549, r26548, r26547, MPFR_RNDN); mpfr_add(r26550, r26542, r26549, MPFR_RNDN); mpfr_set_d(r26551, phi2, MPFR_RNDN); mpfr_cos(r26552, r26551, MPFR_RNDN); mpfr_mul(r26553, r26550, r26552, MPFR_RNDN); mpfr_set_d(r26554, phi1, MPFR_RNDN); mpfr_cos(r26555, r26554, MPFR_RNDN); mpfr_sin(r26556, r26551, MPFR_RNDN); mpfr_mul(r26557, r26555, r26556, MPFR_RNDN); mpfr_sin(r26558, r26554, MPFR_RNDN); mpfr_mul(r26559, r26558, r26552, MPFR_RNDN); mpfr_mul(r26560, r26543, r26541, MPFR_RNDN); mpfr_mul(r26561, r26539, r26545, MPFR_RNDN); mpfr_sub(r26562, r26560, r26561, MPFR_RNDN); mpfr_mul(r26563, r26559, r26562, MPFR_RNDN); mpfr_sub(r26564, r26557, r26563, MPFR_RNDN); mpfr_atan2(r26565, r26553, r26564, MPFR_RNDN); return mpfr_get_d(r26565, MPFR_RNDN); } static mpfr_t r26566, r26567, r26568, r26569, r26570, r26571, r26572, r26573, r26574, r26575, r26576, r26577, r26578, r26579, r26580, r26581, r26582, r26583, r26584, r26585, r26586, r26587, r26588, r26589, r26590, r26591, r26592, r26593; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r26566); mpfr_init(r26567); mpfr_init(r26568); mpfr_init(r26569); mpfr_init(r26570); mpfr_init(r26571); mpfr_init(r26572); mpfr_init(r26573); mpfr_init(r26574); mpfr_init(r26575); mpfr_init(r26576); mpfr_init(r26577); mpfr_init(r26578); mpfr_init(r26579); mpfr_init(r26580); mpfr_init(r26581); mpfr_init(r26582); mpfr_init(r26583); mpfr_init(r26584); mpfr_init(r26585); mpfr_init(r26586); mpfr_init(r26587); mpfr_init(r26588); mpfr_init(r26589); mpfr_init(r26590); mpfr_init(r26591); mpfr_init(r26592); mpfr_init(r26593); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r26566, lambda1, MPFR_RNDN); mpfr_sin(r26567, r26566, MPFR_RNDN); mpfr_set_d(r26568, lambda2, MPFR_RNDN); mpfr_cos(r26569, r26568, MPFR_RNDN); mpfr_mul(r26570, r26567, r26569, MPFR_RNDN); mpfr_cos(r26571, r26566, MPFR_RNDN); mpfr_neg(r26572, r26568, MPFR_RNDN); mpfr_sin(r26573, r26572, MPFR_RNDN); mpfr_mul(r26574, r26571, r26573, MPFR_RNDN); mpfr_cbrt(r26575, r26574, MPFR_RNDN); mpfr_mul(r26576, r26575, r26575, MPFR_RNDN); mpfr_mul(r26577, r26576, r26575, MPFR_RNDN); mpfr_add(r26578, r26570, r26577, MPFR_RNDN); mpfr_set_d(r26579, phi2, MPFR_RNDN); mpfr_cos(r26580, r26579, MPFR_RNDN); mpfr_mul(r26581, r26578, r26580, MPFR_RNDN); mpfr_set_d(r26582, phi1, MPFR_RNDN); mpfr_cos(r26583, r26582, MPFR_RNDN); mpfr_sin(r26584, r26579, MPFR_RNDN); mpfr_mul(r26585, r26583, r26584, MPFR_RNDN); mpfr_sin(r26586, r26582, MPFR_RNDN); mpfr_mul(r26587, r26586, r26580, MPFR_RNDN); mpfr_mul(r26588, r26571, r26569, MPFR_RNDN); mpfr_mul(r26589, r26567, r26573, MPFR_RNDN); mpfr_sub(r26590, r26588, r26589, MPFR_RNDN); mpfr_mul(r26591, r26587, r26590, MPFR_RNDN); mpfr_sub(r26592, r26585, r26591, MPFR_RNDN); mpfr_atan2(r26593, r26581, r26592, MPFR_RNDN); return mpfr_get_d(r26593, MPFR_RNDN); }