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