#include #include #include #include #include char *name = "Toniolo and Linder, Equation (10-)"; double f_if(float t, float l, float k) { float r25056 = 2; float r25057 = t; float r25058 = 3; float r25059 = pow(r25057, r25058); float r25060 = l; float r25061 = r25060 * r25060; float r25062 = r25059 / r25061; float r25063 = k; float r25064 = sin(r25063); float r25065 = r25062 * r25064; float r25066 = tan(r25063); float r25067 = r25065 * r25066; float r25068 = 1; float r25069 = r25063 / r25057; float r25070 = pow(r25069, r25056); float r25071 = r25068 + r25070; float r25072 = r25071 - r25068; float r25073 = r25067 * r25072; float r25074 = r25056 / r25073; return r25074; } double f_id(double t, double l, double k) { double r25075 = 2; double r25076 = t; double r25077 = 3; double r25078 = pow(r25076, r25077); double r25079 = l; double r25080 = r25079 * r25079; double r25081 = r25078 / r25080; double r25082 = k; double r25083 = sin(r25082); double r25084 = r25081 * r25083; double r25085 = tan(r25082); double r25086 = r25084 * r25085; double r25087 = 1; double r25088 = r25082 / r25076; double r25089 = pow(r25088, r25075); double r25090 = r25087 + r25089; double r25091 = r25090 - r25087; double r25092 = r25086 * r25091; double r25093 = r25075 / r25092; return r25093; } double f_of(float t, float l, float k) { float r25094 = 2; float r25095 = k; float r25096 = t; float r25097 = r25095 / r25096; float r25098 = r25097 * r25097; float r25099 = cbrt(r25098); float r25100 = l; float r25101 = r25096 / r25100; float r25102 = r25099 * r25101; float r25103 = r25102 * r25096; float r25104 = tan(r25095); float r25105 = sin(r25095); float r25106 = r25104 * r25105; float r25107 = r25101 * r25106; float r25108 = r25107 * r25099; float r25109 = r25103 * r25108; float r25110 = fabs(r25097); float r25111 = cbrt(r25110); float r25112 = r25111 * r25111; float r25113 = r25109 * r25112; float r25114 = r25094 / r25113; return r25114; } double f_od(double t, double l, double k) { double r25115 = 2; double r25116 = k; double r25117 = t; double r25118 = r25116 / r25117; double r25119 = r25118 * r25118; double r25120 = cbrt(r25119); double r25121 = l; double r25122 = r25117 / r25121; double r25123 = r25120 * r25122; double r25124 = r25123 * r25117; double r25125 = tan(r25116); double r25126 = sin(r25116); double r25127 = r25125 * r25126; double r25128 = r25122 * r25127; double r25129 = r25128 * r25120; double r25130 = r25124 * r25129; double r25131 = fabs(r25118); double r25132 = cbrt(r25131); double r25133 = r25132 * r25132; double r25134 = r25130 * r25133; double r25135 = r25115 / r25134; return r25135; } 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 r25136, r25137, r25138, r25139, r25140, r25141, r25142, r25143, r25144, r25145, r25146, r25147, r25148, r25149, r25150, r25151, r25152, r25153, r25154; void setup_mpfr_f_im() { mpfr_set_default_prec(4432); mpfr_init_set_str(r25136, "2", 10, MPFR_RNDN); mpfr_init(r25137); mpfr_init_set_str(r25138, "3", 10, MPFR_RNDN); mpfr_init(r25139); mpfr_init(r25140); mpfr_init(r25141); mpfr_init(r25142); mpfr_init(r25143); mpfr_init(r25144); mpfr_init(r25145); mpfr_init(r25146); mpfr_init(r25147); mpfr_init_set_str(r25148, "1", 10, MPFR_RNDN); mpfr_init(r25149); mpfr_init(r25150); mpfr_init(r25151); mpfr_init(r25152); mpfr_init(r25153); mpfr_init(r25154); } double f_im(double t, double l, double k) { ; mpfr_set_d(r25137, t, MPFR_RNDN); ; mpfr_pow(r25139, r25137, r25138, MPFR_RNDN); mpfr_set_d(r25140, l, MPFR_RNDN); mpfr_mul(r25141, r25140, r25140, MPFR_RNDN); mpfr_div(r25142, r25139, r25141, MPFR_RNDN); mpfr_set_d(r25143, k, MPFR_RNDN); mpfr_sin(r25144, r25143, MPFR_RNDN); mpfr_mul(r25145, r25142, r25144, MPFR_RNDN); mpfr_tan(r25146, r25143, MPFR_RNDN); mpfr_mul(r25147, r25145, r25146, MPFR_RNDN); ; mpfr_div(r25149, r25143, r25137, MPFR_RNDN); mpfr_pow(r25150, r25149, r25136, MPFR_RNDN); mpfr_add(r25151, r25148, r25150, MPFR_RNDN); mpfr_sub(r25152, r25151, r25148, MPFR_RNDN); mpfr_mul(r25153, r25147, r25152, MPFR_RNDN); mpfr_div(r25154, r25136, r25153, MPFR_RNDN); return mpfr_get_d(r25154, MPFR_RNDN); } static mpfr_t r25155, r25156, r25157, r25158, r25159, r25160, r25161, r25162, r25163, r25164, r25165, r25166, r25167, r25168, r25169, r25170, r25171, r25172, r25173, r25174, r25175; void setup_mpfr_f_fm() { mpfr_set_default_prec(4432); mpfr_init_set_str(r25155, "2", 10, MPFR_RNDN); mpfr_init(r25156); mpfr_init(r25157); mpfr_init(r25158); mpfr_init(r25159); 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); } double f_fm(double t, double l, double k) { ; mpfr_set_d(r25156, k, MPFR_RNDN); mpfr_set_d(r25157, t, MPFR_RNDN); mpfr_div(r25158, r25156, r25157, MPFR_RNDN); mpfr_mul(r25159, r25158, r25158, MPFR_RNDN); mpfr_cbrt(r25160, r25159, MPFR_RNDN); mpfr_set_d(r25161, l, MPFR_RNDN); mpfr_div(r25162, r25157, r25161, MPFR_RNDN); mpfr_mul(r25163, r25160, r25162, MPFR_RNDN); mpfr_mul(r25164, r25163, r25157, MPFR_RNDN); mpfr_tan(r25165, r25156, MPFR_RNDN); mpfr_sin(r25166, r25156, MPFR_RNDN); mpfr_mul(r25167, r25165, r25166, MPFR_RNDN); mpfr_mul(r25168, r25162, r25167, MPFR_RNDN); mpfr_mul(r25169, r25168, r25160, MPFR_RNDN); mpfr_mul(r25170, r25164, r25169, MPFR_RNDN); mpfr_abs(r25171, r25158, MPFR_RNDN); mpfr_cbrt(r25172, r25171, MPFR_RNDN); mpfr_mul(r25173, r25172, r25172, MPFR_RNDN); mpfr_mul(r25174, r25170, r25173, MPFR_RNDN); mpfr_div(r25175, r25155, r25174, MPFR_RNDN); return mpfr_get_d(r25175, MPFR_RNDN); } static mpfr_t r25176, r25177, r25178, r25179, r25180, r25181, r25182, r25183, r25184, r25185, r25186, r25187, r25188, r25189, r25190, r25191, r25192, r25193, r25194, r25195, r25196; void setup_mpfr_f_dm() { mpfr_set_default_prec(4432); mpfr_init_set_str(r25176, "2", 10, MPFR_RNDN); mpfr_init(r25177); mpfr_init(r25178); mpfr_init(r25179); mpfr_init(r25180); mpfr_init(r25181); 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); } double f_dm(double t, double l, double k) { ; mpfr_set_d(r25177, k, MPFR_RNDN); mpfr_set_d(r25178, t, MPFR_RNDN); mpfr_div(r25179, r25177, r25178, MPFR_RNDN); mpfr_mul(r25180, r25179, r25179, MPFR_RNDN); mpfr_cbrt(r25181, r25180, MPFR_RNDN); mpfr_set_d(r25182, l, MPFR_RNDN); mpfr_div(r25183, r25178, r25182, MPFR_RNDN); mpfr_mul(r25184, r25181, r25183, MPFR_RNDN); mpfr_mul(r25185, r25184, r25178, MPFR_RNDN); mpfr_tan(r25186, r25177, MPFR_RNDN); mpfr_sin(r25187, r25177, MPFR_RNDN); mpfr_mul(r25188, r25186, r25187, MPFR_RNDN); mpfr_mul(r25189, r25183, r25188, MPFR_RNDN); mpfr_mul(r25190, r25189, r25181, MPFR_RNDN); mpfr_mul(r25191, r25185, r25190, MPFR_RNDN); mpfr_abs(r25192, r25179, MPFR_RNDN); mpfr_cbrt(r25193, r25192, MPFR_RNDN); mpfr_mul(r25194, r25193, r25193, MPFR_RNDN); mpfr_mul(r25195, r25191, r25194, MPFR_RNDN); mpfr_div(r25196, r25176, r25195, MPFR_RNDN); return mpfr_get_d(r25196, MPFR_RNDN); }