#include #include #include #include #include char *name = "Midpoint on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r24877 = lambda1; float r24878 = phi2; float r24879 = cos(r24878); float r24880 = lambda2; float r24881 = r24877 - r24880; float r24882 = sin(r24881); float r24883 = r24879 * r24882; float r24884 = phi1; float r24885 = cos(r24884); float r24886 = cos(r24881); float r24887 = r24879 * r24886; float r24888 = r24885 + r24887; float r24889 = atan2(r24883, r24888); float r24890 = r24877 + r24889; return r24890; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r24891 = lambda1; double r24892 = phi2; double r24893 = cos(r24892); double r24894 = lambda2; double r24895 = r24891 - r24894; double r24896 = sin(r24895); double r24897 = r24893 * r24896; double r24898 = phi1; double r24899 = cos(r24898); double r24900 = cos(r24895); double r24901 = r24893 * r24900; double r24902 = r24899 + r24901; double r24903 = atan2(r24897, r24902); double r24904 = r24891 + r24903; return r24904; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r24905 = lambda1; float r24906 = phi2; float r24907 = cos(r24906); float r24908 = sin(r24905); float r24909 = lambda2; float r24910 = cos(r24909); float r24911 = r24908 * r24910; float r24912 = cos(r24905); float r24913 = -r24909; float r24914 = sin(r24913); float r24915 = r24912 * r24914; float r24916 = r24911 + r24915; float r24917 = r24907 * r24916; float r24918 = r24912 * r24910; float r24919 = phi1; float r24920 = cos(r24919); float r24921 = fma(r24907, r24918, r24920); float r24922 = sin(r24909); float r24923 = r24908 * r24922; float r24924 = r24907 * r24923; float r24925 = r24921 + r24924; float r24926 = atan2(r24917, r24925); float r24927 = r24905 + r24926; return r24927; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r24928 = lambda1; double r24929 = phi2; double r24930 = cos(r24929); double r24931 = sin(r24928); double r24932 = lambda2; double r24933 = cos(r24932); double r24934 = r24931 * r24933; double r24935 = cos(r24928); double r24936 = -r24932; double r24937 = sin(r24936); double r24938 = r24935 * r24937; double r24939 = r24934 + r24938; double r24940 = r24930 * r24939; double r24941 = r24935 * r24933; double r24942 = phi1; double r24943 = cos(r24942); double r24944 = fma(r24930, r24941, r24943); double r24945 = sin(r24932); double r24946 = r24931 * r24945; double r24947 = r24930 * r24946; double r24948 = r24944 + r24947; double r24949 = atan2(r24940, r24948); double r24950 = r24928 + r24949; return r24950; } 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 r24951, r24952, r24953, r24954, r24955, r24956, r24957, r24958, r24959, r24960, r24961, r24962, r24963, r24964; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r24951); mpfr_init(r24952); mpfr_init(r24953); mpfr_init(r24954); mpfr_init(r24955); mpfr_init(r24956); mpfr_init(r24957); mpfr_init(r24958); mpfr_init(r24959); mpfr_init(r24960); mpfr_init(r24961); mpfr_init(r24962); mpfr_init(r24963); mpfr_init(r24964); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24951, lambda1, MPFR_RNDN); mpfr_set_d(r24952, phi2, MPFR_RNDN); mpfr_cos(r24953, r24952, MPFR_RNDN); mpfr_set_d(r24954, lambda2, MPFR_RNDN); mpfr_sub(r24955, r24951, r24954, MPFR_RNDN); mpfr_sin(r24956, r24955, MPFR_RNDN); mpfr_mul(r24957, r24953, r24956, MPFR_RNDN); mpfr_set_d(r24958, phi1, MPFR_RNDN); mpfr_cos(r24959, r24958, MPFR_RNDN); mpfr_cos(r24960, r24955, MPFR_RNDN); mpfr_mul(r24961, r24953, r24960, MPFR_RNDN); mpfr_add(r24962, r24959, r24961, MPFR_RNDN); mpfr_atan2(r24963, r24957, r24962, MPFR_RNDN); mpfr_add(r24964, r24951, r24963, MPFR_RNDN); return mpfr_get_d(r24964, MPFR_RNDN); } static mpfr_t r24965, r24966, r24967, r24968, r24969, r24970, r24971, r24972, r24973, r24974, r24975, r24976, r24977, r24978, r24979, r24980, r24981, r24982, r24983, r24984, r24985, r24986, r24987; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r24965); mpfr_init(r24966); mpfr_init(r24967); mpfr_init(r24968); mpfr_init(r24969); mpfr_init(r24970); mpfr_init(r24971); mpfr_init(r24972); mpfr_init(r24973); mpfr_init(r24974); mpfr_init(r24975); mpfr_init(r24976); mpfr_init(r24977); mpfr_init(r24978); mpfr_init(r24979); mpfr_init(r24980); mpfr_init(r24981); mpfr_init(r24982); mpfr_init(r24983); mpfr_init(r24984); mpfr_init(r24985); mpfr_init(r24986); mpfr_init(r24987); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24965, lambda1, MPFR_RNDN); mpfr_set_d(r24966, phi2, MPFR_RNDN); mpfr_cos(r24967, r24966, MPFR_RNDN); mpfr_sin(r24968, r24965, MPFR_RNDN); mpfr_set_d(r24969, lambda2, MPFR_RNDN); mpfr_cos(r24970, r24969, MPFR_RNDN); mpfr_mul(r24971, r24968, r24970, MPFR_RNDN); mpfr_cos(r24972, r24965, MPFR_RNDN); mpfr_neg(r24973, r24969, MPFR_RNDN); mpfr_sin(r24974, r24973, MPFR_RNDN); mpfr_mul(r24975, r24972, r24974, MPFR_RNDN); mpfr_add(r24976, r24971, r24975, MPFR_RNDN); mpfr_mul(r24977, r24967, r24976, MPFR_RNDN); mpfr_mul(r24978, r24972, r24970, MPFR_RNDN); mpfr_set_d(r24979, phi1, MPFR_RNDN); mpfr_cos(r24980, r24979, MPFR_RNDN); mpfr_fma(r24981, r24967, r24978, r24980, MPFR_RNDN); mpfr_sin(r24982, r24969, MPFR_RNDN); mpfr_mul(r24983, r24968, r24982, MPFR_RNDN); mpfr_mul(r24984, r24967, r24983, MPFR_RNDN); mpfr_add(r24985, r24981, r24984, MPFR_RNDN); mpfr_atan2(r24986, r24977, r24985, MPFR_RNDN); mpfr_add(r24987, r24965, r24986, MPFR_RNDN); return mpfr_get_d(r24987, MPFR_RNDN); } static mpfr_t r24988, r24989, r24990, r24991, r24992, r24993, r24994, r24995, r24996, r24997, r24998, r24999, r25000, r25001, r25002, r25003, r25004, r25005, r25006, r25007, r25008, r25009, r25010; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r24988); mpfr_init(r24989); mpfr_init(r24990); mpfr_init(r24991); mpfr_init(r24992); mpfr_init(r24993); mpfr_init(r24994); mpfr_init(r24995); mpfr_init(r24996); mpfr_init(r24997); mpfr_init(r24998); mpfr_init(r24999); mpfr_init(r25000); mpfr_init(r25001); mpfr_init(r25002); mpfr_init(r25003); mpfr_init(r25004); mpfr_init(r25005); mpfr_init(r25006); mpfr_init(r25007); mpfr_init(r25008); mpfr_init(r25009); mpfr_init(r25010); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24988, lambda1, MPFR_RNDN); mpfr_set_d(r24989, phi2, MPFR_RNDN); mpfr_cos(r24990, r24989, MPFR_RNDN); mpfr_sin(r24991, r24988, MPFR_RNDN); mpfr_set_d(r24992, lambda2, MPFR_RNDN); mpfr_cos(r24993, r24992, MPFR_RNDN); mpfr_mul(r24994, r24991, r24993, MPFR_RNDN); mpfr_cos(r24995, r24988, MPFR_RNDN); mpfr_neg(r24996, r24992, MPFR_RNDN); mpfr_sin(r24997, r24996, MPFR_RNDN); mpfr_mul(r24998, r24995, r24997, MPFR_RNDN); mpfr_add(r24999, r24994, r24998, MPFR_RNDN); mpfr_mul(r25000, r24990, r24999, MPFR_RNDN); mpfr_mul(r25001, r24995, r24993, MPFR_RNDN); mpfr_set_d(r25002, phi1, MPFR_RNDN); mpfr_cos(r25003, r25002, MPFR_RNDN); mpfr_fma(r25004, r24990, r25001, r25003, MPFR_RNDN); mpfr_sin(r25005, r24992, MPFR_RNDN); mpfr_mul(r25006, r24991, r25005, MPFR_RNDN); mpfr_mul(r25007, r24990, r25006, MPFR_RNDN); mpfr_add(r25008, r25004, r25007, MPFR_RNDN); mpfr_atan2(r25009, r25000, r25008, MPFR_RNDN); mpfr_add(r25010, r24988, r25009, MPFR_RNDN); return mpfr_get_d(r25010, MPFR_RNDN); }