#include #include #include #include #include char *name = "Midpoint on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r24742 = lambda1; float r24743 = phi2; float r24744 = cos(r24743); float r24745 = lambda2; float r24746 = r24742 - r24745; float r24747 = sin(r24746); float r24748 = r24744 * r24747; float r24749 = phi1; float r24750 = cos(r24749); float r24751 = cos(r24746); float r24752 = r24744 * r24751; float r24753 = r24750 + r24752; float r24754 = atan2(r24748, r24753); float r24755 = r24742 + r24754; return r24755; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r24756 = lambda1; double r24757 = phi2; double r24758 = cos(r24757); double r24759 = lambda2; double r24760 = r24756 - r24759; double r24761 = sin(r24760); double r24762 = r24758 * r24761; double r24763 = phi1; double r24764 = cos(r24763); double r24765 = cos(r24760); double r24766 = r24758 * r24765; double r24767 = r24764 + r24766; double r24768 = atan2(r24762, r24767); double r24769 = r24756 + r24768; return r24769; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r24770 = phi2; float r24771 = cos(r24770); float r24772 = lambda1; float r24773 = lambda2; float r24774 = r24772 - r24773; float r24775 = sin(r24774); float r24776 = r24771 * r24775; float r24777 = cos(r24774); float r24778 = log1p(r24777); float r24779 = expm1(r24778); float r24780 = phi1; float r24781 = cos(r24780); float r24782 = fma(r24771, r24779, r24781); float r24783 = atan2(r24776, r24782); float r24784 = r24783 + r24772; return r24784; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r24785 = phi2; double r24786 = cos(r24785); double r24787 = lambda1; double r24788 = lambda2; double r24789 = r24787 - r24788; double r24790 = sin(r24789); double r24791 = r24786 * r24790; double r24792 = cos(r24789); double r24793 = log1p(r24792); double r24794 = expm1(r24793); double r24795 = phi1; double r24796 = cos(r24795); double r24797 = fma(r24786, r24794, r24796); double r24798 = atan2(r24791, r24797); double r24799 = r24798 + r24787; return r24799; } 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 r24800, r24801, r24802, r24803, r24804, r24805, r24806, r24807, r24808, r24809, r24810, r24811, r24812, r24813; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r24800); mpfr_init(r24801); mpfr_init(r24802); mpfr_init(r24803); mpfr_init(r24804); mpfr_init(r24805); mpfr_init(r24806); mpfr_init(r24807); mpfr_init(r24808); mpfr_init(r24809); mpfr_init(r24810); mpfr_init(r24811); mpfr_init(r24812); mpfr_init(r24813); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24800, lambda1, MPFR_RNDN); mpfr_set_d(r24801, phi2, MPFR_RNDN); mpfr_cos(r24802, r24801, MPFR_RNDN); mpfr_set_d(r24803, lambda2, MPFR_RNDN); mpfr_sub(r24804, r24800, r24803, MPFR_RNDN); mpfr_sin(r24805, r24804, MPFR_RNDN); mpfr_mul(r24806, r24802, r24805, MPFR_RNDN); mpfr_set_d(r24807, phi1, MPFR_RNDN); mpfr_cos(r24808, r24807, MPFR_RNDN); mpfr_cos(r24809, r24804, MPFR_RNDN); mpfr_mul(r24810, r24802, r24809, MPFR_RNDN); mpfr_add(r24811, r24808, r24810, MPFR_RNDN); mpfr_atan2(r24812, r24806, r24811, MPFR_RNDN); mpfr_add(r24813, r24800, r24812, MPFR_RNDN); return mpfr_get_d(r24813, MPFR_RNDN); } static mpfr_t r24814, r24815, r24816, r24817, r24818, r24819, r24820, r24821, r24822, r24823, r24824, r24825, r24826, r24827, r24828; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r24814); mpfr_init(r24815); mpfr_init(r24816); mpfr_init(r24817); mpfr_init(r24818); mpfr_init(r24819); mpfr_init(r24820); mpfr_init(r24821); mpfr_init(r24822); mpfr_init(r24823); mpfr_init(r24824); mpfr_init(r24825); mpfr_init(r24826); mpfr_init(r24827); mpfr_init(r24828); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24814, phi2, MPFR_RNDN); mpfr_cos(r24815, r24814, MPFR_RNDN); mpfr_set_d(r24816, lambda1, MPFR_RNDN); mpfr_set_d(r24817, lambda2, MPFR_RNDN); mpfr_sub(r24818, r24816, r24817, MPFR_RNDN); mpfr_sin(r24819, r24818, MPFR_RNDN); mpfr_mul(r24820, r24815, r24819, MPFR_RNDN); mpfr_cos(r24821, r24818, MPFR_RNDN); mpfr_log1p(r24822, r24821, MPFR_RNDN); mpfr_expm1(r24823, r24822, MPFR_RNDN); mpfr_set_d(r24824, phi1, MPFR_RNDN); mpfr_cos(r24825, r24824, MPFR_RNDN); mpfr_fma(r24826, r24815, r24823, r24825, MPFR_RNDN); mpfr_atan2(r24827, r24820, r24826, MPFR_RNDN); mpfr_add(r24828, r24827, r24816, MPFR_RNDN); return mpfr_get_d(r24828, MPFR_RNDN); } static mpfr_t r24829, r24830, r24831, r24832, r24833, r24834, r24835, r24836, r24837, r24838, r24839, r24840, r24841, r24842, r24843; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r24829); mpfr_init(r24830); mpfr_init(r24831); mpfr_init(r24832); mpfr_init(r24833); mpfr_init(r24834); mpfr_init(r24835); mpfr_init(r24836); mpfr_init(r24837); mpfr_init(r24838); mpfr_init(r24839); mpfr_init(r24840); mpfr_init(r24841); mpfr_init(r24842); mpfr_init(r24843); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24829, phi2, MPFR_RNDN); mpfr_cos(r24830, r24829, MPFR_RNDN); mpfr_set_d(r24831, lambda1, MPFR_RNDN); mpfr_set_d(r24832, lambda2, MPFR_RNDN); mpfr_sub(r24833, r24831, r24832, MPFR_RNDN); mpfr_sin(r24834, r24833, MPFR_RNDN); mpfr_mul(r24835, r24830, r24834, MPFR_RNDN); mpfr_cos(r24836, r24833, MPFR_RNDN); mpfr_log1p(r24837, r24836, MPFR_RNDN); mpfr_expm1(r24838, r24837, MPFR_RNDN); mpfr_set_d(r24839, phi1, MPFR_RNDN); mpfr_cos(r24840, r24839, MPFR_RNDN); mpfr_fma(r24841, r24830, r24838, r24840, MPFR_RNDN); mpfr_atan2(r24842, r24835, r24841, MPFR_RNDN); mpfr_add(r24843, r24842, r24831, MPFR_RNDN); return mpfr_get_d(r24843, MPFR_RNDN); }