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