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