#include #include #include #include #include char *name = "Midpoint on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r8372 = lambda1; float r8373 = phi2; float r8374 = cos(r8373); float r8375 = lambda2; float r8376 = r8372 - r8375; float r8377 = sin(r8376); float r8378 = r8374 * r8377; float r8379 = phi1; float r8380 = cos(r8379); float r8381 = cos(r8376); float r8382 = r8374 * r8381; float r8383 = r8380 + r8382; float r8384 = atan2(r8378, r8383); float r8385 = r8372 + r8384; return r8385; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r8386 = lambda1; double r8387 = phi2; double r8388 = cos(r8387); double r8389 = lambda2; double r8390 = r8386 - r8389; double r8391 = sin(r8390); double r8392 = r8388 * r8391; double r8393 = phi1; double r8394 = cos(r8393); double r8395 = cos(r8390); double r8396 = r8388 * r8395; double r8397 = r8394 + r8396; double r8398 = atan2(r8392, r8397); double r8399 = r8386 + r8398; return r8399; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r8400 = phi2; float r8401 = cos(r8400); float r8402 = lambda1; float r8403 = sin(r8402); float r8404 = lambda2; float r8405 = cos(r8404); float r8406 = r8403 * r8405; float r8407 = cos(r8402); float r8408 = sin(r8404); float r8409 = r8407 * r8408; float r8410 = r8406 - r8409; float r8411 = r8401 * r8410; float r8412 = r8408 * r8403; float r8413 = r8412 * r8401; float r8414 = r8407 * r8405; float r8415 = r8401 * r8414; float r8416 = phi1; float r8417 = cos(r8416); float r8418 = r8415 + r8417; float r8419 = r8413 + r8418; float r8420 = atan2(r8411, r8419); float r8421 = r8420 + r8402; return r8421; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r8422 = phi2; double r8423 = cos(r8422); double r8424 = lambda1; double r8425 = sin(r8424); double r8426 = lambda2; double r8427 = cos(r8426); double r8428 = r8425 * r8427; double r8429 = cos(r8424); double r8430 = sin(r8426); double r8431 = r8429 * r8430; double r8432 = r8428 - r8431; double r8433 = r8423 * r8432; double r8434 = r8430 * r8425; double r8435 = r8434 * r8423; double r8436 = r8429 * r8427; double r8437 = r8423 * r8436; double r8438 = phi1; double r8439 = cos(r8438); double r8440 = r8437 + r8439; double r8441 = r8435 + r8440; double r8442 = atan2(r8433, r8441); double r8443 = r8442 + r8424; return r8443; } 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 r8444, r8445, r8446, r8447, r8448, r8449, r8450, r8451, r8452, r8453, r8454, r8455, r8456, r8457; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r8444); mpfr_init(r8445); mpfr_init(r8446); mpfr_init(r8447); mpfr_init(r8448); mpfr_init(r8449); mpfr_init(r8450); mpfr_init(r8451); mpfr_init(r8452); mpfr_init(r8453); mpfr_init(r8454); mpfr_init(r8455); mpfr_init(r8456); mpfr_init(r8457); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r8444, lambda1, MPFR_RNDN); mpfr_set_d(r8445, phi2, MPFR_RNDN); mpfr_cos(r8446, r8445, MPFR_RNDN); mpfr_set_d(r8447, lambda2, MPFR_RNDN); mpfr_sub(r8448, r8444, r8447, MPFR_RNDN); mpfr_sin(r8449, r8448, MPFR_RNDN); mpfr_mul(r8450, r8446, r8449, MPFR_RNDN); mpfr_set_d(r8451, phi1, MPFR_RNDN); mpfr_cos(r8452, r8451, MPFR_RNDN); mpfr_cos(r8453, r8448, MPFR_RNDN); mpfr_mul(r8454, r8446, r8453, MPFR_RNDN); mpfr_add(r8455, r8452, r8454, MPFR_RNDN); mpfr_atan2(r8456, r8450, r8455, MPFR_RNDN); mpfr_add(r8457, r8444, r8456, MPFR_RNDN); return mpfr_get_d(r8457, MPFR_RNDN); } static mpfr_t r8458, r8459, r8460, r8461, r8462, r8463, r8464, r8465, r8466, r8467, r8468, r8469, r8470, r8471, r8472, r8473, r8474, r8475, r8476, r8477, r8478, r8479; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r8458); mpfr_init(r8459); mpfr_init(r8460); mpfr_init(r8461); mpfr_init(r8462); mpfr_init(r8463); mpfr_init(r8464); mpfr_init(r8465); mpfr_init(r8466); mpfr_init(r8467); mpfr_init(r8468); mpfr_init(r8469); mpfr_init(r8470); mpfr_init(r8471); mpfr_init(r8472); mpfr_init(r8473); mpfr_init(r8474); mpfr_init(r8475); mpfr_init(r8476); mpfr_init(r8477); mpfr_init(r8478); mpfr_init(r8479); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r8458, phi2, MPFR_RNDN); mpfr_cos(r8459, r8458, MPFR_RNDN); mpfr_set_d(r8460, lambda1, MPFR_RNDN); mpfr_sin(r8461, r8460, MPFR_RNDN); mpfr_set_d(r8462, lambda2, MPFR_RNDN); mpfr_cos(r8463, r8462, MPFR_RNDN); mpfr_mul(r8464, r8461, r8463, MPFR_RNDN); mpfr_cos(r8465, r8460, MPFR_RNDN); mpfr_sin(r8466, r8462, MPFR_RNDN); mpfr_mul(r8467, r8465, r8466, MPFR_RNDN); mpfr_sub(r8468, r8464, r8467, MPFR_RNDN); mpfr_mul(r8469, r8459, r8468, MPFR_RNDN); mpfr_mul(r8470, r8466, r8461, MPFR_RNDN); mpfr_mul(r8471, r8470, r8459, MPFR_RNDN); mpfr_mul(r8472, r8465, r8463, MPFR_RNDN); mpfr_mul(r8473, r8459, r8472, MPFR_RNDN); mpfr_set_d(r8474, phi1, MPFR_RNDN); mpfr_cos(r8475, r8474, MPFR_RNDN); mpfr_add(r8476, r8473, r8475, MPFR_RNDN); mpfr_add(r8477, r8471, r8476, MPFR_RNDN); mpfr_atan2(r8478, r8469, r8477, MPFR_RNDN); mpfr_add(r8479, r8478, r8460, MPFR_RNDN); return mpfr_get_d(r8479, MPFR_RNDN); } static mpfr_t r8480, r8481, r8482, r8483, r8484, r8485, r8486, r8487, r8488, r8489, r8490, r8491, r8492, r8493, r8494, r8495, r8496, r8497, r8498, r8499, r8500, r8501; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r8480); mpfr_init(r8481); mpfr_init(r8482); mpfr_init(r8483); mpfr_init(r8484); mpfr_init(r8485); mpfr_init(r8486); mpfr_init(r8487); mpfr_init(r8488); mpfr_init(r8489); mpfr_init(r8490); mpfr_init(r8491); mpfr_init(r8492); mpfr_init(r8493); mpfr_init(r8494); mpfr_init(r8495); mpfr_init(r8496); mpfr_init(r8497); mpfr_init(r8498); mpfr_init(r8499); mpfr_init(r8500); mpfr_init(r8501); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r8480, phi2, MPFR_RNDN); mpfr_cos(r8481, r8480, MPFR_RNDN); mpfr_set_d(r8482, lambda1, MPFR_RNDN); mpfr_sin(r8483, r8482, MPFR_RNDN); mpfr_set_d(r8484, lambda2, MPFR_RNDN); mpfr_cos(r8485, r8484, MPFR_RNDN); mpfr_mul(r8486, r8483, r8485, MPFR_RNDN); mpfr_cos(r8487, r8482, MPFR_RNDN); mpfr_sin(r8488, r8484, MPFR_RNDN); mpfr_mul(r8489, r8487, r8488, MPFR_RNDN); mpfr_sub(r8490, r8486, r8489, MPFR_RNDN); mpfr_mul(r8491, r8481, r8490, MPFR_RNDN); mpfr_mul(r8492, r8488, r8483, MPFR_RNDN); mpfr_mul(r8493, r8492, r8481, MPFR_RNDN); mpfr_mul(r8494, r8487, r8485, MPFR_RNDN); mpfr_mul(r8495, r8481, r8494, MPFR_RNDN); mpfr_set_d(r8496, phi1, MPFR_RNDN); mpfr_cos(r8497, r8496, MPFR_RNDN); mpfr_add(r8498, r8495, r8497, MPFR_RNDN); mpfr_add(r8499, r8493, r8498, MPFR_RNDN); mpfr_atan2(r8500, r8491, r8499, MPFR_RNDN); mpfr_add(r8501, r8500, r8482, MPFR_RNDN); return mpfr_get_d(r8501, MPFR_RNDN); }