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