#include #include #include #include #include char *name = "Spherical law of cosines"; double f_if(float R, float lambda1, float lambda2, float phi1, float phi2) { float r25445 = phi1; float r25446 = sin(r25445); float r25447 = phi2; float r25448 = sin(r25447); float r25449 = r25446 * r25448; float r25450 = cos(r25445); float r25451 = cos(r25447); float r25452 = r25450 * r25451; float r25453 = lambda1; float r25454 = lambda2; float r25455 = r25453 - r25454; float r25456 = cos(r25455); float r25457 = r25452 * r25456; float r25458 = r25449 + r25457; float r25459 = acos(r25458); float r25460 = R; float r25461 = r25459 * r25460; return r25461; } double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) { double r25462 = phi1; double r25463 = sin(r25462); double r25464 = phi2; double r25465 = sin(r25464); double r25466 = r25463 * r25465; double r25467 = cos(r25462); double r25468 = cos(r25464); double r25469 = r25467 * r25468; double r25470 = lambda1; double r25471 = lambda2; double r25472 = r25470 - r25471; double r25473 = cos(r25472); double r25474 = r25469 * r25473; double r25475 = r25466 + r25474; double r25476 = acos(r25475); double r25477 = R; double r25478 = r25476 * r25477; return r25478; } double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) { float r25479 = phi1; float r25480 = sin(r25479); float r25481 = phi2; float r25482 = sin(r25481); float r25483 = r25480 * r25482; float r25484 = cos(r25479); float r25485 = cos(r25481); float r25486 = r25484 * r25485; float r25487 = lambda1; float r25488 = cos(r25487); float r25489 = lambda2; float r25490 = cos(r25489); float r25491 = r25488 * r25490; float r25492 = sin(r25489); float r25493 = -r25492; float r25494 = sin(r25487); float r25495 = r25492 * r25494; float r25496 = r25493 * r25495; float r25497 = r25494 * r25494; float r25498 = r25497 * r25492; float r25499 = r25496 * r25498; float r25500 = cbrt(r25499); float r25501 = r25491 - r25500; float r25502 = r25486 * r25501; float r25503 = r25483 + r25502; float r25504 = acos(r25503); float r25505 = R; float r25506 = r25504 * r25505; return r25506; } double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) { double r25507 = phi1; double r25508 = sin(r25507); double r25509 = phi2; double r25510 = sin(r25509); double r25511 = r25508 * r25510; double r25512 = cos(r25507); double r25513 = cos(r25509); double r25514 = r25512 * r25513; double r25515 = lambda1; double r25516 = cos(r25515); double r25517 = lambda2; double r25518 = cos(r25517); double r25519 = r25516 * r25518; double r25520 = sin(r25517); double r25521 = -r25520; double r25522 = sin(r25515); double r25523 = r25520 * r25522; double r25524 = r25521 * r25523; double r25525 = r25522 * r25522; double r25526 = r25525 * r25520; double r25527 = r25524 * r25526; double r25528 = cbrt(r25527); double r25529 = r25519 - r25528; double r25530 = r25514 * r25529; double r25531 = r25511 + r25530; double r25532 = acos(r25531); double r25533 = R; double r25534 = r25532 * r25533; return r25534; } 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 r25535, r25536, r25537, r25538, r25539, r25540, r25541, r25542, r25543, r25544, r25545, r25546, r25547, r25548, r25549, r25550, r25551; void setup_mpfr_f_im() { mpfr_set_default_prec(1936); mpfr_init(r25535); mpfr_init(r25536); mpfr_init(r25537); mpfr_init(r25538); mpfr_init(r25539); mpfr_init(r25540); mpfr_init(r25541); mpfr_init(r25542); mpfr_init(r25543); mpfr_init(r25544); mpfr_init(r25545); mpfr_init(r25546); mpfr_init(r25547); mpfr_init(r25548); mpfr_init(r25549); mpfr_init(r25550); mpfr_init(r25551); } double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r25535, phi1, MPFR_RNDN); mpfr_sin(r25536, r25535, MPFR_RNDN); mpfr_set_d(r25537, phi2, MPFR_RNDN); mpfr_sin(r25538, r25537, MPFR_RNDN); mpfr_mul(r25539, r25536, r25538, MPFR_RNDN); mpfr_cos(r25540, r25535, MPFR_RNDN); mpfr_cos(r25541, r25537, MPFR_RNDN); mpfr_mul(r25542, r25540, r25541, MPFR_RNDN); mpfr_set_d(r25543, lambda1, MPFR_RNDN); mpfr_set_d(r25544, lambda2, MPFR_RNDN); mpfr_sub(r25545, r25543, r25544, MPFR_RNDN); mpfr_cos(r25546, r25545, MPFR_RNDN); mpfr_mul(r25547, r25542, r25546, MPFR_RNDN); mpfr_add(r25548, r25539, r25547, MPFR_RNDN); mpfr_acos(r25549, r25548, MPFR_RNDN); mpfr_set_d(r25550, R, MPFR_RNDN); mpfr_mul(r25551, r25549, r25550, MPFR_RNDN); return mpfr_get_d(r25551, MPFR_RNDN); } static mpfr_t r25552, r25553, r25554, r25555, r25556, r25557, r25558, r25559, r25560, r25561, r25562, r25563, r25564, r25565, r25566, r25567, r25568, r25569, r25570, r25571, r25572, r25573, r25574, r25575, r25576, r25577, r25578, r25579; void setup_mpfr_f_fm() { mpfr_set_default_prec(1936); mpfr_init(r25552); mpfr_init(r25553); mpfr_init(r25554); mpfr_init(r25555); mpfr_init(r25556); mpfr_init(r25557); mpfr_init(r25558); mpfr_init(r25559); mpfr_init(r25560); mpfr_init(r25561); mpfr_init(r25562); mpfr_init(r25563); mpfr_init(r25564); mpfr_init(r25565); mpfr_init(r25566); mpfr_init(r25567); mpfr_init(r25568); mpfr_init(r25569); mpfr_init(r25570); mpfr_init(r25571); mpfr_init(r25572); mpfr_init(r25573); mpfr_init(r25574); mpfr_init(r25575); mpfr_init(r25576); mpfr_init(r25577); mpfr_init(r25578); mpfr_init(r25579); } double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r25552, phi1, MPFR_RNDN); mpfr_sin(r25553, r25552, MPFR_RNDN); mpfr_set_d(r25554, phi2, MPFR_RNDN); mpfr_sin(r25555, r25554, MPFR_RNDN); mpfr_mul(r25556, r25553, r25555, MPFR_RNDN); mpfr_cos(r25557, r25552, MPFR_RNDN); mpfr_cos(r25558, r25554, MPFR_RNDN); mpfr_mul(r25559, r25557, r25558, MPFR_RNDN); mpfr_set_d(r25560, lambda1, MPFR_RNDN); mpfr_cos(r25561, r25560, MPFR_RNDN); mpfr_set_d(r25562, lambda2, MPFR_RNDN); mpfr_cos(r25563, r25562, MPFR_RNDN); mpfr_mul(r25564, r25561, r25563, MPFR_RNDN); mpfr_sin(r25565, r25562, MPFR_RNDN); mpfr_neg(r25566, r25565, MPFR_RNDN); mpfr_sin(r25567, r25560, MPFR_RNDN); mpfr_mul(r25568, r25565, r25567, MPFR_RNDN); mpfr_mul(r25569, r25566, r25568, MPFR_RNDN); mpfr_mul(r25570, r25567, r25567, MPFR_RNDN); mpfr_mul(r25571, r25570, r25565, MPFR_RNDN); mpfr_mul(r25572, r25569, r25571, MPFR_RNDN); mpfr_cbrt(r25573, r25572, MPFR_RNDN); mpfr_sub(r25574, r25564, r25573, MPFR_RNDN); mpfr_mul(r25575, r25559, r25574, MPFR_RNDN); mpfr_add(r25576, r25556, r25575, MPFR_RNDN); mpfr_acos(r25577, r25576, MPFR_RNDN); mpfr_set_d(r25578, R, MPFR_RNDN); mpfr_mul(r25579, r25577, r25578, MPFR_RNDN); return mpfr_get_d(r25579, MPFR_RNDN); } static mpfr_t r25580, r25581, r25582, r25583, r25584, r25585, r25586, r25587, r25588, r25589, r25590, r25591, r25592, r25593, r25594, r25595, r25596, r25597, r25598, r25599, r25600, r25601, r25602, r25603, r25604, r25605, r25606, r25607; void setup_mpfr_f_dm() { mpfr_set_default_prec(1936); mpfr_init(r25580); mpfr_init(r25581); mpfr_init(r25582); mpfr_init(r25583); mpfr_init(r25584); mpfr_init(r25585); mpfr_init(r25586); mpfr_init(r25587); mpfr_init(r25588); mpfr_init(r25589); mpfr_init(r25590); mpfr_init(r25591); mpfr_init(r25592); mpfr_init(r25593); mpfr_init(r25594); mpfr_init(r25595); mpfr_init(r25596); mpfr_init(r25597); mpfr_init(r25598); mpfr_init(r25599); mpfr_init(r25600); mpfr_init(r25601); mpfr_init(r25602); mpfr_init(r25603); mpfr_init(r25604); mpfr_init(r25605); mpfr_init(r25606); mpfr_init(r25607); } double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r25580, phi1, MPFR_RNDN); mpfr_sin(r25581, r25580, MPFR_RNDN); mpfr_set_d(r25582, phi2, MPFR_RNDN); mpfr_sin(r25583, r25582, MPFR_RNDN); mpfr_mul(r25584, r25581, r25583, MPFR_RNDN); mpfr_cos(r25585, r25580, MPFR_RNDN); mpfr_cos(r25586, r25582, MPFR_RNDN); mpfr_mul(r25587, r25585, r25586, MPFR_RNDN); mpfr_set_d(r25588, lambda1, MPFR_RNDN); mpfr_cos(r25589, r25588, MPFR_RNDN); mpfr_set_d(r25590, lambda2, MPFR_RNDN); mpfr_cos(r25591, r25590, MPFR_RNDN); mpfr_mul(r25592, r25589, r25591, MPFR_RNDN); mpfr_sin(r25593, r25590, MPFR_RNDN); mpfr_neg(r25594, r25593, MPFR_RNDN); mpfr_sin(r25595, r25588, MPFR_RNDN); mpfr_mul(r25596, r25593, r25595, MPFR_RNDN); mpfr_mul(r25597, r25594, r25596, MPFR_RNDN); mpfr_mul(r25598, r25595, r25595, MPFR_RNDN); mpfr_mul(r25599, r25598, r25593, MPFR_RNDN); mpfr_mul(r25600, r25597, r25599, MPFR_RNDN); mpfr_cbrt(r25601, r25600, MPFR_RNDN); mpfr_sub(r25602, r25592, r25601, MPFR_RNDN); mpfr_mul(r25603, r25587, r25602, MPFR_RNDN); mpfr_add(r25604, r25584, r25603, MPFR_RNDN); mpfr_acos(r25605, r25604, MPFR_RNDN); mpfr_set_d(r25606, R, MPFR_RNDN); mpfr_mul(r25607, r25605, r25606, MPFR_RNDN); return mpfr_get_d(r25607, MPFR_RNDN); }