#include #include #include #include #include char *name = "Midpoint on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r18325 = lambda1; float r18326 = phi2; float r18327 = cos(r18326); float r18328 = lambda2; float r18329 = r18325 - r18328; float r18330 = sin(r18329); float r18331 = r18327 * r18330; float r18332 = phi1; float r18333 = cos(r18332); float r18334 = cos(r18329); float r18335 = r18327 * r18334; float r18336 = r18333 + r18335; float r18337 = atan2(r18331, r18336); float r18338 = r18325 + r18337; return r18338; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r18339 = lambda1; double r18340 = phi2; double r18341 = cos(r18340); double r18342 = lambda2; double r18343 = r18339 - r18342; double r18344 = sin(r18343); double r18345 = r18341 * r18344; double r18346 = phi1; double r18347 = cos(r18346); double r18348 = cos(r18343); double r18349 = r18341 * r18348; double r18350 = r18347 + r18349; double r18351 = atan2(r18345, r18350); double r18352 = r18339 + r18351; return r18352; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r18353 = lambda1; float r18354 = phi2; float r18355 = cos(r18354); float r18356 = sin(r18353); float r18357 = lambda2; float r18358 = cos(r18357); float r18359 = r18356 * r18358; float r18360 = cos(r18353); float r18361 = -r18357; float r18362 = sin(r18361); float r18363 = r18360 * r18362; float r18364 = r18359 + r18363; float r18365 = r18355 * r18364; float r18366 = phi1; float r18367 = cos(r18366); float r18368 = r18360 * r18358; float r18369 = r18356 * r18362; float r18370 = r18368 - r18369; float r18371 = r18355 * r18370; float r18372 = r18367 + r18371; float r18373 = atan2(r18365, r18372); float r18374 = r18353 + r18373; return r18374; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r18375 = lambda1; double r18376 = phi2; double r18377 = cos(r18376); double r18378 = sin(r18375); double r18379 = lambda2; double r18380 = cos(r18379); double r18381 = r18378 * r18380; double r18382 = cos(r18375); double r18383 = -r18379; double r18384 = sin(r18383); double r18385 = r18382 * r18384; double r18386 = r18381 + r18385; double r18387 = r18377 * r18386; double r18388 = phi1; double r18389 = cos(r18388); double r18390 = r18382 * r18380; double r18391 = r18378 * r18384; double r18392 = r18390 - r18391; double r18393 = r18377 * r18392; double r18394 = r18389 + r18393; double r18395 = atan2(r18387, r18394); double r18396 = r18375 + r18395; return r18396; } 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 r18397, r18398, r18399, r18400, r18401, r18402, r18403, r18404, r18405, r18406, r18407, r18408, r18409, r18410; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r18397); mpfr_init(r18398); mpfr_init(r18399); mpfr_init(r18400); mpfr_init(r18401); mpfr_init(r18402); mpfr_init(r18403); mpfr_init(r18404); mpfr_init(r18405); mpfr_init(r18406); mpfr_init(r18407); mpfr_init(r18408); mpfr_init(r18409); mpfr_init(r18410); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r18397, lambda1, MPFR_RNDN); mpfr_set_d(r18398, phi2, MPFR_RNDN); mpfr_cos(r18399, r18398, MPFR_RNDN); mpfr_set_d(r18400, lambda2, MPFR_RNDN); mpfr_sub(r18401, r18397, r18400, MPFR_RNDN); mpfr_sin(r18402, r18401, MPFR_RNDN); mpfr_mul(r18403, r18399, r18402, MPFR_RNDN); mpfr_set_d(r18404, phi1, MPFR_RNDN); mpfr_cos(r18405, r18404, MPFR_RNDN); mpfr_cos(r18406, r18401, MPFR_RNDN); mpfr_mul(r18407, r18399, r18406, MPFR_RNDN); mpfr_add(r18408, r18405, r18407, MPFR_RNDN); mpfr_atan2(r18409, r18403, r18408, MPFR_RNDN); mpfr_add(r18410, r18397, r18409, MPFR_RNDN); return mpfr_get_d(r18410, MPFR_RNDN); } static mpfr_t r18411, r18412, r18413, r18414, r18415, r18416, r18417, r18418, r18419, r18420, r18421, r18422, r18423, r18424, r18425, r18426, r18427, r18428, r18429, r18430, r18431, r18432; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r18411); mpfr_init(r18412); mpfr_init(r18413); mpfr_init(r18414); mpfr_init(r18415); mpfr_init(r18416); mpfr_init(r18417); mpfr_init(r18418); mpfr_init(r18419); mpfr_init(r18420); mpfr_init(r18421); mpfr_init(r18422); mpfr_init(r18423); mpfr_init(r18424); mpfr_init(r18425); mpfr_init(r18426); mpfr_init(r18427); mpfr_init(r18428); mpfr_init(r18429); mpfr_init(r18430); mpfr_init(r18431); mpfr_init(r18432); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r18411, lambda1, MPFR_RNDN); mpfr_set_d(r18412, phi2, MPFR_RNDN); mpfr_cos(r18413, r18412, MPFR_RNDN); mpfr_sin(r18414, r18411, MPFR_RNDN); mpfr_set_d(r18415, lambda2, MPFR_RNDN); mpfr_cos(r18416, r18415, MPFR_RNDN); mpfr_mul(r18417, r18414, r18416, MPFR_RNDN); mpfr_cos(r18418, r18411, MPFR_RNDN); mpfr_neg(r18419, r18415, MPFR_RNDN); mpfr_sin(r18420, r18419, MPFR_RNDN); mpfr_mul(r18421, r18418, r18420, MPFR_RNDN); mpfr_add(r18422, r18417, r18421, MPFR_RNDN); mpfr_mul(r18423, r18413, r18422, MPFR_RNDN); mpfr_set_d(r18424, phi1, MPFR_RNDN); mpfr_cos(r18425, r18424, MPFR_RNDN); mpfr_mul(r18426, r18418, r18416, MPFR_RNDN); mpfr_mul(r18427, r18414, r18420, MPFR_RNDN); mpfr_sub(r18428, r18426, r18427, MPFR_RNDN); mpfr_mul(r18429, r18413, r18428, MPFR_RNDN); mpfr_add(r18430, r18425, r18429, MPFR_RNDN); mpfr_atan2(r18431, r18423, r18430, MPFR_RNDN); mpfr_add(r18432, r18411, r18431, MPFR_RNDN); return mpfr_get_d(r18432, MPFR_RNDN); } static mpfr_t r18433, r18434, r18435, r18436, r18437, r18438, r18439, r18440, r18441, r18442, r18443, r18444, r18445, r18446, r18447, r18448, r18449, r18450, r18451, r18452, r18453, r18454; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r18433); mpfr_init(r18434); mpfr_init(r18435); mpfr_init(r18436); mpfr_init(r18437); mpfr_init(r18438); mpfr_init(r18439); mpfr_init(r18440); mpfr_init(r18441); mpfr_init(r18442); mpfr_init(r18443); mpfr_init(r18444); mpfr_init(r18445); mpfr_init(r18446); mpfr_init(r18447); mpfr_init(r18448); mpfr_init(r18449); mpfr_init(r18450); mpfr_init(r18451); mpfr_init(r18452); mpfr_init(r18453); mpfr_init(r18454); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r18433, lambda1, MPFR_RNDN); mpfr_set_d(r18434, phi2, MPFR_RNDN); mpfr_cos(r18435, r18434, MPFR_RNDN); mpfr_sin(r18436, r18433, MPFR_RNDN); mpfr_set_d(r18437, lambda2, MPFR_RNDN); mpfr_cos(r18438, r18437, MPFR_RNDN); mpfr_mul(r18439, r18436, r18438, MPFR_RNDN); mpfr_cos(r18440, r18433, MPFR_RNDN); mpfr_neg(r18441, r18437, MPFR_RNDN); mpfr_sin(r18442, r18441, MPFR_RNDN); mpfr_mul(r18443, r18440, r18442, MPFR_RNDN); mpfr_add(r18444, r18439, r18443, MPFR_RNDN); mpfr_mul(r18445, r18435, r18444, MPFR_RNDN); mpfr_set_d(r18446, phi1, MPFR_RNDN); mpfr_cos(r18447, r18446, MPFR_RNDN); mpfr_mul(r18448, r18440, r18438, MPFR_RNDN); mpfr_mul(r18449, r18436, r18442, MPFR_RNDN); mpfr_sub(r18450, r18448, r18449, MPFR_RNDN); mpfr_mul(r18451, r18435, r18450, MPFR_RNDN); mpfr_add(r18452, r18447, r18451, MPFR_RNDN); mpfr_atan2(r18453, r18445, r18452, MPFR_RNDN); mpfr_add(r18454, r18433, r18453, MPFR_RNDN); return mpfr_get_d(r18454, MPFR_RNDN); }