#include #include #include #include #include char *name = "Midpoint on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r8318 = lambda1; float r8319 = phi2; float r8320 = cos(r8319); float r8321 = lambda2; float r8322 = r8318 - r8321; float r8323 = sin(r8322); float r8324 = r8320 * r8323; float r8325 = phi1; float r8326 = cos(r8325); float r8327 = cos(r8322); float r8328 = r8320 * r8327; float r8329 = r8326 + r8328; float r8330 = atan2(r8324, r8329); float r8331 = r8318 + r8330; return r8331; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r8332 = lambda1; double r8333 = phi2; double r8334 = cos(r8333); double r8335 = lambda2; double r8336 = r8332 - r8335; double r8337 = sin(r8336); double r8338 = r8334 * r8337; double r8339 = phi1; double r8340 = cos(r8339); double r8341 = cos(r8336); double r8342 = r8334 * r8341; double r8343 = r8340 + r8342; double r8344 = atan2(r8338, r8343); double r8345 = r8332 + r8344; return r8345; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r8346 = phi2; float r8347 = cos(r8346); float r8348 = lambda1; float r8349 = lambda2; float r8350 = r8348 - r8349; float r8351 = sin(r8350); float r8352 = log1p(r8351); float r8353 = expm1(r8352); float r8354 = r8347 * r8353; float r8355 = cos(r8350); float r8356 = phi1; float r8357 = cos(r8356); float r8358 = fma(r8347, r8355, r8357); float r8359 = atan2(r8354, r8358); float r8360 = r8359 + r8348; return r8360; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r8361 = phi2; double r8362 = cos(r8361); double r8363 = lambda1; double r8364 = lambda2; double r8365 = r8363 - r8364; double r8366 = sin(r8365); double r8367 = log1p(r8366); double r8368 = expm1(r8367); double r8369 = r8362 * r8368; double r8370 = cos(r8365); double r8371 = phi1; double r8372 = cos(r8371); double r8373 = fma(r8362, r8370, r8372); double r8374 = atan2(r8369, r8373); double r8375 = r8374 + r8363; return r8375; } 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 r8376, r8377, r8378, r8379, r8380, r8381, r8382, r8383, r8384, r8385, r8386, r8387, r8388, r8389; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r8376); mpfr_init(r8377); mpfr_init(r8378); mpfr_init(r8379); mpfr_init(r8380); mpfr_init(r8381); mpfr_init(r8382); mpfr_init(r8383); mpfr_init(r8384); mpfr_init(r8385); mpfr_init(r8386); mpfr_init(r8387); mpfr_init(r8388); mpfr_init(r8389); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r8376, lambda1, MPFR_RNDN); mpfr_set_d(r8377, phi2, MPFR_RNDN); mpfr_cos(r8378, r8377, MPFR_RNDN); mpfr_set_d(r8379, lambda2, MPFR_RNDN); mpfr_sub(r8380, r8376, r8379, MPFR_RNDN); mpfr_sin(r8381, r8380, MPFR_RNDN); mpfr_mul(r8382, r8378, r8381, MPFR_RNDN); mpfr_set_d(r8383, phi1, MPFR_RNDN); mpfr_cos(r8384, r8383, MPFR_RNDN); mpfr_cos(r8385, r8380, MPFR_RNDN); mpfr_mul(r8386, r8378, r8385, MPFR_RNDN); mpfr_add(r8387, r8384, r8386, MPFR_RNDN); mpfr_atan2(r8388, r8382, r8387, MPFR_RNDN); mpfr_add(r8389, r8376, r8388, MPFR_RNDN); return mpfr_get_d(r8389, MPFR_RNDN); } static mpfr_t r8390, r8391, r8392, r8393, r8394, r8395, r8396, r8397, r8398, r8399, r8400, r8401, r8402, r8403, r8404; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r8390); mpfr_init(r8391); mpfr_init(r8392); mpfr_init(r8393); mpfr_init(r8394); mpfr_init(r8395); mpfr_init(r8396); mpfr_init(r8397); mpfr_init(r8398); mpfr_init(r8399); mpfr_init(r8400); mpfr_init(r8401); mpfr_init(r8402); mpfr_init(r8403); mpfr_init(r8404); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r8390, phi2, MPFR_RNDN); mpfr_cos(r8391, r8390, MPFR_RNDN); mpfr_set_d(r8392, lambda1, MPFR_RNDN); mpfr_set_d(r8393, lambda2, MPFR_RNDN); mpfr_sub(r8394, r8392, r8393, MPFR_RNDN); mpfr_sin(r8395, r8394, MPFR_RNDN); mpfr_log1p(r8396, r8395, MPFR_RNDN); mpfr_expm1(r8397, r8396, MPFR_RNDN); mpfr_mul(r8398, r8391, r8397, MPFR_RNDN); mpfr_cos(r8399, r8394, MPFR_RNDN); mpfr_set_d(r8400, phi1, MPFR_RNDN); mpfr_cos(r8401, r8400, MPFR_RNDN); mpfr_fma(r8402, r8391, r8399, r8401, MPFR_RNDN); mpfr_atan2(r8403, r8398, r8402, MPFR_RNDN); mpfr_add(r8404, r8403, r8392, MPFR_RNDN); return mpfr_get_d(r8404, MPFR_RNDN); } static mpfr_t r8405, r8406, r8407, r8408, r8409, r8410, r8411, r8412, r8413, r8414, r8415, r8416, r8417, r8418, r8419; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r8405); mpfr_init(r8406); mpfr_init(r8407); mpfr_init(r8408); mpfr_init(r8409); mpfr_init(r8410); mpfr_init(r8411); mpfr_init(r8412); mpfr_init(r8413); mpfr_init(r8414); mpfr_init(r8415); mpfr_init(r8416); mpfr_init(r8417); mpfr_init(r8418); mpfr_init(r8419); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r8405, phi2, MPFR_RNDN); mpfr_cos(r8406, r8405, MPFR_RNDN); mpfr_set_d(r8407, lambda1, MPFR_RNDN); mpfr_set_d(r8408, lambda2, MPFR_RNDN); mpfr_sub(r8409, r8407, r8408, MPFR_RNDN); mpfr_sin(r8410, r8409, MPFR_RNDN); mpfr_log1p(r8411, r8410, MPFR_RNDN); mpfr_expm1(r8412, r8411, MPFR_RNDN); mpfr_mul(r8413, r8406, r8412, MPFR_RNDN); mpfr_cos(r8414, r8409, MPFR_RNDN); mpfr_set_d(r8415, phi1, MPFR_RNDN); mpfr_cos(r8416, r8415, MPFR_RNDN); mpfr_fma(r8417, r8406, r8414, r8416, MPFR_RNDN); mpfr_atan2(r8418, r8413, r8417, MPFR_RNDN); mpfr_add(r8419, r8418, r8407, MPFR_RNDN); return mpfr_get_d(r8419, MPFR_RNDN); }