#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 r21935 = R; float r21936 = lambda1; float r21937 = lambda2; float r21938 = r21936 - r21937; float r21939 = phi1; float r21940 = phi2; float r21941 = r21939 + r21940; float r21942 = 2; float r21943 = r21941 / r21942; float r21944 = cos(r21943); float r21945 = r21938 * r21944; float r21946 = r21945 * r21945; float r21947 = r21939 - r21940; float r21948 = r21947 * r21947; float r21949 = r21946 + r21948; float r21950 = sqrt(r21949); float r21951 = r21935 * r21950; return r21951; } double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) { double r21952 = R; double r21953 = lambda1; double r21954 = lambda2; double r21955 = r21953 - r21954; double r21956 = phi1; double r21957 = phi2; double r21958 = r21956 + r21957; double r21959 = 2; double r21960 = r21958 / r21959; double r21961 = cos(r21960); double r21962 = r21955 * r21961; double r21963 = r21962 * r21962; double r21964 = r21956 - r21957; double r21965 = r21964 * r21964; double r21966 = r21963 + r21965; double r21967 = sqrt(r21966); double r21968 = r21952 * r21967; return r21968; } double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) { float r21969 = phi1; float r21970 = phi2; float r21971 = r21969 + r21970; float r21972 = 2; float r21973 = r21971 / r21972; float r21974 = cos(r21973); float r21975 = exp(r21974); float r21976 = log1p(r21975); float r21977 = expm1(r21976); float r21978 = log(r21977); float r21979 = lambda1; float r21980 = lambda2; float r21981 = r21979 - r21980; float r21982 = r21978 * r21981; float r21983 = r21969 - r21970; float r21984 = hypot(r21982, r21983); float r21985 = R; float r21986 = r21984 * r21985; return r21986; } double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) { double r21987 = phi1; double r21988 = phi2; double r21989 = r21987 + r21988; double r21990 = 2; double r21991 = r21989 / r21990; double r21992 = cos(r21991); double r21993 = exp(r21992); double r21994 = log1p(r21993); double r21995 = expm1(r21994); double r21996 = log(r21995); double r21997 = lambda1; double r21998 = lambda2; double r21999 = r21997 - r21998; double r22000 = r21996 * r21999; double r22001 = r21987 - r21988; double r22002 = hypot(r22000, r22001); double r22003 = R; double r22004 = r22002 * r22003; return r22004; } 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 r22005, r22006, r22007, r22008, r22009, r22010, r22011, r22012, r22013, r22014, r22015, r22016, r22017, r22018, r22019, r22020, r22021; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r22005); mpfr_init(r22006); mpfr_init(r22007); mpfr_init(r22008); mpfr_init(r22009); mpfr_init(r22010); mpfr_init(r22011); mpfr_init_set_str(r22012, "2", 10, MPFR_RNDN); mpfr_init(r22013); mpfr_init(r22014); mpfr_init(r22015); mpfr_init(r22016); mpfr_init(r22017); mpfr_init(r22018); mpfr_init(r22019); mpfr_init(r22020); mpfr_init(r22021); } double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r22005, R, MPFR_RNDN); mpfr_set_d(r22006, lambda1, MPFR_RNDN); mpfr_set_d(r22007, lambda2, MPFR_RNDN); mpfr_sub(r22008, r22006, r22007, MPFR_RNDN); mpfr_set_d(r22009, phi1, MPFR_RNDN); mpfr_set_d(r22010, phi2, MPFR_RNDN); mpfr_add(r22011, r22009, r22010, MPFR_RNDN); ; mpfr_div(r22013, r22011, r22012, MPFR_RNDN); mpfr_cos(r22014, r22013, MPFR_RNDN); mpfr_mul(r22015, r22008, r22014, MPFR_RNDN); mpfr_mul(r22016, r22015, r22015, MPFR_RNDN); mpfr_sub(r22017, r22009, r22010, MPFR_RNDN); mpfr_mul(r22018, r22017, r22017, MPFR_RNDN); mpfr_add(r22019, r22016, r22018, MPFR_RNDN); mpfr_sqrt(r22020, r22019, MPFR_RNDN); mpfr_mul(r22021, r22005, r22020, MPFR_RNDN); return mpfr_get_d(r22021, MPFR_RNDN); } static mpfr_t r22022, r22023, r22024, r22025, r22026, r22027, r22028, r22029, r22030, r22031, r22032, r22033, r22034, r22035, r22036, r22037, r22038, r22039; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r22022); mpfr_init(r22023); mpfr_init(r22024); mpfr_init_set_str(r22025, "2", 10, MPFR_RNDN); mpfr_init(r22026); mpfr_init(r22027); mpfr_init(r22028); mpfr_init(r22029); mpfr_init(r22030); mpfr_init(r22031); mpfr_init(r22032); mpfr_init(r22033); mpfr_init(r22034); mpfr_init(r22035); mpfr_init(r22036); mpfr_init(r22037); mpfr_init(r22038); mpfr_init(r22039); } double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r22022, phi1, MPFR_RNDN); mpfr_set_d(r22023, phi2, MPFR_RNDN); mpfr_add(r22024, r22022, r22023, MPFR_RNDN); ; mpfr_div(r22026, r22024, r22025, MPFR_RNDN); mpfr_cos(r22027, r22026, MPFR_RNDN); mpfr_exp(r22028, r22027, MPFR_RNDN); mpfr_log1p(r22029, r22028, MPFR_RNDN); mpfr_expm1(r22030, r22029, MPFR_RNDN); mpfr_log(r22031, r22030, MPFR_RNDN); mpfr_set_d(r22032, lambda1, MPFR_RNDN); mpfr_set_d(r22033, lambda2, MPFR_RNDN); mpfr_sub(r22034, r22032, r22033, MPFR_RNDN); mpfr_mul(r22035, r22031, r22034, MPFR_RNDN); mpfr_sub(r22036, r22022, r22023, MPFR_RNDN); mpfr_hypot(r22037, r22035, r22036, MPFR_RNDN); mpfr_set_d(r22038, R, MPFR_RNDN); mpfr_mul(r22039, r22037, r22038, MPFR_RNDN); return mpfr_get_d(r22039, MPFR_RNDN); } static mpfr_t r22040, r22041, r22042, r22043, r22044, r22045, r22046, r22047, r22048, r22049, r22050, r22051, r22052, r22053, r22054, r22055, r22056, r22057; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r22040); mpfr_init(r22041); mpfr_init(r22042); mpfr_init_set_str(r22043, "2", 10, MPFR_RNDN); mpfr_init(r22044); mpfr_init(r22045); mpfr_init(r22046); mpfr_init(r22047); mpfr_init(r22048); mpfr_init(r22049); mpfr_init(r22050); mpfr_init(r22051); mpfr_init(r22052); mpfr_init(r22053); mpfr_init(r22054); mpfr_init(r22055); mpfr_init(r22056); mpfr_init(r22057); } double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r22040, phi1, MPFR_RNDN); mpfr_set_d(r22041, phi2, MPFR_RNDN); mpfr_add(r22042, r22040, r22041, MPFR_RNDN); ; mpfr_div(r22044, r22042, r22043, MPFR_RNDN); mpfr_cos(r22045, r22044, MPFR_RNDN); mpfr_exp(r22046, r22045, MPFR_RNDN); mpfr_log1p(r22047, r22046, MPFR_RNDN); mpfr_expm1(r22048, r22047, MPFR_RNDN); mpfr_log(r22049, r22048, MPFR_RNDN); mpfr_set_d(r22050, lambda1, MPFR_RNDN); mpfr_set_d(r22051, lambda2, MPFR_RNDN); mpfr_sub(r22052, r22050, r22051, MPFR_RNDN); mpfr_mul(r22053, r22049, r22052, MPFR_RNDN); mpfr_sub(r22054, r22040, r22041, MPFR_RNDN); mpfr_hypot(r22055, r22053, r22054, MPFR_RNDN); mpfr_set_d(r22056, R, MPFR_RNDN); mpfr_mul(r22057, r22055, r22056, MPFR_RNDN); return mpfr_get_d(r22057, MPFR_RNDN); }