#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3a)"; double f_if(float l, float Om, float kx, float ky) { float r22048 = 1; float r22049 = 2; float r22050 = r22048 / r22049; float r22051 = l; float r22052 = r22049 * r22051; float r22053 = Om; float r22054 = r22052 / r22053; float r22055 = pow(r22054, r22049); float r22056 = kx; float r22057 = sin(r22056); float r22058 = pow(r22057, r22049); float r22059 = ky; float r22060 = sin(r22059); float r22061 = pow(r22060, r22049); float r22062 = r22058 + r22061; float r22063 = r22055 * r22062; float r22064 = r22048 + r22063; float r22065 = sqrt(r22064); float r22066 = r22048 / r22065; float r22067 = r22048 + r22066; float r22068 = r22050 * r22067; float r22069 = sqrt(r22068); return r22069; } double f_id(double l, double Om, double kx, double ky) { double r22070 = 1; double r22071 = 2; double r22072 = r22070 / r22071; double r22073 = l; double r22074 = r22071 * r22073; double r22075 = Om; double r22076 = r22074 / r22075; double r22077 = pow(r22076, r22071); double r22078 = kx; double r22079 = sin(r22078); double r22080 = pow(r22079, r22071); double r22081 = ky; double r22082 = sin(r22081); double r22083 = pow(r22082, r22071); double r22084 = r22080 + r22083; double r22085 = r22077 * r22084; double r22086 = r22070 + r22085; double r22087 = sqrt(r22086); double r22088 = r22070 / r22087; double r22089 = r22070 + r22088; double r22090 = r22072 * r22089; double r22091 = sqrt(r22090); return r22091; } double f_of(float l, float Om, float kx, float ky) { float r22092 = 1; float r22093 = 2; float r22094 = r22092 / r22093; float r22095 = l; float r22096 = r22093 * r22095; float r22097 = Om; float r22098 = r22096 / r22097; float r22099 = kx; float r22100 = sin(r22099); float r22101 = pow(r22100, r22093); float r22102 = ky; float r22103 = sin(r22102); float r22104 = pow(r22103, r22093); float r22105 = r22101 + r22104; float r22106 = r22098 * r22105; float r22107 = r22098 * r22106; float r22108 = r22092 + r22107; float r22109 = sqrt(r22108); float r22110 = r22092 / r22109; float r22111 = r22092 + r22110; float r22112 = r22094 * r22111; float r22113 = sqrt(r22112); return r22113; } double f_od(double l, double Om, double kx, double ky) { double r22114 = 1; double r22115 = 2; double r22116 = r22114 / r22115; double r22117 = l; double r22118 = r22115 * r22117; double r22119 = Om; double r22120 = r22118 / r22119; double r22121 = kx; double r22122 = sin(r22121); double r22123 = pow(r22122, r22115); double r22124 = ky; double r22125 = sin(r22124); double r22126 = pow(r22125, r22115); double r22127 = r22123 + r22126; double r22128 = r22120 * r22127; double r22129 = r22120 * r22128; double r22130 = r22114 + r22129; double r22131 = sqrt(r22130); double r22132 = r22114 / r22131; double r22133 = r22114 + r22132; double r22134 = r22116 * r22133; double r22135 = sqrt(r22134); return r22135; } 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 r22136, r22137, r22138, r22139, r22140, r22141, r22142, r22143, r22144, r22145, r22146, r22147, r22148, r22149, r22150, r22151, r22152, r22153, r22154, r22155, r22156, r22157; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r22136, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22137, "2", 10, MPFR_RNDN); mpfr_init(r22138); mpfr_init(r22139); mpfr_init(r22140); mpfr_init(r22141); mpfr_init(r22142); mpfr_init(r22143); mpfr_init(r22144); mpfr_init(r22145); mpfr_init(r22146); mpfr_init(r22147); mpfr_init(r22148); mpfr_init(r22149); mpfr_init(r22150); mpfr_init(r22151); mpfr_init(r22152); mpfr_init(r22153); mpfr_init(r22154); mpfr_init(r22155); mpfr_init(r22156); mpfr_init(r22157); } double f_im(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22138, r22136, r22137, MPFR_RNDN); mpfr_set_d(r22139, l, MPFR_RNDN); mpfr_mul(r22140, r22137, r22139, MPFR_RNDN); mpfr_set_d(r22141, Om, MPFR_RNDN); mpfr_div(r22142, r22140, r22141, MPFR_RNDN); mpfr_pow(r22143, r22142, r22137, MPFR_RNDN); mpfr_set_d(r22144, kx, MPFR_RNDN); mpfr_sin(r22145, r22144, MPFR_RNDN); mpfr_pow(r22146, r22145, r22137, MPFR_RNDN); mpfr_set_d(r22147, ky, MPFR_RNDN); mpfr_sin(r22148, r22147, MPFR_RNDN); mpfr_pow(r22149, r22148, r22137, MPFR_RNDN); mpfr_add(r22150, r22146, r22149, MPFR_RNDN); mpfr_mul(r22151, r22143, r22150, MPFR_RNDN); mpfr_add(r22152, r22136, r22151, MPFR_RNDN); mpfr_sqrt(r22153, r22152, MPFR_RNDN); mpfr_div(r22154, r22136, r22153, MPFR_RNDN); mpfr_add(r22155, r22136, r22154, MPFR_RNDN); mpfr_mul(r22156, r22138, r22155, MPFR_RNDN); mpfr_sqrt(r22157, r22156, MPFR_RNDN); return mpfr_get_d(r22157, MPFR_RNDN); } static mpfr_t r22158, r22159, r22160, r22161, r22162, r22163, r22164, r22165, r22166, r22167, r22168, r22169, r22170, r22171, r22172, r22173, r22174, r22175, r22176, r22177, r22178, r22179; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r22158, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22159, "2", 10, MPFR_RNDN); mpfr_init(r22160); mpfr_init(r22161); mpfr_init(r22162); mpfr_init(r22163); mpfr_init(r22164); mpfr_init(r22165); mpfr_init(r22166); mpfr_init(r22167); mpfr_init(r22168); mpfr_init(r22169); mpfr_init(r22170); mpfr_init(r22171); mpfr_init(r22172); mpfr_init(r22173); mpfr_init(r22174); mpfr_init(r22175); mpfr_init(r22176); mpfr_init(r22177); mpfr_init(r22178); mpfr_init(r22179); } double f_fm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22160, r22158, r22159, MPFR_RNDN); mpfr_set_d(r22161, l, MPFR_RNDN); mpfr_mul(r22162, r22159, r22161, MPFR_RNDN); mpfr_set_d(r22163, Om, MPFR_RNDN); mpfr_div(r22164, r22162, r22163, MPFR_RNDN); mpfr_set_d(r22165, kx, MPFR_RNDN); mpfr_sin(r22166, r22165, MPFR_RNDN); mpfr_pow(r22167, r22166, r22159, MPFR_RNDN); mpfr_set_d(r22168, ky, MPFR_RNDN); mpfr_sin(r22169, r22168, MPFR_RNDN); mpfr_pow(r22170, r22169, r22159, MPFR_RNDN); mpfr_add(r22171, r22167, r22170, MPFR_RNDN); mpfr_mul(r22172, r22164, r22171, MPFR_RNDN); mpfr_mul(r22173, r22164, r22172, MPFR_RNDN); mpfr_add(r22174, r22158, r22173, MPFR_RNDN); mpfr_sqrt(r22175, r22174, MPFR_RNDN); mpfr_div(r22176, r22158, r22175, MPFR_RNDN); mpfr_add(r22177, r22158, r22176, MPFR_RNDN); mpfr_mul(r22178, r22160, r22177, MPFR_RNDN); mpfr_sqrt(r22179, r22178, MPFR_RNDN); return mpfr_get_d(r22179, MPFR_RNDN); } static mpfr_t r22180, r22181, r22182, r22183, r22184, r22185, r22186, r22187, r22188, r22189, r22190, r22191, r22192, r22193, r22194, r22195, r22196, r22197, r22198, r22199, r22200, r22201; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r22180, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22181, "2", 10, MPFR_RNDN); mpfr_init(r22182); mpfr_init(r22183); mpfr_init(r22184); mpfr_init(r22185); mpfr_init(r22186); mpfr_init(r22187); mpfr_init(r22188); mpfr_init(r22189); mpfr_init(r22190); mpfr_init(r22191); mpfr_init(r22192); mpfr_init(r22193); mpfr_init(r22194); mpfr_init(r22195); mpfr_init(r22196); mpfr_init(r22197); mpfr_init(r22198); mpfr_init(r22199); mpfr_init(r22200); mpfr_init(r22201); } double f_dm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22182, r22180, r22181, MPFR_RNDN); mpfr_set_d(r22183, l, MPFR_RNDN); mpfr_mul(r22184, r22181, r22183, MPFR_RNDN); mpfr_set_d(r22185, Om, MPFR_RNDN); mpfr_div(r22186, r22184, r22185, MPFR_RNDN); mpfr_set_d(r22187, kx, MPFR_RNDN); mpfr_sin(r22188, r22187, MPFR_RNDN); mpfr_pow(r22189, r22188, r22181, MPFR_RNDN); mpfr_set_d(r22190, ky, MPFR_RNDN); mpfr_sin(r22191, r22190, MPFR_RNDN); mpfr_pow(r22192, r22191, r22181, MPFR_RNDN); mpfr_add(r22193, r22189, r22192, MPFR_RNDN); mpfr_mul(r22194, r22186, r22193, MPFR_RNDN); mpfr_mul(r22195, r22186, r22194, MPFR_RNDN); mpfr_add(r22196, r22180, r22195, MPFR_RNDN); mpfr_sqrt(r22197, r22196, MPFR_RNDN); mpfr_div(r22198, r22180, r22197, MPFR_RNDN); mpfr_add(r22199, r22180, r22198, MPFR_RNDN); mpfr_mul(r22200, r22182, r22199, MPFR_RNDN); mpfr_sqrt(r22201, r22200, MPFR_RNDN); return mpfr_get_d(r22201, MPFR_RNDN); }