#include #include #include #include #include char *name = "Midpoint on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r22272 = lambda1; float r22273 = phi2; float r22274 = cos(r22273); float r22275 = lambda2; float r22276 = r22272 - r22275; float r22277 = sin(r22276); float r22278 = r22274 * r22277; float r22279 = phi1; float r22280 = cos(r22279); float r22281 = cos(r22276); float r22282 = r22274 * r22281; float r22283 = r22280 + r22282; float r22284 = atan2(r22278, r22283); float r22285 = r22272 + r22284; return r22285; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r22286 = lambda1; double r22287 = phi2; double r22288 = cos(r22287); double r22289 = lambda2; double r22290 = r22286 - r22289; double r22291 = sin(r22290); double r22292 = r22288 * r22291; double r22293 = phi1; double r22294 = cos(r22293); double r22295 = cos(r22290); double r22296 = r22288 * r22295; double r22297 = r22294 + r22296; double r22298 = atan2(r22292, r22297); double r22299 = r22286 + r22298; return r22299; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r22300 = lambda1; float r22301 = phi2; float r22302 = cos(r22301); float r22303 = lambda2; float r22304 = r22300 - r22303; float r22305 = sin(r22304); float r22306 = r22302 * r22305; float r22307 = phi1; float r22308 = cos(r22307); float r22309 = 3; float r22310 = pow(r22308, r22309); float r22311 = cos(r22304); float r22312 = r22302 * r22311; float r22313 = pow(r22312, r22309); float r22314 = r22310 + r22313; float r22315 = r22308 * r22308; float r22316 = r22312 * r22312; float r22317 = r22308 * r22312; float r22318 = r22316 - r22317; float r22319 = r22315 + r22318; float r22320 = r22314 / r22319; float r22321 = atan2(r22306, r22320); float r22322 = r22300 + r22321; return r22322; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r22323 = lambda1; double r22324 = phi2; double r22325 = cos(r22324); double r22326 = lambda2; double r22327 = r22323 - r22326; double r22328 = sin(r22327); double r22329 = r22325 * r22328; double r22330 = phi1; double r22331 = cos(r22330); double r22332 = 3; double r22333 = pow(r22331, r22332); double r22334 = cos(r22327); double r22335 = r22325 * r22334; double r22336 = pow(r22335, r22332); double r22337 = r22333 + r22336; double r22338 = r22331 * r22331; double r22339 = r22335 * r22335; double r22340 = r22331 * r22335; double r22341 = r22339 - r22340; double r22342 = r22338 + r22341; double r22343 = r22337 / r22342; double r22344 = atan2(r22329, r22343); double r22345 = r22323 + r22344; return r22345; } 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 r22346, r22347, r22348, r22349, r22350, r22351, r22352, r22353, r22354, r22355, r22356, r22357, r22358, r22359; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r22346); mpfr_init(r22347); mpfr_init(r22348); mpfr_init(r22349); mpfr_init(r22350); mpfr_init(r22351); mpfr_init(r22352); mpfr_init(r22353); mpfr_init(r22354); mpfr_init(r22355); mpfr_init(r22356); mpfr_init(r22357); mpfr_init(r22358); mpfr_init(r22359); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r22346, lambda1, MPFR_RNDN); mpfr_set_d(r22347, phi2, MPFR_RNDN); mpfr_cos(r22348, r22347, MPFR_RNDN); mpfr_set_d(r22349, lambda2, MPFR_RNDN); mpfr_sub(r22350, r22346, r22349, MPFR_RNDN); mpfr_sin(r22351, r22350, MPFR_RNDN); mpfr_mul(r22352, r22348, r22351, MPFR_RNDN); mpfr_set_d(r22353, phi1, MPFR_RNDN); mpfr_cos(r22354, r22353, MPFR_RNDN); mpfr_cos(r22355, r22350, MPFR_RNDN); mpfr_mul(r22356, r22348, r22355, MPFR_RNDN); mpfr_add(r22357, r22354, r22356, MPFR_RNDN); mpfr_atan2(r22358, r22352, r22357, MPFR_RNDN); mpfr_add(r22359, r22346, r22358, MPFR_RNDN); return mpfr_get_d(r22359, MPFR_RNDN); } static mpfr_t r22360, r22361, r22362, r22363, r22364, r22365, r22366, r22367, r22368, r22369, r22370, r22371, r22372, r22373, r22374, r22375, r22376, r22377, r22378, r22379, r22380, r22381, r22382; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r22360); mpfr_init(r22361); mpfr_init(r22362); mpfr_init(r22363); mpfr_init(r22364); mpfr_init(r22365); mpfr_init(r22366); mpfr_init(r22367); mpfr_init(r22368); mpfr_init_set_str(r22369, "3", 10, MPFR_RNDN); mpfr_init(r22370); mpfr_init(r22371); mpfr_init(r22372); mpfr_init(r22373); mpfr_init(r22374); mpfr_init(r22375); mpfr_init(r22376); mpfr_init(r22377); mpfr_init(r22378); mpfr_init(r22379); mpfr_init(r22380); mpfr_init(r22381); mpfr_init(r22382); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r22360, lambda1, MPFR_RNDN); mpfr_set_d(r22361, phi2, MPFR_RNDN); mpfr_cos(r22362, r22361, MPFR_RNDN); mpfr_set_d(r22363, lambda2, MPFR_RNDN); mpfr_sub(r22364, r22360, r22363, MPFR_RNDN); mpfr_sin(r22365, r22364, MPFR_RNDN); mpfr_mul(r22366, r22362, r22365, MPFR_RNDN); mpfr_set_d(r22367, phi1, MPFR_RNDN); mpfr_cos(r22368, r22367, MPFR_RNDN); ; mpfr_pow(r22370, r22368, r22369, MPFR_RNDN); mpfr_cos(r22371, r22364, MPFR_RNDN); mpfr_mul(r22372, r22362, r22371, MPFR_RNDN); mpfr_pow(r22373, r22372, r22369, MPFR_RNDN); mpfr_add(r22374, r22370, r22373, MPFR_RNDN); mpfr_mul(r22375, r22368, r22368, MPFR_RNDN); mpfr_mul(r22376, r22372, r22372, MPFR_RNDN); mpfr_mul(r22377, r22368, r22372, MPFR_RNDN); mpfr_sub(r22378, r22376, r22377, MPFR_RNDN); mpfr_add(r22379, r22375, r22378, MPFR_RNDN); mpfr_div(r22380, r22374, r22379, MPFR_RNDN); mpfr_atan2(r22381, r22366, r22380, MPFR_RNDN); mpfr_add(r22382, r22360, r22381, MPFR_RNDN); return mpfr_get_d(r22382, MPFR_RNDN); } static mpfr_t r22383, r22384, r22385, r22386, r22387, r22388, r22389, r22390, r22391, r22392, r22393, r22394, r22395, r22396, r22397, r22398, r22399, r22400, r22401, r22402, r22403, r22404, r22405; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r22383); mpfr_init(r22384); mpfr_init(r22385); mpfr_init(r22386); mpfr_init(r22387); mpfr_init(r22388); mpfr_init(r22389); mpfr_init(r22390); mpfr_init(r22391); mpfr_init_set_str(r22392, "3", 10, MPFR_RNDN); mpfr_init(r22393); mpfr_init(r22394); mpfr_init(r22395); mpfr_init(r22396); mpfr_init(r22397); mpfr_init(r22398); mpfr_init(r22399); mpfr_init(r22400); mpfr_init(r22401); mpfr_init(r22402); mpfr_init(r22403); mpfr_init(r22404); mpfr_init(r22405); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r22383, lambda1, MPFR_RNDN); mpfr_set_d(r22384, phi2, MPFR_RNDN); mpfr_cos(r22385, r22384, MPFR_RNDN); mpfr_set_d(r22386, lambda2, MPFR_RNDN); mpfr_sub(r22387, r22383, r22386, MPFR_RNDN); mpfr_sin(r22388, r22387, MPFR_RNDN); mpfr_mul(r22389, r22385, r22388, MPFR_RNDN); mpfr_set_d(r22390, phi1, MPFR_RNDN); mpfr_cos(r22391, r22390, MPFR_RNDN); ; mpfr_pow(r22393, r22391, r22392, MPFR_RNDN); mpfr_cos(r22394, r22387, MPFR_RNDN); mpfr_mul(r22395, r22385, r22394, MPFR_RNDN); mpfr_pow(r22396, r22395, r22392, MPFR_RNDN); mpfr_add(r22397, r22393, r22396, MPFR_RNDN); mpfr_mul(r22398, r22391, r22391, MPFR_RNDN); mpfr_mul(r22399, r22395, r22395, MPFR_RNDN); mpfr_mul(r22400, r22391, r22395, MPFR_RNDN); mpfr_sub(r22401, r22399, r22400, MPFR_RNDN); mpfr_add(r22402, r22398, r22401, MPFR_RNDN); mpfr_div(r22403, r22397, r22402, MPFR_RNDN); mpfr_atan2(r22404, r22389, r22403, MPFR_RNDN); mpfr_add(r22405, r22383, r22404, MPFR_RNDN); return mpfr_get_d(r22405, MPFR_RNDN); }