#include #include #include #include #include char *name = "Asymptote C"; double f_if(float x) { float r23958 = x; float r23959 = 1; float r23960 = r23958 + r23959; float r23961 = r23958 / r23960; float r23962 = r23958 - r23959; float r23963 = r23960 / r23962; float r23964 = r23961 - r23963; return r23964; } double f_id(double x) { double r23965 = x; double r23966 = 1; double r23967 = r23965 + r23966; double r23968 = r23965 / r23967; double r23969 = r23965 - r23966; double r23970 = r23967 / r23969; double r23971 = r23968 - r23970; return r23971; } double f_of(float x) { float r23972 = x; float r23973 = 1; float r23974 = r23972 + r23973; float r23975 = r23972 / r23974; float r23976 = r23972 - r23973; float r23977 = r23974 / r23976; float r23978 = r23975 - r23977; float r23979 = 6.100063898450969e-09; bool r23980 = r23978 <= r23979; float r23981 = 3; float r23982 = r23981 / r23972; float r23983 = r23973 + r23982; float r23984 = -r23973; float r23985 = r23972 * r23972; float r23986 = r23984 / r23985; float r23987 = -r23982; float r23988 = fma(r23983, r23986, r23987); float r23989 = r23973 + r23972; float r23990 = r23972 / r23989; float r23991 = r23989 / r23976; float r23992 = r23990 - r23991; float r23993 = pow(r23992, r23981); float r23994 = cbrt(r23993); float r23995 = r23980 ? r23988 : r23994; return r23995; } double f_od(double x) { double r23996 = x; double r23997 = 1; double r23998 = r23996 + r23997; double r23999 = r23996 / r23998; double r24000 = r23996 - r23997; double r24001 = r23998 / r24000; double r24002 = r23999 - r24001; double r24003 = 6.100063898450969e-09; bool r24004 = r24002 <= r24003; double r24005 = 3; double r24006 = r24005 / r23996; double r24007 = r23997 + r24006; double r24008 = -r23997; double r24009 = r23996 * r23996; double r24010 = r24008 / r24009; double r24011 = -r24006; double r24012 = fma(r24007, r24010, r24011); double r24013 = r23997 + r23996; double r24014 = r23996 / r24013; double r24015 = r24013 / r24000; double r24016 = r24014 - r24015; double r24017 = pow(r24016, r24005); double r24018 = cbrt(r24017); double r24019 = r24004 ? r24012 : r24018; return r24019; } 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 r24020, r24021, r24022, r24023, r24024, r24025, r24026; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r24020); mpfr_init_set_str(r24021, "1", 10, MPFR_RNDN); mpfr_init(r24022); mpfr_init(r24023); mpfr_init(r24024); mpfr_init(r24025); mpfr_init(r24026); } double f_im(double x) { mpfr_set_d(r24020, x, MPFR_RNDN); ; mpfr_add(r24022, r24020, r24021, MPFR_RNDN); mpfr_div(r24023, r24020, r24022, MPFR_RNDN); mpfr_sub(r24024, r24020, r24021, MPFR_RNDN); mpfr_div(r24025, r24022, r24024, MPFR_RNDN); mpfr_sub(r24026, r24023, r24025, MPFR_RNDN); return mpfr_get_d(r24026, MPFR_RNDN); } static mpfr_t r24027, r24028, r24029, r24030, r24031, r24032, r24033, r24034, r24035, r24036, r24037, r24038, r24039, r24040, r24041, r24042, r24043, r24044, r24045, r24046, r24047, r24048, r24049, r24050; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r24027); mpfr_init_set_str(r24028, "1", 10, MPFR_RNDN); mpfr_init(r24029); mpfr_init(r24030); mpfr_init(r24031); mpfr_init(r24032); mpfr_init(r24033); mpfr_init_set_str(r24034, "6.100063898450969e-09", 10, MPFR_RNDN); mpfr_init(r24035); mpfr_init_set_str(r24036, "3", 10, MPFR_RNDN); mpfr_init(r24037); mpfr_init(r24038); mpfr_init(r24039); mpfr_init(r24040); mpfr_init(r24041); mpfr_init(r24042); mpfr_init(r24043); mpfr_init(r24044); mpfr_init(r24045); mpfr_init(r24046); mpfr_init(r24047); mpfr_init(r24048); mpfr_init(r24049); mpfr_init(r24050); } double f_fm(double x) { mpfr_set_d(r24027, x, MPFR_RNDN); ; mpfr_add(r24029, r24027, r24028, MPFR_RNDN); mpfr_div(r24030, r24027, r24029, MPFR_RNDN); mpfr_sub(r24031, r24027, r24028, MPFR_RNDN); mpfr_div(r24032, r24029, r24031, MPFR_RNDN); mpfr_sub(r24033, r24030, r24032, MPFR_RNDN); ; mpfr_set_si(r24035, mpfr_cmp(r24033, r24034) <= 0, MPFR_RNDN); ; mpfr_div(r24037, r24036, r24027, MPFR_RNDN); mpfr_add(r24038, r24028, r24037, MPFR_RNDN); mpfr_neg(r24039, r24028, MPFR_RNDN); mpfr_mul(r24040, r24027, r24027, MPFR_RNDN); mpfr_div(r24041, r24039, r24040, MPFR_RNDN); mpfr_neg(r24042, r24037, MPFR_RNDN); mpfr_fma(r24043, r24038, r24041, r24042, MPFR_RNDN); mpfr_add(r24044, r24028, r24027, MPFR_RNDN); mpfr_div(r24045, r24027, r24044, MPFR_RNDN); mpfr_div(r24046, r24044, r24031, MPFR_RNDN); mpfr_sub(r24047, r24045, r24046, MPFR_RNDN); mpfr_pow(r24048, r24047, r24036, MPFR_RNDN); mpfr_cbrt(r24049, r24048, MPFR_RNDN); if (mpfr_get_si(r24035, MPFR_RNDN)) { mpfr_set(r24050, r24043, MPFR_RNDN); } else { mpfr_set(r24050, r24049, MPFR_RNDN); }; return mpfr_get_d(r24050, MPFR_RNDN); } static mpfr_t r24051, r24052, r24053, r24054, r24055, r24056, r24057, r24058, r24059, r24060, r24061, r24062, r24063, r24064, r24065, r24066, r24067, r24068, r24069, r24070, r24071, r24072, r24073, r24074; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r24051); mpfr_init_set_str(r24052, "1", 10, MPFR_RNDN); mpfr_init(r24053); mpfr_init(r24054); mpfr_init(r24055); mpfr_init(r24056); mpfr_init(r24057); mpfr_init_set_str(r24058, "6.100063898450969e-09", 10, MPFR_RNDN); mpfr_init(r24059); mpfr_init_set_str(r24060, "3", 10, MPFR_RNDN); mpfr_init(r24061); mpfr_init(r24062); mpfr_init(r24063); mpfr_init(r24064); mpfr_init(r24065); mpfr_init(r24066); mpfr_init(r24067); mpfr_init(r24068); mpfr_init(r24069); mpfr_init(r24070); mpfr_init(r24071); mpfr_init(r24072); mpfr_init(r24073); mpfr_init(r24074); } double f_dm(double x) { mpfr_set_d(r24051, x, MPFR_RNDN); ; mpfr_add(r24053, r24051, r24052, MPFR_RNDN); mpfr_div(r24054, r24051, r24053, MPFR_RNDN); mpfr_sub(r24055, r24051, r24052, MPFR_RNDN); mpfr_div(r24056, r24053, r24055, MPFR_RNDN); mpfr_sub(r24057, r24054, r24056, MPFR_RNDN); ; mpfr_set_si(r24059, mpfr_cmp(r24057, r24058) <= 0, MPFR_RNDN); ; mpfr_div(r24061, r24060, r24051, MPFR_RNDN); mpfr_add(r24062, r24052, r24061, MPFR_RNDN); mpfr_neg(r24063, r24052, MPFR_RNDN); mpfr_mul(r24064, r24051, r24051, MPFR_RNDN); mpfr_div(r24065, r24063, r24064, MPFR_RNDN); mpfr_neg(r24066, r24061, MPFR_RNDN); mpfr_fma(r24067, r24062, r24065, r24066, MPFR_RNDN); mpfr_add(r24068, r24052, r24051, MPFR_RNDN); mpfr_div(r24069, r24051, r24068, MPFR_RNDN); mpfr_div(r24070, r24068, r24055, MPFR_RNDN); mpfr_sub(r24071, r24069, r24070, MPFR_RNDN); mpfr_pow(r24072, r24071, r24060, MPFR_RNDN); mpfr_cbrt(r24073, r24072, MPFR_RNDN); if (mpfr_get_si(r24059, MPFR_RNDN)) { mpfr_set(r24074, r24067, MPFR_RNDN); } else { mpfr_set(r24074, r24073, MPFR_RNDN); }; return mpfr_get_d(r24074, MPFR_RNDN); }