#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3a)"; double f_if(float l, float Om, float kx, float ky) { float r25074 = 1; float r25075 = 2; float r25076 = r25074 / r25075; float r25077 = l; float r25078 = r25075 * r25077; float r25079 = Om; float r25080 = r25078 / r25079; float r25081 = pow(r25080, r25075); float r25082 = kx; float r25083 = sin(r25082); float r25084 = pow(r25083, r25075); float r25085 = ky; float r25086 = sin(r25085); float r25087 = pow(r25086, r25075); float r25088 = r25084 + r25087; float r25089 = r25081 * r25088; float r25090 = r25074 + r25089; float r25091 = sqrt(r25090); float r25092 = r25074 / r25091; float r25093 = r25074 + r25092; float r25094 = r25076 * r25093; float r25095 = sqrt(r25094); return r25095; } double f_id(double l, double Om, double kx, double ky) { double r25096 = 1; double r25097 = 2; double r25098 = r25096 / r25097; double r25099 = l; double r25100 = r25097 * r25099; double r25101 = Om; double r25102 = r25100 / r25101; double r25103 = pow(r25102, r25097); double r25104 = kx; double r25105 = sin(r25104); double r25106 = pow(r25105, r25097); double r25107 = ky; double r25108 = sin(r25107); double r25109 = pow(r25108, r25097); double r25110 = r25106 + r25109; double r25111 = r25103 * r25110; double r25112 = r25096 + r25111; double r25113 = sqrt(r25112); double r25114 = r25096 / r25113; double r25115 = r25096 + r25114; double r25116 = r25098 * r25115; double r25117 = sqrt(r25116); return r25117; } double f_of(float l, float Om, float kx, float ky) { float r25118 = 1; float r25119 = 2; float r25120 = r25118 / r25119; float r25121 = l; float r25122 = Om; float r25123 = r25122 / r25119; float r25124 = r25121 / r25123; float r25125 = kx; float r25126 = sin(r25125); float r25127 = ky; float r25128 = sin(r25127); float r25129 = hypot(r25126, r25128); float r25130 = r25124 * r25129; float r25131 = hypot(r25118, r25130); float r25132 = expm1(r25131); float r25133 = log1p(r25132); float r25134 = r25118 / r25133; float r25135 = r25118 + r25134; float r25136 = r25120 * r25135; float r25137 = sqrt(r25136); return r25137; } double f_od(double l, double Om, double kx, double ky) { double r25138 = 1; double r25139 = 2; double r25140 = r25138 / r25139; double r25141 = l; double r25142 = Om; double r25143 = r25142 / r25139; double r25144 = r25141 / r25143; double r25145 = kx; double r25146 = sin(r25145); double r25147 = ky; double r25148 = sin(r25147); double r25149 = hypot(r25146, r25148); double r25150 = r25144 * r25149; double r25151 = hypot(r25138, r25150); double r25152 = expm1(r25151); double r25153 = log1p(r25152); double r25154 = r25138 / r25153; double r25155 = r25138 + r25154; double r25156 = r25140 * r25155; double r25157 = sqrt(r25156); return r25157; } 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 r25158, r25159, r25160, r25161, r25162, r25163, r25164, r25165, r25166, r25167, r25168, r25169, r25170, r25171, r25172, r25173, r25174, r25175, r25176, r25177, r25178, r25179; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r25158, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25159, "2", 10, MPFR_RNDN); mpfr_init(r25160); mpfr_init(r25161); mpfr_init(r25162); mpfr_init(r25163); mpfr_init(r25164); mpfr_init(r25165); mpfr_init(r25166); mpfr_init(r25167); mpfr_init(r25168); mpfr_init(r25169); mpfr_init(r25170); mpfr_init(r25171); mpfr_init(r25172); mpfr_init(r25173); mpfr_init(r25174); mpfr_init(r25175); mpfr_init(r25176); mpfr_init(r25177); mpfr_init(r25178); mpfr_init(r25179); } double f_im(double l, double Om, double kx, double ky) { ; ; mpfr_div(r25160, r25158, r25159, MPFR_RNDN); mpfr_set_d(r25161, l, MPFR_RNDN); mpfr_mul(r25162, r25159, r25161, MPFR_RNDN); mpfr_set_d(r25163, Om, MPFR_RNDN); mpfr_div(r25164, r25162, r25163, MPFR_RNDN); mpfr_pow(r25165, r25164, r25159, MPFR_RNDN); mpfr_set_d(r25166, kx, MPFR_RNDN); mpfr_sin(r25167, r25166, MPFR_RNDN); mpfr_pow(r25168, r25167, r25159, MPFR_RNDN); mpfr_set_d(r25169, ky, MPFR_RNDN); mpfr_sin(r25170, r25169, MPFR_RNDN); mpfr_pow(r25171, r25170, r25159, MPFR_RNDN); mpfr_add(r25172, r25168, r25171, MPFR_RNDN); mpfr_mul(r25173, r25165, r25172, MPFR_RNDN); mpfr_add(r25174, r25158, r25173, MPFR_RNDN); mpfr_sqrt(r25175, r25174, MPFR_RNDN); mpfr_div(r25176, r25158, r25175, MPFR_RNDN); mpfr_add(r25177, r25158, r25176, MPFR_RNDN); mpfr_mul(r25178, r25160, r25177, MPFR_RNDN); mpfr_sqrt(r25179, r25178, MPFR_RNDN); return mpfr_get_d(r25179, MPFR_RNDN); } static mpfr_t r25180, r25181, r25182, r25183, r25184, r25185, r25186, r25187, r25188, r25189, r25190, r25191, r25192, r25193, r25194, r25195, r25196, r25197, r25198, r25199; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r25180, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25181, "2", 10, MPFR_RNDN); mpfr_init(r25182); mpfr_init(r25183); mpfr_init(r25184); mpfr_init(r25185); mpfr_init(r25186); mpfr_init(r25187); mpfr_init(r25188); mpfr_init(r25189); mpfr_init(r25190); mpfr_init(r25191); mpfr_init(r25192); mpfr_init(r25193); mpfr_init(r25194); mpfr_init(r25195); mpfr_init(r25196); mpfr_init(r25197); mpfr_init(r25198); mpfr_init(r25199); } double f_fm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r25182, r25180, r25181, MPFR_RNDN); mpfr_set_d(r25183, l, MPFR_RNDN); mpfr_set_d(r25184, Om, MPFR_RNDN); mpfr_div(r25185, r25184, r25181, MPFR_RNDN); mpfr_div(r25186, r25183, r25185, MPFR_RNDN); mpfr_set_d(r25187, kx, MPFR_RNDN); mpfr_sin(r25188, r25187, MPFR_RNDN); mpfr_set_d(r25189, ky, MPFR_RNDN); mpfr_sin(r25190, r25189, MPFR_RNDN); mpfr_hypot(r25191, r25188, r25190, MPFR_RNDN); mpfr_mul(r25192, r25186, r25191, MPFR_RNDN); mpfr_hypot(r25193, r25180, r25192, MPFR_RNDN); mpfr_expm1(r25194, r25193, MPFR_RNDN); mpfr_log1p(r25195, r25194, MPFR_RNDN); mpfr_div(r25196, r25180, r25195, MPFR_RNDN); mpfr_add(r25197, r25180, r25196, MPFR_RNDN); mpfr_mul(r25198, r25182, r25197, MPFR_RNDN); mpfr_sqrt(r25199, r25198, MPFR_RNDN); return mpfr_get_d(r25199, MPFR_RNDN); } static mpfr_t r25200, r25201, r25202, r25203, r25204, r25205, r25206, r25207, r25208, r25209, r25210, r25211, r25212, r25213, r25214, r25215, r25216, r25217, r25218, r25219; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r25200, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25201, "2", 10, MPFR_RNDN); mpfr_init(r25202); mpfr_init(r25203); mpfr_init(r25204); mpfr_init(r25205); mpfr_init(r25206); mpfr_init(r25207); mpfr_init(r25208); mpfr_init(r25209); mpfr_init(r25210); mpfr_init(r25211); mpfr_init(r25212); mpfr_init(r25213); mpfr_init(r25214); mpfr_init(r25215); mpfr_init(r25216); mpfr_init(r25217); mpfr_init(r25218); mpfr_init(r25219); } double f_dm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r25202, r25200, r25201, MPFR_RNDN); mpfr_set_d(r25203, l, MPFR_RNDN); mpfr_set_d(r25204, Om, MPFR_RNDN); mpfr_div(r25205, r25204, r25201, MPFR_RNDN); mpfr_div(r25206, r25203, r25205, MPFR_RNDN); mpfr_set_d(r25207, kx, MPFR_RNDN); mpfr_sin(r25208, r25207, MPFR_RNDN); mpfr_set_d(r25209, ky, MPFR_RNDN); mpfr_sin(r25210, r25209, MPFR_RNDN); mpfr_hypot(r25211, r25208, r25210, MPFR_RNDN); mpfr_mul(r25212, r25206, r25211, MPFR_RNDN); mpfr_hypot(r25213, r25200, r25212, MPFR_RNDN); mpfr_expm1(r25214, r25213, MPFR_RNDN); mpfr_log1p(r25215, r25214, MPFR_RNDN); mpfr_div(r25216, r25200, r25215, MPFR_RNDN); mpfr_add(r25217, r25200, r25216, MPFR_RNDN); mpfr_mul(r25218, r25202, r25217, MPFR_RNDN); mpfr_sqrt(r25219, r25218, MPFR_RNDN); return mpfr_get_d(r25219, MPFR_RNDN); }