#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3b), real"; double f_if(float kx, float ky, float th) { float r24134 = ky; float r24135 = sin(r24134); float r24136 = kx; float r24137 = sin(r24136); float r24138 = 2; float r24139 = pow(r24137, r24138); float r24140 = pow(r24135, r24138); float r24141 = r24139 + r24140; float r24142 = sqrt(r24141); float r24143 = r24135 / r24142; float r24144 = th; float r24145 = sin(r24144); float r24146 = r24143 * r24145; return r24146; } double f_id(double kx, double ky, double th) { double r24147 = ky; double r24148 = sin(r24147); double r24149 = kx; double r24150 = sin(r24149); double r24151 = 2; double r24152 = pow(r24150, r24151); double r24153 = pow(r24148, r24151); double r24154 = r24152 + r24153; double r24155 = sqrt(r24154); double r24156 = r24148 / r24155; double r24157 = th; double r24158 = sin(r24157); double r24159 = r24156 * r24158; return r24159; } double f_of(float kx, float ky, float th) { float r24160 = ky; float r24161 = sin(r24160); float r24162 = kx; float r24163 = sin(r24162); float r24164 = 2; float r24165 = pow(r24163, r24164); float r24166 = cbrt(r24165); float r24167 = r24166 * r24166; float r24168 = r24167 * r24166; float r24169 = pow(r24161, r24164); float r24170 = r24168 + r24169; float r24171 = sqrt(r24170); float r24172 = r24161 / r24171; float r24173 = th; float r24174 = sin(r24173); float r24175 = r24172 * r24174; return r24175; } double f_od(double kx, double ky, double th) { double r24176 = ky; double r24177 = sin(r24176); double r24178 = kx; double r24179 = sin(r24178); double r24180 = 2; double r24181 = pow(r24179, r24180); double r24182 = cbrt(r24181); double r24183 = r24182 * r24182; double r24184 = r24183 * r24182; double r24185 = pow(r24177, r24180); double r24186 = r24184 + r24185; double r24187 = sqrt(r24186); double r24188 = r24177 / r24187; double r24189 = th; double r24190 = sin(r24189); double r24191 = r24188 * r24190; return r24191; } 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 r24192, r24193, r24194, r24195, r24196, r24197, r24198, r24199, r24200, r24201, r24202, r24203, r24204; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r24192); mpfr_init(r24193); mpfr_init(r24194); mpfr_init(r24195); mpfr_init_set_str(r24196, "2", 10, MPFR_RNDN); mpfr_init(r24197); mpfr_init(r24198); mpfr_init(r24199); mpfr_init(r24200); mpfr_init(r24201); mpfr_init(r24202); mpfr_init(r24203); mpfr_init(r24204); } double f_im(double kx, double ky, double th) { mpfr_set_d(r24192, ky, MPFR_RNDN); mpfr_sin(r24193, r24192, MPFR_RNDN); mpfr_set_d(r24194, kx, MPFR_RNDN); mpfr_sin(r24195, r24194, MPFR_RNDN); ; mpfr_pow(r24197, r24195, r24196, MPFR_RNDN); mpfr_pow(r24198, r24193, r24196, MPFR_RNDN); mpfr_add(r24199, r24197, r24198, MPFR_RNDN); mpfr_sqrt(r24200, r24199, MPFR_RNDN); mpfr_div(r24201, r24193, r24200, MPFR_RNDN); mpfr_set_d(r24202, th, MPFR_RNDN); mpfr_sin(r24203, r24202, MPFR_RNDN); mpfr_mul(r24204, r24201, r24203, MPFR_RNDN); return mpfr_get_d(r24204, MPFR_RNDN); } static mpfr_t r24205, r24206, r24207, r24208, r24209, r24210, r24211, r24212, r24213, r24214, r24215, r24216, r24217, r24218, r24219, r24220; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r24205); mpfr_init(r24206); mpfr_init(r24207); mpfr_init(r24208); mpfr_init_set_str(r24209, "2", 10, MPFR_RNDN); mpfr_init(r24210); mpfr_init(r24211); mpfr_init(r24212); mpfr_init(r24213); mpfr_init(r24214); mpfr_init(r24215); mpfr_init(r24216); mpfr_init(r24217); mpfr_init(r24218); mpfr_init(r24219); mpfr_init(r24220); } double f_fm(double kx, double ky, double th) { mpfr_set_d(r24205, ky, MPFR_RNDN); mpfr_sin(r24206, r24205, MPFR_RNDN); mpfr_set_d(r24207, kx, MPFR_RNDN); mpfr_sin(r24208, r24207, MPFR_RNDN); ; mpfr_pow(r24210, r24208, r24209, MPFR_RNDN); mpfr_cbrt(r24211, r24210, MPFR_RNDN); mpfr_mul(r24212, r24211, r24211, MPFR_RNDN); mpfr_mul(r24213, r24212, r24211, MPFR_RNDN); mpfr_pow(r24214, r24206, r24209, MPFR_RNDN); mpfr_add(r24215, r24213, r24214, MPFR_RNDN); mpfr_sqrt(r24216, r24215, MPFR_RNDN); mpfr_div(r24217, r24206, r24216, MPFR_RNDN); mpfr_set_d(r24218, th, MPFR_RNDN); mpfr_sin(r24219, r24218, MPFR_RNDN); mpfr_mul(r24220, r24217, r24219, MPFR_RNDN); return mpfr_get_d(r24220, MPFR_RNDN); } static mpfr_t r24221, r24222, r24223, r24224, r24225, r24226, r24227, r24228, r24229, r24230, r24231, r24232, r24233, r24234, r24235, r24236; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r24221); mpfr_init(r24222); mpfr_init(r24223); mpfr_init(r24224); mpfr_init_set_str(r24225, "2", 10, MPFR_RNDN); mpfr_init(r24226); mpfr_init(r24227); mpfr_init(r24228); mpfr_init(r24229); mpfr_init(r24230); mpfr_init(r24231); mpfr_init(r24232); mpfr_init(r24233); mpfr_init(r24234); mpfr_init(r24235); mpfr_init(r24236); } double f_dm(double kx, double ky, double th) { mpfr_set_d(r24221, ky, MPFR_RNDN); mpfr_sin(r24222, r24221, MPFR_RNDN); mpfr_set_d(r24223, kx, MPFR_RNDN); mpfr_sin(r24224, r24223, MPFR_RNDN); ; mpfr_pow(r24226, r24224, r24225, MPFR_RNDN); mpfr_cbrt(r24227, r24226, MPFR_RNDN); mpfr_mul(r24228, r24227, r24227, MPFR_RNDN); mpfr_mul(r24229, r24228, r24227, MPFR_RNDN); mpfr_pow(r24230, r24222, r24225, MPFR_RNDN); mpfr_add(r24231, r24229, r24230, MPFR_RNDN); mpfr_sqrt(r24232, r24231, MPFR_RNDN); mpfr_div(r24233, r24222, r24232, MPFR_RNDN); mpfr_set_d(r24234, th, MPFR_RNDN); mpfr_sin(r24235, r24234, MPFR_RNDN); mpfr_mul(r24236, r24233, r24235, MPFR_RNDN); return mpfr_get_d(r24236, MPFR_RNDN); }