#include #include #include #include #include char *name = "Toniolo and Linder, Equation (10-)"; double f_if(float t, float l, float k) { float r25245 = 2; float r25246 = t; float r25247 = 3; float r25248 = pow(r25246, r25247); float r25249 = l; float r25250 = r25249 * r25249; float r25251 = r25248 / r25250; float r25252 = k; float r25253 = sin(r25252); float r25254 = r25251 * r25253; float r25255 = tan(r25252); float r25256 = r25254 * r25255; float r25257 = 1; float r25258 = r25252 / r25246; float r25259 = pow(r25258, r25245); float r25260 = r25257 + r25259; float r25261 = r25260 - r25257; float r25262 = r25256 * r25261; float r25263 = r25245 / r25262; return r25263; } double f_id(double t, double l, double k) { double r25264 = 2; double r25265 = t; double r25266 = 3; double r25267 = pow(r25265, r25266); double r25268 = l; double r25269 = r25268 * r25268; double r25270 = r25267 / r25269; double r25271 = k; double r25272 = sin(r25271); double r25273 = r25270 * r25272; double r25274 = tan(r25271); double r25275 = r25273 * r25274; double r25276 = 1; double r25277 = r25271 / r25265; double r25278 = pow(r25277, r25264); double r25279 = r25276 + r25278; double r25280 = r25279 - r25276; double r25281 = r25275 * r25280; double r25282 = r25264 / r25281; return r25282; } double f_of(float t, float l, float k) { float r25283 = 2; float r25284 = k; float r25285 = t; float r25286 = r25284 / r25285; float r25287 = r25286 * r25286; float r25288 = cbrt(r25287); float r25289 = l; float r25290 = r25285 / r25289; float r25291 = r25288 * r25290; float r25292 = r25291 * r25285; float r25293 = tan(r25284); float r25294 = sin(r25284); float r25295 = r25293 * r25294; float r25296 = r25290 * r25295; float r25297 = r25296 * r25288; float r25298 = r25292 * r25297; float r25299 = 1; float r25300 = cbrt(r25299); float r25301 = r25300 * r25288; float r25302 = r25298 * r25301; float r25303 = r25283 / r25302; return r25303; } double f_od(double t, double l, double k) { double r25304 = 2; double r25305 = k; double r25306 = t; double r25307 = r25305 / r25306; double r25308 = r25307 * r25307; double r25309 = cbrt(r25308); double r25310 = l; double r25311 = r25306 / r25310; double r25312 = r25309 * r25311; double r25313 = r25312 * r25306; double r25314 = tan(r25305); double r25315 = sin(r25305); double r25316 = r25314 * r25315; double r25317 = r25311 * r25316; double r25318 = r25317 * r25309; double r25319 = r25313 * r25318; double r25320 = 1; double r25321 = cbrt(r25320); double r25322 = r25321 * r25309; double r25323 = r25319 * r25322; double r25324 = r25304 / r25323; return r25324; } 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 r25325, r25326, r25327, r25328, r25329, r25330, r25331, r25332, r25333, r25334, r25335, r25336, r25337, r25338, r25339, r25340, r25341, r25342, r25343; void setup_mpfr_f_im() { mpfr_set_default_prec(4176); mpfr_init_set_str(r25325, "2", 10, MPFR_RNDN); mpfr_init(r25326); mpfr_init_set_str(r25327, "3", 10, MPFR_RNDN); mpfr_init(r25328); mpfr_init(r25329); mpfr_init(r25330); mpfr_init(r25331); mpfr_init(r25332); mpfr_init(r25333); mpfr_init(r25334); mpfr_init(r25335); mpfr_init(r25336); mpfr_init_set_str(r25337, "1", 10, MPFR_RNDN); mpfr_init(r25338); mpfr_init(r25339); mpfr_init(r25340); mpfr_init(r25341); mpfr_init(r25342); mpfr_init(r25343); } double f_im(double t, double l, double k) { ; mpfr_set_d(r25326, t, MPFR_RNDN); ; mpfr_pow(r25328, r25326, r25327, MPFR_RNDN); mpfr_set_d(r25329, l, MPFR_RNDN); mpfr_mul(r25330, r25329, r25329, MPFR_RNDN); mpfr_div(r25331, r25328, r25330, MPFR_RNDN); mpfr_set_d(r25332, k, MPFR_RNDN); mpfr_sin(r25333, r25332, MPFR_RNDN); mpfr_mul(r25334, r25331, r25333, MPFR_RNDN); mpfr_tan(r25335, r25332, MPFR_RNDN); mpfr_mul(r25336, r25334, r25335, MPFR_RNDN); ; mpfr_div(r25338, r25332, r25326, MPFR_RNDN); mpfr_pow(r25339, r25338, r25325, MPFR_RNDN); mpfr_add(r25340, r25337, r25339, MPFR_RNDN); mpfr_sub(r25341, r25340, r25337, MPFR_RNDN); mpfr_mul(r25342, r25336, r25341, MPFR_RNDN); mpfr_div(r25343, r25325, r25342, MPFR_RNDN); return mpfr_get_d(r25343, MPFR_RNDN); } static mpfr_t r25344, r25345, r25346, r25347, r25348, r25349, r25350, r25351, r25352, r25353, r25354, r25355, r25356, r25357, r25358, r25359, r25360, r25361, r25362, r25363, r25364; void setup_mpfr_f_fm() { mpfr_set_default_prec(4176); mpfr_init_set_str(r25344, "2", 10, MPFR_RNDN); mpfr_init(r25345); mpfr_init(r25346); mpfr_init(r25347); mpfr_init(r25348); mpfr_init(r25349); mpfr_init(r25350); mpfr_init(r25351); mpfr_init(r25352); mpfr_init(r25353); mpfr_init(r25354); mpfr_init(r25355); mpfr_init(r25356); mpfr_init(r25357); mpfr_init(r25358); mpfr_init(r25359); mpfr_init_set_str(r25360, "1", 10, MPFR_RNDN); mpfr_init(r25361); mpfr_init(r25362); mpfr_init(r25363); mpfr_init(r25364); } double f_fm(double t, double l, double k) { ; mpfr_set_d(r25345, k, MPFR_RNDN); mpfr_set_d(r25346, t, MPFR_RNDN); mpfr_div(r25347, r25345, r25346, MPFR_RNDN); mpfr_mul(r25348, r25347, r25347, MPFR_RNDN); mpfr_cbrt(r25349, r25348, MPFR_RNDN); mpfr_set_d(r25350, l, MPFR_RNDN); mpfr_div(r25351, r25346, r25350, MPFR_RNDN); mpfr_mul(r25352, r25349, r25351, MPFR_RNDN); mpfr_mul(r25353, r25352, r25346, MPFR_RNDN); mpfr_tan(r25354, r25345, MPFR_RNDN); mpfr_sin(r25355, r25345, MPFR_RNDN); mpfr_mul(r25356, r25354, r25355, MPFR_RNDN); mpfr_mul(r25357, r25351, r25356, MPFR_RNDN); mpfr_mul(r25358, r25357, r25349, MPFR_RNDN); mpfr_mul(r25359, r25353, r25358, MPFR_RNDN); ; mpfr_cbrt(r25361, r25360, MPFR_RNDN); mpfr_mul(r25362, r25361, r25349, MPFR_RNDN); mpfr_mul(r25363, r25359, r25362, MPFR_RNDN); mpfr_div(r25364, r25344, r25363, MPFR_RNDN); return mpfr_get_d(r25364, MPFR_RNDN); } static mpfr_t r25365, r25366, r25367, r25368, r25369, r25370, r25371, r25372, r25373, r25374, r25375, r25376, r25377, r25378, r25379, r25380, r25381, r25382, r25383, r25384, r25385; void setup_mpfr_f_dm() { mpfr_set_default_prec(4176); mpfr_init_set_str(r25365, "2", 10, MPFR_RNDN); mpfr_init(r25366); mpfr_init(r25367); mpfr_init(r25368); mpfr_init(r25369); mpfr_init(r25370); mpfr_init(r25371); mpfr_init(r25372); mpfr_init(r25373); mpfr_init(r25374); mpfr_init(r25375); mpfr_init(r25376); mpfr_init(r25377); mpfr_init(r25378); mpfr_init(r25379); mpfr_init(r25380); mpfr_init_set_str(r25381, "1", 10, MPFR_RNDN); mpfr_init(r25382); mpfr_init(r25383); mpfr_init(r25384); mpfr_init(r25385); } double f_dm(double t, double l, double k) { ; mpfr_set_d(r25366, k, MPFR_RNDN); mpfr_set_d(r25367, t, MPFR_RNDN); mpfr_div(r25368, r25366, r25367, MPFR_RNDN); mpfr_mul(r25369, r25368, r25368, MPFR_RNDN); mpfr_cbrt(r25370, r25369, MPFR_RNDN); mpfr_set_d(r25371, l, MPFR_RNDN); mpfr_div(r25372, r25367, r25371, MPFR_RNDN); mpfr_mul(r25373, r25370, r25372, MPFR_RNDN); mpfr_mul(r25374, r25373, r25367, MPFR_RNDN); mpfr_tan(r25375, r25366, MPFR_RNDN); mpfr_sin(r25376, r25366, MPFR_RNDN); mpfr_mul(r25377, r25375, r25376, MPFR_RNDN); mpfr_mul(r25378, r25372, r25377, MPFR_RNDN); mpfr_mul(r25379, r25378, r25370, MPFR_RNDN); mpfr_mul(r25380, r25374, r25379, MPFR_RNDN); ; mpfr_cbrt(r25382, r25381, MPFR_RNDN); mpfr_mul(r25383, r25382, r25370, MPFR_RNDN); mpfr_mul(r25384, r25380, r25383, MPFR_RNDN); mpfr_div(r25385, r25365, r25384, MPFR_RNDN); return mpfr_get_d(r25385, MPFR_RNDN); }