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