#include #include #include #include #include char *name = "Midpoint on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r24758 = lambda1; float r24759 = phi2; float r24760 = cos(r24759); float r24761 = lambda2; float r24762 = r24758 - r24761; float r24763 = sin(r24762); float r24764 = r24760 * r24763; float r24765 = phi1; float r24766 = cos(r24765); float r24767 = cos(r24762); float r24768 = r24760 * r24767; float r24769 = r24766 + r24768; float r24770 = atan2(r24764, r24769); float r24771 = r24758 + r24770; return r24771; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r24772 = lambda1; double r24773 = phi2; double r24774 = cos(r24773); double r24775 = lambda2; double r24776 = r24772 - r24775; double r24777 = sin(r24776); double r24778 = r24774 * r24777; double r24779 = phi1; double r24780 = cos(r24779); double r24781 = cos(r24776); double r24782 = r24774 * r24781; double r24783 = r24780 + r24782; double r24784 = atan2(r24778, r24783); double r24785 = r24772 + r24784; return r24785; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r24786 = lambda1; float r24787 = phi2; float r24788 = cos(r24787); float r24789 = sin(r24786); float r24790 = lambda2; float r24791 = cos(r24790); float r24792 = r24789 * r24791; float r24793 = cos(r24786); float r24794 = -r24790; float r24795 = sin(r24794); float r24796 = r24793 * r24795; float r24797 = r24792 + r24796; float r24798 = r24788 * r24797; float r24799 = phi1; float r24800 = cos(r24799); float r24801 = r24793 * r24791; float r24802 = r24788 * r24801; float r24803 = r24800 + r24802; float r24804 = sin(r24790); float r24805 = r24789 * r24804; float r24806 = r24788 * r24805; float r24807 = r24803 + r24806; float r24808 = atan2(r24798, r24807); float r24809 = r24786 + r24808; return r24809; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r24810 = lambda1; double r24811 = phi2; double r24812 = cos(r24811); double r24813 = sin(r24810); double r24814 = lambda2; double r24815 = cos(r24814); double r24816 = r24813 * r24815; double r24817 = cos(r24810); double r24818 = -r24814; double r24819 = sin(r24818); double r24820 = r24817 * r24819; double r24821 = r24816 + r24820; double r24822 = r24812 * r24821; double r24823 = phi1; double r24824 = cos(r24823); double r24825 = r24817 * r24815; double r24826 = r24812 * r24825; double r24827 = r24824 + r24826; double r24828 = sin(r24814); double r24829 = r24813 * r24828; double r24830 = r24812 * r24829; double r24831 = r24827 + r24830; double r24832 = atan2(r24822, r24831); double r24833 = r24810 + r24832; return r24833; } 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 r24834, r24835, r24836, r24837, r24838, r24839, r24840, r24841, r24842, r24843, r24844, r24845, r24846, r24847; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); 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); mpfr_init(r24844); mpfr_init(r24845); mpfr_init(r24846); mpfr_init(r24847); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24834, lambda1, MPFR_RNDN); mpfr_set_d(r24835, phi2, MPFR_RNDN); mpfr_cos(r24836, r24835, MPFR_RNDN); mpfr_set_d(r24837, lambda2, MPFR_RNDN); mpfr_sub(r24838, r24834, r24837, MPFR_RNDN); mpfr_sin(r24839, r24838, MPFR_RNDN); mpfr_mul(r24840, r24836, r24839, MPFR_RNDN); mpfr_set_d(r24841, phi1, MPFR_RNDN); mpfr_cos(r24842, r24841, MPFR_RNDN); mpfr_cos(r24843, r24838, MPFR_RNDN); mpfr_mul(r24844, r24836, r24843, MPFR_RNDN); mpfr_add(r24845, r24842, r24844, MPFR_RNDN); mpfr_atan2(r24846, r24840, r24845, MPFR_RNDN); mpfr_add(r24847, r24834, r24846, MPFR_RNDN); return mpfr_get_d(r24847, MPFR_RNDN); } static mpfr_t r24848, r24849, r24850, r24851, r24852, r24853, r24854, r24855, r24856, r24857, r24858, r24859, r24860, r24861, r24862, r24863, r24864, r24865, r24866, r24867, r24868, r24869, r24870, r24871; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r24848); mpfr_init(r24849); mpfr_init(r24850); mpfr_init(r24851); mpfr_init(r24852); mpfr_init(r24853); mpfr_init(r24854); mpfr_init(r24855); mpfr_init(r24856); mpfr_init(r24857); mpfr_init(r24858); mpfr_init(r24859); mpfr_init(r24860); mpfr_init(r24861); mpfr_init(r24862); mpfr_init(r24863); mpfr_init(r24864); mpfr_init(r24865); mpfr_init(r24866); mpfr_init(r24867); mpfr_init(r24868); mpfr_init(r24869); mpfr_init(r24870); mpfr_init(r24871); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24848, lambda1, MPFR_RNDN); mpfr_set_d(r24849, phi2, MPFR_RNDN); mpfr_cos(r24850, r24849, MPFR_RNDN); mpfr_sin(r24851, r24848, MPFR_RNDN); mpfr_set_d(r24852, lambda2, MPFR_RNDN); mpfr_cos(r24853, r24852, MPFR_RNDN); mpfr_mul(r24854, r24851, r24853, MPFR_RNDN); mpfr_cos(r24855, r24848, MPFR_RNDN); mpfr_neg(r24856, r24852, MPFR_RNDN); mpfr_sin(r24857, r24856, MPFR_RNDN); mpfr_mul(r24858, r24855, r24857, MPFR_RNDN); mpfr_add(r24859, r24854, r24858, MPFR_RNDN); mpfr_mul(r24860, r24850, r24859, MPFR_RNDN); mpfr_set_d(r24861, phi1, MPFR_RNDN); mpfr_cos(r24862, r24861, MPFR_RNDN); mpfr_mul(r24863, r24855, r24853, MPFR_RNDN); mpfr_mul(r24864, r24850, r24863, MPFR_RNDN); mpfr_add(r24865, r24862, r24864, MPFR_RNDN); mpfr_sin(r24866, r24852, MPFR_RNDN); mpfr_mul(r24867, r24851, r24866, MPFR_RNDN); mpfr_mul(r24868, r24850, r24867, MPFR_RNDN); mpfr_add(r24869, r24865, r24868, MPFR_RNDN); mpfr_atan2(r24870, r24860, r24869, MPFR_RNDN); mpfr_add(r24871, r24848, r24870, MPFR_RNDN); return mpfr_get_d(r24871, MPFR_RNDN); } static mpfr_t r24872, r24873, r24874, r24875, r24876, r24877, r24878, r24879, r24880, r24881, r24882, r24883, r24884, r24885, r24886, r24887, r24888, r24889, r24890, r24891, r24892, r24893, r24894, r24895; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r24872); mpfr_init(r24873); mpfr_init(r24874); mpfr_init(r24875); mpfr_init(r24876); mpfr_init(r24877); mpfr_init(r24878); mpfr_init(r24879); mpfr_init(r24880); mpfr_init(r24881); mpfr_init(r24882); mpfr_init(r24883); mpfr_init(r24884); mpfr_init(r24885); mpfr_init(r24886); mpfr_init(r24887); mpfr_init(r24888); mpfr_init(r24889); mpfr_init(r24890); mpfr_init(r24891); mpfr_init(r24892); mpfr_init(r24893); mpfr_init(r24894); mpfr_init(r24895); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24872, lambda1, MPFR_RNDN); mpfr_set_d(r24873, phi2, MPFR_RNDN); mpfr_cos(r24874, r24873, MPFR_RNDN); mpfr_sin(r24875, r24872, MPFR_RNDN); mpfr_set_d(r24876, lambda2, MPFR_RNDN); mpfr_cos(r24877, r24876, MPFR_RNDN); mpfr_mul(r24878, r24875, r24877, MPFR_RNDN); mpfr_cos(r24879, r24872, MPFR_RNDN); mpfr_neg(r24880, r24876, MPFR_RNDN); mpfr_sin(r24881, r24880, MPFR_RNDN); mpfr_mul(r24882, r24879, r24881, MPFR_RNDN); mpfr_add(r24883, r24878, r24882, MPFR_RNDN); mpfr_mul(r24884, r24874, r24883, MPFR_RNDN); mpfr_set_d(r24885, phi1, MPFR_RNDN); mpfr_cos(r24886, r24885, MPFR_RNDN); mpfr_mul(r24887, r24879, r24877, MPFR_RNDN); mpfr_mul(r24888, r24874, r24887, MPFR_RNDN); mpfr_add(r24889, r24886, r24888, MPFR_RNDN); mpfr_sin(r24890, r24876, MPFR_RNDN); mpfr_mul(r24891, r24875, r24890, MPFR_RNDN); mpfr_mul(r24892, r24874, r24891, MPFR_RNDN); mpfr_add(r24893, r24889, r24892, MPFR_RNDN); mpfr_atan2(r24894, r24884, r24893, MPFR_RNDN); mpfr_add(r24895, r24872, r24894, MPFR_RNDN); return mpfr_get_d(r24895, MPFR_RNDN); }