#include #include #include #include #include char *name = "Equirectangular approximation to distance on a great circle"; double f_if(float R, float lambda1, float lambda2, float phi1, float phi2) { float r22507 = R; float r22508 = lambda1; float r22509 = lambda2; float r22510 = r22508 - r22509; float r22511 = phi1; float r22512 = phi2; float r22513 = r22511 + r22512; float r22514 = 2; float r22515 = r22513 / r22514; float r22516 = cos(r22515); float r22517 = r22510 * r22516; float r22518 = r22517 * r22517; float r22519 = r22511 - r22512; float r22520 = r22519 * r22519; float r22521 = r22518 + r22520; float r22522 = sqrt(r22521); float r22523 = r22507 * r22522; return r22523; } double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) { double r22524 = R; double r22525 = lambda1; double r22526 = lambda2; double r22527 = r22525 - r22526; double r22528 = phi1; double r22529 = phi2; double r22530 = r22528 + r22529; double r22531 = 2; double r22532 = r22530 / r22531; double r22533 = cos(r22532); double r22534 = r22527 * r22533; double r22535 = r22534 * r22534; double r22536 = r22528 - r22529; double r22537 = r22536 * r22536; double r22538 = r22535 + r22537; double r22539 = sqrt(r22538); double r22540 = r22524 * r22539; return r22540; } double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) { float r22541 = phi1; float r22542 = phi2; float r22543 = r22541 + r22542; float r22544 = 2; float r22545 = r22543 / r22544; float r22546 = cos(r22545); float r22547 = lambda1; float r22548 = lambda2; float r22549 = r22547 - r22548; float r22550 = r22546 * r22549; float r22551 = r22541 - r22542; float r22552 = hypot(r22550, r22551); float r22553 = R; float r22554 = r22552 * r22553; return r22554; } double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) { double r22555 = phi1; double r22556 = phi2; double r22557 = r22555 + r22556; double r22558 = 2; double r22559 = r22557 / r22558; double r22560 = cos(r22559); double r22561 = lambda1; double r22562 = lambda2; double r22563 = r22561 - r22562; double r22564 = r22560 * r22563; double r22565 = r22555 - r22556; double r22566 = hypot(r22564, r22565); double r22567 = R; double r22568 = r22566 * r22567; return r22568; } 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 r22569, r22570, r22571, r22572, r22573, r22574, r22575, r22576, r22577, r22578, r22579, r22580, r22581, r22582, r22583, r22584, r22585; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r22569); mpfr_init(r22570); mpfr_init(r22571); mpfr_init(r22572); mpfr_init(r22573); mpfr_init(r22574); mpfr_init(r22575); mpfr_init_set_str(r22576, "2", 10, MPFR_RNDN); 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); } double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r22569, R, MPFR_RNDN); mpfr_set_d(r22570, lambda1, MPFR_RNDN); mpfr_set_d(r22571, lambda2, MPFR_RNDN); mpfr_sub(r22572, r22570, r22571, MPFR_RNDN); mpfr_set_d(r22573, phi1, MPFR_RNDN); mpfr_set_d(r22574, phi2, MPFR_RNDN); mpfr_add(r22575, r22573, r22574, MPFR_RNDN); ; mpfr_div(r22577, r22575, r22576, MPFR_RNDN); mpfr_cos(r22578, r22577, MPFR_RNDN); mpfr_mul(r22579, r22572, r22578, MPFR_RNDN); mpfr_mul(r22580, r22579, r22579, MPFR_RNDN); mpfr_sub(r22581, r22573, r22574, MPFR_RNDN); mpfr_mul(r22582, r22581, r22581, MPFR_RNDN); mpfr_add(r22583, r22580, r22582, MPFR_RNDN); mpfr_sqrt(r22584, r22583, MPFR_RNDN); mpfr_mul(r22585, r22569, r22584, MPFR_RNDN); return mpfr_get_d(r22585, MPFR_RNDN); } static mpfr_t r22586, r22587, r22588, r22589, r22590, r22591, r22592, r22593, r22594, r22595, r22596, r22597, r22598, r22599; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r22586); mpfr_init(r22587); mpfr_init(r22588); mpfr_init_set_str(r22589, "2", 10, MPFR_RNDN); mpfr_init(r22590); mpfr_init(r22591); mpfr_init(r22592); mpfr_init(r22593); mpfr_init(r22594); mpfr_init(r22595); mpfr_init(r22596); mpfr_init(r22597); mpfr_init(r22598); mpfr_init(r22599); } double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r22586, phi1, MPFR_RNDN); mpfr_set_d(r22587, phi2, MPFR_RNDN); mpfr_add(r22588, r22586, r22587, MPFR_RNDN); ; mpfr_div(r22590, r22588, r22589, MPFR_RNDN); mpfr_cos(r22591, r22590, MPFR_RNDN); mpfr_set_d(r22592, lambda1, MPFR_RNDN); mpfr_set_d(r22593, lambda2, MPFR_RNDN); mpfr_sub(r22594, r22592, r22593, MPFR_RNDN); mpfr_mul(r22595, r22591, r22594, MPFR_RNDN); mpfr_sub(r22596, r22586, r22587, MPFR_RNDN); mpfr_hypot(r22597, r22595, r22596, MPFR_RNDN); mpfr_set_d(r22598, R, MPFR_RNDN); mpfr_mul(r22599, r22597, r22598, MPFR_RNDN); return mpfr_get_d(r22599, MPFR_RNDN); } static mpfr_t r22600, r22601, r22602, r22603, r22604, r22605, r22606, r22607, r22608, r22609, r22610, r22611, r22612, r22613; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r22600); mpfr_init(r22601); mpfr_init(r22602); mpfr_init_set_str(r22603, "2", 10, MPFR_RNDN); 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); } double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r22600, phi1, MPFR_RNDN); mpfr_set_d(r22601, phi2, MPFR_RNDN); mpfr_add(r22602, r22600, r22601, MPFR_RNDN); ; mpfr_div(r22604, r22602, r22603, MPFR_RNDN); mpfr_cos(r22605, r22604, MPFR_RNDN); mpfr_set_d(r22606, lambda1, MPFR_RNDN); mpfr_set_d(r22607, lambda2, MPFR_RNDN); mpfr_sub(r22608, r22606, r22607, MPFR_RNDN); mpfr_mul(r22609, r22605, r22608, MPFR_RNDN); mpfr_sub(r22610, r22600, r22601, MPFR_RNDN); mpfr_hypot(r22611, r22609, r22610, MPFR_RNDN); mpfr_set_d(r22612, R, MPFR_RNDN); mpfr_mul(r22613, r22611, r22612, MPFR_RNDN); return mpfr_get_d(r22613, MPFR_RNDN); }