#include #include #include #include #include char *name = "Midpoint on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r21394 = lambda1; float r21395 = phi2; float r21396 = cos(r21395); float r21397 = lambda2; float r21398 = r21394 - r21397; float r21399 = sin(r21398); float r21400 = r21396 * r21399; float r21401 = phi1; float r21402 = cos(r21401); float r21403 = cos(r21398); float r21404 = r21396 * r21403; float r21405 = r21402 + r21404; float r21406 = atan2(r21400, r21405); float r21407 = r21394 + r21406; return r21407; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r21408 = lambda1; double r21409 = phi2; double r21410 = cos(r21409); double r21411 = lambda2; double r21412 = r21408 - r21411; double r21413 = sin(r21412); double r21414 = r21410 * r21413; double r21415 = phi1; double r21416 = cos(r21415); double r21417 = cos(r21412); double r21418 = r21410 * r21417; double r21419 = r21416 + r21418; double r21420 = atan2(r21414, r21419); double r21421 = r21408 + r21420; return r21421; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r21422 = lambda1; float r21423 = phi2; float r21424 = cos(r21423); float r21425 = sin(r21422); float r21426 = lambda2; float r21427 = cos(r21426); float r21428 = r21425 * r21427; float r21429 = cos(r21422); float r21430 = -r21426; float r21431 = sin(r21430); float r21432 = r21429 * r21431; float r21433 = r21428 + r21432; float r21434 = r21424 * r21433; float r21435 = phi1; float r21436 = cos(r21435); float r21437 = exp(r21436); float r21438 = r21429 * r21427; float r21439 = r21425 * r21431; float r21440 = r21438 - r21439; float r21441 = r21424 * r21440; float r21442 = exp(r21441); float r21443 = r21437 * r21442; float r21444 = log(r21443); float r21445 = atan2(r21434, r21444); float r21446 = r21422 + r21445; return r21446; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r21447 = lambda1; double r21448 = phi2; double r21449 = cos(r21448); double r21450 = sin(r21447); double r21451 = lambda2; double r21452 = cos(r21451); double r21453 = r21450 * r21452; double r21454 = cos(r21447); double r21455 = -r21451; double r21456 = sin(r21455); double r21457 = r21454 * r21456; double r21458 = r21453 + r21457; double r21459 = r21449 * r21458; double r21460 = phi1; double r21461 = cos(r21460); double r21462 = exp(r21461); double r21463 = r21454 * r21452; double r21464 = r21450 * r21456; double r21465 = r21463 - r21464; double r21466 = r21449 * r21465; double r21467 = exp(r21466); double r21468 = r21462 * r21467; double r21469 = log(r21468); double r21470 = atan2(r21459, r21469); double r21471 = r21447 + r21470; return r21471; } 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 r21472, r21473, r21474, r21475, r21476, r21477, r21478, r21479, r21480, r21481, r21482, r21483, r21484, r21485; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r21472); mpfr_init(r21473); mpfr_init(r21474); mpfr_init(r21475); mpfr_init(r21476); mpfr_init(r21477); mpfr_init(r21478); mpfr_init(r21479); mpfr_init(r21480); mpfr_init(r21481); mpfr_init(r21482); mpfr_init(r21483); mpfr_init(r21484); mpfr_init(r21485); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r21472, lambda1, MPFR_RNDN); mpfr_set_d(r21473, phi2, MPFR_RNDN); mpfr_cos(r21474, r21473, MPFR_RNDN); mpfr_set_d(r21475, lambda2, MPFR_RNDN); mpfr_sub(r21476, r21472, r21475, MPFR_RNDN); mpfr_sin(r21477, r21476, MPFR_RNDN); mpfr_mul(r21478, r21474, r21477, MPFR_RNDN); mpfr_set_d(r21479, phi1, MPFR_RNDN); mpfr_cos(r21480, r21479, MPFR_RNDN); mpfr_cos(r21481, r21476, MPFR_RNDN); mpfr_mul(r21482, r21474, r21481, MPFR_RNDN); mpfr_add(r21483, r21480, r21482, MPFR_RNDN); mpfr_atan2(r21484, r21478, r21483, MPFR_RNDN); mpfr_add(r21485, r21472, r21484, MPFR_RNDN); return mpfr_get_d(r21485, MPFR_RNDN); } static mpfr_t r21486, r21487, r21488, r21489, r21490, r21491, r21492, r21493, r21494, r21495, r21496, r21497, r21498, r21499, r21500, r21501, r21502, r21503, r21504, r21505, r21506, r21507, r21508, r21509, r21510; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r21486); mpfr_init(r21487); mpfr_init(r21488); mpfr_init(r21489); mpfr_init(r21490); mpfr_init(r21491); mpfr_init(r21492); mpfr_init(r21493); mpfr_init(r21494); mpfr_init(r21495); mpfr_init(r21496); mpfr_init(r21497); mpfr_init(r21498); mpfr_init(r21499); mpfr_init(r21500); mpfr_init(r21501); mpfr_init(r21502); mpfr_init(r21503); mpfr_init(r21504); mpfr_init(r21505); mpfr_init(r21506); mpfr_init(r21507); mpfr_init(r21508); mpfr_init(r21509); mpfr_init(r21510); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r21486, lambda1, MPFR_RNDN); mpfr_set_d(r21487, phi2, MPFR_RNDN); mpfr_cos(r21488, r21487, MPFR_RNDN); mpfr_sin(r21489, r21486, MPFR_RNDN); mpfr_set_d(r21490, lambda2, MPFR_RNDN); mpfr_cos(r21491, r21490, MPFR_RNDN); mpfr_mul(r21492, r21489, r21491, MPFR_RNDN); mpfr_cos(r21493, r21486, MPFR_RNDN); mpfr_neg(r21494, r21490, MPFR_RNDN); mpfr_sin(r21495, r21494, MPFR_RNDN); mpfr_mul(r21496, r21493, r21495, MPFR_RNDN); mpfr_add(r21497, r21492, r21496, MPFR_RNDN); mpfr_mul(r21498, r21488, r21497, MPFR_RNDN); mpfr_set_d(r21499, phi1, MPFR_RNDN); mpfr_cos(r21500, r21499, MPFR_RNDN); mpfr_exp(r21501, r21500, MPFR_RNDN); mpfr_mul(r21502, r21493, r21491, MPFR_RNDN); mpfr_mul(r21503, r21489, r21495, MPFR_RNDN); mpfr_sub(r21504, r21502, r21503, MPFR_RNDN); mpfr_mul(r21505, r21488, r21504, MPFR_RNDN); mpfr_exp(r21506, r21505, MPFR_RNDN); mpfr_mul(r21507, r21501, r21506, MPFR_RNDN); mpfr_log(r21508, r21507, MPFR_RNDN); mpfr_atan2(r21509, r21498, r21508, MPFR_RNDN); mpfr_add(r21510, r21486, r21509, MPFR_RNDN); return mpfr_get_d(r21510, MPFR_RNDN); } static mpfr_t r21511, r21512, r21513, r21514, r21515, r21516, r21517, r21518, r21519, r21520, r21521, r21522, r21523, r21524, r21525, r21526, r21527, r21528, r21529, r21530, r21531, r21532, r21533, r21534, r21535; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r21511); mpfr_init(r21512); mpfr_init(r21513); mpfr_init(r21514); mpfr_init(r21515); mpfr_init(r21516); mpfr_init(r21517); mpfr_init(r21518); mpfr_init(r21519); mpfr_init(r21520); mpfr_init(r21521); mpfr_init(r21522); mpfr_init(r21523); mpfr_init(r21524); mpfr_init(r21525); mpfr_init(r21526); mpfr_init(r21527); mpfr_init(r21528); mpfr_init(r21529); mpfr_init(r21530); mpfr_init(r21531); mpfr_init(r21532); mpfr_init(r21533); mpfr_init(r21534); mpfr_init(r21535); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r21511, lambda1, MPFR_RNDN); mpfr_set_d(r21512, phi2, MPFR_RNDN); mpfr_cos(r21513, r21512, MPFR_RNDN); mpfr_sin(r21514, r21511, MPFR_RNDN); mpfr_set_d(r21515, lambda2, MPFR_RNDN); mpfr_cos(r21516, r21515, MPFR_RNDN); mpfr_mul(r21517, r21514, r21516, MPFR_RNDN); mpfr_cos(r21518, r21511, MPFR_RNDN); mpfr_neg(r21519, r21515, MPFR_RNDN); mpfr_sin(r21520, r21519, MPFR_RNDN); mpfr_mul(r21521, r21518, r21520, MPFR_RNDN); mpfr_add(r21522, r21517, r21521, MPFR_RNDN); mpfr_mul(r21523, r21513, r21522, MPFR_RNDN); mpfr_set_d(r21524, phi1, MPFR_RNDN); mpfr_cos(r21525, r21524, MPFR_RNDN); mpfr_exp(r21526, r21525, MPFR_RNDN); mpfr_mul(r21527, r21518, r21516, MPFR_RNDN); mpfr_mul(r21528, r21514, r21520, MPFR_RNDN); mpfr_sub(r21529, r21527, r21528, MPFR_RNDN); mpfr_mul(r21530, r21513, r21529, MPFR_RNDN); mpfr_exp(r21531, r21530, MPFR_RNDN); mpfr_mul(r21532, r21526, r21531, MPFR_RNDN); mpfr_log(r21533, r21532, MPFR_RNDN); mpfr_atan2(r21534, r21523, r21533, MPFR_RNDN); mpfr_add(r21535, r21511, r21534, MPFR_RNDN); return mpfr_get_d(r21535, MPFR_RNDN); }