#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3a)"; double f_if(float l, float Om, float kx, float ky) { float r21944 = 1; float r21945 = 2; float r21946 = r21944 / r21945; float r21947 = l; float r21948 = r21945 * r21947; float r21949 = Om; float r21950 = r21948 / r21949; float r21951 = pow(r21950, r21945); float r21952 = kx; float r21953 = sin(r21952); float r21954 = pow(r21953, r21945); float r21955 = ky; float r21956 = sin(r21955); float r21957 = pow(r21956, r21945); float r21958 = r21954 + r21957; float r21959 = r21951 * r21958; float r21960 = r21944 + r21959; float r21961 = sqrt(r21960); float r21962 = r21944 / r21961; float r21963 = r21944 + r21962; float r21964 = r21946 * r21963; float r21965 = sqrt(r21964); return r21965; } double f_id(double l, double Om, double kx, double ky) { double r21966 = 1; double r21967 = 2; double r21968 = r21966 / r21967; double r21969 = l; double r21970 = r21967 * r21969; double r21971 = Om; double r21972 = r21970 / r21971; double r21973 = pow(r21972, r21967); double r21974 = kx; double r21975 = sin(r21974); double r21976 = pow(r21975, r21967); double r21977 = ky; double r21978 = sin(r21977); double r21979 = pow(r21978, r21967); double r21980 = r21976 + r21979; double r21981 = r21973 * r21980; double r21982 = r21966 + r21981; double r21983 = sqrt(r21982); double r21984 = r21966 / r21983; double r21985 = r21966 + r21984; double r21986 = r21968 * r21985; double r21987 = sqrt(r21986); return r21987; } double f_of(float l, float Om, float kx, float ky) { float r21988 = 1; float r21989 = 2; float r21990 = r21988 / r21989; float r21991 = l; float r21992 = r21989 * r21991; float r21993 = Om; float r21994 = r21992 / r21993; float r21995 = pow(r21994, r21989); float r21996 = kx; float r21997 = sin(r21996); float r21998 = pow(r21997, r21989); float r21999 = ky; float r22000 = sin(r21999); float r22001 = pow(r22000, r21989); float r22002 = r21998 + r22001; float r22003 = r21995 * r22002; float r22004 = r21988 + r22003; float r22005 = sqrt(r22004); float r22006 = sqrt(r22005); float r22007 = r22006 * r22006; float r22008 = r21988 / r22007; float r22009 = r21988 + r22008; float r22010 = r21990 * r22009; float r22011 = sqrt(r22010); return r22011; } double f_od(double l, double Om, double kx, double ky) { double r22012 = 1; double r22013 = 2; double r22014 = r22012 / r22013; double r22015 = l; double r22016 = r22013 * r22015; double r22017 = Om; double r22018 = r22016 / r22017; double r22019 = pow(r22018, r22013); double r22020 = kx; double r22021 = sin(r22020); double r22022 = pow(r22021, r22013); double r22023 = ky; double r22024 = sin(r22023); double r22025 = pow(r22024, r22013); double r22026 = r22022 + r22025; double r22027 = r22019 * r22026; double r22028 = r22012 + r22027; double r22029 = sqrt(r22028); double r22030 = sqrt(r22029); double r22031 = r22030 * r22030; double r22032 = r22012 / r22031; double r22033 = r22012 + r22032; double r22034 = r22014 * r22033; double r22035 = sqrt(r22034); return r22035; } 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 r22036, r22037, r22038, r22039, r22040, r22041, r22042, r22043, r22044, r22045, r22046, r22047, r22048, r22049, r22050, r22051, r22052, r22053, r22054, r22055, r22056, r22057; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r22036, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22037, "2", 10, MPFR_RNDN); 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); mpfr_init(r22052); mpfr_init(r22053); mpfr_init(r22054); mpfr_init(r22055); mpfr_init(r22056); mpfr_init(r22057); } double f_im(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22038, r22036, r22037, MPFR_RNDN); mpfr_set_d(r22039, l, MPFR_RNDN); mpfr_mul(r22040, r22037, r22039, MPFR_RNDN); mpfr_set_d(r22041, Om, MPFR_RNDN); mpfr_div(r22042, r22040, r22041, MPFR_RNDN); mpfr_pow(r22043, r22042, r22037, MPFR_RNDN); mpfr_set_d(r22044, kx, MPFR_RNDN); mpfr_sin(r22045, r22044, MPFR_RNDN); mpfr_pow(r22046, r22045, r22037, MPFR_RNDN); mpfr_set_d(r22047, ky, MPFR_RNDN); mpfr_sin(r22048, r22047, MPFR_RNDN); mpfr_pow(r22049, r22048, r22037, MPFR_RNDN); mpfr_add(r22050, r22046, r22049, MPFR_RNDN); mpfr_mul(r22051, r22043, r22050, MPFR_RNDN); mpfr_add(r22052, r22036, r22051, MPFR_RNDN); mpfr_sqrt(r22053, r22052, MPFR_RNDN); mpfr_div(r22054, r22036, r22053, MPFR_RNDN); mpfr_add(r22055, r22036, r22054, MPFR_RNDN); mpfr_mul(r22056, r22038, r22055, MPFR_RNDN); mpfr_sqrt(r22057, r22056, MPFR_RNDN); return mpfr_get_d(r22057, MPFR_RNDN); } static mpfr_t r22058, r22059, r22060, r22061, r22062, r22063, r22064, r22065, r22066, r22067, r22068, r22069, r22070, r22071, r22072, r22073, r22074, r22075, r22076, r22077, r22078, r22079, r22080, r22081; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22058, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22059, "2", 10, MPFR_RNDN); 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); mpfr_init(r22076); mpfr_init(r22077); mpfr_init(r22078); mpfr_init(r22079); mpfr_init(r22080); mpfr_init(r22081); } double f_fm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22060, r22058, r22059, MPFR_RNDN); mpfr_set_d(r22061, l, MPFR_RNDN); mpfr_mul(r22062, r22059, r22061, MPFR_RNDN); mpfr_set_d(r22063, Om, MPFR_RNDN); mpfr_div(r22064, r22062, r22063, MPFR_RNDN); mpfr_pow(r22065, r22064, r22059, MPFR_RNDN); mpfr_set_d(r22066, kx, MPFR_RNDN); mpfr_sin(r22067, r22066, MPFR_RNDN); mpfr_pow(r22068, r22067, r22059, MPFR_RNDN); mpfr_set_d(r22069, ky, MPFR_RNDN); mpfr_sin(r22070, r22069, MPFR_RNDN); mpfr_pow(r22071, r22070, r22059, MPFR_RNDN); mpfr_add(r22072, r22068, r22071, MPFR_RNDN); mpfr_mul(r22073, r22065, r22072, MPFR_RNDN); mpfr_add(r22074, r22058, r22073, MPFR_RNDN); mpfr_sqrt(r22075, r22074, MPFR_RNDN); mpfr_sqrt(r22076, r22075, MPFR_RNDN); mpfr_mul(r22077, r22076, r22076, MPFR_RNDN); mpfr_div(r22078, r22058, r22077, MPFR_RNDN); mpfr_add(r22079, r22058, r22078, MPFR_RNDN); mpfr_mul(r22080, r22060, r22079, MPFR_RNDN); mpfr_sqrt(r22081, r22080, MPFR_RNDN); return mpfr_get_d(r22081, MPFR_RNDN); } static mpfr_t r22082, r22083, r22084, r22085, r22086, r22087, r22088, r22089, r22090, r22091, r22092, r22093, r22094, r22095, r22096, r22097, r22098, r22099, r22100, r22101, r22102, r22103, r22104, r22105; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22082, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22083, "2", 10, MPFR_RNDN); mpfr_init(r22084); mpfr_init(r22085); mpfr_init(r22086); mpfr_init(r22087); mpfr_init(r22088); mpfr_init(r22089); mpfr_init(r22090); mpfr_init(r22091); mpfr_init(r22092); mpfr_init(r22093); mpfr_init(r22094); mpfr_init(r22095); mpfr_init(r22096); mpfr_init(r22097); mpfr_init(r22098); mpfr_init(r22099); mpfr_init(r22100); mpfr_init(r22101); mpfr_init(r22102); mpfr_init(r22103); mpfr_init(r22104); mpfr_init(r22105); } double f_dm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22084, r22082, r22083, MPFR_RNDN); mpfr_set_d(r22085, l, MPFR_RNDN); mpfr_mul(r22086, r22083, r22085, MPFR_RNDN); mpfr_set_d(r22087, Om, MPFR_RNDN); mpfr_div(r22088, r22086, r22087, MPFR_RNDN); mpfr_pow(r22089, r22088, r22083, MPFR_RNDN); mpfr_set_d(r22090, kx, MPFR_RNDN); mpfr_sin(r22091, r22090, MPFR_RNDN); mpfr_pow(r22092, r22091, r22083, MPFR_RNDN); mpfr_set_d(r22093, ky, MPFR_RNDN); mpfr_sin(r22094, r22093, MPFR_RNDN); mpfr_pow(r22095, r22094, r22083, MPFR_RNDN); mpfr_add(r22096, r22092, r22095, MPFR_RNDN); mpfr_mul(r22097, r22089, r22096, MPFR_RNDN); mpfr_add(r22098, r22082, r22097, MPFR_RNDN); mpfr_sqrt(r22099, r22098, MPFR_RNDN); mpfr_sqrt(r22100, r22099, MPFR_RNDN); mpfr_mul(r22101, r22100, r22100, MPFR_RNDN); mpfr_div(r22102, r22082, r22101, MPFR_RNDN); mpfr_add(r22103, r22082, r22102, MPFR_RNDN); mpfr_mul(r22104, r22084, r22103, MPFR_RNDN); mpfr_sqrt(r22105, r22104, MPFR_RNDN); return mpfr_get_d(r22105, MPFR_RNDN); }