#include #include #include #include #include char *name = "Midpoint on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r22374 = lambda1; float r22375 = phi2; float r22376 = cos(r22375); float r22377 = lambda2; float r22378 = r22374 - r22377; float r22379 = sin(r22378); float r22380 = r22376 * r22379; float r22381 = phi1; float r22382 = cos(r22381); float r22383 = cos(r22378); float r22384 = r22376 * r22383; float r22385 = r22382 + r22384; float r22386 = atan2(r22380, r22385); float r22387 = r22374 + r22386; return r22387; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r22388 = lambda1; double r22389 = phi2; double r22390 = cos(r22389); double r22391 = lambda2; double r22392 = r22388 - r22391; double r22393 = sin(r22392); double r22394 = r22390 * r22393; double r22395 = phi1; double r22396 = cos(r22395); double r22397 = cos(r22392); double r22398 = r22390 * r22397; double r22399 = r22396 + r22398; double r22400 = atan2(r22394, r22399); double r22401 = r22388 + r22400; return r22401; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r22402 = lambda1; float r22403 = phi2; float r22404 = cos(r22403); float r22405 = sin(r22402); float r22406 = lambda2; float r22407 = cos(r22406); float r22408 = r22405 * r22407; float r22409 = cos(r22402); float r22410 = -r22406; float r22411 = sin(r22410); float r22412 = r22409 * r22411; float r22413 = r22408 + r22412; float r22414 = r22404 * r22413; float r22415 = phi1; float r22416 = cos(r22415); float r22417 = exp(r22416); float r22418 = exp(r22404); float r22419 = sin(r22406); float r22420 = r22419 * r22405; float r22421 = r22407 * r22409; float r22422 = r22420 + r22421; float r22423 = pow(r22418, r22422); float r22424 = r22417 * r22423; float r22425 = log(r22424); float r22426 = atan2(r22414, r22425); float r22427 = r22402 + r22426; return r22427; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r22428 = lambda1; double r22429 = phi2; double r22430 = cos(r22429); double r22431 = sin(r22428); double r22432 = lambda2; double r22433 = cos(r22432); double r22434 = r22431 * r22433; double r22435 = cos(r22428); double r22436 = -r22432; double r22437 = sin(r22436); double r22438 = r22435 * r22437; double r22439 = r22434 + r22438; double r22440 = r22430 * r22439; double r22441 = phi1; double r22442 = cos(r22441); double r22443 = exp(r22442); double r22444 = exp(r22430); double r22445 = sin(r22432); double r22446 = r22445 * r22431; double r22447 = r22433 * r22435; double r22448 = r22446 + r22447; double r22449 = pow(r22444, r22448); double r22450 = r22443 * r22449; double r22451 = log(r22450); double r22452 = atan2(r22440, r22451); double r22453 = r22428 + r22452; return r22453; } 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 r22454, r22455, r22456, r22457, r22458, r22459, r22460, r22461, r22462, r22463, r22464, r22465, r22466, r22467; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r22454); mpfr_init(r22455); mpfr_init(r22456); mpfr_init(r22457); mpfr_init(r22458); mpfr_init(r22459); mpfr_init(r22460); mpfr_init(r22461); mpfr_init(r22462); mpfr_init(r22463); mpfr_init(r22464); mpfr_init(r22465); mpfr_init(r22466); mpfr_init(r22467); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r22454, lambda1, MPFR_RNDN); mpfr_set_d(r22455, phi2, MPFR_RNDN); mpfr_cos(r22456, r22455, MPFR_RNDN); mpfr_set_d(r22457, lambda2, MPFR_RNDN); mpfr_sub(r22458, r22454, r22457, MPFR_RNDN); mpfr_sin(r22459, r22458, MPFR_RNDN); mpfr_mul(r22460, r22456, r22459, MPFR_RNDN); mpfr_set_d(r22461, phi1, MPFR_RNDN); mpfr_cos(r22462, r22461, MPFR_RNDN); mpfr_cos(r22463, r22458, MPFR_RNDN); mpfr_mul(r22464, r22456, r22463, MPFR_RNDN); mpfr_add(r22465, r22462, r22464, MPFR_RNDN); mpfr_atan2(r22466, r22460, r22465, MPFR_RNDN); mpfr_add(r22467, r22454, r22466, MPFR_RNDN); return mpfr_get_d(r22467, MPFR_RNDN); } static mpfr_t r22468, r22469, r22470, r22471, r22472, r22473, r22474, r22475, r22476, r22477, r22478, r22479, r22480, r22481, r22482, r22483, r22484, r22485, r22486, r22487, r22488, r22489, r22490, r22491, r22492, r22493; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r22468); mpfr_init(r22469); mpfr_init(r22470); mpfr_init(r22471); mpfr_init(r22472); mpfr_init(r22473); mpfr_init(r22474); mpfr_init(r22475); mpfr_init(r22476); mpfr_init(r22477); mpfr_init(r22478); mpfr_init(r22479); mpfr_init(r22480); mpfr_init(r22481); mpfr_init(r22482); mpfr_init(r22483); mpfr_init(r22484); mpfr_init(r22485); mpfr_init(r22486); mpfr_init(r22487); mpfr_init(r22488); mpfr_init(r22489); mpfr_init(r22490); mpfr_init(r22491); mpfr_init(r22492); mpfr_init(r22493); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r22468, lambda1, MPFR_RNDN); mpfr_set_d(r22469, phi2, MPFR_RNDN); mpfr_cos(r22470, r22469, MPFR_RNDN); mpfr_sin(r22471, r22468, MPFR_RNDN); mpfr_set_d(r22472, lambda2, MPFR_RNDN); mpfr_cos(r22473, r22472, MPFR_RNDN); mpfr_mul(r22474, r22471, r22473, MPFR_RNDN); mpfr_cos(r22475, r22468, MPFR_RNDN); mpfr_neg(r22476, r22472, MPFR_RNDN); mpfr_sin(r22477, r22476, MPFR_RNDN); mpfr_mul(r22478, r22475, r22477, MPFR_RNDN); mpfr_add(r22479, r22474, r22478, MPFR_RNDN); mpfr_mul(r22480, r22470, r22479, MPFR_RNDN); mpfr_set_d(r22481, phi1, MPFR_RNDN); mpfr_cos(r22482, r22481, MPFR_RNDN); mpfr_exp(r22483, r22482, MPFR_RNDN); mpfr_exp(r22484, r22470, MPFR_RNDN); mpfr_sin(r22485, r22472, MPFR_RNDN); mpfr_mul(r22486, r22485, r22471, MPFR_RNDN); mpfr_mul(r22487, r22473, r22475, MPFR_RNDN); mpfr_add(r22488, r22486, r22487, MPFR_RNDN); mpfr_pow(r22489, r22484, r22488, MPFR_RNDN); mpfr_mul(r22490, r22483, r22489, MPFR_RNDN); mpfr_log(r22491, r22490, MPFR_RNDN); mpfr_atan2(r22492, r22480, r22491, MPFR_RNDN); mpfr_add(r22493, r22468, r22492, MPFR_RNDN); return mpfr_get_d(r22493, MPFR_RNDN); } static mpfr_t r22494, r22495, r22496, r22497, r22498, r22499, r22500, r22501, r22502, r22503, r22504, r22505, r22506, r22507, r22508, r22509, r22510, r22511, r22512, r22513, r22514, r22515, r22516, r22517, r22518, r22519; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r22494); mpfr_init(r22495); mpfr_init(r22496); mpfr_init(r22497); mpfr_init(r22498); mpfr_init(r22499); mpfr_init(r22500); mpfr_init(r22501); mpfr_init(r22502); mpfr_init(r22503); mpfr_init(r22504); mpfr_init(r22505); mpfr_init(r22506); mpfr_init(r22507); mpfr_init(r22508); mpfr_init(r22509); mpfr_init(r22510); mpfr_init(r22511); mpfr_init(r22512); mpfr_init(r22513); mpfr_init(r22514); mpfr_init(r22515); mpfr_init(r22516); mpfr_init(r22517); mpfr_init(r22518); mpfr_init(r22519); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r22494, lambda1, MPFR_RNDN); mpfr_set_d(r22495, phi2, MPFR_RNDN); mpfr_cos(r22496, r22495, MPFR_RNDN); mpfr_sin(r22497, r22494, MPFR_RNDN); mpfr_set_d(r22498, lambda2, MPFR_RNDN); mpfr_cos(r22499, r22498, MPFR_RNDN); mpfr_mul(r22500, r22497, r22499, MPFR_RNDN); mpfr_cos(r22501, r22494, MPFR_RNDN); mpfr_neg(r22502, r22498, MPFR_RNDN); mpfr_sin(r22503, r22502, MPFR_RNDN); mpfr_mul(r22504, r22501, r22503, MPFR_RNDN); mpfr_add(r22505, r22500, r22504, MPFR_RNDN); mpfr_mul(r22506, r22496, r22505, MPFR_RNDN); mpfr_set_d(r22507, phi1, MPFR_RNDN); mpfr_cos(r22508, r22507, MPFR_RNDN); mpfr_exp(r22509, r22508, MPFR_RNDN); mpfr_exp(r22510, r22496, MPFR_RNDN); mpfr_sin(r22511, r22498, MPFR_RNDN); mpfr_mul(r22512, r22511, r22497, MPFR_RNDN); mpfr_mul(r22513, r22499, r22501, MPFR_RNDN); mpfr_add(r22514, r22512, r22513, MPFR_RNDN); mpfr_pow(r22515, r22510, r22514, MPFR_RNDN); mpfr_mul(r22516, r22509, r22515, MPFR_RNDN); mpfr_log(r22517, r22516, MPFR_RNDN); mpfr_atan2(r22518, r22506, r22517, MPFR_RNDN); mpfr_add(r22519, r22494, r22518, MPFR_RNDN); return mpfr_get_d(r22519, MPFR_RNDN); }