#include #include #include #include #include char *name = "Midpoint on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r28266 = lambda1; float r28267 = phi2; float r28268 = cos(r28267); float r28269 = lambda2; float r28270 = r28266 - r28269; float r28271 = sin(r28270); float r28272 = r28268 * r28271; float r28273 = phi1; float r28274 = cos(r28273); float r28275 = cos(r28270); float r28276 = r28268 * r28275; float r28277 = r28274 + r28276; float r28278 = atan2(r28272, r28277); float r28279 = r28266 + r28278; return r28279; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r28280 = lambda1; double r28281 = phi2; double r28282 = cos(r28281); double r28283 = lambda2; double r28284 = r28280 - r28283; double r28285 = sin(r28284); double r28286 = r28282 * r28285; double r28287 = phi1; double r28288 = cos(r28287); double r28289 = cos(r28284); double r28290 = r28282 * r28289; double r28291 = r28288 + r28290; double r28292 = atan2(r28286, r28291); double r28293 = r28280 + r28292; return r28293; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r28294 = lambda1; float r28295 = phi2; float r28296 = cos(r28295); float r28297 = lambda2; float r28298 = r28294 - r28297; float r28299 = sin(r28298); float r28300 = r28296 * r28299; float r28301 = phi1; float r28302 = cos(r28301); float r28303 = 3; float r28304 = pow(r28302, r28303); float r28305 = cos(r28298); float r28306 = r28296 * r28305; float r28307 = pow(r28306, r28303); float r28308 = r28304 + r28307; float r28309 = r28302 * r28302; float r28310 = r28306 * r28306; float r28311 = r28302 * r28306; float r28312 = r28310 - r28311; float r28313 = r28309 + r28312; float r28314 = r28308 / r28313; float r28315 = atan2(r28300, r28314); float r28316 = r28294 + r28315; return r28316; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r28317 = lambda1; double r28318 = phi2; double r28319 = cos(r28318); double r28320 = lambda2; double r28321 = r28317 - r28320; double r28322 = sin(r28321); double r28323 = r28319 * r28322; double r28324 = phi1; double r28325 = cos(r28324); double r28326 = 3; double r28327 = pow(r28325, r28326); double r28328 = cos(r28321); double r28329 = r28319 * r28328; double r28330 = pow(r28329, r28326); double r28331 = r28327 + r28330; double r28332 = r28325 * r28325; double r28333 = r28329 * r28329; double r28334 = r28325 * r28329; double r28335 = r28333 - r28334; double r28336 = r28332 + r28335; double r28337 = r28331 / r28336; double r28338 = atan2(r28323, r28337); double r28339 = r28317 + r28338; return r28339; } 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 r28340, r28341, r28342, r28343, r28344, r28345, r28346, r28347, r28348, r28349, r28350, r28351, r28352, r28353; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r28340); mpfr_init(r28341); mpfr_init(r28342); mpfr_init(r28343); mpfr_init(r28344); mpfr_init(r28345); mpfr_init(r28346); mpfr_init(r28347); mpfr_init(r28348); mpfr_init(r28349); mpfr_init(r28350); mpfr_init(r28351); mpfr_init(r28352); mpfr_init(r28353); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r28340, lambda1, MPFR_RNDN); mpfr_set_d(r28341, phi2, MPFR_RNDN); mpfr_cos(r28342, r28341, MPFR_RNDN); mpfr_set_d(r28343, lambda2, MPFR_RNDN); mpfr_sub(r28344, r28340, r28343, MPFR_RNDN); mpfr_sin(r28345, r28344, MPFR_RNDN); mpfr_mul(r28346, r28342, r28345, MPFR_RNDN); mpfr_set_d(r28347, phi1, MPFR_RNDN); mpfr_cos(r28348, r28347, MPFR_RNDN); mpfr_cos(r28349, r28344, MPFR_RNDN); mpfr_mul(r28350, r28342, r28349, MPFR_RNDN); mpfr_add(r28351, r28348, r28350, MPFR_RNDN); mpfr_atan2(r28352, r28346, r28351, MPFR_RNDN); mpfr_add(r28353, r28340, r28352, MPFR_RNDN); return mpfr_get_d(r28353, MPFR_RNDN); } static mpfr_t r28354, r28355, r28356, r28357, r28358, r28359, r28360, r28361, r28362, r28363, r28364, r28365, r28366, r28367, r28368, r28369, r28370, r28371, r28372, r28373, r28374, r28375, r28376; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r28354); mpfr_init(r28355); mpfr_init(r28356); mpfr_init(r28357); mpfr_init(r28358); mpfr_init(r28359); mpfr_init(r28360); mpfr_init(r28361); mpfr_init(r28362); mpfr_init_set_str(r28363, "3", 10, MPFR_RNDN); mpfr_init(r28364); mpfr_init(r28365); mpfr_init(r28366); mpfr_init(r28367); mpfr_init(r28368); mpfr_init(r28369); mpfr_init(r28370); mpfr_init(r28371); mpfr_init(r28372); mpfr_init(r28373); mpfr_init(r28374); mpfr_init(r28375); mpfr_init(r28376); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r28354, lambda1, MPFR_RNDN); mpfr_set_d(r28355, phi2, MPFR_RNDN); mpfr_cos(r28356, r28355, MPFR_RNDN); mpfr_set_d(r28357, lambda2, MPFR_RNDN); mpfr_sub(r28358, r28354, r28357, MPFR_RNDN); mpfr_sin(r28359, r28358, MPFR_RNDN); mpfr_mul(r28360, r28356, r28359, MPFR_RNDN); mpfr_set_d(r28361, phi1, MPFR_RNDN); mpfr_cos(r28362, r28361, MPFR_RNDN); ; mpfr_pow(r28364, r28362, r28363, MPFR_RNDN); mpfr_cos(r28365, r28358, MPFR_RNDN); mpfr_mul(r28366, r28356, r28365, MPFR_RNDN); mpfr_pow(r28367, r28366, r28363, MPFR_RNDN); mpfr_add(r28368, r28364, r28367, MPFR_RNDN); mpfr_mul(r28369, r28362, r28362, MPFR_RNDN); mpfr_mul(r28370, r28366, r28366, MPFR_RNDN); mpfr_mul(r28371, r28362, r28366, MPFR_RNDN); mpfr_sub(r28372, r28370, r28371, MPFR_RNDN); mpfr_add(r28373, r28369, r28372, MPFR_RNDN); mpfr_div(r28374, r28368, r28373, MPFR_RNDN); mpfr_atan2(r28375, r28360, r28374, MPFR_RNDN); mpfr_add(r28376, r28354, r28375, MPFR_RNDN); return mpfr_get_d(r28376, MPFR_RNDN); } static mpfr_t r28377, r28378, r28379, r28380, r28381, r28382, r28383, r28384, r28385, r28386, r28387, r28388, r28389, r28390, r28391, r28392, r28393, r28394, r28395, r28396, r28397, r28398, r28399; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r28377); mpfr_init(r28378); mpfr_init(r28379); mpfr_init(r28380); mpfr_init(r28381); mpfr_init(r28382); mpfr_init(r28383); mpfr_init(r28384); mpfr_init(r28385); mpfr_init_set_str(r28386, "3", 10, MPFR_RNDN); mpfr_init(r28387); mpfr_init(r28388); mpfr_init(r28389); mpfr_init(r28390); mpfr_init(r28391); mpfr_init(r28392); mpfr_init(r28393); mpfr_init(r28394); mpfr_init(r28395); mpfr_init(r28396); mpfr_init(r28397); mpfr_init(r28398); mpfr_init(r28399); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r28377, lambda1, MPFR_RNDN); mpfr_set_d(r28378, phi2, MPFR_RNDN); mpfr_cos(r28379, r28378, MPFR_RNDN); mpfr_set_d(r28380, lambda2, MPFR_RNDN); mpfr_sub(r28381, r28377, r28380, MPFR_RNDN); mpfr_sin(r28382, r28381, MPFR_RNDN); mpfr_mul(r28383, r28379, r28382, MPFR_RNDN); mpfr_set_d(r28384, phi1, MPFR_RNDN); mpfr_cos(r28385, r28384, MPFR_RNDN); ; mpfr_pow(r28387, r28385, r28386, MPFR_RNDN); mpfr_cos(r28388, r28381, MPFR_RNDN); mpfr_mul(r28389, r28379, r28388, MPFR_RNDN); mpfr_pow(r28390, r28389, r28386, MPFR_RNDN); mpfr_add(r28391, r28387, r28390, MPFR_RNDN); mpfr_mul(r28392, r28385, r28385, MPFR_RNDN); mpfr_mul(r28393, r28389, r28389, MPFR_RNDN); mpfr_mul(r28394, r28385, r28389, MPFR_RNDN); mpfr_sub(r28395, r28393, r28394, MPFR_RNDN); mpfr_add(r28396, r28392, r28395, MPFR_RNDN); mpfr_div(r28397, r28391, r28396, MPFR_RNDN); mpfr_atan2(r28398, r28383, r28397, MPFR_RNDN); mpfr_add(r28399, r28377, r28398, MPFR_RNDN); return mpfr_get_d(r28399, MPFR_RNDN); }