#include #include #include #include #include char *name = "Midpoint on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r19825 = lambda1; float r19826 = phi2; float r19827 = cos(r19826); float r19828 = lambda2; float r19829 = r19825 - r19828; float r19830 = sin(r19829); float r19831 = r19827 * r19830; float r19832 = phi1; float r19833 = cos(r19832); float r19834 = cos(r19829); float r19835 = r19827 * r19834; float r19836 = r19833 + r19835; float r19837 = atan2(r19831, r19836); float r19838 = r19825 + r19837; return r19838; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r19839 = lambda1; double r19840 = phi2; double r19841 = cos(r19840); double r19842 = lambda2; double r19843 = r19839 - r19842; double r19844 = sin(r19843); double r19845 = r19841 * r19844; double r19846 = phi1; double r19847 = cos(r19846); double r19848 = cos(r19843); double r19849 = r19841 * r19848; double r19850 = r19847 + r19849; double r19851 = atan2(r19845, r19850); double r19852 = r19839 + r19851; return r19852; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r19853 = lambda1; float r19854 = phi2; float r19855 = cos(r19854); float r19856 = sin(r19853); float r19857 = lambda2; float r19858 = cos(r19857); float r19859 = r19856 * r19858; float r19860 = cos(r19853); float r19861 = -r19857; float r19862 = sin(r19861); float r19863 = r19860 * r19862; float r19864 = r19859 + r19863; float r19865 = r19855 * r19864; float r19866 = phi1; float r19867 = cos(r19866); float r19868 = exp(r19867); float r19869 = r19860 * r19858; float r19870 = r19856 * r19862; float r19871 = r19869 - r19870; float r19872 = r19855 * r19871; float r19873 = exp(r19872); float r19874 = r19868 * r19873; float r19875 = log(r19874); float r19876 = atan2(r19865, r19875); float r19877 = r19853 + r19876; return r19877; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r19878 = lambda1; double r19879 = phi2; double r19880 = cos(r19879); double r19881 = sin(r19878); double r19882 = lambda2; double r19883 = cos(r19882); double r19884 = r19881 * r19883; double r19885 = cos(r19878); double r19886 = -r19882; double r19887 = sin(r19886); double r19888 = r19885 * r19887; double r19889 = r19884 + r19888; double r19890 = r19880 * r19889; double r19891 = phi1; double r19892 = cos(r19891); double r19893 = exp(r19892); double r19894 = r19885 * r19883; double r19895 = r19881 * r19887; double r19896 = r19894 - r19895; double r19897 = r19880 * r19896; double r19898 = exp(r19897); double r19899 = r19893 * r19898; double r19900 = log(r19899); double r19901 = atan2(r19890, r19900); double r19902 = r19878 + r19901; return r19902; } 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 r19903, r19904, r19905, r19906, r19907, r19908, r19909, r19910, r19911, r19912, r19913, r19914, r19915, r19916; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r19903); mpfr_init(r19904); mpfr_init(r19905); mpfr_init(r19906); mpfr_init(r19907); mpfr_init(r19908); mpfr_init(r19909); mpfr_init(r19910); mpfr_init(r19911); mpfr_init(r19912); mpfr_init(r19913); mpfr_init(r19914); mpfr_init(r19915); mpfr_init(r19916); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r19903, lambda1, MPFR_RNDN); mpfr_set_d(r19904, phi2, MPFR_RNDN); mpfr_cos(r19905, r19904, MPFR_RNDN); mpfr_set_d(r19906, lambda2, MPFR_RNDN); mpfr_sub(r19907, r19903, r19906, MPFR_RNDN); mpfr_sin(r19908, r19907, MPFR_RNDN); mpfr_mul(r19909, r19905, r19908, MPFR_RNDN); mpfr_set_d(r19910, phi1, MPFR_RNDN); mpfr_cos(r19911, r19910, MPFR_RNDN); mpfr_cos(r19912, r19907, MPFR_RNDN); mpfr_mul(r19913, r19905, r19912, MPFR_RNDN); mpfr_add(r19914, r19911, r19913, MPFR_RNDN); mpfr_atan2(r19915, r19909, r19914, MPFR_RNDN); mpfr_add(r19916, r19903, r19915, MPFR_RNDN); return mpfr_get_d(r19916, MPFR_RNDN); } static mpfr_t r19917, r19918, r19919, r19920, r19921, r19922, r19923, r19924, r19925, r19926, r19927, r19928, r19929, r19930, r19931, r19932, r19933, r19934, r19935, r19936, r19937, r19938, r19939, r19940, r19941; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r19917); mpfr_init(r19918); mpfr_init(r19919); mpfr_init(r19920); mpfr_init(r19921); mpfr_init(r19922); mpfr_init(r19923); mpfr_init(r19924); mpfr_init(r19925); mpfr_init(r19926); mpfr_init(r19927); mpfr_init(r19928); mpfr_init(r19929); mpfr_init(r19930); mpfr_init(r19931); mpfr_init(r19932); mpfr_init(r19933); mpfr_init(r19934); mpfr_init(r19935); mpfr_init(r19936); mpfr_init(r19937); mpfr_init(r19938); mpfr_init(r19939); mpfr_init(r19940); mpfr_init(r19941); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r19917, lambda1, MPFR_RNDN); mpfr_set_d(r19918, phi2, MPFR_RNDN); mpfr_cos(r19919, r19918, MPFR_RNDN); mpfr_sin(r19920, r19917, MPFR_RNDN); mpfr_set_d(r19921, lambda2, MPFR_RNDN); mpfr_cos(r19922, r19921, MPFR_RNDN); mpfr_mul(r19923, r19920, r19922, MPFR_RNDN); mpfr_cos(r19924, r19917, MPFR_RNDN); mpfr_neg(r19925, r19921, MPFR_RNDN); mpfr_sin(r19926, r19925, MPFR_RNDN); mpfr_mul(r19927, r19924, r19926, MPFR_RNDN); mpfr_add(r19928, r19923, r19927, MPFR_RNDN); mpfr_mul(r19929, r19919, r19928, MPFR_RNDN); mpfr_set_d(r19930, phi1, MPFR_RNDN); mpfr_cos(r19931, r19930, MPFR_RNDN); mpfr_exp(r19932, r19931, MPFR_RNDN); mpfr_mul(r19933, r19924, r19922, MPFR_RNDN); mpfr_mul(r19934, r19920, r19926, MPFR_RNDN); mpfr_sub(r19935, r19933, r19934, MPFR_RNDN); mpfr_mul(r19936, r19919, r19935, MPFR_RNDN); mpfr_exp(r19937, r19936, MPFR_RNDN); mpfr_mul(r19938, r19932, r19937, MPFR_RNDN); mpfr_log(r19939, r19938, MPFR_RNDN); mpfr_atan2(r19940, r19929, r19939, MPFR_RNDN); mpfr_add(r19941, r19917, r19940, MPFR_RNDN); return mpfr_get_d(r19941, MPFR_RNDN); } static mpfr_t r19942, r19943, r19944, r19945, r19946, r19947, r19948, r19949, r19950, r19951, r19952, r19953, r19954, r19955, r19956, r19957, r19958, r19959, r19960, r19961, r19962, r19963, r19964, r19965, r19966; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r19942); mpfr_init(r19943); mpfr_init(r19944); mpfr_init(r19945); mpfr_init(r19946); mpfr_init(r19947); mpfr_init(r19948); mpfr_init(r19949); mpfr_init(r19950); mpfr_init(r19951); mpfr_init(r19952); mpfr_init(r19953); mpfr_init(r19954); mpfr_init(r19955); mpfr_init(r19956); mpfr_init(r19957); mpfr_init(r19958); mpfr_init(r19959); mpfr_init(r19960); mpfr_init(r19961); mpfr_init(r19962); mpfr_init(r19963); mpfr_init(r19964); mpfr_init(r19965); mpfr_init(r19966); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r19942, lambda1, MPFR_RNDN); mpfr_set_d(r19943, phi2, MPFR_RNDN); mpfr_cos(r19944, r19943, MPFR_RNDN); mpfr_sin(r19945, r19942, MPFR_RNDN); mpfr_set_d(r19946, lambda2, MPFR_RNDN); mpfr_cos(r19947, r19946, MPFR_RNDN); mpfr_mul(r19948, r19945, r19947, MPFR_RNDN); mpfr_cos(r19949, r19942, MPFR_RNDN); mpfr_neg(r19950, r19946, MPFR_RNDN); mpfr_sin(r19951, r19950, MPFR_RNDN); mpfr_mul(r19952, r19949, r19951, MPFR_RNDN); mpfr_add(r19953, r19948, r19952, MPFR_RNDN); mpfr_mul(r19954, r19944, r19953, MPFR_RNDN); mpfr_set_d(r19955, phi1, MPFR_RNDN); mpfr_cos(r19956, r19955, MPFR_RNDN); mpfr_exp(r19957, r19956, MPFR_RNDN); mpfr_mul(r19958, r19949, r19947, MPFR_RNDN); mpfr_mul(r19959, r19945, r19951, MPFR_RNDN); mpfr_sub(r19960, r19958, r19959, MPFR_RNDN); mpfr_mul(r19961, r19944, r19960, MPFR_RNDN); mpfr_exp(r19962, r19961, MPFR_RNDN); mpfr_mul(r19963, r19957, r19962, MPFR_RNDN); mpfr_log(r19964, r19963, MPFR_RNDN); mpfr_atan2(r19965, r19954, r19964, MPFR_RNDN); mpfr_add(r19966, r19942, r19965, MPFR_RNDN); return mpfr_get_d(r19966, MPFR_RNDN); }