#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 r23074 = R; float r23075 = lambda1; float r23076 = lambda2; float r23077 = r23075 - r23076; float r23078 = phi1; float r23079 = phi2; float r23080 = r23078 + r23079; float r23081 = 2; float r23082 = r23080 / r23081; float r23083 = cos(r23082); float r23084 = r23077 * r23083; float r23085 = r23084 * r23084; float r23086 = r23078 - r23079; float r23087 = r23086 * r23086; float r23088 = r23085 + r23087; float r23089 = sqrt(r23088); float r23090 = r23074 * r23089; return r23090; } double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) { double r23091 = R; double r23092 = lambda1; double r23093 = lambda2; double r23094 = r23092 - r23093; double r23095 = phi1; double r23096 = phi2; double r23097 = r23095 + r23096; double r23098 = 2; double r23099 = r23097 / r23098; double r23100 = cos(r23099); double r23101 = r23094 * r23100; double r23102 = r23101 * r23101; double r23103 = r23095 - r23096; double r23104 = r23103 * r23103; double r23105 = r23102 + r23104; double r23106 = sqrt(r23105); double r23107 = r23091 * r23106; return r23107; } double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) { float r23108 = phi1; float r23109 = phi2; float r23110 = r23108 + r23109; float r23111 = 2; float r23112 = r23110 / r23111; float r23113 = cos(r23112); float r23114 = lambda1; float r23115 = lambda2; float r23116 = r23114 - r23115; float r23117 = r23113 * r23116; float r23118 = r23108 - r23109; float r23119 = hypot(r23117, r23118); float r23120 = R; float r23121 = r23119 * r23120; return r23121; } double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) { double r23122 = phi1; double r23123 = phi2; double r23124 = r23122 + r23123; double r23125 = 2; double r23126 = r23124 / r23125; double r23127 = cos(r23126); double r23128 = lambda1; double r23129 = lambda2; double r23130 = r23128 - r23129; double r23131 = r23127 * r23130; double r23132 = r23122 - r23123; double r23133 = hypot(r23131, r23132); double r23134 = R; double r23135 = r23133 * r23134; return r23135; } 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 r23136, r23137, r23138, r23139, r23140, r23141, r23142, r23143, r23144, r23145, r23146, r23147, r23148, r23149, r23150, r23151, r23152; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r23136); mpfr_init(r23137); mpfr_init(r23138); mpfr_init(r23139); mpfr_init(r23140); mpfr_init(r23141); mpfr_init(r23142); mpfr_init_set_str(r23143, "2", 10, MPFR_RNDN); mpfr_init(r23144); mpfr_init(r23145); mpfr_init(r23146); mpfr_init(r23147); mpfr_init(r23148); mpfr_init(r23149); mpfr_init(r23150); mpfr_init(r23151); mpfr_init(r23152); } double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r23136, R, MPFR_RNDN); mpfr_set_d(r23137, lambda1, MPFR_RNDN); mpfr_set_d(r23138, lambda2, MPFR_RNDN); mpfr_sub(r23139, r23137, r23138, MPFR_RNDN); mpfr_set_d(r23140, phi1, MPFR_RNDN); mpfr_set_d(r23141, phi2, MPFR_RNDN); mpfr_add(r23142, r23140, r23141, MPFR_RNDN); ; mpfr_div(r23144, r23142, r23143, MPFR_RNDN); mpfr_cos(r23145, r23144, MPFR_RNDN); mpfr_mul(r23146, r23139, r23145, MPFR_RNDN); mpfr_mul(r23147, r23146, r23146, MPFR_RNDN); mpfr_sub(r23148, r23140, r23141, MPFR_RNDN); mpfr_mul(r23149, r23148, r23148, MPFR_RNDN); mpfr_add(r23150, r23147, r23149, MPFR_RNDN); mpfr_sqrt(r23151, r23150, MPFR_RNDN); mpfr_mul(r23152, r23136, r23151, MPFR_RNDN); return mpfr_get_d(r23152, MPFR_RNDN); } static mpfr_t r23153, r23154, r23155, r23156, r23157, r23158, r23159, r23160, r23161, r23162, r23163, r23164, r23165, r23166; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r23153); mpfr_init(r23154); mpfr_init(r23155); mpfr_init_set_str(r23156, "2", 10, MPFR_RNDN); mpfr_init(r23157); mpfr_init(r23158); mpfr_init(r23159); mpfr_init(r23160); mpfr_init(r23161); mpfr_init(r23162); mpfr_init(r23163); mpfr_init(r23164); mpfr_init(r23165); mpfr_init(r23166); } double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r23153, phi1, MPFR_RNDN); mpfr_set_d(r23154, phi2, MPFR_RNDN); mpfr_add(r23155, r23153, r23154, MPFR_RNDN); ; mpfr_div(r23157, r23155, r23156, MPFR_RNDN); mpfr_cos(r23158, r23157, MPFR_RNDN); mpfr_set_d(r23159, lambda1, MPFR_RNDN); mpfr_set_d(r23160, lambda2, MPFR_RNDN); mpfr_sub(r23161, r23159, r23160, MPFR_RNDN); mpfr_mul(r23162, r23158, r23161, MPFR_RNDN); mpfr_sub(r23163, r23153, r23154, MPFR_RNDN); mpfr_hypot(r23164, r23162, r23163, MPFR_RNDN); mpfr_set_d(r23165, R, MPFR_RNDN); mpfr_mul(r23166, r23164, r23165, MPFR_RNDN); return mpfr_get_d(r23166, MPFR_RNDN); } static mpfr_t r23167, r23168, r23169, r23170, r23171, r23172, r23173, r23174, r23175, r23176, r23177, r23178, r23179, r23180; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r23167); mpfr_init(r23168); mpfr_init(r23169); mpfr_init_set_str(r23170, "2", 10, MPFR_RNDN); mpfr_init(r23171); mpfr_init(r23172); mpfr_init(r23173); mpfr_init(r23174); mpfr_init(r23175); mpfr_init(r23176); mpfr_init(r23177); mpfr_init(r23178); mpfr_init(r23179); mpfr_init(r23180); } double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r23167, phi1, MPFR_RNDN); mpfr_set_d(r23168, phi2, MPFR_RNDN); mpfr_add(r23169, r23167, r23168, MPFR_RNDN); ; mpfr_div(r23171, r23169, r23170, MPFR_RNDN); mpfr_cos(r23172, r23171, MPFR_RNDN); mpfr_set_d(r23173, lambda1, MPFR_RNDN); mpfr_set_d(r23174, lambda2, MPFR_RNDN); mpfr_sub(r23175, r23173, r23174, MPFR_RNDN); mpfr_mul(r23176, r23172, r23175, MPFR_RNDN); mpfr_sub(r23177, r23167, r23168, MPFR_RNDN); mpfr_hypot(r23178, r23176, r23177, MPFR_RNDN); mpfr_set_d(r23179, R, MPFR_RNDN); mpfr_mul(r23180, r23178, r23179, MPFR_RNDN); return mpfr_get_d(r23180, MPFR_RNDN); }