#include #include #include #include #include char *name = "Midpoint on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r26543 = lambda1; float r26544 = phi2; float r26545 = cos(r26544); float r26546 = lambda2; float r26547 = r26543 - r26546; float r26548 = sin(r26547); float r26549 = r26545 * r26548; float r26550 = phi1; float r26551 = cos(r26550); float r26552 = cos(r26547); float r26553 = r26545 * r26552; float r26554 = r26551 + r26553; float r26555 = atan2(r26549, r26554); float r26556 = r26543 + r26555; return r26556; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r26557 = lambda1; double r26558 = phi2; double r26559 = cos(r26558); double r26560 = lambda2; double r26561 = r26557 - r26560; double r26562 = sin(r26561); double r26563 = r26559 * r26562; double r26564 = phi1; double r26565 = cos(r26564); double r26566 = cos(r26561); double r26567 = r26559 * r26566; double r26568 = r26565 + r26567; double r26569 = atan2(r26563, r26568); double r26570 = r26557 + r26569; return r26570; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r26571 = lambda1; float r26572 = phi2; float r26573 = cos(r26572); float r26574 = lambda2; float r26575 = r26571 - r26574; float r26576 = sin(r26575); float r26577 = r26573 * r26576; float r26578 = cos(r26575); float r26579 = r26578 * r26573; float r26580 = r26579 * (r26579 * r26579); float r26581 = phi1; float r26582 = cos(r26581); float r26583 = r26582 * (r26582 * r26582); float r26584 = r26580 + r26583; float r26585 = r26582 * r26582; float r26586 = r26573 * r26578; float r26587 = r26586 * r26586; float r26588 = r26582 * r26586; float r26589 = r26587 - r26588; float r26590 = r26585 + r26589; float r26591 = r26584 / r26590; float r26592 = atan2(r26577, r26591); float r26593 = r26571 + r26592; return r26593; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r26594 = lambda1; double r26595 = phi2; double r26596 = cos(r26595); double r26597 = lambda2; double r26598 = r26594 - r26597; double r26599 = sin(r26598); double r26600 = r26596 * r26599; double r26601 = cos(r26598); double r26602 = r26601 * r26596; double r26603 = r26602 * (r26602 * r26602); double r26604 = phi1; double r26605 = cos(r26604); double r26606 = r26605 * (r26605 * r26605); double r26607 = r26603 + r26606; double r26608 = r26605 * r26605; double r26609 = r26596 * r26601; double r26610 = r26609 * r26609; double r26611 = r26605 * r26609; double r26612 = r26610 - r26611; double r26613 = r26608 + r26612; double r26614 = r26607 / r26613; double r26615 = atan2(r26600, r26614); double r26616 = r26594 + r26615; return r26616; } 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 r26617, r26618, r26619, r26620, r26621, r26622, r26623, r26624, r26625, r26626, r26627, r26628, r26629, r26630; void setup_mpfr_f_im() { mpfr_set_default_prec(144); 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); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r26617, lambda1, MPFR_RNDN); mpfr_set_d(r26618, phi2, MPFR_RNDN); mpfr_cos(r26619, r26618, MPFR_RNDN); mpfr_set_d(r26620, lambda2, MPFR_RNDN); mpfr_sub(r26621, r26617, r26620, MPFR_RNDN); mpfr_sin(r26622, r26621, MPFR_RNDN); mpfr_mul(r26623, r26619, r26622, MPFR_RNDN); mpfr_set_d(r26624, phi1, MPFR_RNDN); mpfr_cos(r26625, r26624, MPFR_RNDN); mpfr_cos(r26626, r26621, MPFR_RNDN); mpfr_mul(r26627, r26619, r26626, MPFR_RNDN); mpfr_add(r26628, r26625, r26627, MPFR_RNDN); mpfr_atan2(r26629, r26623, r26628, MPFR_RNDN); mpfr_add(r26630, r26617, r26629, MPFR_RNDN); return mpfr_get_d(r26630, MPFR_RNDN); } static mpfr_t r26631, r26632, r26633, r26634, r26635, r26636, r26637, r26638, r26639, r26640, r26641, r26642, r26643, r26644, r26645, r26646, r26647, r26648, r26649, r26650, r26651, r26652, r26653; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r26631); mpfr_init(r26632); mpfr_init(r26633); mpfr_init(r26634); mpfr_init(r26635); mpfr_init(r26636); mpfr_init(r26637); mpfr_init(r26638); mpfr_init(r26639); mpfr_init(r26640); mpfr_init(r26641); mpfr_init(r26642); mpfr_init(r26643); mpfr_init(r26644); mpfr_init(r26645); mpfr_init(r26646); mpfr_init(r26647); mpfr_init(r26648); mpfr_init(r26649); mpfr_init(r26650); mpfr_init(r26651); mpfr_init(r26652); mpfr_init(r26653); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r26631, lambda1, MPFR_RNDN); mpfr_set_d(r26632, phi2, MPFR_RNDN); mpfr_cos(r26633, r26632, MPFR_RNDN); mpfr_set_d(r26634, lambda2, MPFR_RNDN); mpfr_sub(r26635, r26631, r26634, MPFR_RNDN); mpfr_sin(r26636, r26635, MPFR_RNDN); mpfr_mul(r26637, r26633, r26636, MPFR_RNDN); mpfr_cos(r26638, r26635, MPFR_RNDN); mpfr_mul(r26639, r26638, r26633, MPFR_RNDN); mpfr_mul(r26640, r26639, r26639, MPFR_RNDN); mpfr_mul(r26640, r26640, r26639, MPFR_RNDN); mpfr_set_d(r26641, phi1, MPFR_RNDN); mpfr_cos(r26642, r26641, MPFR_RNDN); mpfr_mul(r26643, r26642, r26642, MPFR_RNDN); mpfr_mul(r26643, r26643, r26642, MPFR_RNDN); mpfr_add(r26644, r26640, r26643, MPFR_RNDN); mpfr_sqr(r26645, r26642, MPFR_RNDN); mpfr_mul(r26646, r26633, r26638, MPFR_RNDN); mpfr_sqr(r26647, r26646, MPFR_RNDN); mpfr_mul(r26648, r26642, r26646, MPFR_RNDN); mpfr_sub(r26649, r26647, r26648, MPFR_RNDN); mpfr_add(r26650, r26645, r26649, MPFR_RNDN); mpfr_div(r26651, r26644, r26650, MPFR_RNDN); mpfr_atan2(r26652, r26637, r26651, MPFR_RNDN); mpfr_add(r26653, r26631, r26652, MPFR_RNDN); return mpfr_get_d(r26653, MPFR_RNDN); } static mpfr_t r26654, r26655, r26656, r26657, r26658, r26659, r26660, r26661, r26662, r26663, r26664, r26665, r26666, r26667, r26668, r26669, r26670, r26671, r26672, r26673, r26674, r26675, r26676; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r26654); mpfr_init(r26655); mpfr_init(r26656); mpfr_init(r26657); mpfr_init(r26658); mpfr_init(r26659); mpfr_init(r26660); mpfr_init(r26661); mpfr_init(r26662); mpfr_init(r26663); mpfr_init(r26664); mpfr_init(r26665); mpfr_init(r26666); mpfr_init(r26667); mpfr_init(r26668); mpfr_init(r26669); mpfr_init(r26670); mpfr_init(r26671); mpfr_init(r26672); mpfr_init(r26673); mpfr_init(r26674); mpfr_init(r26675); mpfr_init(r26676); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r26654, lambda1, MPFR_RNDN); mpfr_set_d(r26655, phi2, MPFR_RNDN); mpfr_cos(r26656, r26655, MPFR_RNDN); mpfr_set_d(r26657, lambda2, MPFR_RNDN); mpfr_sub(r26658, r26654, r26657, MPFR_RNDN); mpfr_sin(r26659, r26658, MPFR_RNDN); mpfr_mul(r26660, r26656, r26659, MPFR_RNDN); mpfr_cos(r26661, r26658, MPFR_RNDN); mpfr_mul(r26662, r26661, r26656, MPFR_RNDN); mpfr_mul(r26663, r26662, r26662, MPFR_RNDN); mpfr_mul(r26663, r26663, r26662, MPFR_RNDN); mpfr_set_d(r26664, phi1, MPFR_RNDN); mpfr_cos(r26665, r26664, MPFR_RNDN); mpfr_mul(r26666, r26665, r26665, MPFR_RNDN); mpfr_mul(r26666, r26666, r26665, MPFR_RNDN); mpfr_add(r26667, r26663, r26666, MPFR_RNDN); mpfr_sqr(r26668, r26665, MPFR_RNDN); mpfr_mul(r26669, r26656, r26661, MPFR_RNDN); mpfr_sqr(r26670, r26669, MPFR_RNDN); mpfr_mul(r26671, r26665, r26669, MPFR_RNDN); mpfr_sub(r26672, r26670, r26671, MPFR_RNDN); mpfr_add(r26673, r26668, r26672, MPFR_RNDN); mpfr_div(r26674, r26667, r26673, MPFR_RNDN); mpfr_atan2(r26675, r26660, r26674, MPFR_RNDN); mpfr_add(r26676, r26654, r26675, MPFR_RNDN); return mpfr_get_d(r26676, MPFR_RNDN); }