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