#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 r22658 = R; float r22659 = lambda1; float r22660 = lambda2; float r22661 = r22659 - r22660; float r22662 = phi1; float r22663 = phi2; float r22664 = r22662 + r22663; float r22665 = 2; float r22666 = r22664 / r22665; float r22667 = cos(r22666); float r22668 = r22661 * r22667; float r22669 = r22668 * r22668; float r22670 = r22662 - r22663; float r22671 = r22670 * r22670; float r22672 = r22669 + r22671; float r22673 = sqrt(r22672); float r22674 = r22658 * r22673; return r22674; } double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) { double r22675 = R; double r22676 = lambda1; double r22677 = lambda2; double r22678 = r22676 - r22677; double r22679 = phi1; double r22680 = phi2; double r22681 = r22679 + r22680; double r22682 = 2; double r22683 = r22681 / r22682; double r22684 = cos(r22683); double r22685 = r22678 * r22684; double r22686 = r22685 * r22685; double r22687 = r22679 - r22680; double r22688 = r22687 * r22687; double r22689 = r22686 + r22688; double r22690 = sqrt(r22689); double r22691 = r22675 * r22690; return r22691; } double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) { float r22692 = phi1; float r22693 = phi2; float r22694 = r22692 + r22693; float r22695 = 2; float r22696 = r22694 / r22695; float r22697 = cos(r22696); float r22698 = lambda1; float r22699 = lambda2; float r22700 = r22698 - r22699; float r22701 = r22697 * r22700; float r22702 = r22692 - r22693; float r22703 = hypot(r22701, r22702); float r22704 = R; float r22705 = r22703 * r22704; return r22705; } double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) { double r22706 = phi1; double r22707 = phi2; double r22708 = r22706 + r22707; double r22709 = 2; double r22710 = r22708 / r22709; double r22711 = cos(r22710); double r22712 = lambda1; double r22713 = lambda2; double r22714 = r22712 - r22713; double r22715 = r22711 * r22714; double r22716 = r22706 - r22707; double r22717 = hypot(r22715, r22716); double r22718 = R; double r22719 = r22717 * r22718; return r22719; } 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 r22720, r22721, r22722, r22723, r22724, r22725, r22726, r22727, r22728, r22729, r22730, r22731, r22732, r22733, r22734, r22735, r22736; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r22720); mpfr_init(r22721); mpfr_init(r22722); mpfr_init(r22723); mpfr_init(r22724); mpfr_init(r22725); mpfr_init(r22726); mpfr_init_set_str(r22727, "2", 10, MPFR_RNDN); mpfr_init(r22728); mpfr_init(r22729); mpfr_init(r22730); mpfr_init(r22731); mpfr_init(r22732); mpfr_init(r22733); mpfr_init(r22734); mpfr_init(r22735); mpfr_init(r22736); } double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r22720, R, MPFR_RNDN); mpfr_set_d(r22721, lambda1, MPFR_RNDN); mpfr_set_d(r22722, lambda2, MPFR_RNDN); mpfr_sub(r22723, r22721, r22722, MPFR_RNDN); mpfr_set_d(r22724, phi1, MPFR_RNDN); mpfr_set_d(r22725, phi2, MPFR_RNDN); mpfr_add(r22726, r22724, r22725, MPFR_RNDN); ; mpfr_div(r22728, r22726, r22727, MPFR_RNDN); mpfr_cos(r22729, r22728, MPFR_RNDN); mpfr_mul(r22730, r22723, r22729, MPFR_RNDN); mpfr_mul(r22731, r22730, r22730, MPFR_RNDN); mpfr_sub(r22732, r22724, r22725, MPFR_RNDN); mpfr_mul(r22733, r22732, r22732, MPFR_RNDN); mpfr_add(r22734, r22731, r22733, MPFR_RNDN); mpfr_sqrt(r22735, r22734, MPFR_RNDN); mpfr_mul(r22736, r22720, r22735, MPFR_RNDN); return mpfr_get_d(r22736, MPFR_RNDN); } static mpfr_t r22737, r22738, r22739, r22740, r22741, r22742, r22743, r22744, r22745, r22746, r22747, r22748, r22749, r22750; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r22737); mpfr_init(r22738); mpfr_init(r22739); mpfr_init_set_str(r22740, "2", 10, MPFR_RNDN); mpfr_init(r22741); mpfr_init(r22742); mpfr_init(r22743); mpfr_init(r22744); mpfr_init(r22745); mpfr_init(r22746); mpfr_init(r22747); mpfr_init(r22748); mpfr_init(r22749); mpfr_init(r22750); } double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r22737, phi1, MPFR_RNDN); mpfr_set_d(r22738, phi2, MPFR_RNDN); mpfr_add(r22739, r22737, r22738, MPFR_RNDN); ; mpfr_div(r22741, r22739, r22740, MPFR_RNDN); mpfr_cos(r22742, r22741, MPFR_RNDN); mpfr_set_d(r22743, lambda1, MPFR_RNDN); mpfr_set_d(r22744, lambda2, MPFR_RNDN); mpfr_sub(r22745, r22743, r22744, MPFR_RNDN); mpfr_mul(r22746, r22742, r22745, MPFR_RNDN); mpfr_sub(r22747, r22737, r22738, MPFR_RNDN); mpfr_hypot(r22748, r22746, r22747, MPFR_RNDN); mpfr_set_d(r22749, R, MPFR_RNDN); mpfr_mul(r22750, r22748, r22749, MPFR_RNDN); return mpfr_get_d(r22750, MPFR_RNDN); } static mpfr_t r22751, r22752, r22753, r22754, r22755, r22756, r22757, r22758, r22759, r22760, r22761, r22762, r22763, r22764; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r22751); mpfr_init(r22752); mpfr_init(r22753); mpfr_init_set_str(r22754, "2", 10, MPFR_RNDN); mpfr_init(r22755); mpfr_init(r22756); mpfr_init(r22757); mpfr_init(r22758); mpfr_init(r22759); mpfr_init(r22760); mpfr_init(r22761); mpfr_init(r22762); mpfr_init(r22763); mpfr_init(r22764); } double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r22751, phi1, MPFR_RNDN); mpfr_set_d(r22752, phi2, MPFR_RNDN); mpfr_add(r22753, r22751, r22752, MPFR_RNDN); ; mpfr_div(r22755, r22753, r22754, MPFR_RNDN); mpfr_cos(r22756, r22755, MPFR_RNDN); mpfr_set_d(r22757, lambda1, MPFR_RNDN); mpfr_set_d(r22758, lambda2, MPFR_RNDN); mpfr_sub(r22759, r22757, r22758, MPFR_RNDN); mpfr_mul(r22760, r22756, r22759, MPFR_RNDN); mpfr_sub(r22761, r22751, r22752, MPFR_RNDN); mpfr_hypot(r22762, r22760, r22761, MPFR_RNDN); mpfr_set_d(r22763, R, MPFR_RNDN); mpfr_mul(r22764, r22762, r22763, MPFR_RNDN); return mpfr_get_d(r22764, MPFR_RNDN); }