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