#include #include #include #include #include char *name = "Bearing on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r23836 = lambda1; float r23837 = lambda2; float r23838 = r23836 - r23837; float r23839 = sin(r23838); float r23840 = phi2; float r23841 = cos(r23840); float r23842 = r23839 * r23841; float r23843 = phi1; float r23844 = cos(r23843); float r23845 = sin(r23840); float r23846 = r23844 * r23845; float r23847 = sin(r23843); float r23848 = r23847 * r23841; float r23849 = cos(r23838); float r23850 = r23848 * r23849; float r23851 = r23846 - r23850; float r23852 = atan2(r23842, r23851); return r23852; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r23853 = lambda1; double r23854 = lambda2; double r23855 = r23853 - r23854; double r23856 = sin(r23855); double r23857 = phi2; double r23858 = cos(r23857); double r23859 = r23856 * r23858; double r23860 = phi1; double r23861 = cos(r23860); double r23862 = sin(r23857); double r23863 = r23861 * r23862; double r23864 = sin(r23860); double r23865 = r23864 * r23858; double r23866 = cos(r23855); double r23867 = r23865 * r23866; double r23868 = r23863 - r23867; double r23869 = atan2(r23859, r23868); return r23869; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r23870 = lambda1; float r23871 = sin(r23870); float r23872 = lambda2; float r23873 = cos(r23872); float r23874 = r23871 * r23873; float r23875 = cos(r23870); float r23876 = -r23872; float r23877 = sin(r23876); float r23878 = r23875 * r23877; float r23879 = r23874 + r23878; float r23880 = phi2; float r23881 = cos(r23880); float r23882 = r23879 * r23881; float r23883 = phi1; float r23884 = cos(r23883); float r23885 = cbrt(r23884); float r23886 = r23885 * r23885; float r23887 = sin(r23880); float r23888 = r23885 * r23887; float r23889 = r23886 * r23888; float r23890 = sin(r23883); float r23891 = r23890 * r23881; float r23892 = r23875 * r23873; float r23893 = r23871 * r23877; float r23894 = r23892 - r23893; float r23895 = r23891 * r23894; float r23896 = r23889 - r23895; float r23897 = atan2(r23882, r23896); return r23897; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r23898 = lambda1; double r23899 = sin(r23898); double r23900 = lambda2; double r23901 = cos(r23900); double r23902 = r23899 * r23901; double r23903 = cos(r23898); double r23904 = -r23900; double r23905 = sin(r23904); double r23906 = r23903 * r23905; double r23907 = r23902 + r23906; double r23908 = phi2; double r23909 = cos(r23908); double r23910 = r23907 * r23909; double r23911 = phi1; double r23912 = cos(r23911); double r23913 = cbrt(r23912); double r23914 = r23913 * r23913; double r23915 = sin(r23908); double r23916 = r23913 * r23915; double r23917 = r23914 * r23916; double r23918 = sin(r23911); double r23919 = r23918 * r23909; double r23920 = r23903 * r23901; double r23921 = r23899 * r23905; double r23922 = r23920 - r23921; double r23923 = r23919 * r23922; double r23924 = r23917 - r23923; double r23925 = atan2(r23910, r23924); return r23925; } 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 r23926, r23927, r23928, r23929, r23930, r23931, r23932, r23933, r23934, r23935, r23936, r23937, r23938, r23939, r23940, r23941, r23942; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r23926); mpfr_init(r23927); mpfr_init(r23928); mpfr_init(r23929); mpfr_init(r23930); mpfr_init(r23931); mpfr_init(r23932); mpfr_init(r23933); mpfr_init(r23934); mpfr_init(r23935); mpfr_init(r23936); mpfr_init(r23937); mpfr_init(r23938); mpfr_init(r23939); mpfr_init(r23940); mpfr_init(r23941); mpfr_init(r23942); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r23926, lambda1, MPFR_RNDN); mpfr_set_d(r23927, lambda2, MPFR_RNDN); mpfr_sub(r23928, r23926, r23927, MPFR_RNDN); mpfr_sin(r23929, r23928, MPFR_RNDN); mpfr_set_d(r23930, phi2, MPFR_RNDN); mpfr_cos(r23931, r23930, MPFR_RNDN); mpfr_mul(r23932, r23929, r23931, MPFR_RNDN); mpfr_set_d(r23933, phi1, MPFR_RNDN); mpfr_cos(r23934, r23933, MPFR_RNDN); mpfr_sin(r23935, r23930, MPFR_RNDN); mpfr_mul(r23936, r23934, r23935, MPFR_RNDN); mpfr_sin(r23937, r23933, MPFR_RNDN); mpfr_mul(r23938, r23937, r23931, MPFR_RNDN); mpfr_cos(r23939, r23928, MPFR_RNDN); mpfr_mul(r23940, r23938, r23939, MPFR_RNDN); mpfr_sub(r23941, r23936, r23940, MPFR_RNDN); mpfr_atan2(r23942, r23932, r23941, MPFR_RNDN); return mpfr_get_d(r23942, MPFR_RNDN); } static mpfr_t r23943, r23944, r23945, r23946, r23947, r23948, r23949, r23950, r23951, r23952, r23953, r23954, r23955, r23956, r23957, r23958, r23959, r23960, r23961, r23962, r23963, r23964, r23965, r23966, r23967, r23968, r23969, r23970; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r23943); mpfr_init(r23944); mpfr_init(r23945); mpfr_init(r23946); mpfr_init(r23947); mpfr_init(r23948); mpfr_init(r23949); mpfr_init(r23950); mpfr_init(r23951); mpfr_init(r23952); mpfr_init(r23953); mpfr_init(r23954); mpfr_init(r23955); mpfr_init(r23956); mpfr_init(r23957); mpfr_init(r23958); mpfr_init(r23959); mpfr_init(r23960); mpfr_init(r23961); mpfr_init(r23962); mpfr_init(r23963); mpfr_init(r23964); mpfr_init(r23965); mpfr_init(r23966); mpfr_init(r23967); mpfr_init(r23968); mpfr_init(r23969); mpfr_init(r23970); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r23943, lambda1, MPFR_RNDN); mpfr_sin(r23944, r23943, MPFR_RNDN); mpfr_set_d(r23945, lambda2, MPFR_RNDN); mpfr_cos(r23946, r23945, MPFR_RNDN); mpfr_mul(r23947, r23944, r23946, MPFR_RNDN); mpfr_cos(r23948, r23943, MPFR_RNDN); mpfr_neg(r23949, r23945, MPFR_RNDN); mpfr_sin(r23950, r23949, MPFR_RNDN); mpfr_mul(r23951, r23948, r23950, MPFR_RNDN); mpfr_add(r23952, r23947, r23951, MPFR_RNDN); mpfr_set_d(r23953, phi2, MPFR_RNDN); mpfr_cos(r23954, r23953, MPFR_RNDN); mpfr_mul(r23955, r23952, r23954, MPFR_RNDN); mpfr_set_d(r23956, phi1, MPFR_RNDN); mpfr_cos(r23957, r23956, MPFR_RNDN); mpfr_cbrt(r23958, r23957, MPFR_RNDN); mpfr_mul(r23959, r23958, r23958, MPFR_RNDN); mpfr_sin(r23960, r23953, MPFR_RNDN); mpfr_mul(r23961, r23958, r23960, MPFR_RNDN); mpfr_mul(r23962, r23959, r23961, MPFR_RNDN); mpfr_sin(r23963, r23956, MPFR_RNDN); mpfr_mul(r23964, r23963, r23954, MPFR_RNDN); mpfr_mul(r23965, r23948, r23946, MPFR_RNDN); mpfr_mul(r23966, r23944, r23950, MPFR_RNDN); mpfr_sub(r23967, r23965, r23966, MPFR_RNDN); mpfr_mul(r23968, r23964, r23967, MPFR_RNDN); mpfr_sub(r23969, r23962, r23968, MPFR_RNDN); mpfr_atan2(r23970, r23955, r23969, MPFR_RNDN); return mpfr_get_d(r23970, MPFR_RNDN); } static mpfr_t r23971, r23972, r23973, r23974, r23975, r23976, r23977, r23978, r23979, r23980, r23981, r23982, r23983, r23984, r23985, r23986, r23987, r23988, r23989, r23990, r23991, r23992, r23993, r23994, r23995, r23996, r23997, r23998; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r23971); mpfr_init(r23972); mpfr_init(r23973); mpfr_init(r23974); mpfr_init(r23975); mpfr_init(r23976); mpfr_init(r23977); mpfr_init(r23978); mpfr_init(r23979); mpfr_init(r23980); mpfr_init(r23981); mpfr_init(r23982); mpfr_init(r23983); mpfr_init(r23984); mpfr_init(r23985); mpfr_init(r23986); mpfr_init(r23987); mpfr_init(r23988); mpfr_init(r23989); mpfr_init(r23990); mpfr_init(r23991); mpfr_init(r23992); mpfr_init(r23993); mpfr_init(r23994); mpfr_init(r23995); mpfr_init(r23996); mpfr_init(r23997); mpfr_init(r23998); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r23971, lambda1, MPFR_RNDN); mpfr_sin(r23972, r23971, MPFR_RNDN); mpfr_set_d(r23973, lambda2, MPFR_RNDN); mpfr_cos(r23974, r23973, MPFR_RNDN); mpfr_mul(r23975, r23972, r23974, MPFR_RNDN); mpfr_cos(r23976, r23971, MPFR_RNDN); mpfr_neg(r23977, r23973, MPFR_RNDN); mpfr_sin(r23978, r23977, MPFR_RNDN); mpfr_mul(r23979, r23976, r23978, MPFR_RNDN); mpfr_add(r23980, r23975, r23979, MPFR_RNDN); mpfr_set_d(r23981, phi2, MPFR_RNDN); mpfr_cos(r23982, r23981, MPFR_RNDN); mpfr_mul(r23983, r23980, r23982, MPFR_RNDN); mpfr_set_d(r23984, phi1, MPFR_RNDN); mpfr_cos(r23985, r23984, MPFR_RNDN); mpfr_cbrt(r23986, r23985, MPFR_RNDN); mpfr_mul(r23987, r23986, r23986, MPFR_RNDN); mpfr_sin(r23988, r23981, MPFR_RNDN); mpfr_mul(r23989, r23986, r23988, MPFR_RNDN); mpfr_mul(r23990, r23987, r23989, MPFR_RNDN); mpfr_sin(r23991, r23984, MPFR_RNDN); mpfr_mul(r23992, r23991, r23982, MPFR_RNDN); mpfr_mul(r23993, r23976, r23974, MPFR_RNDN); mpfr_mul(r23994, r23972, r23978, MPFR_RNDN); mpfr_sub(r23995, r23993, r23994, MPFR_RNDN); mpfr_mul(r23996, r23992, r23995, MPFR_RNDN); mpfr_sub(r23997, r23990, r23996, MPFR_RNDN); mpfr_atan2(r23998, r23983, r23997, MPFR_RNDN); return mpfr_get_d(r23998, MPFR_RNDN); }