#include #include #include #include #include char *name = "Bearing on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r26478 = lambda1; float r26479 = lambda2; float r26480 = r26478 - r26479; float r26481 = sin(r26480); float r26482 = phi2; float r26483 = cos(r26482); float r26484 = r26481 * r26483; float r26485 = phi1; float r26486 = cos(r26485); float r26487 = sin(r26482); float r26488 = r26486 * r26487; float r26489 = sin(r26485); float r26490 = r26489 * r26483; float r26491 = cos(r26480); float r26492 = r26490 * r26491; float r26493 = r26488 - r26492; float r26494 = atan2(r26484, r26493); return r26494; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r26495 = lambda1; double r26496 = lambda2; double r26497 = r26495 - r26496; double r26498 = sin(r26497); double r26499 = phi2; double r26500 = cos(r26499); double r26501 = r26498 * r26500; double r26502 = phi1; double r26503 = cos(r26502); double r26504 = sin(r26499); double r26505 = r26503 * r26504; double r26506 = sin(r26502); double r26507 = r26506 * r26500; double r26508 = cos(r26497); double r26509 = r26507 * r26508; double r26510 = r26505 - r26509; double r26511 = atan2(r26501, r26510); return r26511; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r26512 = lambda1; float r26513 = sin(r26512); float r26514 = lambda2; float r26515 = cos(r26514); float r26516 = r26513 * r26515; float r26517 = cos(r26512); float r26518 = -r26514; float r26519 = sin(r26518); float r26520 = r26517 * r26519; float r26521 = r26516 + r26520; float r26522 = phi2; float r26523 = cos(r26522); float r26524 = r26521 * r26523; float r26525 = phi1; float r26526 = cos(r26525); float r26527 = sin(r26522); float r26528 = r26526 * r26527; float r26529 = sin(r26525); float r26530 = r26529 * r26523; float r26531 = sin(r26514); float r26532 = r26513 * r26531; float r26533 = fma(r26515, r26517, r26532); float r26534 = 3; float r26535 = pow(r26533, r26534); float r26536 = cbrt(r26535); float r26537 = r26530 * r26536; float r26538 = r26528 - r26537; float r26539 = atan2(r26524, r26538); return r26539; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r26540 = lambda1; double r26541 = sin(r26540); double r26542 = lambda2; double r26543 = cos(r26542); double r26544 = r26541 * r26543; double r26545 = cos(r26540); double r26546 = -r26542; double r26547 = sin(r26546); double r26548 = r26545 * r26547; double r26549 = r26544 + r26548; double r26550 = phi2; double r26551 = cos(r26550); double r26552 = r26549 * r26551; double r26553 = phi1; double r26554 = cos(r26553); double r26555 = sin(r26550); double r26556 = r26554 * r26555; double r26557 = sin(r26553); double r26558 = r26557 * r26551; double r26559 = sin(r26542); double r26560 = r26541 * r26559; double r26561 = fma(r26543, r26545, r26560); double r26562 = 3; double r26563 = pow(r26561, r26562); double r26564 = cbrt(r26563); double r26565 = r26558 * r26564; double r26566 = r26556 - r26565; double r26567 = atan2(r26552, r26566); return r26567; } 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 r26568, r26569, r26570, r26571, r26572, r26573, r26574, r26575, r26576, r26577, r26578, r26579, r26580, r26581, r26582, r26583, r26584; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); 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); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r26568, lambda1, MPFR_RNDN); mpfr_set_d(r26569, lambda2, MPFR_RNDN); mpfr_sub(r26570, r26568, r26569, MPFR_RNDN); mpfr_sin(r26571, r26570, MPFR_RNDN); mpfr_set_d(r26572, phi2, MPFR_RNDN); mpfr_cos(r26573, r26572, MPFR_RNDN); mpfr_mul(r26574, r26571, r26573, MPFR_RNDN); mpfr_set_d(r26575, phi1, MPFR_RNDN); mpfr_cos(r26576, r26575, MPFR_RNDN); mpfr_sin(r26577, r26572, MPFR_RNDN); mpfr_mul(r26578, r26576, r26577, MPFR_RNDN); mpfr_sin(r26579, r26575, MPFR_RNDN); mpfr_mul(r26580, r26579, r26573, MPFR_RNDN); mpfr_cos(r26581, r26570, MPFR_RNDN); mpfr_mul(r26582, r26580, r26581, MPFR_RNDN); mpfr_sub(r26583, r26578, r26582, MPFR_RNDN); mpfr_atan2(r26584, r26574, r26583, MPFR_RNDN); return mpfr_get_d(r26584, MPFR_RNDN); } static mpfr_t r26585, r26586, r26587, r26588, r26589, r26590, r26591, r26592, r26593, r26594, r26595, r26596, r26597, r26598, r26599, r26600, r26601, r26602, r26603, r26604, r26605, r26606, r26607, r26608, r26609, r26610, r26611, r26612; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); 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); mpfr_init(r26594); mpfr_init(r26595); mpfr_init(r26596); mpfr_init(r26597); mpfr_init(r26598); mpfr_init(r26599); mpfr_init(r26600); mpfr_init(r26601); mpfr_init(r26602); mpfr_init(r26603); mpfr_init(r26604); mpfr_init(r26605); mpfr_init(r26606); mpfr_init_set_str(r26607, "3", 10, MPFR_RNDN); mpfr_init(r26608); mpfr_init(r26609); mpfr_init(r26610); mpfr_init(r26611); mpfr_init(r26612); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r26585, lambda1, MPFR_RNDN); mpfr_sin(r26586, r26585, MPFR_RNDN); mpfr_set_d(r26587, lambda2, MPFR_RNDN); mpfr_cos(r26588, r26587, MPFR_RNDN); mpfr_mul(r26589, r26586, r26588, MPFR_RNDN); mpfr_cos(r26590, r26585, MPFR_RNDN); mpfr_neg(r26591, r26587, MPFR_RNDN); mpfr_sin(r26592, r26591, MPFR_RNDN); mpfr_mul(r26593, r26590, r26592, MPFR_RNDN); mpfr_add(r26594, r26589, r26593, MPFR_RNDN); mpfr_set_d(r26595, phi2, MPFR_RNDN); mpfr_cos(r26596, r26595, MPFR_RNDN); mpfr_mul(r26597, r26594, r26596, MPFR_RNDN); mpfr_set_d(r26598, phi1, MPFR_RNDN); mpfr_cos(r26599, r26598, MPFR_RNDN); mpfr_sin(r26600, r26595, MPFR_RNDN); mpfr_mul(r26601, r26599, r26600, MPFR_RNDN); mpfr_sin(r26602, r26598, MPFR_RNDN); mpfr_mul(r26603, r26602, r26596, MPFR_RNDN); mpfr_sin(r26604, r26587, MPFR_RNDN); mpfr_mul(r26605, r26586, r26604, MPFR_RNDN); mpfr_fma(r26606, r26588, r26590, r26605, MPFR_RNDN); ; mpfr_pow(r26608, r26606, r26607, MPFR_RNDN); mpfr_cbrt(r26609, r26608, MPFR_RNDN); mpfr_mul(r26610, r26603, r26609, MPFR_RNDN); mpfr_sub(r26611, r26601, r26610, MPFR_RNDN); mpfr_atan2(r26612, r26597, r26611, MPFR_RNDN); return mpfr_get_d(r26612, MPFR_RNDN); } static mpfr_t r26613, r26614, r26615, r26616, r26617, r26618, r26619, r26620, r26621, r26622, r26623, r26624, r26625, r26626, r26627, r26628, r26629, r26630, r26631, r26632, r26633, r26634, r26635, r26636, r26637, r26638, r26639, r26640; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r26613); mpfr_init(r26614); mpfr_init(r26615); mpfr_init(r26616); mpfr_init(r26617); mpfr_init(r26618); mpfr_init(r26619); mpfr_init(r26620); mpfr_init(r26621); mpfr_init(r26622); mpfr_init(r26623); mpfr_init(r26624); mpfr_init(r26625); mpfr_init(r26626); mpfr_init(r26627); mpfr_init(r26628); mpfr_init(r26629); mpfr_init(r26630); mpfr_init(r26631); mpfr_init(r26632); mpfr_init(r26633); mpfr_init(r26634); mpfr_init_set_str(r26635, "3", 10, MPFR_RNDN); mpfr_init(r26636); mpfr_init(r26637); mpfr_init(r26638); mpfr_init(r26639); mpfr_init(r26640); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r26613, lambda1, MPFR_RNDN); mpfr_sin(r26614, r26613, MPFR_RNDN); mpfr_set_d(r26615, lambda2, MPFR_RNDN); mpfr_cos(r26616, r26615, MPFR_RNDN); mpfr_mul(r26617, r26614, r26616, MPFR_RNDN); mpfr_cos(r26618, r26613, MPFR_RNDN); mpfr_neg(r26619, r26615, MPFR_RNDN); mpfr_sin(r26620, r26619, MPFR_RNDN); mpfr_mul(r26621, r26618, r26620, MPFR_RNDN); mpfr_add(r26622, r26617, r26621, MPFR_RNDN); mpfr_set_d(r26623, phi2, MPFR_RNDN); mpfr_cos(r26624, r26623, MPFR_RNDN); mpfr_mul(r26625, r26622, r26624, MPFR_RNDN); mpfr_set_d(r26626, phi1, MPFR_RNDN); mpfr_cos(r26627, r26626, MPFR_RNDN); mpfr_sin(r26628, r26623, MPFR_RNDN); mpfr_mul(r26629, r26627, r26628, MPFR_RNDN); mpfr_sin(r26630, r26626, MPFR_RNDN); mpfr_mul(r26631, r26630, r26624, MPFR_RNDN); mpfr_sin(r26632, r26615, MPFR_RNDN); mpfr_mul(r26633, r26614, r26632, MPFR_RNDN); mpfr_fma(r26634, r26616, r26618, r26633, MPFR_RNDN); ; mpfr_pow(r26636, r26634, r26635, MPFR_RNDN); mpfr_cbrt(r26637, r26636, MPFR_RNDN); mpfr_mul(r26638, r26631, r26637, MPFR_RNDN); mpfr_sub(r26639, r26629, r26638, MPFR_RNDN); mpfr_atan2(r26640, r26625, r26639, MPFR_RNDN); return mpfr_get_d(r26640, MPFR_RNDN); }