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