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