#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3b), real"; double f_if(float kx, float ky, float th) { float r24086 = ky; float r24087 = sin(r24086); float r24088 = kx; float r24089 = sin(r24088); float r24090 = 2; float r24091 = pow(r24089, r24090); float r24092 = pow(r24087, r24090); float r24093 = r24091 + r24092; float r24094 = sqrt(r24093); float r24095 = r24087 / r24094; float r24096 = th; float r24097 = sin(r24096); float r24098 = r24095 * r24097; return r24098; } double f_id(double kx, double ky, double th) { double r24099 = ky; double r24100 = sin(r24099); double r24101 = kx; double r24102 = sin(r24101); double r24103 = 2; double r24104 = pow(r24102, r24103); double r24105 = pow(r24100, r24103); double r24106 = r24104 + r24105; double r24107 = sqrt(r24106); double r24108 = r24100 / r24107; double r24109 = th; double r24110 = sin(r24109); double r24111 = r24108 * r24110; return r24111; } double f_of(float kx, float ky, float th) { float r24112 = ky; float r24113 = sin(r24112); float r24114 = kx; float r24115 = sin(r24114); float r24116 = 2; float r24117 = pow(r24115, r24116); float r24118 = cbrt(r24117); float r24119 = r24118 * r24118; float r24120 = r24119 * r24118; float r24121 = pow(r24113, r24116); float r24122 = r24120 + r24121; float r24123 = sqrt(r24122); float r24124 = r24113 / r24123; float r24125 = th; float r24126 = sin(r24125); float r24127 = r24124 * r24126; return r24127; } double f_od(double kx, double ky, double th) { double r24128 = ky; double r24129 = sin(r24128); double r24130 = kx; double r24131 = sin(r24130); double r24132 = 2; double r24133 = pow(r24131, r24132); double r24134 = cbrt(r24133); double r24135 = r24134 * r24134; double r24136 = r24135 * r24134; double r24137 = pow(r24129, r24132); double r24138 = r24136 + r24137; double r24139 = sqrt(r24138); double r24140 = r24129 / r24139; double r24141 = th; double r24142 = sin(r24141); double r24143 = r24140 * r24142; return r24143; } 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 r24144, r24145, r24146, r24147, r24148, r24149, r24150, r24151, r24152, r24153, r24154, r24155, r24156; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r24144); mpfr_init(r24145); mpfr_init(r24146); mpfr_init(r24147); mpfr_init_set_str(r24148, "2", 10, MPFR_RNDN); mpfr_init(r24149); mpfr_init(r24150); mpfr_init(r24151); mpfr_init(r24152); mpfr_init(r24153); mpfr_init(r24154); mpfr_init(r24155); mpfr_init(r24156); } double f_im(double kx, double ky, double th) { mpfr_set_d(r24144, ky, MPFR_RNDN); mpfr_sin(r24145, r24144, MPFR_RNDN); mpfr_set_d(r24146, kx, MPFR_RNDN); mpfr_sin(r24147, r24146, MPFR_RNDN); ; mpfr_pow(r24149, r24147, r24148, MPFR_RNDN); mpfr_pow(r24150, r24145, r24148, MPFR_RNDN); mpfr_add(r24151, r24149, r24150, MPFR_RNDN); mpfr_sqrt(r24152, r24151, MPFR_RNDN); mpfr_div(r24153, r24145, r24152, MPFR_RNDN); mpfr_set_d(r24154, th, MPFR_RNDN); mpfr_sin(r24155, r24154, MPFR_RNDN); mpfr_mul(r24156, r24153, r24155, MPFR_RNDN); return mpfr_get_d(r24156, MPFR_RNDN); } static mpfr_t r24157, r24158, r24159, r24160, r24161, r24162, r24163, r24164, r24165, r24166, r24167, r24168, r24169, r24170, r24171, r24172; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r24157); mpfr_init(r24158); mpfr_init(r24159); mpfr_init(r24160); mpfr_init_set_str(r24161, "2", 10, MPFR_RNDN); mpfr_init(r24162); mpfr_init(r24163); mpfr_init(r24164); mpfr_init(r24165); mpfr_init(r24166); mpfr_init(r24167); mpfr_init(r24168); mpfr_init(r24169); mpfr_init(r24170); mpfr_init(r24171); mpfr_init(r24172); } double f_fm(double kx, double ky, double th) { mpfr_set_d(r24157, ky, MPFR_RNDN); mpfr_sin(r24158, r24157, MPFR_RNDN); mpfr_set_d(r24159, kx, MPFR_RNDN); mpfr_sin(r24160, r24159, MPFR_RNDN); ; mpfr_pow(r24162, r24160, r24161, MPFR_RNDN); mpfr_cbrt(r24163, r24162, MPFR_RNDN); mpfr_mul(r24164, r24163, r24163, MPFR_RNDN); mpfr_mul(r24165, r24164, r24163, MPFR_RNDN); mpfr_pow(r24166, r24158, r24161, MPFR_RNDN); mpfr_add(r24167, r24165, r24166, MPFR_RNDN); mpfr_sqrt(r24168, r24167, MPFR_RNDN); mpfr_div(r24169, r24158, r24168, MPFR_RNDN); mpfr_set_d(r24170, th, MPFR_RNDN); mpfr_sin(r24171, r24170, MPFR_RNDN); mpfr_mul(r24172, r24169, r24171, MPFR_RNDN); return mpfr_get_d(r24172, MPFR_RNDN); } static mpfr_t r24173, r24174, r24175, r24176, r24177, r24178, r24179, r24180, r24181, r24182, r24183, r24184, r24185, r24186, r24187, r24188; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r24173); mpfr_init(r24174); mpfr_init(r24175); mpfr_init(r24176); mpfr_init_set_str(r24177, "2", 10, MPFR_RNDN); mpfr_init(r24178); mpfr_init(r24179); mpfr_init(r24180); mpfr_init(r24181); mpfr_init(r24182); mpfr_init(r24183); mpfr_init(r24184); mpfr_init(r24185); mpfr_init(r24186); mpfr_init(r24187); mpfr_init(r24188); } double f_dm(double kx, double ky, double th) { mpfr_set_d(r24173, ky, MPFR_RNDN); mpfr_sin(r24174, r24173, MPFR_RNDN); mpfr_set_d(r24175, kx, MPFR_RNDN); mpfr_sin(r24176, r24175, MPFR_RNDN); ; mpfr_pow(r24178, r24176, r24177, MPFR_RNDN); mpfr_cbrt(r24179, r24178, MPFR_RNDN); mpfr_mul(r24180, r24179, r24179, MPFR_RNDN); mpfr_mul(r24181, r24180, r24179, MPFR_RNDN); mpfr_pow(r24182, r24174, r24177, MPFR_RNDN); mpfr_add(r24183, r24181, r24182, MPFR_RNDN); mpfr_sqrt(r24184, r24183, MPFR_RNDN); mpfr_div(r24185, r24174, r24184, MPFR_RNDN); mpfr_set_d(r24186, th, MPFR_RNDN); mpfr_sin(r24187, r24186, MPFR_RNDN); mpfr_mul(r24188, r24185, r24187, MPFR_RNDN); return mpfr_get_d(r24188, MPFR_RNDN); }