#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3b), real"; double f_if(float kx, float ky, float th) { float r22090 = ky; float r22091 = sin(r22090); float r22092 = kx; float r22093 = sin(r22092); float r22094 = 2; float r22095 = pow(r22093, r22094); float r22096 = pow(r22091, r22094); float r22097 = r22095 + r22096; float r22098 = sqrt(r22097); float r22099 = r22091 / r22098; float r22100 = th; float r22101 = sin(r22100); float r22102 = r22099 * r22101; return r22102; } double f_id(double kx, double ky, double th) { double r22103 = ky; double r22104 = sin(r22103); double r22105 = kx; double r22106 = sin(r22105); double r22107 = 2; double r22108 = pow(r22106, r22107); double r22109 = pow(r22104, r22107); double r22110 = r22108 + r22109; double r22111 = sqrt(r22110); double r22112 = r22104 / r22111; double r22113 = th; double r22114 = sin(r22113); double r22115 = r22112 * r22114; return r22115; } double f_of(float kx, float ky, float th) { float r22116 = 1; float r22117 = kx; float r22118 = sin(r22117); float r22119 = 2; float r22120 = pow(r22118, r22119); float r22121 = ky; float r22122 = sin(r22121); float r22123 = pow(r22122, r22119); float r22124 = r22120 + r22123; float r22125 = sqrt(r22124); float r22126 = r22125 / r22122; float r22127 = r22116 / r22126; float r22128 = th; float r22129 = sin(r22128); float r22130 = r22127 * r22129; return r22130; } double f_od(double kx, double ky, double th) { double r22131 = 1; double r22132 = kx; double r22133 = sin(r22132); double r22134 = 2; double r22135 = pow(r22133, r22134); double r22136 = ky; double r22137 = sin(r22136); double r22138 = pow(r22137, r22134); double r22139 = r22135 + r22138; double r22140 = sqrt(r22139); double r22141 = r22140 / r22137; double r22142 = r22131 / r22141; double r22143 = th; double r22144 = sin(r22143); double r22145 = r22142 * r22144; return r22145; } 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 r22146, r22147, r22148, r22149, r22150, r22151, r22152, r22153, r22154, r22155, r22156, r22157, r22158; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22146); mpfr_init(r22147); mpfr_init(r22148); mpfr_init(r22149); mpfr_init_set_str(r22150, "2", 10, MPFR_RNDN); mpfr_init(r22151); mpfr_init(r22152); mpfr_init(r22153); mpfr_init(r22154); mpfr_init(r22155); mpfr_init(r22156); mpfr_init(r22157); mpfr_init(r22158); } double f_im(double kx, double ky, double th) { mpfr_set_d(r22146, ky, MPFR_RNDN); mpfr_sin(r22147, r22146, MPFR_RNDN); mpfr_set_d(r22148, kx, MPFR_RNDN); mpfr_sin(r22149, r22148, MPFR_RNDN); ; mpfr_pow(r22151, r22149, r22150, MPFR_RNDN); mpfr_pow(r22152, r22147, r22150, MPFR_RNDN); mpfr_add(r22153, r22151, r22152, MPFR_RNDN); mpfr_sqrt(r22154, r22153, MPFR_RNDN); mpfr_div(r22155, r22147, r22154, MPFR_RNDN); mpfr_set_d(r22156, th, MPFR_RNDN); mpfr_sin(r22157, r22156, MPFR_RNDN); mpfr_mul(r22158, r22155, r22157, MPFR_RNDN); return mpfr_get_d(r22158, MPFR_RNDN); } static mpfr_t r22159, r22160, r22161, r22162, r22163, r22164, r22165, r22166, r22167, r22168, r22169, r22170, r22171, r22172, r22173; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22159, "1", 10, MPFR_RNDN); mpfr_init(r22160); mpfr_init(r22161); mpfr_init_set_str(r22162, "2", 10, MPFR_RNDN); 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); } double f_fm(double kx, double ky, double th) { ; mpfr_set_d(r22160, kx, MPFR_RNDN); mpfr_sin(r22161, r22160, MPFR_RNDN); ; mpfr_pow(r22163, r22161, r22162, MPFR_RNDN); mpfr_set_d(r22164, ky, MPFR_RNDN); mpfr_sin(r22165, r22164, MPFR_RNDN); mpfr_pow(r22166, r22165, r22162, MPFR_RNDN); mpfr_add(r22167, r22163, r22166, MPFR_RNDN); mpfr_sqrt(r22168, r22167, MPFR_RNDN); mpfr_div(r22169, r22168, r22165, MPFR_RNDN); mpfr_div(r22170, r22159, r22169, MPFR_RNDN); mpfr_set_d(r22171, th, MPFR_RNDN); mpfr_sin(r22172, r22171, MPFR_RNDN); mpfr_mul(r22173, r22170, r22172, MPFR_RNDN); return mpfr_get_d(r22173, MPFR_RNDN); } static mpfr_t r22174, r22175, r22176, r22177, r22178, r22179, r22180, r22181, r22182, r22183, r22184, r22185, r22186, r22187, r22188; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22174, "1", 10, MPFR_RNDN); mpfr_init(r22175); mpfr_init(r22176); mpfr_init_set_str(r22177, "2", 10, MPFR_RNDN); mpfr_init(r22178); mpfr_init(r22179); mpfr_init(r22180); mpfr_init(r22181); mpfr_init(r22182); mpfr_init(r22183); mpfr_init(r22184); mpfr_init(r22185); mpfr_init(r22186); mpfr_init(r22187); mpfr_init(r22188); } double f_dm(double kx, double ky, double th) { ; mpfr_set_d(r22175, kx, MPFR_RNDN); mpfr_sin(r22176, r22175, MPFR_RNDN); ; mpfr_pow(r22178, r22176, r22177, MPFR_RNDN); mpfr_set_d(r22179, ky, MPFR_RNDN); mpfr_sin(r22180, r22179, MPFR_RNDN); mpfr_pow(r22181, r22180, r22177, MPFR_RNDN); mpfr_add(r22182, r22178, r22181, MPFR_RNDN); mpfr_sqrt(r22183, r22182, MPFR_RNDN); mpfr_div(r22184, r22183, r22180, MPFR_RNDN); mpfr_div(r22185, r22174, r22184, MPFR_RNDN); mpfr_set_d(r22186, th, MPFR_RNDN); mpfr_sin(r22187, r22186, MPFR_RNDN); mpfr_mul(r22188, r22185, r22187, MPFR_RNDN); return mpfr_get_d(r22188, MPFR_RNDN); }