#include #include #include #include #include char *name = "Toniolo and Linder, Equation (10-)"; double f_if(float t, float l, float k) { float r25028 = 2; float r25029 = t; float r25030 = 3; float r25031 = pow(r25029, r25030); float r25032 = l; float r25033 = r25032 * r25032; float r25034 = r25031 / r25033; float r25035 = k; float r25036 = sin(r25035); float r25037 = r25034 * r25036; float r25038 = tan(r25035); float r25039 = r25037 * r25038; float r25040 = 1; float r25041 = r25035 / r25029; float r25042 = pow(r25041, r25028); float r25043 = r25040 + r25042; float r25044 = r25043 - r25040; float r25045 = r25039 * r25044; float r25046 = r25028 / r25045; return r25046; } double f_id(double t, double l, double k) { double r25047 = 2; double r25048 = t; double r25049 = 3; double r25050 = pow(r25048, r25049); double r25051 = l; double r25052 = r25051 * r25051; double r25053 = r25050 / r25052; double r25054 = k; double r25055 = sin(r25054); double r25056 = r25053 * r25055; double r25057 = tan(r25054); double r25058 = r25056 * r25057; double r25059 = 1; double r25060 = r25054 / r25048; double r25061 = pow(r25060, r25047); double r25062 = r25059 + r25061; double r25063 = r25062 - r25059; double r25064 = r25058 * r25063; double r25065 = r25047 / r25064; return r25065; } double f_of(float t, float l, float k) { float r25066 = 2; float r25067 = 1; float r25068 = r25066 / r25067; float r25069 = k; float r25070 = l; float r25071 = r25069 / r25070; float r25072 = sin(r25069); float r25073 = r25071 * r25072; float r25074 = t; float r25075 = cos(r25069); float r25076 = r25074 / r25075; float r25077 = r25073 * r25076; float r25078 = r25073 * r25077; float r25079 = r25068 / r25078; return r25079; } double f_od(double t, double l, double k) { double r25080 = 2; double r25081 = 1; double r25082 = r25080 / r25081; double r25083 = k; double r25084 = l; double r25085 = r25083 / r25084; double r25086 = sin(r25083); double r25087 = r25085 * r25086; double r25088 = t; double r25089 = cos(r25083); double r25090 = r25088 / r25089; double r25091 = r25087 * r25090; double r25092 = r25087 * r25091; double r25093 = r25082 / r25092; return r25093; } 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 r25094, r25095, r25096, r25097, r25098, r25099, r25100, r25101, r25102, r25103, r25104, r25105, r25106, r25107, r25108, r25109, r25110, r25111, r25112; void setup_mpfr_f_im() { mpfr_set_default_prec(4432); mpfr_init_set_str(r25094, "2", 10, MPFR_RNDN); mpfr_init(r25095); mpfr_init_set_str(r25096, "3", 10, MPFR_RNDN); mpfr_init(r25097); mpfr_init(r25098); mpfr_init(r25099); mpfr_init(r25100); mpfr_init(r25101); mpfr_init(r25102); mpfr_init(r25103); mpfr_init(r25104); mpfr_init(r25105); mpfr_init_set_str(r25106, "1", 10, MPFR_RNDN); mpfr_init(r25107); mpfr_init(r25108); mpfr_init(r25109); mpfr_init(r25110); mpfr_init(r25111); mpfr_init(r25112); } double f_im(double t, double l, double k) { ; mpfr_set_d(r25095, t, MPFR_RNDN); ; mpfr_pow(r25097, r25095, r25096, MPFR_RNDN); mpfr_set_d(r25098, l, MPFR_RNDN); mpfr_mul(r25099, r25098, r25098, MPFR_RNDN); mpfr_div(r25100, r25097, r25099, MPFR_RNDN); mpfr_set_d(r25101, k, MPFR_RNDN); mpfr_sin(r25102, r25101, MPFR_RNDN); mpfr_mul(r25103, r25100, r25102, MPFR_RNDN); mpfr_tan(r25104, r25101, MPFR_RNDN); mpfr_mul(r25105, r25103, r25104, MPFR_RNDN); ; mpfr_div(r25107, r25101, r25095, MPFR_RNDN); mpfr_pow(r25108, r25107, r25094, MPFR_RNDN); mpfr_add(r25109, r25106, r25108, MPFR_RNDN); mpfr_sub(r25110, r25109, r25106, MPFR_RNDN); mpfr_mul(r25111, r25105, r25110, MPFR_RNDN); mpfr_div(r25112, r25094, r25111, MPFR_RNDN); return mpfr_get_d(r25112, MPFR_RNDN); } static mpfr_t r25113, r25114, r25115, r25116, r25117, r25118, r25119, r25120, r25121, r25122, r25123, r25124, r25125, r25126; void setup_mpfr_f_fm() { mpfr_set_default_prec(4432); mpfr_init_set_str(r25113, "2", 10, MPFR_RNDN); mpfr_init_set_str(r25114, "1", 10, MPFR_RNDN); mpfr_init(r25115); mpfr_init(r25116); mpfr_init(r25117); mpfr_init(r25118); mpfr_init(r25119); mpfr_init(r25120); mpfr_init(r25121); mpfr_init(r25122); mpfr_init(r25123); mpfr_init(r25124); mpfr_init(r25125); mpfr_init(r25126); } double f_fm(double t, double l, double k) { ; ; mpfr_div(r25115, r25113, r25114, MPFR_RNDN); mpfr_set_d(r25116, k, MPFR_RNDN); mpfr_set_d(r25117, l, MPFR_RNDN); mpfr_div(r25118, r25116, r25117, MPFR_RNDN); mpfr_sin(r25119, r25116, MPFR_RNDN); mpfr_mul(r25120, r25118, r25119, MPFR_RNDN); mpfr_set_d(r25121, t, MPFR_RNDN); mpfr_cos(r25122, r25116, MPFR_RNDN); mpfr_div(r25123, r25121, r25122, MPFR_RNDN); mpfr_mul(r25124, r25120, r25123, MPFR_RNDN); mpfr_mul(r25125, r25120, r25124, MPFR_RNDN); mpfr_div(r25126, r25115, r25125, MPFR_RNDN); return mpfr_get_d(r25126, MPFR_RNDN); } static mpfr_t r25127, r25128, r25129, r25130, r25131, r25132, r25133, r25134, r25135, r25136, r25137, r25138, r25139, r25140; void setup_mpfr_f_dm() { mpfr_set_default_prec(4432); mpfr_init_set_str(r25127, "2", 10, MPFR_RNDN); mpfr_init_set_str(r25128, "1", 10, MPFR_RNDN); mpfr_init(r25129); mpfr_init(r25130); mpfr_init(r25131); mpfr_init(r25132); mpfr_init(r25133); mpfr_init(r25134); mpfr_init(r25135); mpfr_init(r25136); mpfr_init(r25137); mpfr_init(r25138); mpfr_init(r25139); mpfr_init(r25140); } double f_dm(double t, double l, double k) { ; ; mpfr_div(r25129, r25127, r25128, MPFR_RNDN); mpfr_set_d(r25130, k, MPFR_RNDN); mpfr_set_d(r25131, l, MPFR_RNDN); mpfr_div(r25132, r25130, r25131, MPFR_RNDN); mpfr_sin(r25133, r25130, MPFR_RNDN); mpfr_mul(r25134, r25132, r25133, MPFR_RNDN); mpfr_set_d(r25135, t, MPFR_RNDN); mpfr_cos(r25136, r25130, MPFR_RNDN); mpfr_div(r25137, r25135, r25136, MPFR_RNDN); mpfr_mul(r25138, r25134, r25137, MPFR_RNDN); mpfr_mul(r25139, r25134, r25138, MPFR_RNDN); mpfr_div(r25140, r25129, r25139, MPFR_RNDN); return mpfr_get_d(r25140, MPFR_RNDN); }