#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3a)"; double f_if(float l, float Om, float kx, float ky) { float r21914 = 1; float r21915 = 2; float r21916 = r21914 / r21915; float r21917 = l; float r21918 = r21915 * r21917; float r21919 = Om; float r21920 = r21918 / r21919; float r21921 = pow(r21920, r21915); float r21922 = kx; float r21923 = sin(r21922); float r21924 = pow(r21923, r21915); float r21925 = ky; float r21926 = sin(r21925); float r21927 = pow(r21926, r21915); float r21928 = r21924 + r21927; float r21929 = r21921 * r21928; float r21930 = r21914 + r21929; float r21931 = sqrt(r21930); float r21932 = r21914 / r21931; float r21933 = r21914 + r21932; float r21934 = r21916 * r21933; float r21935 = sqrt(r21934); return r21935; } double f_id(double l, double Om, double kx, double ky) { double r21936 = 1; double r21937 = 2; double r21938 = r21936 / r21937; double r21939 = l; double r21940 = r21937 * r21939; double r21941 = Om; double r21942 = r21940 / r21941; double r21943 = pow(r21942, r21937); double r21944 = kx; double r21945 = sin(r21944); double r21946 = pow(r21945, r21937); double r21947 = ky; double r21948 = sin(r21947); double r21949 = pow(r21948, r21937); double r21950 = r21946 + r21949; double r21951 = r21943 * r21950; double r21952 = r21936 + r21951; double r21953 = sqrt(r21952); double r21954 = r21936 / r21953; double r21955 = r21936 + r21954; double r21956 = r21938 * r21955; double r21957 = sqrt(r21956); return r21957; } double f_of(float l, float Om, float kx, float ky) { float r21958 = 1; float r21959 = 2; float r21960 = r21958 / r21959; float r21961 = l; float r21962 = r21959 * r21961; float r21963 = Om; float r21964 = r21962 / r21963; float r21965 = pow(r21964, r21959); float r21966 = kx; float r21967 = sin(r21966); float r21968 = pow(r21967, r21959); float r21969 = ky; float r21970 = sin(r21969); float r21971 = pow(r21970, r21959); float r21972 = r21968 + r21971; float r21973 = r21965 * r21972; float r21974 = r21958 + r21973; float r21975 = sqrt(r21974); float r21976 = r21958 / r21975; float r21977 = exp(r21976); float r21978 = log(r21977); float r21979 = r21958 + r21978; float r21980 = r21960 * r21979; float r21981 = sqrt(r21980); return r21981; } double f_od(double l, double Om, double kx, double ky) { double r21982 = 1; double r21983 = 2; double r21984 = r21982 / r21983; double r21985 = l; double r21986 = r21983 * r21985; double r21987 = Om; double r21988 = r21986 / r21987; double r21989 = pow(r21988, r21983); double r21990 = kx; double r21991 = sin(r21990); double r21992 = pow(r21991, r21983); double r21993 = ky; double r21994 = sin(r21993); double r21995 = pow(r21994, r21983); double r21996 = r21992 + r21995; double r21997 = r21989 * r21996; double r21998 = r21982 + r21997; double r21999 = sqrt(r21998); double r22000 = r21982 / r21999; double r22001 = exp(r22000); double r22002 = log(r22001); double r22003 = r21982 + r22002; double r22004 = r21984 * r22003; double r22005 = sqrt(r22004); return r22005; } 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 r22006, r22007, r22008, r22009, r22010, r22011, r22012, r22013, r22014, r22015, r22016, r22017, r22018, r22019, r22020, r22021, r22022, r22023, r22024, r22025, r22026, r22027; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r22006, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22007, "2", 10, MPFR_RNDN); mpfr_init(r22008); mpfr_init(r22009); mpfr_init(r22010); mpfr_init(r22011); mpfr_init(r22012); mpfr_init(r22013); mpfr_init(r22014); mpfr_init(r22015); mpfr_init(r22016); mpfr_init(r22017); mpfr_init(r22018); mpfr_init(r22019); mpfr_init(r22020); mpfr_init(r22021); mpfr_init(r22022); mpfr_init(r22023); mpfr_init(r22024); mpfr_init(r22025); mpfr_init(r22026); mpfr_init(r22027); } double f_im(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22008, r22006, r22007, MPFR_RNDN); mpfr_set_d(r22009, l, MPFR_RNDN); mpfr_mul(r22010, r22007, r22009, MPFR_RNDN); mpfr_set_d(r22011, Om, MPFR_RNDN); mpfr_div(r22012, r22010, r22011, MPFR_RNDN); mpfr_pow(r22013, r22012, r22007, MPFR_RNDN); mpfr_set_d(r22014, kx, MPFR_RNDN); mpfr_sin(r22015, r22014, MPFR_RNDN); mpfr_pow(r22016, r22015, r22007, MPFR_RNDN); mpfr_set_d(r22017, ky, MPFR_RNDN); mpfr_sin(r22018, r22017, MPFR_RNDN); mpfr_pow(r22019, r22018, r22007, MPFR_RNDN); mpfr_add(r22020, r22016, r22019, MPFR_RNDN); mpfr_mul(r22021, r22013, r22020, MPFR_RNDN); mpfr_add(r22022, r22006, r22021, MPFR_RNDN); mpfr_sqrt(r22023, r22022, MPFR_RNDN); mpfr_div(r22024, r22006, r22023, MPFR_RNDN); mpfr_add(r22025, r22006, r22024, MPFR_RNDN); mpfr_mul(r22026, r22008, r22025, MPFR_RNDN); mpfr_sqrt(r22027, r22026, MPFR_RNDN); return mpfr_get_d(r22027, MPFR_RNDN); } static mpfr_t r22028, r22029, r22030, r22031, r22032, r22033, r22034, r22035, r22036, r22037, r22038, r22039, r22040, r22041, r22042, r22043, r22044, r22045, r22046, r22047, r22048, r22049, r22050, r22051; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r22028, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22029, "2", 10, MPFR_RNDN); mpfr_init(r22030); mpfr_init(r22031); mpfr_init(r22032); mpfr_init(r22033); mpfr_init(r22034); mpfr_init(r22035); mpfr_init(r22036); mpfr_init(r22037); mpfr_init(r22038); mpfr_init(r22039); mpfr_init(r22040); mpfr_init(r22041); mpfr_init(r22042); mpfr_init(r22043); mpfr_init(r22044); mpfr_init(r22045); mpfr_init(r22046); mpfr_init(r22047); mpfr_init(r22048); mpfr_init(r22049); mpfr_init(r22050); mpfr_init(r22051); } double f_fm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22030, r22028, r22029, MPFR_RNDN); mpfr_set_d(r22031, l, MPFR_RNDN); mpfr_mul(r22032, r22029, r22031, MPFR_RNDN); mpfr_set_d(r22033, Om, MPFR_RNDN); mpfr_div(r22034, r22032, r22033, MPFR_RNDN); mpfr_pow(r22035, r22034, r22029, MPFR_RNDN); mpfr_set_d(r22036, kx, MPFR_RNDN); mpfr_sin(r22037, r22036, MPFR_RNDN); mpfr_pow(r22038, r22037, r22029, MPFR_RNDN); mpfr_set_d(r22039, ky, MPFR_RNDN); mpfr_sin(r22040, r22039, MPFR_RNDN); mpfr_pow(r22041, r22040, r22029, MPFR_RNDN); mpfr_add(r22042, r22038, r22041, MPFR_RNDN); mpfr_mul(r22043, r22035, r22042, MPFR_RNDN); mpfr_add(r22044, r22028, r22043, MPFR_RNDN); mpfr_sqrt(r22045, r22044, MPFR_RNDN); mpfr_div(r22046, r22028, r22045, MPFR_RNDN); mpfr_exp(r22047, r22046, MPFR_RNDN); mpfr_log(r22048, r22047, MPFR_RNDN); mpfr_add(r22049, r22028, r22048, MPFR_RNDN); mpfr_mul(r22050, r22030, r22049, MPFR_RNDN); mpfr_sqrt(r22051, r22050, MPFR_RNDN); return mpfr_get_d(r22051, MPFR_RNDN); } static mpfr_t r22052, r22053, r22054, r22055, r22056, r22057, r22058, r22059, r22060, r22061, r22062, r22063, r22064, r22065, r22066, r22067, r22068, r22069, r22070, r22071, r22072, r22073, r22074, r22075; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r22052, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22053, "2", 10, MPFR_RNDN); mpfr_init(r22054); mpfr_init(r22055); mpfr_init(r22056); mpfr_init(r22057); mpfr_init(r22058); mpfr_init(r22059); mpfr_init(r22060); mpfr_init(r22061); mpfr_init(r22062); mpfr_init(r22063); mpfr_init(r22064); mpfr_init(r22065); mpfr_init(r22066); mpfr_init(r22067); mpfr_init(r22068); mpfr_init(r22069); mpfr_init(r22070); mpfr_init(r22071); mpfr_init(r22072); mpfr_init(r22073); mpfr_init(r22074); mpfr_init(r22075); } double f_dm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22054, r22052, r22053, MPFR_RNDN); mpfr_set_d(r22055, l, MPFR_RNDN); mpfr_mul(r22056, r22053, r22055, MPFR_RNDN); mpfr_set_d(r22057, Om, MPFR_RNDN); mpfr_div(r22058, r22056, r22057, MPFR_RNDN); mpfr_pow(r22059, r22058, r22053, MPFR_RNDN); mpfr_set_d(r22060, kx, MPFR_RNDN); mpfr_sin(r22061, r22060, MPFR_RNDN); mpfr_pow(r22062, r22061, r22053, MPFR_RNDN); mpfr_set_d(r22063, ky, MPFR_RNDN); mpfr_sin(r22064, r22063, MPFR_RNDN); mpfr_pow(r22065, r22064, r22053, MPFR_RNDN); mpfr_add(r22066, r22062, r22065, MPFR_RNDN); mpfr_mul(r22067, r22059, r22066, MPFR_RNDN); mpfr_add(r22068, r22052, r22067, MPFR_RNDN); mpfr_sqrt(r22069, r22068, MPFR_RNDN); mpfr_div(r22070, r22052, r22069, MPFR_RNDN); mpfr_exp(r22071, r22070, MPFR_RNDN); mpfr_log(r22072, r22071, MPFR_RNDN); mpfr_add(r22073, r22052, r22072, MPFR_RNDN); mpfr_mul(r22074, r22054, r22073, MPFR_RNDN); mpfr_sqrt(r22075, r22074, MPFR_RNDN); return mpfr_get_d(r22075, MPFR_RNDN); }