#include #include #include #include #include char *name = "Destination given bearing on a great circle"; double f_if(float lambda1, float phi1, float __attribute__((unused)) phi2, float delta, float theta) { float r23618 = lambda1; float r23619 = theta; float r23620 = sin(r23619); float r23621 = delta; float r23622 = sin(r23621); float r23623 = r23620 * r23622; float r23624 = phi1; float r23625 = cos(r23624); float r23626 = r23623 * r23625; float r23627 = cos(r23621); float r23628 = sin(r23624); float r23629 = r23628 * r23627; float r23630 = r23625 * r23622; float r23631 = cos(r23619); float r23632 = r23630 * r23631; float r23633 = r23629 + r23632; float r23634 = asin(r23633); float r23635 = sin(r23634); float r23636 = r23628 * r23635; float r23637 = r23627 - r23636; float r23638 = atan2(r23626, r23637); float r23639 = r23618 + r23638; return r23639; } double f_id(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { double r23640 = lambda1; double r23641 = theta; double r23642 = sin(r23641); double r23643 = delta; double r23644 = sin(r23643); double r23645 = r23642 * r23644; double r23646 = phi1; double r23647 = cos(r23646); double r23648 = r23645 * r23647; double r23649 = cos(r23643); double r23650 = sin(r23646); double r23651 = r23650 * r23649; double r23652 = r23647 * r23644; double r23653 = cos(r23641); double r23654 = r23652 * r23653; double r23655 = r23651 + r23654; double r23656 = asin(r23655); double r23657 = sin(r23656); double r23658 = r23650 * r23657; double r23659 = r23649 - r23658; double r23660 = atan2(r23648, r23659); double r23661 = r23640 + r23660; return r23661; } double f_of(float lambda1, float phi1, float __attribute__((unused)) phi2, float delta, float theta) { float r23662 = lambda1; float r23663 = theta; float r23664 = sin(r23663); float r23665 = delta; float r23666 = sin(r23665); float r23667 = r23664 * r23666; float r23668 = phi1; float r23669 = cos(r23668); float r23670 = r23667 * r23669; float r23671 = cos(r23665); float r23672 = sin(r23668); float r23673 = cos(r23663); float r23674 = r23673 * r23669; float r23675 = r23666 * r23674; float r23676 = r23672 * r23675; float r23677 = 2; float r23678 = pow(r23672, r23677); float r23679 = exp(r23678); float r23680 = log(r23679); float r23681 = r23671 * r23680; float r23682 = r23676 + r23681; float r23683 = r23671 - r23682; float r23684 = atan2(r23670, r23683); float r23685 = r23662 + r23684; return r23685; } double f_od(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { double r23686 = lambda1; double r23687 = theta; double r23688 = sin(r23687); double r23689 = delta; double r23690 = sin(r23689); double r23691 = r23688 * r23690; double r23692 = phi1; double r23693 = cos(r23692); double r23694 = r23691 * r23693; double r23695 = cos(r23689); double r23696 = sin(r23692); double r23697 = cos(r23687); double r23698 = r23697 * r23693; double r23699 = r23690 * r23698; double r23700 = r23696 * r23699; double r23701 = 2; double r23702 = pow(r23696, r23701); double r23703 = exp(r23702); double r23704 = log(r23703); double r23705 = r23695 * r23704; double r23706 = r23700 + r23705; double r23707 = r23695 - r23706; double r23708 = atan2(r23694, r23707); double r23709 = r23686 + r23708; return r23709; } 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 r23710, r23711, r23712, r23713, r23714, r23715, r23716, r23717, r23718, r23719, r23720, r23721, r23722, r23723, r23724, r23725, r23726, r23727, r23728, r23729, r23730, r23731; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r23710); mpfr_init(r23711); mpfr_init(r23712); mpfr_init(r23713); mpfr_init(r23714); mpfr_init(r23715); mpfr_init(r23716); mpfr_init(r23717); mpfr_init(r23718); mpfr_init(r23719); mpfr_init(r23720); mpfr_init(r23721); mpfr_init(r23722); mpfr_init(r23723); mpfr_init(r23724); mpfr_init(r23725); mpfr_init(r23726); mpfr_init(r23727); mpfr_init(r23728); mpfr_init(r23729); mpfr_init(r23730); mpfr_init(r23731); } double f_im(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { mpfr_set_d(r23710, lambda1, MPFR_RNDN); mpfr_set_d(r23711, theta, MPFR_RNDN); mpfr_sin(r23712, r23711, MPFR_RNDN); mpfr_set_d(r23713, delta, MPFR_RNDN); mpfr_sin(r23714, r23713, MPFR_RNDN); mpfr_mul(r23715, r23712, r23714, MPFR_RNDN); mpfr_set_d(r23716, phi1, MPFR_RNDN); mpfr_cos(r23717, r23716, MPFR_RNDN); mpfr_mul(r23718, r23715, r23717, MPFR_RNDN); mpfr_cos(r23719, r23713, MPFR_RNDN); mpfr_sin(r23720, r23716, MPFR_RNDN); mpfr_mul(r23721, r23720, r23719, MPFR_RNDN); mpfr_mul(r23722, r23717, r23714, MPFR_RNDN); mpfr_cos(r23723, r23711, MPFR_RNDN); mpfr_mul(r23724, r23722, r23723, MPFR_RNDN); mpfr_add(r23725, r23721, r23724, MPFR_RNDN); mpfr_asin(r23726, r23725, MPFR_RNDN); mpfr_sin(r23727, r23726, MPFR_RNDN); mpfr_mul(r23728, r23720, r23727, MPFR_RNDN); mpfr_sub(r23729, r23719, r23728, MPFR_RNDN); mpfr_atan2(r23730, r23718, r23729, MPFR_RNDN); mpfr_add(r23731, r23710, r23730, MPFR_RNDN); return mpfr_get_d(r23731, MPFR_RNDN); } static mpfr_t r23732, r23733, r23734, r23735, r23736, r23737, r23738, r23739, r23740, r23741, r23742, r23743, r23744, r23745, r23746, r23747, r23748, r23749, r23750, r23751, r23752, r23753, r23754, r23755; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r23732); mpfr_init(r23733); mpfr_init(r23734); mpfr_init(r23735); mpfr_init(r23736); mpfr_init(r23737); mpfr_init(r23738); mpfr_init(r23739); mpfr_init(r23740); mpfr_init(r23741); mpfr_init(r23742); mpfr_init(r23743); mpfr_init(r23744); mpfr_init(r23745); mpfr_init(r23746); mpfr_init_set_str(r23747, "2", 10, MPFR_RNDN); mpfr_init(r23748); mpfr_init(r23749); mpfr_init(r23750); mpfr_init(r23751); mpfr_init(r23752); mpfr_init(r23753); mpfr_init(r23754); mpfr_init(r23755); } double f_fm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { mpfr_set_d(r23732, lambda1, MPFR_RNDN); mpfr_set_d(r23733, theta, MPFR_RNDN); mpfr_sin(r23734, r23733, MPFR_RNDN); mpfr_set_d(r23735, delta, MPFR_RNDN); mpfr_sin(r23736, r23735, MPFR_RNDN); mpfr_mul(r23737, r23734, r23736, MPFR_RNDN); mpfr_set_d(r23738, phi1, MPFR_RNDN); mpfr_cos(r23739, r23738, MPFR_RNDN); mpfr_mul(r23740, r23737, r23739, MPFR_RNDN); mpfr_cos(r23741, r23735, MPFR_RNDN); mpfr_sin(r23742, r23738, MPFR_RNDN); mpfr_cos(r23743, r23733, MPFR_RNDN); mpfr_mul(r23744, r23743, r23739, MPFR_RNDN); mpfr_mul(r23745, r23736, r23744, MPFR_RNDN); mpfr_mul(r23746, r23742, r23745, MPFR_RNDN); ; mpfr_pow(r23748, r23742, r23747, MPFR_RNDN); mpfr_exp(r23749, r23748, MPFR_RNDN); mpfr_log(r23750, r23749, MPFR_RNDN); mpfr_mul(r23751, r23741, r23750, MPFR_RNDN); mpfr_add(r23752, r23746, r23751, MPFR_RNDN); mpfr_sub(r23753, r23741, r23752, MPFR_RNDN); mpfr_atan2(r23754, r23740, r23753, MPFR_RNDN); mpfr_add(r23755, r23732, r23754, MPFR_RNDN); return mpfr_get_d(r23755, MPFR_RNDN); } static mpfr_t r23756, r23757, r23758, r23759, r23760, r23761, r23762, r23763, r23764, r23765, r23766, r23767, r23768, r23769, r23770, r23771, r23772, r23773, r23774, r23775, r23776, r23777, r23778, r23779; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r23756); mpfr_init(r23757); mpfr_init(r23758); mpfr_init(r23759); mpfr_init(r23760); mpfr_init(r23761); mpfr_init(r23762); mpfr_init(r23763); mpfr_init(r23764); mpfr_init(r23765); mpfr_init(r23766); mpfr_init(r23767); mpfr_init(r23768); mpfr_init(r23769); mpfr_init(r23770); mpfr_init_set_str(r23771, "2", 10, MPFR_RNDN); mpfr_init(r23772); mpfr_init(r23773); mpfr_init(r23774); mpfr_init(r23775); mpfr_init(r23776); mpfr_init(r23777); mpfr_init(r23778); mpfr_init(r23779); } double f_dm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { mpfr_set_d(r23756, lambda1, MPFR_RNDN); mpfr_set_d(r23757, theta, MPFR_RNDN); mpfr_sin(r23758, r23757, MPFR_RNDN); mpfr_set_d(r23759, delta, MPFR_RNDN); mpfr_sin(r23760, r23759, MPFR_RNDN); mpfr_mul(r23761, r23758, r23760, MPFR_RNDN); mpfr_set_d(r23762, phi1, MPFR_RNDN); mpfr_cos(r23763, r23762, MPFR_RNDN); mpfr_mul(r23764, r23761, r23763, MPFR_RNDN); mpfr_cos(r23765, r23759, MPFR_RNDN); mpfr_sin(r23766, r23762, MPFR_RNDN); mpfr_cos(r23767, r23757, MPFR_RNDN); mpfr_mul(r23768, r23767, r23763, MPFR_RNDN); mpfr_mul(r23769, r23760, r23768, MPFR_RNDN); mpfr_mul(r23770, r23766, r23769, MPFR_RNDN); ; mpfr_pow(r23772, r23766, r23771, MPFR_RNDN); mpfr_exp(r23773, r23772, MPFR_RNDN); mpfr_log(r23774, r23773, MPFR_RNDN); mpfr_mul(r23775, r23765, r23774, MPFR_RNDN); mpfr_add(r23776, r23770, r23775, MPFR_RNDN); mpfr_sub(r23777, r23765, r23776, MPFR_RNDN); mpfr_atan2(r23778, r23764, r23777, MPFR_RNDN); mpfr_add(r23779, r23756, r23778, MPFR_RNDN); return mpfr_get_d(r23779, MPFR_RNDN); }