#include #include #include #include #include char *name = "Midpoint on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r24777 = lambda1; float r24778 = phi2; float r24779 = cos(r24778); float r24780 = lambda2; float r24781 = r24777 - r24780; float r24782 = sin(r24781); float r24783 = r24779 * r24782; float r24784 = phi1; float r24785 = cos(r24784); float r24786 = cos(r24781); float r24787 = r24779 * r24786; float r24788 = r24785 + r24787; float r24789 = atan2(r24783, r24788); float r24790 = r24777 + r24789; return r24790; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r24791 = lambda1; double r24792 = phi2; double r24793 = cos(r24792); double r24794 = lambda2; double r24795 = r24791 - r24794; double r24796 = sin(r24795); double r24797 = r24793 * r24796; double r24798 = phi1; double r24799 = cos(r24798); double r24800 = cos(r24795); double r24801 = r24793 * r24800; double r24802 = r24799 + r24801; double r24803 = atan2(r24797, r24802); double r24804 = r24791 + r24803; return r24804; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r24805 = lambda1; float r24806 = sin(r24805); float r24807 = lambda2; float r24808 = cos(r24807); float r24809 = r24806 * r24808; float r24810 = cos(r24805); float r24811 = sin(r24807); float r24812 = r24810 * r24811; float r24813 = r24809 - r24812; float r24814 = phi2; float r24815 = cos(r24814); float r24816 = r24813 * r24815; float r24817 = expm1(r24816); float r24818 = log1p(r24817); float r24819 = phi1; float r24820 = cos(r24819); float r24821 = r24810 * r24808; float r24822 = -r24807; float r24823 = sin(r24822); float r24824 = r24806 * r24823; float r24825 = r24821 - r24824; float r24826 = r24815 * r24825; float r24827 = r24820 + r24826; float r24828 = atan2(r24818, r24827); float r24829 = r24805 + r24828; return r24829; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r24830 = lambda1; double r24831 = sin(r24830); double r24832 = lambda2; double r24833 = cos(r24832); double r24834 = r24831 * r24833; double r24835 = cos(r24830); double r24836 = sin(r24832); double r24837 = r24835 * r24836; double r24838 = r24834 - r24837; double r24839 = phi2; double r24840 = cos(r24839); double r24841 = r24838 * r24840; double r24842 = expm1(r24841); double r24843 = log1p(r24842); double r24844 = phi1; double r24845 = cos(r24844); double r24846 = r24835 * r24833; double r24847 = -r24832; double r24848 = sin(r24847); double r24849 = r24831 * r24848; double r24850 = r24846 - r24849; double r24851 = r24840 * r24850; double r24852 = r24845 + r24851; double r24853 = atan2(r24843, r24852); double r24854 = r24830 + r24853; return r24854; } 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 r24855, r24856, r24857, r24858, r24859, r24860, r24861, r24862, r24863, r24864, r24865, r24866, r24867, r24868; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); 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); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24855, lambda1, MPFR_RNDN); mpfr_set_d(r24856, phi2, MPFR_RNDN); mpfr_cos(r24857, r24856, MPFR_RNDN); mpfr_set_d(r24858, lambda2, MPFR_RNDN); mpfr_sub(r24859, r24855, r24858, MPFR_RNDN); mpfr_sin(r24860, r24859, MPFR_RNDN); mpfr_mul(r24861, r24857, r24860, MPFR_RNDN); mpfr_set_d(r24862, phi1, MPFR_RNDN); mpfr_cos(r24863, r24862, MPFR_RNDN); mpfr_cos(r24864, r24859, MPFR_RNDN); mpfr_mul(r24865, r24857, r24864, MPFR_RNDN); mpfr_add(r24866, r24863, r24865, MPFR_RNDN); mpfr_atan2(r24867, r24861, r24866, MPFR_RNDN); mpfr_add(r24868, r24855, r24867, MPFR_RNDN); return mpfr_get_d(r24868, MPFR_RNDN); } static mpfr_t r24869, r24870, r24871, r24872, r24873, r24874, r24875, r24876, r24877, r24878, r24879, r24880, r24881, r24882, r24883, r24884, r24885, r24886, r24887, r24888, r24889, r24890, r24891, r24892, r24893; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r24869); mpfr_init(r24870); mpfr_init(r24871); 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); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24869, lambda1, MPFR_RNDN); mpfr_sin(r24870, r24869, MPFR_RNDN); mpfr_set_d(r24871, lambda2, MPFR_RNDN); mpfr_cos(r24872, r24871, MPFR_RNDN); mpfr_mul(r24873, r24870, r24872, MPFR_RNDN); mpfr_cos(r24874, r24869, MPFR_RNDN); mpfr_sin(r24875, r24871, MPFR_RNDN); mpfr_mul(r24876, r24874, r24875, MPFR_RNDN); mpfr_sub(r24877, r24873, r24876, MPFR_RNDN); mpfr_set_d(r24878, phi2, MPFR_RNDN); mpfr_cos(r24879, r24878, MPFR_RNDN); mpfr_mul(r24880, r24877, r24879, MPFR_RNDN); mpfr_expm1(r24881, r24880, MPFR_RNDN); mpfr_log1p(r24882, r24881, MPFR_RNDN); mpfr_set_d(r24883, phi1, MPFR_RNDN); mpfr_cos(r24884, r24883, MPFR_RNDN); mpfr_mul(r24885, r24874, r24872, MPFR_RNDN); mpfr_neg(r24886, r24871, MPFR_RNDN); mpfr_sin(r24887, r24886, MPFR_RNDN); mpfr_mul(r24888, r24870, r24887, MPFR_RNDN); mpfr_sub(r24889, r24885, r24888, MPFR_RNDN); mpfr_mul(r24890, r24879, r24889, MPFR_RNDN); mpfr_add(r24891, r24884, r24890, MPFR_RNDN); mpfr_atan2(r24892, r24882, r24891, MPFR_RNDN); mpfr_add(r24893, r24869, r24892, MPFR_RNDN); return mpfr_get_d(r24893, MPFR_RNDN); } static mpfr_t r24894, r24895, r24896, r24897, r24898, r24899, r24900, r24901, r24902, r24903, r24904, r24905, r24906, r24907, r24908, r24909, r24910, r24911, r24912, r24913, r24914, r24915, r24916, r24917, r24918; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r24894); mpfr_init(r24895); mpfr_init(r24896); mpfr_init(r24897); mpfr_init(r24898); mpfr_init(r24899); mpfr_init(r24900); mpfr_init(r24901); mpfr_init(r24902); mpfr_init(r24903); mpfr_init(r24904); mpfr_init(r24905); mpfr_init(r24906); mpfr_init(r24907); mpfr_init(r24908); mpfr_init(r24909); mpfr_init(r24910); mpfr_init(r24911); mpfr_init(r24912); mpfr_init(r24913); mpfr_init(r24914); mpfr_init(r24915); mpfr_init(r24916); mpfr_init(r24917); mpfr_init(r24918); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24894, lambda1, MPFR_RNDN); mpfr_sin(r24895, r24894, MPFR_RNDN); mpfr_set_d(r24896, lambda2, MPFR_RNDN); mpfr_cos(r24897, r24896, MPFR_RNDN); mpfr_mul(r24898, r24895, r24897, MPFR_RNDN); mpfr_cos(r24899, r24894, MPFR_RNDN); mpfr_sin(r24900, r24896, MPFR_RNDN); mpfr_mul(r24901, r24899, r24900, MPFR_RNDN); mpfr_sub(r24902, r24898, r24901, MPFR_RNDN); mpfr_set_d(r24903, phi2, MPFR_RNDN); mpfr_cos(r24904, r24903, MPFR_RNDN); mpfr_mul(r24905, r24902, r24904, MPFR_RNDN); mpfr_expm1(r24906, r24905, MPFR_RNDN); mpfr_log1p(r24907, r24906, MPFR_RNDN); mpfr_set_d(r24908, phi1, MPFR_RNDN); mpfr_cos(r24909, r24908, MPFR_RNDN); mpfr_mul(r24910, r24899, r24897, MPFR_RNDN); mpfr_neg(r24911, r24896, MPFR_RNDN); mpfr_sin(r24912, r24911, MPFR_RNDN); mpfr_mul(r24913, r24895, r24912, MPFR_RNDN); mpfr_sub(r24914, r24910, r24913, MPFR_RNDN); mpfr_mul(r24915, r24904, r24914, MPFR_RNDN); mpfr_add(r24916, r24909, r24915, MPFR_RNDN); mpfr_atan2(r24917, r24907, r24916, MPFR_RNDN); mpfr_add(r24918, r24894, r24917, MPFR_RNDN); return mpfr_get_d(r24918, MPFR_RNDN); }