#include #include #include #include #include char *name = "Asymptote C"; double f_if(float x) { float r23113 = x; float r23114 = 1; float r23115 = r23113 + r23114; float r23116 = r23113 / r23115; float r23117 = r23113 - r23114; float r23118 = r23115 / r23117; float r23119 = r23116 - r23118; return r23119; } double f_id(double x) { double r23120 = x; double r23121 = 1; double r23122 = r23120 + r23121; double r23123 = r23120 / r23122; double r23124 = r23120 - r23121; double r23125 = r23122 / r23124; double r23126 = r23123 - r23125; return r23126; } double f_of(float x) { float r23127 = x; float r23128 = 1; float r23129 = r23127 + r23128; float r23130 = r23127 / r23129; float r23131 = r23127 - r23128; float r23132 = r23129 / r23131; float r23133 = r23130 - r23132; float r23134 = 1.1213252548714081e-14; bool r23135 = r23133 <= r23134; float r23136 = 3; float r23137 = r23136 / r23127; float r23138 = r23128 + r23137; float r23139 = -r23128; float r23140 = r23127 * r23127; float r23141 = r23139 / r23140; float r23142 = -r23136; float r23143 = r23142 / r23127; float r23144 = fma(r23138, r23141, r23143); float r23145 = cbrt(r23129); float r23146 = r23145 * r23145; float r23147 = r23127 / r23146; float r23148 = r23147 / r23145; float r23149 = r23148 - r23132; float r23150 = r23135 ? r23144 : r23149; return r23150; } double f_od(double x) { double r23151 = x; double r23152 = 1; double r23153 = r23151 + r23152; double r23154 = r23151 / r23153; double r23155 = r23151 - r23152; double r23156 = r23153 / r23155; double r23157 = r23154 - r23156; double r23158 = 1.1213252548714081e-14; bool r23159 = r23157 <= r23158; double r23160 = 3; double r23161 = r23160 / r23151; double r23162 = r23152 + r23161; double r23163 = -r23152; double r23164 = r23151 * r23151; double r23165 = r23163 / r23164; double r23166 = -r23160; double r23167 = r23166 / r23151; double r23168 = fma(r23162, r23165, r23167); double r23169 = cbrt(r23153); double r23170 = r23169 * r23169; double r23171 = r23151 / r23170; double r23172 = r23171 / r23169; double r23173 = r23172 - r23156; double r23174 = r23159 ? r23168 : r23173; return r23174; } 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 r23175, r23176, r23177, r23178, r23179, r23180, r23181; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r23175); mpfr_init_set_str(r23176, "1", 10, MPFR_RNDN); mpfr_init(r23177); mpfr_init(r23178); mpfr_init(r23179); mpfr_init(r23180); mpfr_init(r23181); } double f_im(double x) { mpfr_set_d(r23175, x, MPFR_RNDN); ; mpfr_add(r23177, r23175, r23176, MPFR_RNDN); mpfr_div(r23178, r23175, r23177, MPFR_RNDN); mpfr_sub(r23179, r23175, r23176, MPFR_RNDN); mpfr_div(r23180, r23177, r23179, MPFR_RNDN); mpfr_sub(r23181, r23178, r23180, MPFR_RNDN); return mpfr_get_d(r23181, MPFR_RNDN); } static mpfr_t r23182, r23183, r23184, r23185, r23186, r23187, r23188, r23189, r23190, r23191, r23192, r23193, r23194, r23195, r23196, r23197, r23198, r23199, r23200, r23201, r23202, r23203, r23204, r23205; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r23182); mpfr_init_set_str(r23183, "1", 10, MPFR_RNDN); mpfr_init(r23184); mpfr_init(r23185); mpfr_init(r23186); mpfr_init(r23187); mpfr_init(r23188); mpfr_init_set_str(r23189, "1.1213252548714081e-14", 10, MPFR_RNDN); mpfr_init(r23190); mpfr_init_set_str(r23191, "3", 10, MPFR_RNDN); mpfr_init(r23192); mpfr_init(r23193); mpfr_init(r23194); mpfr_init(r23195); mpfr_init(r23196); mpfr_init(r23197); mpfr_init(r23198); mpfr_init(r23199); mpfr_init(r23200); mpfr_init(r23201); mpfr_init(r23202); mpfr_init(r23203); mpfr_init(r23204); mpfr_init(r23205); } double f_fm(double x) { mpfr_set_d(r23182, x, MPFR_RNDN); ; mpfr_add(r23184, r23182, r23183, MPFR_RNDN); mpfr_div(r23185, r23182, r23184, MPFR_RNDN); mpfr_sub(r23186, r23182, r23183, MPFR_RNDN); mpfr_div(r23187, r23184, r23186, MPFR_RNDN); mpfr_sub(r23188, r23185, r23187, MPFR_RNDN); ; mpfr_set_si(r23190, mpfr_cmp(r23188, r23189) <= 0, MPFR_RNDN); ; mpfr_div(r23192, r23191, r23182, MPFR_RNDN); mpfr_add(r23193, r23183, r23192, MPFR_RNDN); mpfr_neg(r23194, r23183, MPFR_RNDN); mpfr_mul(r23195, r23182, r23182, MPFR_RNDN); mpfr_div(r23196, r23194, r23195, MPFR_RNDN); mpfr_neg(r23197, r23191, MPFR_RNDN); mpfr_div(r23198, r23197, r23182, MPFR_RNDN); mpfr_fma(r23199, r23193, r23196, r23198, MPFR_RNDN); mpfr_cbrt(r23200, r23184, MPFR_RNDN); mpfr_mul(r23201, r23200, r23200, MPFR_RNDN); mpfr_div(r23202, r23182, r23201, MPFR_RNDN); mpfr_div(r23203, r23202, r23200, MPFR_RNDN); mpfr_sub(r23204, r23203, r23187, MPFR_RNDN); if (mpfr_get_si(r23190, MPFR_RNDN)) { mpfr_set(r23205, r23199, MPFR_RNDN); } else { mpfr_set(r23205, r23204, MPFR_RNDN); }; return mpfr_get_d(r23205, MPFR_RNDN); } static mpfr_t r23206, r23207, r23208, r23209, r23210, r23211, r23212, r23213, r23214, r23215, r23216, r23217, r23218, r23219, r23220, r23221, r23222, r23223, r23224, r23225, r23226, r23227, r23228, r23229; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r23206); mpfr_init_set_str(r23207, "1", 10, MPFR_RNDN); mpfr_init(r23208); mpfr_init(r23209); mpfr_init(r23210); mpfr_init(r23211); mpfr_init(r23212); mpfr_init_set_str(r23213, "1.1213252548714081e-14", 10, MPFR_RNDN); mpfr_init(r23214); mpfr_init_set_str(r23215, "3", 10, MPFR_RNDN); mpfr_init(r23216); mpfr_init(r23217); mpfr_init(r23218); mpfr_init(r23219); mpfr_init(r23220); mpfr_init(r23221); mpfr_init(r23222); mpfr_init(r23223); mpfr_init(r23224); mpfr_init(r23225); mpfr_init(r23226); mpfr_init(r23227); mpfr_init(r23228); mpfr_init(r23229); } double f_dm(double x) { mpfr_set_d(r23206, x, MPFR_RNDN); ; mpfr_add(r23208, r23206, r23207, MPFR_RNDN); mpfr_div(r23209, r23206, r23208, MPFR_RNDN); mpfr_sub(r23210, r23206, r23207, MPFR_RNDN); mpfr_div(r23211, r23208, r23210, MPFR_RNDN); mpfr_sub(r23212, r23209, r23211, MPFR_RNDN); ; mpfr_set_si(r23214, mpfr_cmp(r23212, r23213) <= 0, MPFR_RNDN); ; mpfr_div(r23216, r23215, r23206, MPFR_RNDN); mpfr_add(r23217, r23207, r23216, MPFR_RNDN); mpfr_neg(r23218, r23207, MPFR_RNDN); mpfr_mul(r23219, r23206, r23206, MPFR_RNDN); mpfr_div(r23220, r23218, r23219, MPFR_RNDN); mpfr_neg(r23221, r23215, MPFR_RNDN); mpfr_div(r23222, r23221, r23206, MPFR_RNDN); mpfr_fma(r23223, r23217, r23220, r23222, MPFR_RNDN); mpfr_cbrt(r23224, r23208, MPFR_RNDN); mpfr_mul(r23225, r23224, r23224, MPFR_RNDN); mpfr_div(r23226, r23206, r23225, MPFR_RNDN); mpfr_div(r23227, r23226, r23224, MPFR_RNDN); mpfr_sub(r23228, r23227, r23211, MPFR_RNDN); if (mpfr_get_si(r23214, MPFR_RNDN)) { mpfr_set(r23229, r23223, MPFR_RNDN); } else { mpfr_set(r23229, r23228, MPFR_RNDN); }; return mpfr_get_d(r23229, MPFR_RNDN); }