#include #include #include #include #include char *name = "Distance on a great circle"; double f_if(float R, float lambda1, float lambda2, float phi1, float phi2) { float r9527 = R; float r9528 = 2; float r9529 = phi1; float r9530 = phi2; float r9531 = r9529 - r9530; float r9532 = r9531 / r9528; float r9533 = sin(r9532); float r9534 = pow(r9533, r9528); float r9535 = cos(r9529); float r9536 = cos(r9530); float r9537 = r9535 * r9536; float r9538 = lambda1; float r9539 = lambda2; float r9540 = r9538 - r9539; float r9541 = r9540 / r9528; float r9542 = sin(r9541); float r9543 = r9537 * r9542; float r9544 = r9543 * r9542; float r9545 = r9534 + r9544; float r9546 = sqrt(r9545); float r9547 = 1; float r9548 = r9547 - r9545; float r9549 = sqrt(r9548); float r9550 = atan2(r9546, r9549); float r9551 = r9528 * r9550; float r9552 = r9527 * r9551; return r9552; } double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) { double r9553 = R; double r9554 = 2; double r9555 = phi1; double r9556 = phi2; double r9557 = r9555 - r9556; double r9558 = r9557 / r9554; double r9559 = sin(r9558); double r9560 = pow(r9559, r9554); double r9561 = cos(r9555); double r9562 = cos(r9556); double r9563 = r9561 * r9562; double r9564 = lambda1; double r9565 = lambda2; double r9566 = r9564 - r9565; double r9567 = r9566 / r9554; double r9568 = sin(r9567); double r9569 = r9563 * r9568; double r9570 = r9569 * r9568; double r9571 = r9560 + r9570; double r9572 = sqrt(r9571); double r9573 = 1; double r9574 = r9573 - r9571; double r9575 = sqrt(r9574); double r9576 = atan2(r9572, r9575); double r9577 = r9554 * r9576; double r9578 = r9553 * r9577; return r9578; } double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) { float r9579 = R; float r9580 = 2; float r9581 = phi1; float r9582 = phi2; float r9583 = r9581 - r9582; float r9584 = r9583 / r9580; float r9585 = sin(r9584); float r9586 = pow(r9585, r9580); float r9587 = cos(r9581); float r9588 = cos(r9582); float r9589 = r9587 * r9588; float r9590 = lambda1; float r9591 = lambda2; float r9592 = r9590 - r9591; float r9593 = r9592 / r9580; float r9594 = sin(r9593); float r9595 = r9589 * r9594; float r9596 = r9595 * r9594; float r9597 = r9586 + r9596; float r9598 = sqrt(r9597); float r9599 = 1; float r9600 = log1p(r9594); float r9601 = expm1(r9600); float r9602 = r9589 * r9601; float r9603 = r9602 * r9601; float r9604 = r9586 + r9603; float r9605 = r9599 - r9604; float r9606 = sqrt(r9605); float r9607 = atan2(r9598, r9606); float r9608 = r9580 * r9607; float r9609 = r9579 * r9608; return r9609; } double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) { double r9610 = R; double r9611 = 2; double r9612 = phi1; double r9613 = phi2; double r9614 = r9612 - r9613; double r9615 = r9614 / r9611; double r9616 = sin(r9615); double r9617 = pow(r9616, r9611); double r9618 = cos(r9612); double r9619 = cos(r9613); double r9620 = r9618 * r9619; double r9621 = lambda1; double r9622 = lambda2; double r9623 = r9621 - r9622; double r9624 = r9623 / r9611; double r9625 = sin(r9624); double r9626 = r9620 * r9625; double r9627 = r9626 * r9625; double r9628 = r9617 + r9627; double r9629 = sqrt(r9628); double r9630 = 1; double r9631 = log1p(r9625); double r9632 = expm1(r9631); double r9633 = r9620 * r9632; double r9634 = r9633 * r9632; double r9635 = r9617 + r9634; double r9636 = r9630 - r9635; double r9637 = sqrt(r9636); double r9638 = atan2(r9629, r9637); double r9639 = r9611 * r9638; double r9640 = r9610 * r9639; return r9640; } 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 r9641, r9642, r9643, r9644, r9645, r9646, r9647, r9648, r9649, r9650, r9651, r9652, r9653, r9654, r9655, r9656, r9657, r9658, r9659, r9660, r9661, r9662, r9663, r9664, r9665, r9666; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r9641); mpfr_init_set_str(r9642, "2", 10, MPFR_RNDN); mpfr_init(r9643); mpfr_init(r9644); mpfr_init(r9645); mpfr_init(r9646); mpfr_init(r9647); mpfr_init(r9648); mpfr_init(r9649); mpfr_init(r9650); mpfr_init(r9651); mpfr_init(r9652); mpfr_init(r9653); mpfr_init(r9654); mpfr_init(r9655); mpfr_init(r9656); mpfr_init(r9657); mpfr_init(r9658); mpfr_init(r9659); mpfr_init(r9660); mpfr_init_set_str(r9661, "1", 10, MPFR_RNDN); mpfr_init(r9662); mpfr_init(r9663); mpfr_init(r9664); mpfr_init(r9665); mpfr_init(r9666); } double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r9641, R, MPFR_RNDN); ; mpfr_set_d(r9643, phi1, MPFR_RNDN); mpfr_set_d(r9644, phi2, MPFR_RNDN); mpfr_sub(r9645, r9643, r9644, MPFR_RNDN); mpfr_div(r9646, r9645, r9642, MPFR_RNDN); mpfr_sin(r9647, r9646, MPFR_RNDN); mpfr_pow(r9648, r9647, r9642, MPFR_RNDN); mpfr_cos(r9649, r9643, MPFR_RNDN); mpfr_cos(r9650, r9644, MPFR_RNDN); mpfr_mul(r9651, r9649, r9650, MPFR_RNDN); mpfr_set_d(r9652, lambda1, MPFR_RNDN); mpfr_set_d(r9653, lambda2, MPFR_RNDN); mpfr_sub(r9654, r9652, r9653, MPFR_RNDN); mpfr_div(r9655, r9654, r9642, MPFR_RNDN); mpfr_sin(r9656, r9655, MPFR_RNDN); mpfr_mul(r9657, r9651, r9656, MPFR_RNDN); mpfr_mul(r9658, r9657, r9656, MPFR_RNDN); mpfr_add(r9659, r9648, r9658, MPFR_RNDN); mpfr_sqrt(r9660, r9659, MPFR_RNDN); ; mpfr_sub(r9662, r9661, r9659, MPFR_RNDN); mpfr_sqrt(r9663, r9662, MPFR_RNDN); mpfr_atan2(r9664, r9660, r9663, MPFR_RNDN); mpfr_mul(r9665, r9642, r9664, MPFR_RNDN); mpfr_mul(r9666, r9641, r9665, MPFR_RNDN); return mpfr_get_d(r9666, MPFR_RNDN); } static mpfr_t r9667, r9668, r9669, r9670, r9671, r9672, r9673, r9674, r9675, r9676, r9677, r9678, r9679, r9680, r9681, r9682, r9683, r9684, r9685, r9686, r9687, r9688, r9689, r9690, r9691, r9692, r9693, r9694, r9695, r9696, r9697; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r9667); mpfr_init_set_str(r9668, "2", 10, MPFR_RNDN); mpfr_init(r9669); mpfr_init(r9670); mpfr_init(r9671); mpfr_init(r9672); mpfr_init(r9673); mpfr_init(r9674); mpfr_init(r9675); mpfr_init(r9676); mpfr_init(r9677); mpfr_init(r9678); mpfr_init(r9679); mpfr_init(r9680); mpfr_init(r9681); mpfr_init(r9682); mpfr_init(r9683); mpfr_init(r9684); mpfr_init(r9685); mpfr_init(r9686); mpfr_init_set_str(r9687, "1", 10, MPFR_RNDN); mpfr_init(r9688); mpfr_init(r9689); mpfr_init(r9690); mpfr_init(r9691); mpfr_init(r9692); mpfr_init(r9693); mpfr_init(r9694); mpfr_init(r9695); mpfr_init(r9696); mpfr_init(r9697); } double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r9667, R, MPFR_RNDN); ; mpfr_set_d(r9669, phi1, MPFR_RNDN); mpfr_set_d(r9670, phi2, MPFR_RNDN); mpfr_sub(r9671, r9669, r9670, MPFR_RNDN); mpfr_div(r9672, r9671, r9668, MPFR_RNDN); mpfr_sin(r9673, r9672, MPFR_RNDN); mpfr_pow(r9674, r9673, r9668, MPFR_RNDN); mpfr_cos(r9675, r9669, MPFR_RNDN); mpfr_cos(r9676, r9670, MPFR_RNDN); mpfr_mul(r9677, r9675, r9676, MPFR_RNDN); mpfr_set_d(r9678, lambda1, MPFR_RNDN); mpfr_set_d(r9679, lambda2, MPFR_RNDN); mpfr_sub(r9680, r9678, r9679, MPFR_RNDN); mpfr_div(r9681, r9680, r9668, MPFR_RNDN); mpfr_sin(r9682, r9681, MPFR_RNDN); mpfr_mul(r9683, r9677, r9682, MPFR_RNDN); mpfr_mul(r9684, r9683, r9682, MPFR_RNDN); mpfr_add(r9685, r9674, r9684, MPFR_RNDN); mpfr_sqrt(r9686, r9685, MPFR_RNDN); ; mpfr_log1p(r9688, r9682, MPFR_RNDN); mpfr_expm1(r9689, r9688, MPFR_RNDN); mpfr_mul(r9690, r9677, r9689, MPFR_RNDN); mpfr_mul(r9691, r9690, r9689, MPFR_RNDN); mpfr_add(r9692, r9674, r9691, MPFR_RNDN); mpfr_sub(r9693, r9687, r9692, MPFR_RNDN); mpfr_sqrt(r9694, r9693, MPFR_RNDN); mpfr_atan2(r9695, r9686, r9694, MPFR_RNDN); mpfr_mul(r9696, r9668, r9695, MPFR_RNDN); mpfr_mul(r9697, r9667, r9696, MPFR_RNDN); return mpfr_get_d(r9697, MPFR_RNDN); } static mpfr_t r9698, r9699, r9700, r9701, r9702, r9703, r9704, r9705, r9706, r9707, r9708, r9709, r9710, r9711, r9712, r9713, r9714, r9715, r9716, r9717, r9718, r9719, r9720, r9721, r9722, r9723, r9724, r9725, r9726, r9727, r9728; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r9698); mpfr_init_set_str(r9699, "2", 10, MPFR_RNDN); mpfr_init(r9700); mpfr_init(r9701); mpfr_init(r9702); mpfr_init(r9703); mpfr_init(r9704); mpfr_init(r9705); mpfr_init(r9706); mpfr_init(r9707); mpfr_init(r9708); mpfr_init(r9709); mpfr_init(r9710); mpfr_init(r9711); mpfr_init(r9712); mpfr_init(r9713); mpfr_init(r9714); mpfr_init(r9715); mpfr_init(r9716); mpfr_init(r9717); mpfr_init_set_str(r9718, "1", 10, MPFR_RNDN); mpfr_init(r9719); mpfr_init(r9720); mpfr_init(r9721); mpfr_init(r9722); mpfr_init(r9723); mpfr_init(r9724); mpfr_init(r9725); mpfr_init(r9726); mpfr_init(r9727); mpfr_init(r9728); } double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r9698, R, MPFR_RNDN); ; mpfr_set_d(r9700, phi1, MPFR_RNDN); mpfr_set_d(r9701, phi2, MPFR_RNDN); mpfr_sub(r9702, r9700, r9701, MPFR_RNDN); mpfr_div(r9703, r9702, r9699, MPFR_RNDN); mpfr_sin(r9704, r9703, MPFR_RNDN); mpfr_pow(r9705, r9704, r9699, MPFR_RNDN); mpfr_cos(r9706, r9700, MPFR_RNDN); mpfr_cos(r9707, r9701, MPFR_RNDN); mpfr_mul(r9708, r9706, r9707, MPFR_RNDN); mpfr_set_d(r9709, lambda1, MPFR_RNDN); mpfr_set_d(r9710, lambda2, MPFR_RNDN); mpfr_sub(r9711, r9709, r9710, MPFR_RNDN); mpfr_div(r9712, r9711, r9699, MPFR_RNDN); mpfr_sin(r9713, r9712, MPFR_RNDN); mpfr_mul(r9714, r9708, r9713, MPFR_RNDN); mpfr_mul(r9715, r9714, r9713, MPFR_RNDN); mpfr_add(r9716, r9705, r9715, MPFR_RNDN); mpfr_sqrt(r9717, r9716, MPFR_RNDN); ; mpfr_log1p(r9719, r9713, MPFR_RNDN); mpfr_expm1(r9720, r9719, MPFR_RNDN); mpfr_mul(r9721, r9708, r9720, MPFR_RNDN); mpfr_mul(r9722, r9721, r9720, MPFR_RNDN); mpfr_add(r9723, r9705, r9722, MPFR_RNDN); mpfr_sub(r9724, r9718, r9723, MPFR_RNDN); mpfr_sqrt(r9725, r9724, MPFR_RNDN); mpfr_atan2(r9726, r9717, r9725, MPFR_RNDN); mpfr_mul(r9727, r9699, r9726, MPFR_RNDN); mpfr_mul(r9728, r9698, r9727, MPFR_RNDN); return mpfr_get_d(r9728, MPFR_RNDN); }