#include #include #include #include #include char *name = "Toniolo and Linder, Equation (10-)"; double f_if(float t, float l, float k) { float r25181 = 2; float r25182 = t; float r25183 = 3; float r25184 = pow(r25182, r25183); float r25185 = l; float r25186 = r25185 * r25185; float r25187 = r25184 / r25186; float r25188 = k; float r25189 = sin(r25188); float r25190 = r25187 * r25189; float r25191 = tan(r25188); float r25192 = r25190 * r25191; float r25193 = 1; float r25194 = r25188 / r25182; float r25195 = pow(r25194, r25181); float r25196 = r25193 + r25195; float r25197 = r25196 - r25193; float r25198 = r25192 * r25197; float r25199 = r25181 / r25198; return r25199; } double f_id(double t, double l, double k) { double r25200 = 2; double r25201 = t; double r25202 = 3; double r25203 = pow(r25201, r25202); double r25204 = l; double r25205 = r25204 * r25204; double r25206 = r25203 / r25205; double r25207 = k; double r25208 = sin(r25207); double r25209 = r25206 * r25208; double r25210 = tan(r25207); double r25211 = r25209 * r25210; double r25212 = 1; double r25213 = r25207 / r25201; double r25214 = pow(r25213, r25200); double r25215 = r25212 + r25214; double r25216 = r25215 - r25212; double r25217 = r25211 * r25216; double r25218 = r25200 / r25217; return r25218; } double f_of(float t, float l, float k) { float r25219 = 2; float r25220 = 1; float r25221 = r25219 / r25220; float r25222 = k; float r25223 = l; float r25224 = r25222 / r25223; float r25225 = sin(r25222); float r25226 = r25224 * r25225; float r25227 = t; float r25228 = cos(r25222); float r25229 = r25227 / r25228; float r25230 = r25226 * r25229; float r25231 = r25226 * r25230; float r25232 = r25221 / r25231; return r25232; } double f_od(double t, double l, double k) { double r25233 = 2; double r25234 = 1; double r25235 = r25233 / r25234; double r25236 = k; double r25237 = l; double r25238 = r25236 / r25237; double r25239 = sin(r25236); double r25240 = r25238 * r25239; double r25241 = t; double r25242 = cos(r25236); double r25243 = r25241 / r25242; double r25244 = r25240 * r25243; double r25245 = r25240 * r25244; double r25246 = r25235 / r25245; return r25246; } 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 r25247, r25248, r25249, r25250, r25251, r25252, r25253, r25254, r25255, r25256, r25257, r25258, r25259, r25260, r25261, r25262, r25263, r25264, r25265; void setup_mpfr_f_im() { mpfr_set_default_prec(4176); mpfr_init_set_str(r25247, "2", 10, MPFR_RNDN); mpfr_init(r25248); mpfr_init_set_str(r25249, "3", 10, MPFR_RNDN); mpfr_init(r25250); mpfr_init(r25251); mpfr_init(r25252); mpfr_init(r25253); mpfr_init(r25254); mpfr_init(r25255); mpfr_init(r25256); mpfr_init(r25257); mpfr_init(r25258); mpfr_init_set_str(r25259, "1", 10, MPFR_RNDN); mpfr_init(r25260); mpfr_init(r25261); mpfr_init(r25262); mpfr_init(r25263); mpfr_init(r25264); mpfr_init(r25265); } double f_im(double t, double l, double k) { ; mpfr_set_d(r25248, t, MPFR_RNDN); ; mpfr_pow(r25250, r25248, r25249, MPFR_RNDN); mpfr_set_d(r25251, l, MPFR_RNDN); mpfr_mul(r25252, r25251, r25251, MPFR_RNDN); mpfr_div(r25253, r25250, r25252, MPFR_RNDN); mpfr_set_d(r25254, k, MPFR_RNDN); mpfr_sin(r25255, r25254, MPFR_RNDN); mpfr_mul(r25256, r25253, r25255, MPFR_RNDN); mpfr_tan(r25257, r25254, MPFR_RNDN); mpfr_mul(r25258, r25256, r25257, MPFR_RNDN); ; mpfr_div(r25260, r25254, r25248, MPFR_RNDN); mpfr_pow(r25261, r25260, r25247, MPFR_RNDN); mpfr_add(r25262, r25259, r25261, MPFR_RNDN); mpfr_sub(r25263, r25262, r25259, MPFR_RNDN); mpfr_mul(r25264, r25258, r25263, MPFR_RNDN); mpfr_div(r25265, r25247, r25264, MPFR_RNDN); return mpfr_get_d(r25265, MPFR_RNDN); } static mpfr_t r25266, r25267, r25268, r25269, r25270, r25271, r25272, r25273, r25274, r25275, r25276, r25277, r25278, r25279; void setup_mpfr_f_fm() { mpfr_set_default_prec(4176); mpfr_init_set_str(r25266, "2", 10, MPFR_RNDN); mpfr_init_set_str(r25267, "1", 10, MPFR_RNDN); mpfr_init(r25268); mpfr_init(r25269); mpfr_init(r25270); mpfr_init(r25271); mpfr_init(r25272); mpfr_init(r25273); mpfr_init(r25274); mpfr_init(r25275); mpfr_init(r25276); mpfr_init(r25277); mpfr_init(r25278); mpfr_init(r25279); } double f_fm(double t, double l, double k) { ; ; mpfr_div(r25268, r25266, r25267, MPFR_RNDN); mpfr_set_d(r25269, k, MPFR_RNDN); mpfr_set_d(r25270, l, MPFR_RNDN); mpfr_div(r25271, r25269, r25270, MPFR_RNDN); mpfr_sin(r25272, r25269, MPFR_RNDN); mpfr_mul(r25273, r25271, r25272, MPFR_RNDN); mpfr_set_d(r25274, t, MPFR_RNDN); mpfr_cos(r25275, r25269, MPFR_RNDN); mpfr_div(r25276, r25274, r25275, MPFR_RNDN); mpfr_mul(r25277, r25273, r25276, MPFR_RNDN); mpfr_mul(r25278, r25273, r25277, MPFR_RNDN); mpfr_div(r25279, r25268, r25278, MPFR_RNDN); return mpfr_get_d(r25279, MPFR_RNDN); } static mpfr_t r25280, r25281, r25282, r25283, r25284, r25285, r25286, r25287, r25288, r25289, r25290, r25291, r25292, r25293; void setup_mpfr_f_dm() { mpfr_set_default_prec(4176); mpfr_init_set_str(r25280, "2", 10, MPFR_RNDN); mpfr_init_set_str(r25281, "1", 10, MPFR_RNDN); mpfr_init(r25282); mpfr_init(r25283); mpfr_init(r25284); mpfr_init(r25285); mpfr_init(r25286); mpfr_init(r25287); mpfr_init(r25288); mpfr_init(r25289); mpfr_init(r25290); mpfr_init(r25291); mpfr_init(r25292); mpfr_init(r25293); } double f_dm(double t, double l, double k) { ; ; mpfr_div(r25282, r25280, r25281, MPFR_RNDN); mpfr_set_d(r25283, k, MPFR_RNDN); mpfr_set_d(r25284, l, MPFR_RNDN); mpfr_div(r25285, r25283, r25284, MPFR_RNDN); mpfr_sin(r25286, r25283, MPFR_RNDN); mpfr_mul(r25287, r25285, r25286, MPFR_RNDN); mpfr_set_d(r25288, t, MPFR_RNDN); mpfr_cos(r25289, r25283, MPFR_RNDN); mpfr_div(r25290, r25288, r25289, MPFR_RNDN); mpfr_mul(r25291, r25287, r25290, MPFR_RNDN); mpfr_mul(r25292, r25287, r25291, MPFR_RNDN); mpfr_div(r25293, r25282, r25292, MPFR_RNDN); return mpfr_get_d(r25293, MPFR_RNDN); }