#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 r21437 = phi1; float r21438 = sin(r21437); float r21439 = phi2; float r21440 = sin(r21439); float r21441 = r21438 * r21440; float r21442 = cos(r21437); float r21443 = cos(r21439); float r21444 = r21442 * r21443; float r21445 = lambda1; float r21446 = lambda2; float r21447 = r21445 - r21446; float r21448 = cos(r21447); float r21449 = r21444 * r21448; float r21450 = r21441 + r21449; float r21451 = acos(r21450); float r21452 = R; float r21453 = r21451 * r21452; return r21453; } double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) { double r21454 = phi1; double r21455 = sin(r21454); double r21456 = phi2; double r21457 = sin(r21456); double r21458 = r21455 * r21457; double r21459 = cos(r21454); double r21460 = cos(r21456); double r21461 = r21459 * r21460; double r21462 = lambda1; double r21463 = lambda2; double r21464 = r21462 - r21463; double r21465 = cos(r21464); double r21466 = r21461 * r21465; double r21467 = r21458 + r21466; double r21468 = acos(r21467); double r21469 = R; double r21470 = r21468 * r21469; return r21470; } double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) { float r21471 = phi1; float r21472 = sin(r21471); float r21473 = phi2; float r21474 = sin(r21473); float r21475 = r21472 * r21474; float r21476 = cos(r21473); float r21477 = cos(r21471); float r21478 = r21476 * r21477; float r21479 = exp(r21478); float r21480 = lambda2; float r21481 = cos(r21480); float r21482 = lambda1; float r21483 = cos(r21482); float r21484 = r21481 * r21483; float r21485 = sin(r21482); float r21486 = sin(r21480); float r21487 = r21485 * r21486; float r21488 = r21484 + r21487; float r21489 = pow(r21479, r21488); float r21490 = log(r21489); float r21491 = r21475 + r21490; float r21492 = acos(r21491); float r21493 = R; float r21494 = r21492 * r21493; return r21494; } double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) { double r21495 = phi1; double r21496 = sin(r21495); double r21497 = phi2; double r21498 = sin(r21497); double r21499 = r21496 * r21498; double r21500 = cos(r21497); double r21501 = cos(r21495); double r21502 = r21500 * r21501; double r21503 = exp(r21502); double r21504 = lambda2; double r21505 = cos(r21504); double r21506 = lambda1; double r21507 = cos(r21506); double r21508 = r21505 * r21507; double r21509 = sin(r21506); double r21510 = sin(r21504); double r21511 = r21509 * r21510; double r21512 = r21508 + r21511; double r21513 = pow(r21503, r21512); double r21514 = log(r21513); double r21515 = r21499 + r21514; double r21516 = acos(r21515); double r21517 = R; double r21518 = r21516 * r21517; return r21518; } 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 r21519, r21520, r21521, r21522, r21523, r21524, r21525, r21526, r21527, r21528, r21529, r21530, r21531, r21532, r21533, r21534, r21535; void setup_mpfr_f_im() { mpfr_set_default_prec(1936); 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_im(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r21519, phi1, MPFR_RNDN); mpfr_sin(r21520, r21519, MPFR_RNDN); mpfr_set_d(r21521, phi2, MPFR_RNDN); mpfr_sin(r21522, r21521, MPFR_RNDN); mpfr_mul(r21523, r21520, r21522, MPFR_RNDN); mpfr_cos(r21524, r21519, MPFR_RNDN); mpfr_cos(r21525, r21521, MPFR_RNDN); mpfr_mul(r21526, r21524, r21525, MPFR_RNDN); mpfr_set_d(r21527, lambda1, MPFR_RNDN); mpfr_set_d(r21528, lambda2, MPFR_RNDN); mpfr_sub(r21529, r21527, r21528, MPFR_RNDN); mpfr_cos(r21530, r21529, MPFR_RNDN); mpfr_mul(r21531, r21526, r21530, MPFR_RNDN); mpfr_add(r21532, r21523, r21531, MPFR_RNDN); mpfr_acos(r21533, r21532, MPFR_RNDN); mpfr_set_d(r21534, R, MPFR_RNDN); mpfr_mul(r21535, r21533, r21534, MPFR_RNDN); return mpfr_get_d(r21535, MPFR_RNDN); } static mpfr_t r21536, r21537, r21538, r21539, r21540, r21541, r21542, r21543, r21544, r21545, r21546, r21547, r21548, r21549, r21550, r21551, r21552, r21553, r21554, r21555, r21556, r21557, r21558, r21559; void setup_mpfr_f_fm() { mpfr_set_default_prec(1936); mpfr_init(r21536); mpfr_init(r21537); mpfr_init(r21538); mpfr_init(r21539); mpfr_init(r21540); mpfr_init(r21541); mpfr_init(r21542); mpfr_init(r21543); mpfr_init(r21544); mpfr_init(r21545); mpfr_init(r21546); mpfr_init(r21547); mpfr_init(r21548); mpfr_init(r21549); mpfr_init(r21550); mpfr_init(r21551); mpfr_init(r21552); mpfr_init(r21553); mpfr_init(r21554); mpfr_init(r21555); mpfr_init(r21556); mpfr_init(r21557); mpfr_init(r21558); mpfr_init(r21559); } double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r21536, phi1, MPFR_RNDN); mpfr_sin(r21537, r21536, MPFR_RNDN); mpfr_set_d(r21538, phi2, MPFR_RNDN); mpfr_sin(r21539, r21538, MPFR_RNDN); mpfr_mul(r21540, r21537, r21539, MPFR_RNDN); mpfr_cos(r21541, r21538, MPFR_RNDN); mpfr_cos(r21542, r21536, MPFR_RNDN); mpfr_mul(r21543, r21541, r21542, MPFR_RNDN); mpfr_exp(r21544, r21543, MPFR_RNDN); mpfr_set_d(r21545, lambda2, MPFR_RNDN); mpfr_cos(r21546, r21545, MPFR_RNDN); mpfr_set_d(r21547, lambda1, MPFR_RNDN); mpfr_cos(r21548, r21547, MPFR_RNDN); mpfr_mul(r21549, r21546, r21548, MPFR_RNDN); mpfr_sin(r21550, r21547, MPFR_RNDN); mpfr_sin(r21551, r21545, MPFR_RNDN); mpfr_mul(r21552, r21550, r21551, MPFR_RNDN); mpfr_add(r21553, r21549, r21552, MPFR_RNDN); mpfr_pow(r21554, r21544, r21553, MPFR_RNDN); mpfr_log(r21555, r21554, MPFR_RNDN); mpfr_add(r21556, r21540, r21555, MPFR_RNDN); mpfr_acos(r21557, r21556, MPFR_RNDN); mpfr_set_d(r21558, R, MPFR_RNDN); mpfr_mul(r21559, r21557, r21558, MPFR_RNDN); return mpfr_get_d(r21559, MPFR_RNDN); } static mpfr_t r21560, r21561, r21562, r21563, r21564, r21565, r21566, r21567, r21568, r21569, r21570, r21571, r21572, r21573, r21574, r21575, r21576, r21577, r21578, r21579, r21580, r21581, r21582, r21583; void setup_mpfr_f_dm() { mpfr_set_default_prec(1936); mpfr_init(r21560); mpfr_init(r21561); mpfr_init(r21562); mpfr_init(r21563); mpfr_init(r21564); mpfr_init(r21565); mpfr_init(r21566); mpfr_init(r21567); mpfr_init(r21568); mpfr_init(r21569); mpfr_init(r21570); mpfr_init(r21571); mpfr_init(r21572); mpfr_init(r21573); mpfr_init(r21574); mpfr_init(r21575); mpfr_init(r21576); mpfr_init(r21577); mpfr_init(r21578); mpfr_init(r21579); mpfr_init(r21580); mpfr_init(r21581); mpfr_init(r21582); mpfr_init(r21583); } double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r21560, phi1, MPFR_RNDN); mpfr_sin(r21561, r21560, MPFR_RNDN); mpfr_set_d(r21562, phi2, MPFR_RNDN); mpfr_sin(r21563, r21562, MPFR_RNDN); mpfr_mul(r21564, r21561, r21563, MPFR_RNDN); mpfr_cos(r21565, r21562, MPFR_RNDN); mpfr_cos(r21566, r21560, MPFR_RNDN); mpfr_mul(r21567, r21565, r21566, MPFR_RNDN); mpfr_exp(r21568, r21567, MPFR_RNDN); mpfr_set_d(r21569, lambda2, MPFR_RNDN); mpfr_cos(r21570, r21569, MPFR_RNDN); mpfr_set_d(r21571, lambda1, MPFR_RNDN); mpfr_cos(r21572, r21571, MPFR_RNDN); mpfr_mul(r21573, r21570, r21572, MPFR_RNDN); mpfr_sin(r21574, r21571, MPFR_RNDN); mpfr_sin(r21575, r21569, MPFR_RNDN); mpfr_mul(r21576, r21574, r21575, MPFR_RNDN); mpfr_add(r21577, r21573, r21576, MPFR_RNDN); mpfr_pow(r21578, r21568, r21577, MPFR_RNDN); mpfr_log(r21579, r21578, MPFR_RNDN); mpfr_add(r21580, r21564, r21579, MPFR_RNDN); mpfr_acos(r21581, r21580, MPFR_RNDN); mpfr_set_d(r21582, R, MPFR_RNDN); mpfr_mul(r21583, r21581, r21582, MPFR_RNDN); return mpfr_get_d(r21583, MPFR_RNDN); }