#include #include #include #include #include char *name = "Numeric.SpecFunctions.Extra:bd0 from math-functions-0.1.5.2"; double f_if(float x, float y, float z) { float r45077 = x; float r45078 = y; float r45079 = r45077 / r45078; float r45080 = log(r45079); float r45081 = r45077 * r45080; float r45082 = z; float r45083 = r45081 - r45082; return r45083; } double f_id(double x, double y, double z) { double r45084 = x; double r45085 = y; double r45086 = r45084 / r45085; double r45087 = log(r45086); double r45088 = r45084 * r45087; double r45089 = z; double r45090 = r45088 - r45089; return r45090; } double f_of(float x, float y, float z) { float r45091 = x; float r45092 = -4.16329420604066e-310; bool r45093 = r45091 <= r45092; float r45094 = y; float r45095 = r45091 / r45094; float r45096 = cbrt(r45095); float r45097 = log(r45096); float r45098 = r45091 + r45091; float r45099 = r45097 * r45098; float r45100 = r45091 * r45097; float r45101 = r45099 + r45100; float r45102 = z; float r45103 = r45101 - r45102; float r45104 = log(r45091); float r45105 = log(r45094); float r45106 = r45104 - r45105; float r45107 = r45091 * r45106; float r45108 = r45107 - r45102; float r45109 = r45093 ? r45103 : r45108; return r45109; } double f_od(double x, double y, double z) { double r45110 = x; double r45111 = -4.16329420604066e-310; bool r45112 = r45110 <= r45111; double r45113 = y; double r45114 = r45110 / r45113; double r45115 = cbrt(r45114); double r45116 = log(r45115); double r45117 = r45110 + r45110; double r45118 = r45116 * r45117; double r45119 = r45110 * r45116; double r45120 = r45118 + r45119; double r45121 = z; double r45122 = r45120 - r45121; double r45123 = log(r45110); double r45124 = log(r45113); double r45125 = r45123 - r45124; double r45126 = r45110 * r45125; double r45127 = r45126 - r45121; double r45128 = r45112 ? r45122 : r45127; return r45128; } 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 r45129, r45130, r45131, r45132, r45133, r45134, r45135; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r45129); mpfr_init(r45130); mpfr_init(r45131); mpfr_init(r45132); mpfr_init(r45133); mpfr_init(r45134); mpfr_init(r45135); } double f_im(double x, double y, double z) { mpfr_set_d(r45129, x, MPFR_RNDN); mpfr_set_d(r45130, y, MPFR_RNDN); mpfr_div(r45131, r45129, r45130, MPFR_RNDN); mpfr_log(r45132, r45131, MPFR_RNDN); mpfr_mul(r45133, r45129, r45132, MPFR_RNDN); mpfr_set_d(r45134, z, MPFR_RNDN); mpfr_sub(r45135, r45133, r45134, MPFR_RNDN); return mpfr_get_d(r45135, MPFR_RNDN); } static mpfr_t r45136, r45137, r45138, r45139, r45140, r45141, r45142, r45143, r45144, r45145, r45146, r45147, r45148, r45149, r45150, r45151, r45152, r45153, r45154; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r45136); mpfr_init_set_str(r45137, "-4.16329420604066e-310", 10, MPFR_RNDN); mpfr_init(r45138); mpfr_init(r45139); mpfr_init(r45140); mpfr_init(r45141); mpfr_init(r45142); mpfr_init(r45143); mpfr_init(r45144); mpfr_init(r45145); mpfr_init(r45146); mpfr_init(r45147); mpfr_init(r45148); mpfr_init(r45149); mpfr_init(r45150); mpfr_init(r45151); mpfr_init(r45152); mpfr_init(r45153); mpfr_init(r45154); } double f_fm(double x, double y, double z) { mpfr_set_d(r45136, x, MPFR_RNDN); ; mpfr_set_si(r45138, mpfr_cmp(r45136, r45137) <= 0, MPFR_RNDN); mpfr_set_d(r45139, y, MPFR_RNDN); mpfr_div(r45140, r45136, r45139, MPFR_RNDN); mpfr_cbrt(r45141, r45140, MPFR_RNDN); mpfr_log(r45142, r45141, MPFR_RNDN); mpfr_add(r45143, r45136, r45136, MPFR_RNDN); mpfr_mul(r45144, r45142, r45143, MPFR_RNDN); mpfr_mul(r45145, r45136, r45142, MPFR_RNDN); mpfr_add(r45146, r45144, r45145, MPFR_RNDN); mpfr_set_d(r45147, z, MPFR_RNDN); mpfr_sub(r45148, r45146, r45147, MPFR_RNDN); mpfr_log(r45149, r45136, MPFR_RNDN); mpfr_log(r45150, r45139, MPFR_RNDN); mpfr_sub(r45151, r45149, r45150, MPFR_RNDN); mpfr_mul(r45152, r45136, r45151, MPFR_RNDN); mpfr_sub(r45153, r45152, r45147, MPFR_RNDN); if (mpfr_get_si(r45138, MPFR_RNDN)) { mpfr_set(r45154, r45148, MPFR_RNDN); } else { mpfr_set(r45154, r45153, MPFR_RNDN); }; return mpfr_get_d(r45154, MPFR_RNDN); } static mpfr_t r45155, r45156, r45157, r45158, r45159, r45160, r45161, r45162, r45163, r45164, r45165, r45166, r45167, r45168, r45169, r45170, r45171, r45172, r45173; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r45155); mpfr_init_set_str(r45156, "-4.16329420604066e-310", 10, MPFR_RNDN); mpfr_init(r45157); mpfr_init(r45158); mpfr_init(r45159); mpfr_init(r45160); mpfr_init(r45161); mpfr_init(r45162); mpfr_init(r45163); mpfr_init(r45164); mpfr_init(r45165); mpfr_init(r45166); mpfr_init(r45167); mpfr_init(r45168); mpfr_init(r45169); mpfr_init(r45170); mpfr_init(r45171); mpfr_init(r45172); mpfr_init(r45173); } double f_dm(double x, double y, double z) { mpfr_set_d(r45155, x, MPFR_RNDN); ; mpfr_set_si(r45157, mpfr_cmp(r45155, r45156) <= 0, MPFR_RNDN); mpfr_set_d(r45158, y, MPFR_RNDN); mpfr_div(r45159, r45155, r45158, MPFR_RNDN); mpfr_cbrt(r45160, r45159, MPFR_RNDN); mpfr_log(r45161, r45160, MPFR_RNDN); mpfr_add(r45162, r45155, r45155, MPFR_RNDN); mpfr_mul(r45163, r45161, r45162, MPFR_RNDN); mpfr_mul(r45164, r45155, r45161, MPFR_RNDN); mpfr_add(r45165, r45163, r45164, MPFR_RNDN); mpfr_set_d(r45166, z, MPFR_RNDN); mpfr_sub(r45167, r45165, r45166, MPFR_RNDN); mpfr_log(r45168, r45155, MPFR_RNDN); mpfr_log(r45169, r45158, MPFR_RNDN); mpfr_sub(r45170, r45168, r45169, MPFR_RNDN); mpfr_mul(r45171, r45155, r45170, MPFR_RNDN); mpfr_sub(r45172, r45171, r45166, MPFR_RNDN); if (mpfr_get_si(r45157, MPFR_RNDN)) { mpfr_set(r45173, r45167, MPFR_RNDN); } else { mpfr_set(r45173, r45172, MPFR_RNDN); }; return mpfr_get_d(r45173, MPFR_RNDN); }