#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 r25936 = lambda1; float r25937 = theta; float r25938 = sin(r25937); float r25939 = delta; float r25940 = sin(r25939); float r25941 = r25938 * r25940; float r25942 = phi1; float r25943 = cos(r25942); float r25944 = r25941 * r25943; float r25945 = cos(r25939); float r25946 = sin(r25942); float r25947 = r25946 * r25945; float r25948 = r25943 * r25940; float r25949 = cos(r25937); float r25950 = r25948 * r25949; float r25951 = r25947 + r25950; float r25952 = asin(r25951); float r25953 = sin(r25952); float r25954 = r25946 * r25953; float r25955 = r25945 - r25954; float r25956 = atan2(r25944, r25955); float r25957 = r25936 + r25956; return r25957; } double f_id(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { double r25958 = lambda1; double r25959 = theta; double r25960 = sin(r25959); double r25961 = delta; double r25962 = sin(r25961); double r25963 = r25960 * r25962; double r25964 = phi1; double r25965 = cos(r25964); double r25966 = r25963 * r25965; double r25967 = cos(r25961); double r25968 = sin(r25964); double r25969 = r25968 * r25967; double r25970 = r25965 * r25962; double r25971 = cos(r25959); double r25972 = r25970 * r25971; double r25973 = r25969 + r25972; double r25974 = asin(r25973); double r25975 = sin(r25974); double r25976 = r25968 * r25975; double r25977 = r25967 - r25976; double r25978 = atan2(r25966, r25977); double r25979 = r25958 + r25978; return r25979; } double f_of(float lambda1, float phi1, float __attribute__((unused)) phi2, float delta, float theta) { float r25980 = delta; float r25981 = sin(r25980); float r25982 = phi1; float r25983 = cos(r25982); float r25984 = theta; float r25985 = sin(r25984); float r25986 = r25983 * r25985; float r25987 = r25981 * r25986; float r25988 = sin(r25982); float r25989 = -r25988; float r25990 = cos(r25984); float r25991 = r25983 * r25981; float r25992 = cos(r25980); float r25993 = r25992 * r25988; float r25994 = fma(r25990, r25991, r25993); float r25995 = 3; float r25996 = pow(r25994, r25995); float r25997 = cbrt(r25996); float r25998 = asin(r25997); float r25999 = sin(r25998); float r26000 = fma(r25989, r25999, r25992); float r26001 = atan2(r25987, r26000); float r26002 = lambda1; float r26003 = r26001 + r26002; return r26003; } double f_od(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { double r26004 = delta; double r26005 = sin(r26004); double r26006 = phi1; double r26007 = cos(r26006); double r26008 = theta; double r26009 = sin(r26008); double r26010 = r26007 * r26009; double r26011 = r26005 * r26010; double r26012 = sin(r26006); double r26013 = -r26012; double r26014 = cos(r26008); double r26015 = r26007 * r26005; double r26016 = cos(r26004); double r26017 = r26016 * r26012; double r26018 = fma(r26014, r26015, r26017); double r26019 = 3; double r26020 = pow(r26018, r26019); double r26021 = cbrt(r26020); double r26022 = asin(r26021); double r26023 = sin(r26022); double r26024 = fma(r26013, r26023, r26016); double r26025 = atan2(r26011, r26024); double r26026 = lambda1; double r26027 = r26025 + r26026; return r26027; } 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 r26028, r26029, r26030, r26031, r26032, r26033, r26034, r26035, r26036, r26037, r26038, r26039, r26040, r26041, r26042, r26043, r26044, r26045, r26046, r26047, r26048, r26049; void setup_mpfr_f_im() { mpfr_set_default_prec(592); 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); mpfr_init(r26047); mpfr_init(r26048); mpfr_init(r26049); } double f_im(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { mpfr_set_d(r26028, lambda1, MPFR_RNDN); mpfr_set_d(r26029, theta, MPFR_RNDN); mpfr_sin(r26030, r26029, MPFR_RNDN); mpfr_set_d(r26031, delta, MPFR_RNDN); mpfr_sin(r26032, r26031, MPFR_RNDN); mpfr_mul(r26033, r26030, r26032, MPFR_RNDN); mpfr_set_d(r26034, phi1, MPFR_RNDN); mpfr_cos(r26035, r26034, MPFR_RNDN); mpfr_mul(r26036, r26033, r26035, MPFR_RNDN); mpfr_cos(r26037, r26031, MPFR_RNDN); mpfr_sin(r26038, r26034, MPFR_RNDN); mpfr_mul(r26039, r26038, r26037, MPFR_RNDN); mpfr_mul(r26040, r26035, r26032, MPFR_RNDN); mpfr_cos(r26041, r26029, MPFR_RNDN); mpfr_mul(r26042, r26040, r26041, MPFR_RNDN); mpfr_add(r26043, r26039, r26042, MPFR_RNDN); mpfr_asin(r26044, r26043, MPFR_RNDN); mpfr_sin(r26045, r26044, MPFR_RNDN); mpfr_mul(r26046, r26038, r26045, MPFR_RNDN); mpfr_sub(r26047, r26037, r26046, MPFR_RNDN); mpfr_atan2(r26048, r26036, r26047, MPFR_RNDN); mpfr_add(r26049, r26028, r26048, MPFR_RNDN); return mpfr_get_d(r26049, MPFR_RNDN); } static mpfr_t r26050, r26051, r26052, r26053, r26054, r26055, r26056, r26057, r26058, r26059, r26060, r26061, r26062, r26063, r26064, r26065, r26066, r26067, r26068, r26069, r26070, r26071, r26072, r26073; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r26050); mpfr_init(r26051); mpfr_init(r26052); mpfr_init(r26053); mpfr_init(r26054); mpfr_init(r26055); mpfr_init(r26056); mpfr_init(r26057); mpfr_init(r26058); mpfr_init(r26059); mpfr_init(r26060); mpfr_init(r26061); mpfr_init(r26062); mpfr_init(r26063); mpfr_init(r26064); mpfr_init_set_str(r26065, "3", 10, MPFR_RNDN); mpfr_init(r26066); mpfr_init(r26067); mpfr_init(r26068); mpfr_init(r26069); mpfr_init(r26070); mpfr_init(r26071); mpfr_init(r26072); mpfr_init(r26073); } double f_fm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { mpfr_set_d(r26050, delta, MPFR_RNDN); mpfr_sin(r26051, r26050, MPFR_RNDN); mpfr_set_d(r26052, phi1, MPFR_RNDN); mpfr_cos(r26053, r26052, MPFR_RNDN); mpfr_set_d(r26054, theta, MPFR_RNDN); mpfr_sin(r26055, r26054, MPFR_RNDN); mpfr_mul(r26056, r26053, r26055, MPFR_RNDN); mpfr_mul(r26057, r26051, r26056, MPFR_RNDN); mpfr_sin(r26058, r26052, MPFR_RNDN); mpfr_neg(r26059, r26058, MPFR_RNDN); mpfr_cos(r26060, r26054, MPFR_RNDN); mpfr_mul(r26061, r26053, r26051, MPFR_RNDN); mpfr_cos(r26062, r26050, MPFR_RNDN); mpfr_mul(r26063, r26062, r26058, MPFR_RNDN); mpfr_fma(r26064, r26060, r26061, r26063, MPFR_RNDN); ; mpfr_pow(r26066, r26064, r26065, MPFR_RNDN); mpfr_cbrt(r26067, r26066, MPFR_RNDN); mpfr_asin(r26068, r26067, MPFR_RNDN); mpfr_sin(r26069, r26068, MPFR_RNDN); mpfr_fma(r26070, r26059, r26069, r26062, MPFR_RNDN); mpfr_atan2(r26071, r26057, r26070, MPFR_RNDN); mpfr_set_d(r26072, lambda1, MPFR_RNDN); mpfr_add(r26073, r26071, r26072, MPFR_RNDN); return mpfr_get_d(r26073, MPFR_RNDN); } static mpfr_t r26074, r26075, r26076, r26077, r26078, r26079, r26080, r26081, r26082, r26083, r26084, r26085, r26086, r26087, r26088, r26089, r26090, r26091, r26092, r26093, r26094, r26095, r26096, r26097; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r26074); mpfr_init(r26075); mpfr_init(r26076); mpfr_init(r26077); mpfr_init(r26078); mpfr_init(r26079); mpfr_init(r26080); mpfr_init(r26081); mpfr_init(r26082); mpfr_init(r26083); mpfr_init(r26084); mpfr_init(r26085); mpfr_init(r26086); mpfr_init(r26087); mpfr_init(r26088); mpfr_init_set_str(r26089, "3", 10, MPFR_RNDN); mpfr_init(r26090); mpfr_init(r26091); mpfr_init(r26092); mpfr_init(r26093); mpfr_init(r26094); mpfr_init(r26095); mpfr_init(r26096); mpfr_init(r26097); } double f_dm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) { mpfr_set_d(r26074, delta, MPFR_RNDN); mpfr_sin(r26075, r26074, MPFR_RNDN); mpfr_set_d(r26076, phi1, MPFR_RNDN); mpfr_cos(r26077, r26076, MPFR_RNDN); mpfr_set_d(r26078, theta, MPFR_RNDN); mpfr_sin(r26079, r26078, MPFR_RNDN); mpfr_mul(r26080, r26077, r26079, MPFR_RNDN); mpfr_mul(r26081, r26075, r26080, MPFR_RNDN); mpfr_sin(r26082, r26076, MPFR_RNDN); mpfr_neg(r26083, r26082, MPFR_RNDN); mpfr_cos(r26084, r26078, MPFR_RNDN); mpfr_mul(r26085, r26077, r26075, MPFR_RNDN); mpfr_cos(r26086, r26074, MPFR_RNDN); mpfr_mul(r26087, r26086, r26082, MPFR_RNDN); mpfr_fma(r26088, r26084, r26085, r26087, MPFR_RNDN); ; mpfr_pow(r26090, r26088, r26089, MPFR_RNDN); mpfr_cbrt(r26091, r26090, MPFR_RNDN); mpfr_asin(r26092, r26091, MPFR_RNDN); mpfr_sin(r26093, r26092, MPFR_RNDN); mpfr_fma(r26094, r26083, r26093, r26086, MPFR_RNDN); mpfr_atan2(r26095, r26081, r26094, MPFR_RNDN); mpfr_set_d(r26096, lambda1, MPFR_RNDN); mpfr_add(r26097, r26095, r26096, MPFR_RNDN); return mpfr_get_d(r26097, MPFR_RNDN); }