#include #include #include #include #include char *name = "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2"; double f_if(float x, float y, float z, float t, float a, float b, float c) { float r23073 = x; float r23074 = y; float r23075 = 2.0; float r23076 = z; float r23077 = t; float r23078 = a; float r23079 = r23077 + r23078; float r23080 = sqrt(r23079); float r23081 = r23076 * r23080; float r23082 = r23081 / r23077; float r23083 = b; float r23084 = c; float r23085 = r23083 - r23084; float r23086 = 5.0; float r23087 = 6.0; float r23088 = r23086 / r23087; float r23089 = r23078 + r23088; float r23090 = 3.0; float r23091 = r23077 * r23090; float r23092 = r23075 / r23091; float r23093 = r23089 - r23092; float r23094 = r23085 * r23093; float r23095 = r23082 - r23094; float r23096 = r23075 * r23095; float r23097 = exp(r23096); float r23098 = r23074 * r23097; float r23099 = r23073 + r23098; float r23100 = r23073 / r23099; return r23100; } double f_id(double x, double y, double z, double t, double a, double b, double c) { double r23101 = x; double r23102 = y; double r23103 = 2.0; double r23104 = z; double r23105 = t; double r23106 = a; double r23107 = r23105 + r23106; double r23108 = sqrt(r23107); double r23109 = r23104 * r23108; double r23110 = r23109 / r23105; double r23111 = b; double r23112 = c; double r23113 = r23111 - r23112; double r23114 = 5.0; double r23115 = 6.0; double r23116 = r23114 / r23115; double r23117 = r23106 + r23116; double r23118 = 3.0; double r23119 = r23105 * r23118; double r23120 = r23103 / r23119; double r23121 = r23117 - r23120; double r23122 = r23113 * r23121; double r23123 = r23110 - r23122; double r23124 = r23103 * r23123; double r23125 = exp(r23124); double r23126 = r23102 * r23125; double r23127 = r23101 + r23126; double r23128 = r23101 / r23127; return r23128; } double f_of(float x, float y, float z, float t, float a, float b, float c) { float r23129 = x; float r23130 = y; float r23131 = 2.0; float r23132 = a; float r23133 = 5.0; float r23134 = 6.0; float r23135 = r23133 / r23134; float r23136 = r23132 + r23135; float r23137 = 3.0; float r23138 = t; float r23139 = r23137 * r23138; float r23140 = r23131 / r23139; float r23141 = r23136 - r23140; float r23142 = c; float r23143 = b; float r23144 = r23142 - r23143; float r23145 = r23138 + r23132; float r23146 = sqrt(r23145); float r23147 = z; float r23148 = r23147 / r23138; float r23149 = r23146 * r23148; float r23150 = fma(r23141, r23144, r23149); float r23151 = 0; float r23152 = r23150 + r23151; float r23153 = r23131 * r23152; float r23154 = exp(r23153); float r23155 = r23130 * r23154; float r23156 = r23129 + r23155; float r23157 = r23129 / r23156; return r23157; } double f_od(double x, double y, double z, double t, double a, double b, double c) { double r23158 = x; double r23159 = y; double r23160 = 2.0; double r23161 = a; double r23162 = 5.0; double r23163 = 6.0; double r23164 = r23162 / r23163; double r23165 = r23161 + r23164; double r23166 = 3.0; double r23167 = t; double r23168 = r23166 * r23167; double r23169 = r23160 / r23168; double r23170 = r23165 - r23169; double r23171 = c; double r23172 = b; double r23173 = r23171 - r23172; double r23174 = r23167 + r23161; double r23175 = sqrt(r23174); double r23176 = z; double r23177 = r23176 / r23167; double r23178 = r23175 * r23177; double r23179 = fma(r23170, r23173, r23178); double r23180 = 0; double r23181 = r23179 + r23180; double r23182 = r23160 * r23181; double r23183 = exp(r23182); double r23184 = r23159 * r23183; double r23185 = r23158 + r23184; double r23186 = r23158 / r23185; return r23186; } 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 r23187, r23188, r23189, r23190, r23191, r23192, r23193, r23194, r23195, r23196, r23197, r23198, r23199, r23200, r23201, r23202, r23203, r23204, r23205, r23206, r23207, r23208, r23209, r23210, r23211, r23212, r23213, r23214; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r23187); mpfr_init(r23188); mpfr_init_set_str(r23189, "2.0", 10, MPFR_RNDN); mpfr_init(r23190); mpfr_init(r23191); 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_set_str(r23200, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r23201, "6.0", 10, MPFR_RNDN); mpfr_init(r23202); mpfr_init(r23203); mpfr_init_set_str(r23204, "3.0", 10, MPFR_RNDN); mpfr_init(r23205); mpfr_init(r23206); mpfr_init(r23207); mpfr_init(r23208); mpfr_init(r23209); mpfr_init(r23210); mpfr_init(r23211); mpfr_init(r23212); mpfr_init(r23213); mpfr_init(r23214); } double f_im(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r23187, x, MPFR_RNDN); mpfr_set_d(r23188, y, MPFR_RNDN); ; mpfr_set_d(r23190, z, MPFR_RNDN); mpfr_set_d(r23191, t, MPFR_RNDN); mpfr_set_d(r23192, a, MPFR_RNDN); mpfr_add(r23193, r23191, r23192, MPFR_RNDN); mpfr_sqrt(r23194, r23193, MPFR_RNDN); mpfr_mul(r23195, r23190, r23194, MPFR_RNDN); mpfr_div(r23196, r23195, r23191, MPFR_RNDN); mpfr_set_d(r23197, b, MPFR_RNDN); mpfr_set_d(r23198, c, MPFR_RNDN); mpfr_sub(r23199, r23197, r23198, MPFR_RNDN); ; ; mpfr_div(r23202, r23200, r23201, MPFR_RNDN); mpfr_add(r23203, r23192, r23202, MPFR_RNDN); ; mpfr_mul(r23205, r23191, r23204, MPFR_RNDN); mpfr_div(r23206, r23189, r23205, MPFR_RNDN); mpfr_sub(r23207, r23203, r23206, MPFR_RNDN); mpfr_mul(r23208, r23199, r23207, MPFR_RNDN); mpfr_sub(r23209, r23196, r23208, MPFR_RNDN); mpfr_mul(r23210, r23189, r23209, MPFR_RNDN); mpfr_exp(r23211, r23210, MPFR_RNDN); mpfr_mul(r23212, r23188, r23211, MPFR_RNDN); mpfr_add(r23213, r23187, r23212, MPFR_RNDN); mpfr_div(r23214, r23187, r23213, MPFR_RNDN); return mpfr_get_d(r23214, MPFR_RNDN); } static mpfr_t r23215, r23216, r23217, r23218, r23219, r23220, r23221, r23222, r23223, r23224, r23225, r23226, r23227, r23228, r23229, r23230, r23231, r23232, r23233, r23234, r23235, r23236, r23237, r23238, r23239, r23240, r23241, r23242, r23243; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r23215); mpfr_init(r23216); mpfr_init_set_str(r23217, "2.0", 10, MPFR_RNDN); mpfr_init(r23218); mpfr_init_set_str(r23219, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r23220, "6.0", 10, MPFR_RNDN); mpfr_init(r23221); mpfr_init(r23222); mpfr_init_set_str(r23223, "3.0", 10, MPFR_RNDN); mpfr_init(r23224); mpfr_init(r23225); mpfr_init(r23226); mpfr_init(r23227); mpfr_init(r23228); mpfr_init(r23229); mpfr_init(r23230); mpfr_init(r23231); mpfr_init(r23232); mpfr_init(r23233); mpfr_init(r23234); mpfr_init(r23235); mpfr_init(r23236); mpfr_init_set_str(r23237, "0", 10, MPFR_RNDN); mpfr_init(r23238); mpfr_init(r23239); mpfr_init(r23240); mpfr_init(r23241); mpfr_init(r23242); mpfr_init(r23243); } double f_fm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r23215, x, MPFR_RNDN); mpfr_set_d(r23216, y, MPFR_RNDN); ; mpfr_set_d(r23218, a, MPFR_RNDN); ; ; mpfr_div(r23221, r23219, r23220, MPFR_RNDN); mpfr_add(r23222, r23218, r23221, MPFR_RNDN); ; mpfr_set_d(r23224, t, MPFR_RNDN); mpfr_mul(r23225, r23223, r23224, MPFR_RNDN); mpfr_div(r23226, r23217, r23225, MPFR_RNDN); mpfr_sub(r23227, r23222, r23226, MPFR_RNDN); mpfr_set_d(r23228, c, MPFR_RNDN); mpfr_set_d(r23229, b, MPFR_RNDN); mpfr_sub(r23230, r23228, r23229, MPFR_RNDN); mpfr_add(r23231, r23224, r23218, MPFR_RNDN); mpfr_sqrt(r23232, r23231, MPFR_RNDN); mpfr_set_d(r23233, z, MPFR_RNDN); mpfr_div(r23234, r23233, r23224, MPFR_RNDN); mpfr_mul(r23235, r23232, r23234, MPFR_RNDN); mpfr_fma(r23236, r23227, r23230, r23235, MPFR_RNDN); ; mpfr_add(r23238, r23236, r23237, MPFR_RNDN); mpfr_mul(r23239, r23217, r23238, MPFR_RNDN); mpfr_exp(r23240, r23239, MPFR_RNDN); mpfr_mul(r23241, r23216, r23240, MPFR_RNDN); mpfr_add(r23242, r23215, r23241, MPFR_RNDN); mpfr_div(r23243, r23215, r23242, MPFR_RNDN); return mpfr_get_d(r23243, MPFR_RNDN); } static mpfr_t r23244, r23245, r23246, r23247, r23248, r23249, r23250, r23251, r23252, r23253, r23254, r23255, r23256, r23257, r23258, r23259, r23260, r23261, r23262, r23263, r23264, r23265, r23266, r23267, r23268, r23269, r23270, r23271, r23272; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r23244); mpfr_init(r23245); mpfr_init_set_str(r23246, "2.0", 10, MPFR_RNDN); mpfr_init(r23247); mpfr_init_set_str(r23248, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r23249, "6.0", 10, MPFR_RNDN); mpfr_init(r23250); mpfr_init(r23251); mpfr_init_set_str(r23252, "3.0", 10, MPFR_RNDN); mpfr_init(r23253); mpfr_init(r23254); mpfr_init(r23255); mpfr_init(r23256); mpfr_init(r23257); mpfr_init(r23258); mpfr_init(r23259); mpfr_init(r23260); mpfr_init(r23261); mpfr_init(r23262); mpfr_init(r23263); mpfr_init(r23264); mpfr_init(r23265); mpfr_init_set_str(r23266, "0", 10, MPFR_RNDN); mpfr_init(r23267); mpfr_init(r23268); mpfr_init(r23269); mpfr_init(r23270); mpfr_init(r23271); mpfr_init(r23272); } double f_dm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r23244, x, MPFR_RNDN); mpfr_set_d(r23245, y, MPFR_RNDN); ; mpfr_set_d(r23247, a, MPFR_RNDN); ; ; mpfr_div(r23250, r23248, r23249, MPFR_RNDN); mpfr_add(r23251, r23247, r23250, MPFR_RNDN); ; mpfr_set_d(r23253, t, MPFR_RNDN); mpfr_mul(r23254, r23252, r23253, MPFR_RNDN); mpfr_div(r23255, r23246, r23254, MPFR_RNDN); mpfr_sub(r23256, r23251, r23255, MPFR_RNDN); mpfr_set_d(r23257, c, MPFR_RNDN); mpfr_set_d(r23258, b, MPFR_RNDN); mpfr_sub(r23259, r23257, r23258, MPFR_RNDN); mpfr_add(r23260, r23253, r23247, MPFR_RNDN); mpfr_sqrt(r23261, r23260, MPFR_RNDN); mpfr_set_d(r23262, z, MPFR_RNDN); mpfr_div(r23263, r23262, r23253, MPFR_RNDN); mpfr_mul(r23264, r23261, r23263, MPFR_RNDN); mpfr_fma(r23265, r23256, r23259, r23264, MPFR_RNDN); ; mpfr_add(r23267, r23265, r23266, MPFR_RNDN); mpfr_mul(r23268, r23246, r23267, MPFR_RNDN); mpfr_exp(r23269, r23268, MPFR_RNDN); mpfr_mul(r23270, r23245, r23269, MPFR_RNDN); mpfr_add(r23271, r23244, r23270, MPFR_RNDN); mpfr_div(r23272, r23244, r23271, MPFR_RNDN); return mpfr_get_d(r23272, MPFR_RNDN); }