#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3b), real"; double f_if(float kx, float ky, float th) { float r21950 = ky; float r21951 = sin(r21950); float r21952 = kx; float r21953 = sin(r21952); float r21954 = 2; float r21955 = pow(r21953, r21954); float r21956 = pow(r21951, r21954); float r21957 = r21955 + r21956; float r21958 = sqrt(r21957); float r21959 = r21951 / r21958; float r21960 = th; float r21961 = sin(r21960); float r21962 = r21959 * r21961; return r21962; } double f_id(double kx, double ky, double th) { double r21963 = ky; double r21964 = sin(r21963); double r21965 = kx; double r21966 = sin(r21965); double r21967 = 2; double r21968 = pow(r21966, r21967); double r21969 = pow(r21964, r21967); double r21970 = r21968 + r21969; double r21971 = sqrt(r21970); double r21972 = r21964 / r21971; double r21973 = th; double r21974 = sin(r21973); double r21975 = r21972 * r21974; return r21975; } double f_of(float kx, float ky, float th) { float r21976 = ky; float r21977 = sin(r21976); float r21978 = kx; float r21979 = sin(r21978); float r21980 = 2; float r21981 = pow(r21979, r21980); float r21982 = cbrt(r21981); float r21983 = r21982 * r21982; float r21984 = r21983 * r21982; float r21985 = pow(r21977, r21980); float r21986 = r21984 + r21985; float r21987 = sqrt(r21986); float r21988 = r21977 / r21987; float r21989 = th; float r21990 = sin(r21989); float r21991 = r21988 * r21990; return r21991; } double f_od(double kx, double ky, double th) { double r21992 = ky; double r21993 = sin(r21992); double r21994 = kx; double r21995 = sin(r21994); double r21996 = 2; double r21997 = pow(r21995, r21996); double r21998 = cbrt(r21997); double r21999 = r21998 * r21998; double r22000 = r21999 * r21998; double r22001 = pow(r21993, r21996); double r22002 = r22000 + r22001; double r22003 = sqrt(r22002); double r22004 = r21993 / r22003; double r22005 = th; double r22006 = sin(r22005); double r22007 = r22004 * r22006; return r22007; } 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 r22008, r22009, r22010, r22011, r22012, r22013, r22014, r22015, r22016, r22017, r22018, r22019, r22020; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22008); mpfr_init(r22009); mpfr_init(r22010); mpfr_init(r22011); mpfr_init_set_str(r22012, "2", 10, MPFR_RNDN); mpfr_init(r22013); mpfr_init(r22014); mpfr_init(r22015); mpfr_init(r22016); mpfr_init(r22017); mpfr_init(r22018); mpfr_init(r22019); mpfr_init(r22020); } double f_im(double kx, double ky, double th) { mpfr_set_d(r22008, ky, MPFR_RNDN); mpfr_sin(r22009, r22008, MPFR_RNDN); mpfr_set_d(r22010, kx, MPFR_RNDN); mpfr_sin(r22011, r22010, MPFR_RNDN); ; mpfr_pow(r22013, r22011, r22012, MPFR_RNDN); mpfr_pow(r22014, r22009, r22012, MPFR_RNDN); mpfr_add(r22015, r22013, r22014, MPFR_RNDN); mpfr_sqrt(r22016, r22015, MPFR_RNDN); mpfr_div(r22017, r22009, r22016, MPFR_RNDN); mpfr_set_d(r22018, th, MPFR_RNDN); mpfr_sin(r22019, r22018, MPFR_RNDN); mpfr_mul(r22020, r22017, r22019, MPFR_RNDN); return mpfr_get_d(r22020, MPFR_RNDN); } static mpfr_t r22021, r22022, r22023, r22024, r22025, r22026, r22027, r22028, r22029, r22030, r22031, r22032, r22033, r22034, r22035, r22036; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22021); mpfr_init(r22022); mpfr_init(r22023); mpfr_init(r22024); mpfr_init_set_str(r22025, "2", 10, MPFR_RNDN); mpfr_init(r22026); mpfr_init(r22027); mpfr_init(r22028); mpfr_init(r22029); mpfr_init(r22030); mpfr_init(r22031); mpfr_init(r22032); mpfr_init(r22033); mpfr_init(r22034); mpfr_init(r22035); mpfr_init(r22036); } double f_fm(double kx, double ky, double th) { mpfr_set_d(r22021, ky, MPFR_RNDN); mpfr_sin(r22022, r22021, MPFR_RNDN); mpfr_set_d(r22023, kx, MPFR_RNDN); mpfr_sin(r22024, r22023, MPFR_RNDN); ; mpfr_pow(r22026, r22024, r22025, MPFR_RNDN); mpfr_cbrt(r22027, r22026, MPFR_RNDN); mpfr_mul(r22028, r22027, r22027, MPFR_RNDN); mpfr_mul(r22029, r22028, r22027, MPFR_RNDN); mpfr_pow(r22030, r22022, r22025, MPFR_RNDN); mpfr_add(r22031, r22029, r22030, MPFR_RNDN); mpfr_sqrt(r22032, r22031, MPFR_RNDN); mpfr_div(r22033, r22022, r22032, MPFR_RNDN); mpfr_set_d(r22034, th, MPFR_RNDN); mpfr_sin(r22035, r22034, MPFR_RNDN); mpfr_mul(r22036, r22033, r22035, MPFR_RNDN); return mpfr_get_d(r22036, MPFR_RNDN); } static mpfr_t r22037, r22038, r22039, r22040, r22041, r22042, r22043, r22044, r22045, r22046, r22047, r22048, r22049, r22050, r22051, r22052; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22037); mpfr_init(r22038); mpfr_init(r22039); mpfr_init(r22040); mpfr_init_set_str(r22041, "2", 10, MPFR_RNDN); 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); } double f_dm(double kx, double ky, double th) { mpfr_set_d(r22037, ky, MPFR_RNDN); mpfr_sin(r22038, r22037, MPFR_RNDN); mpfr_set_d(r22039, kx, MPFR_RNDN); mpfr_sin(r22040, r22039, MPFR_RNDN); ; mpfr_pow(r22042, r22040, r22041, MPFR_RNDN); mpfr_cbrt(r22043, r22042, MPFR_RNDN); mpfr_mul(r22044, r22043, r22043, MPFR_RNDN); mpfr_mul(r22045, r22044, r22043, MPFR_RNDN); mpfr_pow(r22046, r22038, r22041, MPFR_RNDN); mpfr_add(r22047, r22045, r22046, MPFR_RNDN); mpfr_sqrt(r22048, r22047, MPFR_RNDN); mpfr_div(r22049, r22038, r22048, MPFR_RNDN); mpfr_set_d(r22050, th, MPFR_RNDN); mpfr_sin(r22051, r22050, MPFR_RNDN); mpfr_mul(r22052, r22049, r22051, MPFR_RNDN); return mpfr_get_d(r22052, MPFR_RNDN); }