#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3b), real"; double f_if(float kx, float ky, float th) { float r21107 = ky; float r21108 = sin(r21107); float r21109 = kx; float r21110 = sin(r21109); float r21111 = 2; float r21112 = pow(r21110, r21111); float r21113 = pow(r21108, r21111); float r21114 = r21112 + r21113; float r21115 = sqrt(r21114); float r21116 = r21108 / r21115; float r21117 = th; float r21118 = sin(r21117); float r21119 = r21116 * r21118; return r21119; } double f_id(double kx, double ky, double th) { double r21120 = ky; double r21121 = sin(r21120); double r21122 = kx; double r21123 = sin(r21122); double r21124 = 2; double r21125 = pow(r21123, r21124); double r21126 = pow(r21121, r21124); double r21127 = r21125 + r21126; double r21128 = sqrt(r21127); double r21129 = r21121 / r21128; double r21130 = th; double r21131 = sin(r21130); double r21132 = r21129 * r21131; return r21132; } double f_of(float kx, float ky, float th) { float r21133 = 1; float r21134 = kx; float r21135 = sin(r21134); float r21136 = 2; float r21137 = pow(r21135, r21136); float r21138 = ky; float r21139 = sin(r21138); float r21140 = pow(r21139, r21136); float r21141 = r21137 + r21140; float r21142 = sqrt(r21141); float r21143 = r21142 / r21139; float r21144 = r21133 / r21143; float r21145 = th; float r21146 = sin(r21145); float r21147 = r21144 * r21146; return r21147; } double f_od(double kx, double ky, double th) { double r21148 = 1; double r21149 = kx; double r21150 = sin(r21149); double r21151 = 2; double r21152 = pow(r21150, r21151); double r21153 = ky; double r21154 = sin(r21153); double r21155 = pow(r21154, r21151); double r21156 = r21152 + r21155; double r21157 = sqrt(r21156); double r21158 = r21157 / r21154; double r21159 = r21148 / r21158; double r21160 = th; double r21161 = sin(r21160); double r21162 = r21159 * r21161; return r21162; } 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 r21163, r21164, r21165, r21166, r21167, r21168, r21169, r21170, r21171, r21172, r21173, r21174, r21175; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r21163); mpfr_init(r21164); mpfr_init(r21165); mpfr_init(r21166); mpfr_init_set_str(r21167, "2", 10, MPFR_RNDN); mpfr_init(r21168); mpfr_init(r21169); mpfr_init(r21170); mpfr_init(r21171); mpfr_init(r21172); mpfr_init(r21173); mpfr_init(r21174); mpfr_init(r21175); } double f_im(double kx, double ky, double th) { mpfr_set_d(r21163, ky, MPFR_RNDN); mpfr_sin(r21164, r21163, MPFR_RNDN); mpfr_set_d(r21165, kx, MPFR_RNDN); mpfr_sin(r21166, r21165, MPFR_RNDN); ; mpfr_pow(r21168, r21166, r21167, MPFR_RNDN); mpfr_pow(r21169, r21164, r21167, MPFR_RNDN); mpfr_add(r21170, r21168, r21169, MPFR_RNDN); mpfr_sqrt(r21171, r21170, MPFR_RNDN); mpfr_div(r21172, r21164, r21171, MPFR_RNDN); mpfr_set_d(r21173, th, MPFR_RNDN); mpfr_sin(r21174, r21173, MPFR_RNDN); mpfr_mul(r21175, r21172, r21174, MPFR_RNDN); return mpfr_get_d(r21175, MPFR_RNDN); } static mpfr_t r21176, r21177, r21178, r21179, r21180, r21181, r21182, r21183, r21184, r21185, r21186, r21187, r21188, r21189, r21190; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21176, "1", 10, MPFR_RNDN); mpfr_init(r21177); mpfr_init(r21178); mpfr_init_set_str(r21179, "2", 10, MPFR_RNDN); mpfr_init(r21180); mpfr_init(r21181); mpfr_init(r21182); mpfr_init(r21183); mpfr_init(r21184); mpfr_init(r21185); mpfr_init(r21186); mpfr_init(r21187); mpfr_init(r21188); mpfr_init(r21189); mpfr_init(r21190); } double f_fm(double kx, double ky, double th) { ; mpfr_set_d(r21177, kx, MPFR_RNDN); mpfr_sin(r21178, r21177, MPFR_RNDN); ; mpfr_pow(r21180, r21178, r21179, MPFR_RNDN); mpfr_set_d(r21181, ky, MPFR_RNDN); mpfr_sin(r21182, r21181, MPFR_RNDN); mpfr_pow(r21183, r21182, r21179, MPFR_RNDN); mpfr_add(r21184, r21180, r21183, MPFR_RNDN); mpfr_sqrt(r21185, r21184, MPFR_RNDN); mpfr_div(r21186, r21185, r21182, MPFR_RNDN); mpfr_div(r21187, r21176, r21186, MPFR_RNDN); mpfr_set_d(r21188, th, MPFR_RNDN); mpfr_sin(r21189, r21188, MPFR_RNDN); mpfr_mul(r21190, r21187, r21189, MPFR_RNDN); return mpfr_get_d(r21190, MPFR_RNDN); } static mpfr_t r21191, r21192, r21193, r21194, r21195, r21196, r21197, r21198, r21199, r21200, r21201, r21202, r21203, r21204, r21205; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21191, "1", 10, MPFR_RNDN); mpfr_init(r21192); mpfr_init(r21193); mpfr_init_set_str(r21194, "2", 10, MPFR_RNDN); mpfr_init(r21195); mpfr_init(r21196); mpfr_init(r21197); mpfr_init(r21198); mpfr_init(r21199); mpfr_init(r21200); mpfr_init(r21201); mpfr_init(r21202); mpfr_init(r21203); mpfr_init(r21204); mpfr_init(r21205); } double f_dm(double kx, double ky, double th) { ; mpfr_set_d(r21192, kx, MPFR_RNDN); mpfr_sin(r21193, r21192, MPFR_RNDN); ; mpfr_pow(r21195, r21193, r21194, MPFR_RNDN); mpfr_set_d(r21196, ky, MPFR_RNDN); mpfr_sin(r21197, r21196, MPFR_RNDN); mpfr_pow(r21198, r21197, r21194, MPFR_RNDN); mpfr_add(r21199, r21195, r21198, MPFR_RNDN); mpfr_sqrt(r21200, r21199, MPFR_RNDN); mpfr_div(r21201, r21200, r21197, MPFR_RNDN); mpfr_div(r21202, r21191, r21201, MPFR_RNDN); mpfr_set_d(r21203, th, MPFR_RNDN); mpfr_sin(r21204, r21203, MPFR_RNDN); mpfr_mul(r21205, r21202, r21204, MPFR_RNDN); return mpfr_get_d(r21205, MPFR_RNDN); }