#include #include #include #include #include char *name = "Destination given bearing on a great circle"; double f_if(float lambda1, float phi1, float __attribute__((unused)) phi2, float delta, float theta) { float r25901 = lambda1; float r25902 = theta; float r25903 = sin(r25902); float r25904 = delta; float r25905 = sin(r25904); float r25906 = r25903 * r25905; float r25907 = phi1; float r25908 = cos(r25907); float r25909 = r25906 * r25908; float r25910 = cos(r25904); float r25911 = sin(r25907); float r25912 = r25911 * r25910; float r25913 = r25908 * r25905; float r25914 = cos(r25902); float r25915 = r25913 * r25914; float r25916 = r25912 + r25915; float r25917 = asin(r25916); float r25918 = sin(r25917); float r25919 = r25911 * r25918; float r25920 = r25910 - r25919; float r25921 = atan2(r25909, r25920); float r25922 = r25901 + r25921; return r25922; } double f_id(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { double r25923 = lambda1; double r25924 = theta; double r25925 = sin(r25924); double r25926 = delta; double r25927 = sin(r25926); double r25928 = r25925 * r25927; double r25929 = phi1; double r25930 = cos(r25929); double r25931 = r25928 * r25930; double r25932 = cos(r25926); double r25933 = sin(r25929); double r25934 = r25933 * r25932; double r25935 = r25930 * r25927; double r25936 = cos(r25924); double r25937 = r25935 * r25936; double r25938 = r25934 + r25937; double r25939 = asin(r25938); double r25940 = sin(r25939); double r25941 = r25933 * r25940; double r25942 = r25932 - r25941; double r25943 = atan2(r25931, r25942); double r25944 = r25923 + r25943; return r25944; } double f_of(float lambda1, float phi1, float __attribute__((unused)) phi2, float delta, float theta) { float r25945 = lambda1; float r25946 = phi1; float r25947 = cos(r25946); float r25948 = delta; float r25949 = sin(r25948); float r25950 = r25947 * r25949; float r25951 = theta; float r25952 = sin(r25951); float r25953 = r25950 * r25952; float r25954 = cos(r25948); float r25955 = sin(r25946); float r25956 = r25954 * r25955; float r25957 = r25955 * r25956; float r25958 = r25954 - r25957; float r25959 = cos(r25951); float r25960 = r25959 * r25955; float r25961 = r25960 * r25950; float r25962 = r25958 - r25961; float r25963 = atan2(r25953, r25962); float r25964 = r25945 + r25963; return r25964; } double f_od(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { double r25965 = lambda1; double r25966 = phi1; double r25967 = cos(r25966); double r25968 = delta; double r25969 = sin(r25968); double r25970 = r25967 * r25969; double r25971 = theta; double r25972 = sin(r25971); double r25973 = r25970 * r25972; double r25974 = cos(r25968); double r25975 = sin(r25966); double r25976 = r25974 * r25975; double r25977 = r25975 * r25976; double r25978 = r25974 - r25977; double r25979 = cos(r25971); double r25980 = r25979 * r25975; double r25981 = r25980 * r25970; double r25982 = r25978 - r25981; double r25983 = atan2(r25973, r25982); double r25984 = r25965 + r25983; return r25984; } 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 r25985, r25986, r25987, r25988, r25989, r25990, r25991, r25992, r25993, r25994, r25995, r25996, r25997, r25998, r25999, r26000, r26001, r26002, r26003, r26004, r26005, r26006; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25985); mpfr_init(r25986); mpfr_init(r25987); mpfr_init(r25988); mpfr_init(r25989); mpfr_init(r25990); mpfr_init(r25991); mpfr_init(r25992); mpfr_init(r25993); mpfr_init(r25994); mpfr_init(r25995); mpfr_init(r25996); mpfr_init(r25997); mpfr_init(r25998); mpfr_init(r25999); mpfr_init(r26000); mpfr_init(r26001); mpfr_init(r26002); mpfr_init(r26003); mpfr_init(r26004); mpfr_init(r26005); mpfr_init(r26006); } double f_im(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { mpfr_set_d(r25985, lambda1, MPFR_RNDN); mpfr_set_d(r25986, theta, MPFR_RNDN); mpfr_sin(r25987, r25986, MPFR_RNDN); mpfr_set_d(r25988, delta, MPFR_RNDN); mpfr_sin(r25989, r25988, MPFR_RNDN); mpfr_mul(r25990, r25987, r25989, MPFR_RNDN); mpfr_set_d(r25991, phi1, MPFR_RNDN); mpfr_cos(r25992, r25991, MPFR_RNDN); mpfr_mul(r25993, r25990, r25992, MPFR_RNDN); mpfr_cos(r25994, r25988, MPFR_RNDN); mpfr_sin(r25995, r25991, MPFR_RNDN); mpfr_mul(r25996, r25995, r25994, MPFR_RNDN); mpfr_mul(r25997, r25992, r25989, MPFR_RNDN); mpfr_cos(r25998, r25986, MPFR_RNDN); mpfr_mul(r25999, r25997, r25998, MPFR_RNDN); mpfr_add(r26000, r25996, r25999, MPFR_RNDN); mpfr_asin(r26001, r26000, MPFR_RNDN); mpfr_sin(r26002, r26001, MPFR_RNDN); mpfr_mul(r26003, r25995, r26002, MPFR_RNDN); mpfr_sub(r26004, r25994, r26003, MPFR_RNDN); mpfr_atan2(r26005, r25993, r26004, MPFR_RNDN); mpfr_add(r26006, r25985, r26005, MPFR_RNDN); return mpfr_get_d(r26006, MPFR_RNDN); } static mpfr_t r26007, r26008, r26009, r26010, r26011, r26012, r26013, r26014, r26015, r26016, r26017, r26018, r26019, r26020, r26021, r26022, r26023, r26024, r26025, r26026; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r26007); mpfr_init(r26008); mpfr_init(r26009); mpfr_init(r26010); mpfr_init(r26011); mpfr_init(r26012); mpfr_init(r26013); mpfr_init(r26014); mpfr_init(r26015); mpfr_init(r26016); mpfr_init(r26017); mpfr_init(r26018); mpfr_init(r26019); mpfr_init(r26020); mpfr_init(r26021); mpfr_init(r26022); mpfr_init(r26023); mpfr_init(r26024); mpfr_init(r26025); mpfr_init(r26026); } double f_fm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { mpfr_set_d(r26007, lambda1, MPFR_RNDN); mpfr_set_d(r26008, phi1, MPFR_RNDN); mpfr_cos(r26009, r26008, MPFR_RNDN); mpfr_set_d(r26010, delta, MPFR_RNDN); mpfr_sin(r26011, r26010, MPFR_RNDN); mpfr_mul(r26012, r26009, r26011, MPFR_RNDN); mpfr_set_d(r26013, theta, MPFR_RNDN); mpfr_sin(r26014, r26013, MPFR_RNDN); mpfr_mul(r26015, r26012, r26014, MPFR_RNDN); mpfr_cos(r26016, r26010, MPFR_RNDN); mpfr_sin(r26017, r26008, MPFR_RNDN); mpfr_mul(r26018, r26016, r26017, MPFR_RNDN); mpfr_mul(r26019, r26017, r26018, MPFR_RNDN); mpfr_sub(r26020, r26016, r26019, MPFR_RNDN); mpfr_cos(r26021, r26013, MPFR_RNDN); mpfr_mul(r26022, r26021, r26017, MPFR_RNDN); mpfr_mul(r26023, r26022, r26012, MPFR_RNDN); mpfr_sub(r26024, r26020, r26023, MPFR_RNDN); mpfr_atan2(r26025, r26015, r26024, MPFR_RNDN); mpfr_add(r26026, r26007, r26025, MPFR_RNDN); return mpfr_get_d(r26026, MPFR_RNDN); } static mpfr_t r26027, r26028, r26029, r26030, r26031, r26032, r26033, r26034, r26035, r26036, r26037, r26038, r26039, r26040, r26041, r26042, r26043, r26044, r26045, r26046; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r26027); mpfr_init(r26028); mpfr_init(r26029); mpfr_init(r26030); mpfr_init(r26031); mpfr_init(r26032); mpfr_init(r26033); mpfr_init(r26034); mpfr_init(r26035); mpfr_init(r26036); mpfr_init(r26037); mpfr_init(r26038); mpfr_init(r26039); mpfr_init(r26040); mpfr_init(r26041); mpfr_init(r26042); mpfr_init(r26043); mpfr_init(r26044); mpfr_init(r26045); mpfr_init(r26046); } double f_dm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { mpfr_set_d(r26027, lambda1, MPFR_RNDN); mpfr_set_d(r26028, phi1, MPFR_RNDN); mpfr_cos(r26029, r26028, MPFR_RNDN); mpfr_set_d(r26030, delta, MPFR_RNDN); mpfr_sin(r26031, r26030, MPFR_RNDN); mpfr_mul(r26032, r26029, r26031, MPFR_RNDN); mpfr_set_d(r26033, theta, MPFR_RNDN); mpfr_sin(r26034, r26033, MPFR_RNDN); mpfr_mul(r26035, r26032, r26034, MPFR_RNDN); mpfr_cos(r26036, r26030, MPFR_RNDN); mpfr_sin(r26037, r26028, MPFR_RNDN); mpfr_mul(r26038, r26036, r26037, MPFR_RNDN); mpfr_mul(r26039, r26037, r26038, MPFR_RNDN); mpfr_sub(r26040, r26036, r26039, MPFR_RNDN); mpfr_cos(r26041, r26033, MPFR_RNDN); mpfr_mul(r26042, r26041, r26037, MPFR_RNDN); mpfr_mul(r26043, r26042, r26032, MPFR_RNDN); mpfr_sub(r26044, r26040, r26043, MPFR_RNDN); mpfr_atan2(r26045, r26035, r26044, MPFR_RNDN); mpfr_add(r26046, r26027, r26045, MPFR_RNDN); return mpfr_get_d(r26046, MPFR_RNDN); }