#include #include #include #include #include char *name = "Toniolo and Linder, Equation (10-)"; double f_if(float t, float l, float k) { float r25152 = 2; float r25153 = t; float r25154 = 3; float r25155 = pow(r25153, r25154); float r25156 = l; float r25157 = r25156 * r25156; float r25158 = r25155 / r25157; float r25159 = k; float r25160 = sin(r25159); float r25161 = r25158 * r25160; float r25162 = tan(r25159); float r25163 = r25161 * r25162; float r25164 = 1; float r25165 = r25159 / r25153; float r25166 = pow(r25165, r25152); float r25167 = r25164 + r25166; float r25168 = r25167 - r25164; float r25169 = r25163 * r25168; float r25170 = r25152 / r25169; return r25170; } double f_id(double t, double l, double k) { double r25171 = 2; double r25172 = t; double r25173 = 3; double r25174 = pow(r25172, r25173); double r25175 = l; double r25176 = r25175 * r25175; double r25177 = r25174 / r25176; double r25178 = k; double r25179 = sin(r25178); double r25180 = r25177 * r25179; double r25181 = tan(r25178); double r25182 = r25180 * r25181; double r25183 = 1; double r25184 = r25178 / r25172; double r25185 = pow(r25184, r25171); double r25186 = r25183 + r25185; double r25187 = r25186 - r25183; double r25188 = r25182 * r25187; double r25189 = r25171 / r25188; return r25189; } double f_of(float t, float l, float k) { float r25190 = 2; float r25191 = k; float r25192 = sin(r25191); float r25193 = l; float r25194 = r25193 / r25191; float r25195 = r25192 / r25194; float r25196 = r25190 / r25195; float r25197 = r25191 / r25193; float r25198 = r25197 * r25192; float r25199 = t; float r25200 = cos(r25191); float r25201 = r25199 / r25200; float r25202 = r25198 * r25201; float r25203 = r25196 / r25202; return r25203; } double f_od(double t, double l, double k) { double r25204 = 2; double r25205 = k; double r25206 = sin(r25205); double r25207 = l; double r25208 = r25207 / r25205; double r25209 = r25206 / r25208; double r25210 = r25204 / r25209; double r25211 = r25205 / r25207; double r25212 = r25211 * r25206; double r25213 = t; double r25214 = cos(r25205); double r25215 = r25213 / r25214; double r25216 = r25212 * r25215; double r25217 = r25210 / r25216; return r25217; } 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 r25218, r25219, r25220, r25221, r25222, r25223, r25224, r25225, r25226, r25227, r25228, r25229, r25230, r25231, r25232, r25233, r25234, r25235, r25236; void setup_mpfr_f_im() { mpfr_set_default_prec(4432); mpfr_init_set_str(r25218, "2", 10, MPFR_RNDN); mpfr_init(r25219); mpfr_init_set_str(r25220, "3", 10, MPFR_RNDN); mpfr_init(r25221); mpfr_init(r25222); mpfr_init(r25223); mpfr_init(r25224); mpfr_init(r25225); mpfr_init(r25226); mpfr_init(r25227); mpfr_init(r25228); mpfr_init(r25229); mpfr_init_set_str(r25230, "1", 10, MPFR_RNDN); mpfr_init(r25231); mpfr_init(r25232); mpfr_init(r25233); mpfr_init(r25234); mpfr_init(r25235); mpfr_init(r25236); } double f_im(double t, double l, double k) { ; mpfr_set_d(r25219, t, MPFR_RNDN); ; mpfr_pow(r25221, r25219, r25220, MPFR_RNDN); mpfr_set_d(r25222, l, MPFR_RNDN); mpfr_mul(r25223, r25222, r25222, MPFR_RNDN); mpfr_div(r25224, r25221, r25223, MPFR_RNDN); mpfr_set_d(r25225, k, MPFR_RNDN); mpfr_sin(r25226, r25225, MPFR_RNDN); mpfr_mul(r25227, r25224, r25226, MPFR_RNDN); mpfr_tan(r25228, r25225, MPFR_RNDN); mpfr_mul(r25229, r25227, r25228, MPFR_RNDN); ; mpfr_div(r25231, r25225, r25219, MPFR_RNDN); mpfr_pow(r25232, r25231, r25218, MPFR_RNDN); mpfr_add(r25233, r25230, r25232, MPFR_RNDN); mpfr_sub(r25234, r25233, r25230, MPFR_RNDN); mpfr_mul(r25235, r25229, r25234, MPFR_RNDN); mpfr_div(r25236, r25218, r25235, MPFR_RNDN); return mpfr_get_d(r25236, MPFR_RNDN); } static mpfr_t r25237, r25238, r25239, r25240, r25241, r25242, r25243, r25244, r25245, r25246, r25247, r25248, r25249, r25250; void setup_mpfr_f_fm() { mpfr_set_default_prec(4432); mpfr_init_set_str(r25237, "2", 10, MPFR_RNDN); mpfr_init(r25238); mpfr_init(r25239); mpfr_init(r25240); mpfr_init(r25241); mpfr_init(r25242); mpfr_init(r25243); mpfr_init(r25244); mpfr_init(r25245); mpfr_init(r25246); mpfr_init(r25247); mpfr_init(r25248); mpfr_init(r25249); mpfr_init(r25250); } double f_fm(double t, double l, double k) { ; mpfr_set_d(r25238, k, MPFR_RNDN); mpfr_sin(r25239, r25238, MPFR_RNDN); mpfr_set_d(r25240, l, MPFR_RNDN); mpfr_div(r25241, r25240, r25238, MPFR_RNDN); mpfr_div(r25242, r25239, r25241, MPFR_RNDN); mpfr_div(r25243, r25237, r25242, MPFR_RNDN); mpfr_div(r25244, r25238, r25240, MPFR_RNDN); mpfr_mul(r25245, r25244, r25239, MPFR_RNDN); mpfr_set_d(r25246, t, MPFR_RNDN); mpfr_cos(r25247, r25238, MPFR_RNDN); mpfr_div(r25248, r25246, r25247, MPFR_RNDN); mpfr_mul(r25249, r25245, r25248, MPFR_RNDN); mpfr_div(r25250, r25243, r25249, MPFR_RNDN); return mpfr_get_d(r25250, MPFR_RNDN); } static mpfr_t r25251, r25252, r25253, r25254, r25255, r25256, r25257, r25258, r25259, r25260, r25261, r25262, r25263, r25264; void setup_mpfr_f_dm() { mpfr_set_default_prec(4432); mpfr_init_set_str(r25251, "2", 10, MPFR_RNDN); mpfr_init(r25252); mpfr_init(r25253); mpfr_init(r25254); mpfr_init(r25255); mpfr_init(r25256); mpfr_init(r25257); mpfr_init(r25258); mpfr_init(r25259); mpfr_init(r25260); mpfr_init(r25261); mpfr_init(r25262); mpfr_init(r25263); mpfr_init(r25264); } double f_dm(double t, double l, double k) { ; mpfr_set_d(r25252, k, MPFR_RNDN); mpfr_sin(r25253, r25252, MPFR_RNDN); mpfr_set_d(r25254, l, MPFR_RNDN); mpfr_div(r25255, r25254, r25252, MPFR_RNDN); mpfr_div(r25256, r25253, r25255, MPFR_RNDN); mpfr_div(r25257, r25251, r25256, MPFR_RNDN); mpfr_div(r25258, r25252, r25254, MPFR_RNDN); mpfr_mul(r25259, r25258, r25253, MPFR_RNDN); mpfr_set_d(r25260, t, MPFR_RNDN); mpfr_cos(r25261, r25252, MPFR_RNDN); mpfr_div(r25262, r25260, r25261, MPFR_RNDN); mpfr_mul(r25263, r25259, r25262, MPFR_RNDN); mpfr_div(r25264, r25257, r25263, MPFR_RNDN); return mpfr_get_d(r25264, MPFR_RNDN); }