#include #include #include #include #include char *name = "r*sin(b)/cos(a+b), A"; double f_if(float r, float a, float b) { float r24021 = r; float r24022 = b; float r24023 = sin(r24022); float r24024 = r24021 * r24023; float r24025 = a; float r24026 = r24025 + r24022; float r24027 = cos(r24026); float r24028 = r24024 / r24027; return r24028; } double f_id(double r, double a, double b) { double r24029 = r; double r24030 = b; double r24031 = sin(r24030); double r24032 = r24029 * r24031; double r24033 = a; double r24034 = r24033 + r24030; double r24035 = cos(r24034); double r24036 = r24032 / r24035; return r24036; } double f_of(float r, float a, float b) { float r24037 = b; float r24038 = sin(r24037); float r24039 = a; float r24040 = sin(r24039); float r24041 = r24038 * r24040; float r24042 = cos(r24037); float r24043 = cos(r24039); float r24044 = r24042 * r24043; float r24045 = r24041 + r24044; float r24046 = r24038 / r24045; float r24047 = r; float r24048 = r24044 - r24041; float r24049 = r24047 / r24048; float r24050 = r24046 * r24049; float r24051 = r24043 * r24042; float r24052 = r24040 * r24038; float r24053 = exp(r24052); float r24054 = log(r24053); float r24055 = r24051 + r24054; float r24056 = r24050 * r24055; return r24056; } double f_od(double r, double a, double b) { double r24057 = b; double r24058 = sin(r24057); double r24059 = a; double r24060 = sin(r24059); double r24061 = r24058 * r24060; double r24062 = cos(r24057); double r24063 = cos(r24059); double r24064 = r24062 * r24063; double r24065 = r24061 + r24064; double r24066 = r24058 / r24065; double r24067 = r; double r24068 = r24064 - r24061; double r24069 = r24067 / r24068; double r24070 = r24066 * r24069; double r24071 = r24063 * r24062; double r24072 = r24060 * r24058; double r24073 = exp(r24072); double r24074 = log(r24073); double r24075 = r24071 + r24074; double r24076 = r24070 * r24075; return r24076; } 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 r24077, r24078, r24079, r24080, r24081, r24082, r24083, r24084; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r24077); mpfr_init(r24078); mpfr_init(r24079); mpfr_init(r24080); mpfr_init(r24081); mpfr_init(r24082); mpfr_init(r24083); mpfr_init(r24084); } double f_im(double r, double a, double b) { mpfr_set_d(r24077, r, MPFR_RNDN); mpfr_set_d(r24078, b, MPFR_RNDN); mpfr_sin(r24079, r24078, MPFR_RNDN); mpfr_mul(r24080, r24077, r24079, MPFR_RNDN); mpfr_set_d(r24081, a, MPFR_RNDN); mpfr_add(r24082, r24081, r24078, MPFR_RNDN); mpfr_cos(r24083, r24082, MPFR_RNDN); mpfr_div(r24084, r24080, r24083, MPFR_RNDN); return mpfr_get_d(r24084, MPFR_RNDN); } static mpfr_t r24085, r24086, r24087, r24088, r24089, r24090, r24091, r24092, r24093, r24094, r24095, r24096, r24097, r24098, r24099, r24100, r24101, r24102, r24103, r24104; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r24085); mpfr_init(r24086); mpfr_init(r24087); mpfr_init(r24088); mpfr_init(r24089); mpfr_init(r24090); mpfr_init(r24091); mpfr_init(r24092); mpfr_init(r24093); mpfr_init(r24094); mpfr_init(r24095); mpfr_init(r24096); mpfr_init(r24097); mpfr_init(r24098); mpfr_init(r24099); mpfr_init(r24100); mpfr_init(r24101); mpfr_init(r24102); mpfr_init(r24103); mpfr_init(r24104); } double f_fm(double r, double a, double b) { mpfr_set_d(r24085, b, MPFR_RNDN); mpfr_sin(r24086, r24085, MPFR_RNDN); mpfr_set_d(r24087, a, MPFR_RNDN); mpfr_sin(r24088, r24087, MPFR_RNDN); mpfr_mul(r24089, r24086, r24088, MPFR_RNDN); mpfr_cos(r24090, r24085, MPFR_RNDN); mpfr_cos(r24091, r24087, MPFR_RNDN); mpfr_mul(r24092, r24090, r24091, MPFR_RNDN); mpfr_add(r24093, r24089, r24092, MPFR_RNDN); mpfr_div(r24094, r24086, r24093, MPFR_RNDN); mpfr_set_d(r24095, r, MPFR_RNDN); mpfr_sub(r24096, r24092, r24089, MPFR_RNDN); mpfr_div(r24097, r24095, r24096, MPFR_RNDN); mpfr_mul(r24098, r24094, r24097, MPFR_RNDN); mpfr_mul(r24099, r24091, r24090, MPFR_RNDN); mpfr_mul(r24100, r24088, r24086, MPFR_RNDN); mpfr_exp(r24101, r24100, MPFR_RNDN); mpfr_log(r24102, r24101, MPFR_RNDN); mpfr_add(r24103, r24099, r24102, MPFR_RNDN); mpfr_mul(r24104, r24098, r24103, MPFR_RNDN); return mpfr_get_d(r24104, MPFR_RNDN); } static mpfr_t r24105, r24106, r24107, r24108, r24109, r24110, r24111, r24112, r24113, r24114, r24115, r24116, r24117, r24118, r24119, r24120, r24121, r24122, r24123, r24124; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r24105); mpfr_init(r24106); mpfr_init(r24107); mpfr_init(r24108); mpfr_init(r24109); mpfr_init(r24110); mpfr_init(r24111); mpfr_init(r24112); mpfr_init(r24113); mpfr_init(r24114); mpfr_init(r24115); mpfr_init(r24116); mpfr_init(r24117); mpfr_init(r24118); mpfr_init(r24119); mpfr_init(r24120); mpfr_init(r24121); mpfr_init(r24122); mpfr_init(r24123); mpfr_init(r24124); } double f_dm(double r, double a, double b) { mpfr_set_d(r24105, b, MPFR_RNDN); mpfr_sin(r24106, r24105, MPFR_RNDN); mpfr_set_d(r24107, a, MPFR_RNDN); mpfr_sin(r24108, r24107, MPFR_RNDN); mpfr_mul(r24109, r24106, r24108, MPFR_RNDN); mpfr_cos(r24110, r24105, MPFR_RNDN); mpfr_cos(r24111, r24107, MPFR_RNDN); mpfr_mul(r24112, r24110, r24111, MPFR_RNDN); mpfr_add(r24113, r24109, r24112, MPFR_RNDN); mpfr_div(r24114, r24106, r24113, MPFR_RNDN); mpfr_set_d(r24115, r, MPFR_RNDN); mpfr_sub(r24116, r24112, r24109, MPFR_RNDN); mpfr_div(r24117, r24115, r24116, MPFR_RNDN); mpfr_mul(r24118, r24114, r24117, MPFR_RNDN); mpfr_mul(r24119, r24111, r24110, MPFR_RNDN); mpfr_mul(r24120, r24108, r24106, MPFR_RNDN); mpfr_exp(r24121, r24120, MPFR_RNDN); mpfr_log(r24122, r24121, MPFR_RNDN); mpfr_add(r24123, r24119, r24122, MPFR_RNDN); mpfr_mul(r24124, r24118, r24123, MPFR_RNDN); return mpfr_get_d(r24124, MPFR_RNDN); }