#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 r23128 = x; float r23129 = y; float r23130 = 2.0; float r23131 = z; float r23132 = t; float r23133 = a; float r23134 = r23132 + r23133; float r23135 = sqrt(r23134); float r23136 = r23131 * r23135; float r23137 = r23136 / r23132; float r23138 = b; float r23139 = c; float r23140 = r23138 - r23139; float r23141 = 5.0; float r23142 = 6.0; float r23143 = r23141 / r23142; float r23144 = r23133 + r23143; float r23145 = 3.0; float r23146 = r23132 * r23145; float r23147 = r23130 / r23146; float r23148 = r23144 - r23147; float r23149 = r23140 * r23148; float r23150 = r23137 - r23149; float r23151 = r23130 * r23150; float r23152 = exp(r23151); float r23153 = r23129 * r23152; float r23154 = r23128 + r23153; float r23155 = r23128 / r23154; return r23155; } double f_id(double x, double y, double z, double t, double a, double b, double c) { double r23156 = x; double r23157 = y; double r23158 = 2.0; double r23159 = z; double r23160 = t; double r23161 = a; double r23162 = r23160 + r23161; double r23163 = sqrt(r23162); double r23164 = r23159 * r23163; double r23165 = r23164 / r23160; double r23166 = b; double r23167 = c; double r23168 = r23166 - r23167; double r23169 = 5.0; double r23170 = 6.0; double r23171 = r23169 / r23170; double r23172 = r23161 + r23171; double r23173 = 3.0; double r23174 = r23160 * r23173; double r23175 = r23158 / r23174; double r23176 = r23172 - r23175; double r23177 = r23168 * r23176; double r23178 = r23165 - r23177; double r23179 = r23158 * r23178; double r23180 = exp(r23179); double r23181 = r23157 * r23180; double r23182 = r23156 + r23181; double r23183 = r23156 / r23182; return r23183; } double f_of(float x, float y, float z, float t, float a, float b, float c) { float r23184 = x; float r23185 = y; float r23186 = 2.0; float r23187 = a; float r23188 = 5.0; float r23189 = 6.0; float r23190 = r23188 / r23189; float r23191 = r23187 + r23190; float r23192 = 3.0; float r23193 = t; float r23194 = r23192 * r23193; float r23195 = r23186 / r23194; float r23196 = r23191 - r23195; float r23197 = c; float r23198 = b; float r23199 = r23197 - r23198; float r23200 = r23193 + r23187; float r23201 = sqrt(r23200); float r23202 = z; float r23203 = r23202 / r23193; float r23204 = r23201 * r23203; float r23205 = fma(r23196, r23199, r23204); float r23206 = 0; float r23207 = r23205 + r23206; float r23208 = r23186 * r23207; float r23209 = exp(r23208); float r23210 = r23185 * r23209; float r23211 = r23184 + r23210; float r23212 = r23184 / r23211; return r23212; } double f_od(double x, double y, double z, double t, double a, double b, double c) { double r23213 = x; double r23214 = y; double r23215 = 2.0; double r23216 = a; double r23217 = 5.0; double r23218 = 6.0; double r23219 = r23217 / r23218; double r23220 = r23216 + r23219; double r23221 = 3.0; double r23222 = t; double r23223 = r23221 * r23222; double r23224 = r23215 / r23223; double r23225 = r23220 - r23224; double r23226 = c; double r23227 = b; double r23228 = r23226 - r23227; double r23229 = r23222 + r23216; double r23230 = sqrt(r23229); double r23231 = z; double r23232 = r23231 / r23222; double r23233 = r23230 * r23232; double r23234 = fma(r23225, r23228, r23233); double r23235 = 0; double r23236 = r23234 + r23235; double r23237 = r23215 * r23236; double r23238 = exp(r23237); double r23239 = r23214 * r23238; double r23240 = r23213 + r23239; double r23241 = r23213 / r23240; return r23241; } 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 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_im() { mpfr_set_default_prec(336); mpfr_init(r23242); mpfr_init(r23243); mpfr_init_set_str(r23244, "2.0", 10, MPFR_RNDN); mpfr_init(r23245); mpfr_init(r23246); mpfr_init(r23247); mpfr_init(r23248); mpfr_init(r23249); mpfr_init(r23250); mpfr_init(r23251); mpfr_init(r23252); mpfr_init(r23253); mpfr_init(r23254); mpfr_init_set_str(r23255, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r23256, "6.0", 10, MPFR_RNDN); mpfr_init(r23257); mpfr_init(r23258); mpfr_init_set_str(r23259, "3.0", 10, MPFR_RNDN); mpfr_init(r23260); mpfr_init(r23261); mpfr_init(r23262); mpfr_init(r23263); mpfr_init(r23264); mpfr_init(r23265); mpfr_init(r23266); mpfr_init(r23267); mpfr_init(r23268); mpfr_init(r23269); } double f_im(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r23242, x, MPFR_RNDN); mpfr_set_d(r23243, y, MPFR_RNDN); ; mpfr_set_d(r23245, z, MPFR_RNDN); mpfr_set_d(r23246, t, MPFR_RNDN); mpfr_set_d(r23247, a, MPFR_RNDN); mpfr_add(r23248, r23246, r23247, MPFR_RNDN); mpfr_sqrt(r23249, r23248, MPFR_RNDN); mpfr_mul(r23250, r23245, r23249, MPFR_RNDN); mpfr_div(r23251, r23250, r23246, MPFR_RNDN); mpfr_set_d(r23252, b, MPFR_RNDN); mpfr_set_d(r23253, c, MPFR_RNDN); mpfr_sub(r23254, r23252, r23253, MPFR_RNDN); ; ; mpfr_div(r23257, r23255, r23256, MPFR_RNDN); mpfr_add(r23258, r23247, r23257, MPFR_RNDN); ; mpfr_mul(r23260, r23246, r23259, MPFR_RNDN); mpfr_div(r23261, r23244, r23260, MPFR_RNDN); mpfr_sub(r23262, r23258, r23261, MPFR_RNDN); mpfr_mul(r23263, r23254, r23262, MPFR_RNDN); mpfr_sub(r23264, r23251, r23263, MPFR_RNDN); mpfr_mul(r23265, r23244, r23264, MPFR_RNDN); mpfr_exp(r23266, r23265, MPFR_RNDN); mpfr_mul(r23267, r23243, r23266, MPFR_RNDN); mpfr_add(r23268, r23242, r23267, MPFR_RNDN); mpfr_div(r23269, r23242, r23268, MPFR_RNDN); return mpfr_get_d(r23269, MPFR_RNDN); } static mpfr_t r23270, r23271, r23272, r23273, r23274, r23275, r23276, r23277, r23278, r23279, r23280, r23281, r23282, r23283, r23284, r23285, r23286, r23287, r23288, r23289, r23290, r23291, r23292, r23293, r23294, r23295, r23296, r23297, r23298; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r23270); mpfr_init(r23271); mpfr_init_set_str(r23272, "2.0", 10, MPFR_RNDN); mpfr_init(r23273); mpfr_init_set_str(r23274, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r23275, "6.0", 10, MPFR_RNDN); mpfr_init(r23276); mpfr_init(r23277); mpfr_init_set_str(r23278, "3.0", 10, MPFR_RNDN); mpfr_init(r23279); mpfr_init(r23280); mpfr_init(r23281); mpfr_init(r23282); mpfr_init(r23283); mpfr_init(r23284); mpfr_init(r23285); mpfr_init(r23286); mpfr_init(r23287); mpfr_init(r23288); mpfr_init(r23289); mpfr_init(r23290); mpfr_init(r23291); mpfr_init_set_str(r23292, "0", 10, MPFR_RNDN); mpfr_init(r23293); mpfr_init(r23294); mpfr_init(r23295); mpfr_init(r23296); mpfr_init(r23297); mpfr_init(r23298); } double f_fm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r23270, x, MPFR_RNDN); mpfr_set_d(r23271, y, MPFR_RNDN); ; mpfr_set_d(r23273, a, MPFR_RNDN); ; ; mpfr_div(r23276, r23274, r23275, MPFR_RNDN); mpfr_add(r23277, r23273, r23276, MPFR_RNDN); ; mpfr_set_d(r23279, t, MPFR_RNDN); mpfr_mul(r23280, r23278, r23279, MPFR_RNDN); mpfr_div(r23281, r23272, r23280, MPFR_RNDN); mpfr_sub(r23282, r23277, r23281, MPFR_RNDN); mpfr_set_d(r23283, c, MPFR_RNDN); mpfr_set_d(r23284, b, MPFR_RNDN); mpfr_sub(r23285, r23283, r23284, MPFR_RNDN); mpfr_add(r23286, r23279, r23273, MPFR_RNDN); mpfr_sqrt(r23287, r23286, MPFR_RNDN); mpfr_set_d(r23288, z, MPFR_RNDN); mpfr_div(r23289, r23288, r23279, MPFR_RNDN); mpfr_mul(r23290, r23287, r23289, MPFR_RNDN); mpfr_fma(r23291, r23282, r23285, r23290, MPFR_RNDN); ; mpfr_add(r23293, r23291, r23292, MPFR_RNDN); mpfr_mul(r23294, r23272, r23293, MPFR_RNDN); mpfr_exp(r23295, r23294, MPFR_RNDN); mpfr_mul(r23296, r23271, r23295, MPFR_RNDN); mpfr_add(r23297, r23270, r23296, MPFR_RNDN); mpfr_div(r23298, r23270, r23297, MPFR_RNDN); return mpfr_get_d(r23298, MPFR_RNDN); } static mpfr_t r23299, r23300, r23301, r23302, r23303, r23304, r23305, r23306, r23307, r23308, r23309, r23310, r23311, r23312, r23313, r23314, r23315, r23316, r23317, r23318, r23319, r23320, r23321, r23322, r23323, r23324, r23325, r23326, r23327; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r23299); mpfr_init(r23300); mpfr_init_set_str(r23301, "2.0", 10, MPFR_RNDN); mpfr_init(r23302); mpfr_init_set_str(r23303, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r23304, "6.0", 10, MPFR_RNDN); mpfr_init(r23305); mpfr_init(r23306); mpfr_init_set_str(r23307, "3.0", 10, MPFR_RNDN); mpfr_init(r23308); mpfr_init(r23309); mpfr_init(r23310); mpfr_init(r23311); mpfr_init(r23312); mpfr_init(r23313); mpfr_init(r23314); mpfr_init(r23315); mpfr_init(r23316); mpfr_init(r23317); mpfr_init(r23318); mpfr_init(r23319); mpfr_init(r23320); mpfr_init_set_str(r23321, "0", 10, MPFR_RNDN); mpfr_init(r23322); mpfr_init(r23323); mpfr_init(r23324); mpfr_init(r23325); mpfr_init(r23326); mpfr_init(r23327); } double f_dm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r23299, x, MPFR_RNDN); mpfr_set_d(r23300, y, MPFR_RNDN); ; mpfr_set_d(r23302, a, MPFR_RNDN); ; ; mpfr_div(r23305, r23303, r23304, MPFR_RNDN); mpfr_add(r23306, r23302, r23305, MPFR_RNDN); ; mpfr_set_d(r23308, t, MPFR_RNDN); mpfr_mul(r23309, r23307, r23308, MPFR_RNDN); mpfr_div(r23310, r23301, r23309, MPFR_RNDN); mpfr_sub(r23311, r23306, r23310, MPFR_RNDN); mpfr_set_d(r23312, c, MPFR_RNDN); mpfr_set_d(r23313, b, MPFR_RNDN); mpfr_sub(r23314, r23312, r23313, MPFR_RNDN); mpfr_add(r23315, r23308, r23302, MPFR_RNDN); mpfr_sqrt(r23316, r23315, MPFR_RNDN); mpfr_set_d(r23317, z, MPFR_RNDN); mpfr_div(r23318, r23317, r23308, MPFR_RNDN); mpfr_mul(r23319, r23316, r23318, MPFR_RNDN); mpfr_fma(r23320, r23311, r23314, r23319, MPFR_RNDN); ; mpfr_add(r23322, r23320, r23321, MPFR_RNDN); mpfr_mul(r23323, r23301, r23322, MPFR_RNDN); mpfr_exp(r23324, r23323, MPFR_RNDN); mpfr_mul(r23325, r23300, r23324, MPFR_RNDN); mpfr_add(r23326, r23299, r23325, MPFR_RNDN); mpfr_div(r23327, r23299, r23326, MPFR_RNDN); return mpfr_get_d(r23327, MPFR_RNDN); }