#include #include #include #include #include char *name = "Midpoint on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r22478 = lambda1; float r22479 = phi2; float r22480 = cos(r22479); float r22481 = lambda2; float r22482 = r22478 - r22481; float r22483 = sin(r22482); float r22484 = r22480 * r22483; float r22485 = phi1; float r22486 = cos(r22485); float r22487 = cos(r22482); float r22488 = r22480 * r22487; float r22489 = r22486 + r22488; float r22490 = atan2(r22484, r22489); float r22491 = r22478 + r22490; return r22491; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r22492 = lambda1; double r22493 = phi2; double r22494 = cos(r22493); double r22495 = lambda2; double r22496 = r22492 - r22495; double r22497 = sin(r22496); double r22498 = r22494 * r22497; double r22499 = phi1; double r22500 = cos(r22499); double r22501 = cos(r22496); double r22502 = r22494 * r22501; double r22503 = r22500 + r22502; double r22504 = atan2(r22498, r22503); double r22505 = r22492 + r22504; return r22505; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r22506 = phi2; float r22507 = cos(r22506); float r22508 = lambda1; float r22509 = sin(r22508); float r22510 = lambda2; float r22511 = cos(r22510); float r22512 = r22509 * r22511; float r22513 = cos(r22508); float r22514 = -r22510; float r22515 = sin(r22514); float r22516 = r22513 * r22515; float r22517 = r22512 + r22516; float r22518 = r22507 * r22517; float r22519 = sin(r22510); float r22520 = r22519 * r22509; float r22521 = fma(r22511, r22513, r22520); float r22522 = phi1; float r22523 = cos(r22522); float r22524 = fma(r22521, r22507, r22523); float r22525 = expm1(r22524); float r22526 = exp(r22525); float r22527 = log(r22526); float r22528 = log1p(r22527); float r22529 = atan2(r22518, r22528); float r22530 = r22529 + r22508; return r22530; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r22531 = phi2; double r22532 = cos(r22531); double r22533 = lambda1; double r22534 = sin(r22533); double r22535 = lambda2; double r22536 = cos(r22535); double r22537 = r22534 * r22536; double r22538 = cos(r22533); double r22539 = -r22535; double r22540 = sin(r22539); double r22541 = r22538 * r22540; double r22542 = r22537 + r22541; double r22543 = r22532 * r22542; double r22544 = sin(r22535); double r22545 = r22544 * r22534; double r22546 = fma(r22536, r22538, r22545); double r22547 = phi1; double r22548 = cos(r22547); double r22549 = fma(r22546, r22532, r22548); double r22550 = expm1(r22549); double r22551 = exp(r22550); double r22552 = log(r22551); double r22553 = log1p(r22552); double r22554 = atan2(r22543, r22553); double r22555 = r22554 + r22533; return r22555; } 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 r22556, r22557, r22558, r22559, r22560, r22561, r22562, r22563, r22564, r22565, r22566, r22567, r22568, r22569; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r22556); mpfr_init(r22557); mpfr_init(r22558); mpfr_init(r22559); mpfr_init(r22560); mpfr_init(r22561); mpfr_init(r22562); mpfr_init(r22563); mpfr_init(r22564); mpfr_init(r22565); mpfr_init(r22566); mpfr_init(r22567); mpfr_init(r22568); mpfr_init(r22569); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r22556, lambda1, MPFR_RNDN); mpfr_set_d(r22557, phi2, MPFR_RNDN); mpfr_cos(r22558, r22557, MPFR_RNDN); mpfr_set_d(r22559, lambda2, MPFR_RNDN); mpfr_sub(r22560, r22556, r22559, MPFR_RNDN); mpfr_sin(r22561, r22560, MPFR_RNDN); mpfr_mul(r22562, r22558, r22561, MPFR_RNDN); mpfr_set_d(r22563, phi1, MPFR_RNDN); mpfr_cos(r22564, r22563, MPFR_RNDN); mpfr_cos(r22565, r22560, MPFR_RNDN); mpfr_mul(r22566, r22558, r22565, MPFR_RNDN); mpfr_add(r22567, r22564, r22566, MPFR_RNDN); mpfr_atan2(r22568, r22562, r22567, MPFR_RNDN); mpfr_add(r22569, r22556, r22568, MPFR_RNDN); return mpfr_get_d(r22569, MPFR_RNDN); } static mpfr_t r22570, r22571, r22572, r22573, r22574, r22575, r22576, r22577, r22578, r22579, r22580, r22581, r22582, r22583, r22584, r22585, r22586, r22587, r22588, r22589, r22590, r22591, r22592, r22593, r22594; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r22570); mpfr_init(r22571); mpfr_init(r22572); mpfr_init(r22573); mpfr_init(r22574); mpfr_init(r22575); mpfr_init(r22576); mpfr_init(r22577); mpfr_init(r22578); mpfr_init(r22579); mpfr_init(r22580); mpfr_init(r22581); mpfr_init(r22582); mpfr_init(r22583); mpfr_init(r22584); mpfr_init(r22585); mpfr_init(r22586); mpfr_init(r22587); mpfr_init(r22588); mpfr_init(r22589); mpfr_init(r22590); mpfr_init(r22591); mpfr_init(r22592); mpfr_init(r22593); mpfr_init(r22594); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r22570, phi2, MPFR_RNDN); mpfr_cos(r22571, r22570, MPFR_RNDN); mpfr_set_d(r22572, lambda1, MPFR_RNDN); mpfr_sin(r22573, r22572, MPFR_RNDN); mpfr_set_d(r22574, lambda2, MPFR_RNDN); mpfr_cos(r22575, r22574, MPFR_RNDN); mpfr_mul(r22576, r22573, r22575, MPFR_RNDN); mpfr_cos(r22577, r22572, MPFR_RNDN); mpfr_neg(r22578, r22574, MPFR_RNDN); mpfr_sin(r22579, r22578, MPFR_RNDN); mpfr_mul(r22580, r22577, r22579, MPFR_RNDN); mpfr_add(r22581, r22576, r22580, MPFR_RNDN); mpfr_mul(r22582, r22571, r22581, MPFR_RNDN); mpfr_sin(r22583, r22574, MPFR_RNDN); mpfr_mul(r22584, r22583, r22573, MPFR_RNDN); mpfr_fma(r22585, r22575, r22577, r22584, MPFR_RNDN); mpfr_set_d(r22586, phi1, MPFR_RNDN); mpfr_cos(r22587, r22586, MPFR_RNDN); mpfr_fma(r22588, r22585, r22571, r22587, MPFR_RNDN); mpfr_expm1(r22589, r22588, MPFR_RNDN); mpfr_exp(r22590, r22589, MPFR_RNDN); mpfr_log(r22591, r22590, MPFR_RNDN); mpfr_log1p(r22592, r22591, MPFR_RNDN); mpfr_atan2(r22593, r22582, r22592, MPFR_RNDN); mpfr_add(r22594, r22593, r22572, MPFR_RNDN); return mpfr_get_d(r22594, MPFR_RNDN); } static mpfr_t r22595, r22596, r22597, r22598, r22599, r22600, r22601, r22602, r22603, r22604, r22605, r22606, r22607, r22608, r22609, r22610, r22611, r22612, r22613, r22614, r22615, r22616, r22617, r22618, r22619; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r22595); mpfr_init(r22596); mpfr_init(r22597); mpfr_init(r22598); mpfr_init(r22599); mpfr_init(r22600); mpfr_init(r22601); mpfr_init(r22602); mpfr_init(r22603); mpfr_init(r22604); mpfr_init(r22605); mpfr_init(r22606); mpfr_init(r22607); mpfr_init(r22608); mpfr_init(r22609); mpfr_init(r22610); mpfr_init(r22611); mpfr_init(r22612); mpfr_init(r22613); mpfr_init(r22614); mpfr_init(r22615); mpfr_init(r22616); mpfr_init(r22617); mpfr_init(r22618); mpfr_init(r22619); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r22595, phi2, MPFR_RNDN); mpfr_cos(r22596, r22595, MPFR_RNDN); mpfr_set_d(r22597, lambda1, MPFR_RNDN); mpfr_sin(r22598, r22597, MPFR_RNDN); mpfr_set_d(r22599, lambda2, MPFR_RNDN); mpfr_cos(r22600, r22599, MPFR_RNDN); mpfr_mul(r22601, r22598, r22600, MPFR_RNDN); mpfr_cos(r22602, r22597, MPFR_RNDN); mpfr_neg(r22603, r22599, MPFR_RNDN); mpfr_sin(r22604, r22603, MPFR_RNDN); mpfr_mul(r22605, r22602, r22604, MPFR_RNDN); mpfr_add(r22606, r22601, r22605, MPFR_RNDN); mpfr_mul(r22607, r22596, r22606, MPFR_RNDN); mpfr_sin(r22608, r22599, MPFR_RNDN); mpfr_mul(r22609, r22608, r22598, MPFR_RNDN); mpfr_fma(r22610, r22600, r22602, r22609, MPFR_RNDN); mpfr_set_d(r22611, phi1, MPFR_RNDN); mpfr_cos(r22612, r22611, MPFR_RNDN); mpfr_fma(r22613, r22610, r22596, r22612, MPFR_RNDN); mpfr_expm1(r22614, r22613, MPFR_RNDN); mpfr_exp(r22615, r22614, MPFR_RNDN); mpfr_log(r22616, r22615, MPFR_RNDN); mpfr_log1p(r22617, r22616, MPFR_RNDN); mpfr_atan2(r22618, r22607, r22617, MPFR_RNDN); mpfr_add(r22619, r22618, r22597, MPFR_RNDN); return mpfr_get_d(r22619, MPFR_RNDN); }