#include #include #include #include #include char *name = "Numeric.SpecFunctions:logGamma 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 i) { float r22141 = x; float r22142 = y; float r22143 = r22141 * r22142; float r22144 = z; float r22145 = r22143 + r22144; float r22146 = r22145 * r22142; float r22147 = 27464.7644705; float r22148 = r22146 + r22147; float r22149 = r22148 * r22142; float r22150 = 230661.510616; float r22151 = r22149 + r22150; float r22152 = r22151 * r22142; float r22153 = t; float r22154 = r22152 + r22153; float r22155 = a; float r22156 = r22142 + r22155; float r22157 = r22156 * r22142; float r22158 = b; float r22159 = r22157 + r22158; float r22160 = r22159 * r22142; float r22161 = c; float r22162 = r22160 + r22161; float r22163 = r22162 * r22142; float r22164 = i; float r22165 = r22163 + r22164; float r22166 = r22154 / r22165; return r22166; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i) { double r22167 = x; double r22168 = y; double r22169 = r22167 * r22168; double r22170 = z; double r22171 = r22169 + r22170; double r22172 = r22171 * r22168; double r22173 = 27464.7644705; double r22174 = r22172 + r22173; double r22175 = r22174 * r22168; double r22176 = 230661.510616; double r22177 = r22175 + r22176; double r22178 = r22177 * r22168; double r22179 = t; double r22180 = r22178 + r22179; double r22181 = a; double r22182 = r22168 + r22181; double r22183 = r22182 * r22168; double r22184 = b; double r22185 = r22183 + r22184; double r22186 = r22185 * r22168; double r22187 = c; double r22188 = r22186 + r22187; double r22189 = r22188 * r22168; double r22190 = i; double r22191 = r22189 + r22190; double r22192 = r22180 / r22191; return r22192; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i) { float r22193 = x; float r22194 = y; float r22195 = r22193 * r22194; float r22196 = z; float r22197 = r22195 + r22196; float r22198 = r22197 * r22194; float r22199 = 27464.7644705; float r22200 = r22198 + r22199; float r22201 = r22200 * r22194; float r22202 = 230661.510616; float r22203 = r22201 + r22202; float r22204 = r22203 * r22194; float r22205 = t; float r22206 = r22204 + r22205; float r22207 = 1; float r22208 = a; float r22209 = r22194 + r22208; float r22210 = r22209 * r22194; float r22211 = b; float r22212 = r22210 + r22211; float r22213 = r22212 * r22194; float r22214 = c; float r22215 = r22213 + r22214; float r22216 = r22215 * r22194; float r22217 = i; float r22218 = r22216 + r22217; float r22219 = r22207 / r22218; float r22220 = r22206 * r22219; return r22220; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i) { double r22221 = x; double r22222 = y; double r22223 = r22221 * r22222; double r22224 = z; double r22225 = r22223 + r22224; double r22226 = r22225 * r22222; double r22227 = 27464.7644705; double r22228 = r22226 + r22227; double r22229 = r22228 * r22222; double r22230 = 230661.510616; double r22231 = r22229 + r22230; double r22232 = r22231 * r22222; double r22233 = t; double r22234 = r22232 + r22233; double r22235 = 1; double r22236 = a; double r22237 = r22222 + r22236; double r22238 = r22237 * r22222; double r22239 = b; double r22240 = r22238 + r22239; double r22241 = r22240 * r22222; double r22242 = c; double r22243 = r22241 + r22242; double r22244 = r22243 * r22222; double r22245 = i; double r22246 = r22244 + r22245; double r22247 = r22235 / r22246; double r22248 = r22234 * r22247; return r22248; } 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 r22249, r22250, r22251, r22252, r22253, r22254, r22255, r22256, r22257, r22258, r22259, r22260, r22261, r22262, r22263, r22264, r22265, r22266, r22267, r22268, r22269, r22270, r22271, r22272, r22273, r22274; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22249); mpfr_init(r22250); mpfr_init(r22251); mpfr_init(r22252); mpfr_init(r22253); mpfr_init(r22254); mpfr_init_set_str(r22255, "27464.7644705", 10, MPFR_RNDN); mpfr_init(r22256); mpfr_init(r22257); mpfr_init_set_str(r22258, "230661.510616", 10, MPFR_RNDN); mpfr_init(r22259); mpfr_init(r22260); mpfr_init(r22261); mpfr_init(r22262); mpfr_init(r22263); mpfr_init(r22264); mpfr_init(r22265); mpfr_init(r22266); mpfr_init(r22267); mpfr_init(r22268); mpfr_init(r22269); mpfr_init(r22270); mpfr_init(r22271); mpfr_init(r22272); mpfr_init(r22273); mpfr_init(r22274); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r22249, x, MPFR_RNDN); mpfr_set_d(r22250, y, MPFR_RNDN); mpfr_mul(r22251, r22249, r22250, MPFR_RNDN); mpfr_set_d(r22252, z, MPFR_RNDN); mpfr_add(r22253, r22251, r22252, MPFR_RNDN); mpfr_mul(r22254, r22253, r22250, MPFR_RNDN); ; mpfr_add(r22256, r22254, r22255, MPFR_RNDN); mpfr_mul(r22257, r22256, r22250, MPFR_RNDN); ; mpfr_add(r22259, r22257, r22258, MPFR_RNDN); mpfr_mul(r22260, r22259, r22250, MPFR_RNDN); mpfr_set_d(r22261, t, MPFR_RNDN); mpfr_add(r22262, r22260, r22261, MPFR_RNDN); mpfr_set_d(r22263, a, MPFR_RNDN); mpfr_add(r22264, r22250, r22263, MPFR_RNDN); mpfr_mul(r22265, r22264, r22250, MPFR_RNDN); mpfr_set_d(r22266, b, MPFR_RNDN); mpfr_add(r22267, r22265, r22266, MPFR_RNDN); mpfr_mul(r22268, r22267, r22250, MPFR_RNDN); mpfr_set_d(r22269, c, MPFR_RNDN); mpfr_add(r22270, r22268, r22269, MPFR_RNDN); mpfr_mul(r22271, r22270, r22250, MPFR_RNDN); mpfr_set_d(r22272, i, MPFR_RNDN); mpfr_add(r22273, r22271, r22272, MPFR_RNDN); mpfr_div(r22274, r22262, r22273, MPFR_RNDN); return mpfr_get_d(r22274, MPFR_RNDN); } static mpfr_t r22275, r22276, r22277, r22278, r22279, r22280, r22281, r22282, r22283, r22284, r22285, r22286, r22287, r22288, r22289, r22290, r22291, r22292, r22293, r22294, r22295, r22296, r22297, r22298, r22299, r22300, r22301, r22302; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22275); mpfr_init(r22276); mpfr_init(r22277); mpfr_init(r22278); mpfr_init(r22279); mpfr_init(r22280); mpfr_init_set_str(r22281, "27464.7644705", 10, MPFR_RNDN); mpfr_init(r22282); mpfr_init(r22283); mpfr_init_set_str(r22284, "230661.510616", 10, MPFR_RNDN); mpfr_init(r22285); mpfr_init(r22286); mpfr_init(r22287); mpfr_init(r22288); mpfr_init_set_str(r22289, "1", 10, MPFR_RNDN); mpfr_init(r22290); mpfr_init(r22291); mpfr_init(r22292); mpfr_init(r22293); mpfr_init(r22294); mpfr_init(r22295); mpfr_init(r22296); mpfr_init(r22297); mpfr_init(r22298); mpfr_init(r22299); mpfr_init(r22300); mpfr_init(r22301); mpfr_init(r22302); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r22275, x, MPFR_RNDN); mpfr_set_d(r22276, y, MPFR_RNDN); mpfr_mul(r22277, r22275, r22276, MPFR_RNDN); mpfr_set_d(r22278, z, MPFR_RNDN); mpfr_add(r22279, r22277, r22278, MPFR_RNDN); mpfr_mul(r22280, r22279, r22276, MPFR_RNDN); ; mpfr_add(r22282, r22280, r22281, MPFR_RNDN); mpfr_mul(r22283, r22282, r22276, MPFR_RNDN); ; mpfr_add(r22285, r22283, r22284, MPFR_RNDN); mpfr_mul(r22286, r22285, r22276, MPFR_RNDN); mpfr_set_d(r22287, t, MPFR_RNDN); mpfr_add(r22288, r22286, r22287, MPFR_RNDN); ; mpfr_set_d(r22290, a, MPFR_RNDN); mpfr_add(r22291, r22276, r22290, MPFR_RNDN); mpfr_mul(r22292, r22291, r22276, MPFR_RNDN); mpfr_set_d(r22293, b, MPFR_RNDN); mpfr_add(r22294, r22292, r22293, MPFR_RNDN); mpfr_mul(r22295, r22294, r22276, MPFR_RNDN); mpfr_set_d(r22296, c, MPFR_RNDN); mpfr_add(r22297, r22295, r22296, MPFR_RNDN); mpfr_mul(r22298, r22297, r22276, MPFR_RNDN); mpfr_set_d(r22299, i, MPFR_RNDN); mpfr_add(r22300, r22298, r22299, MPFR_RNDN); mpfr_div(r22301, r22289, r22300, MPFR_RNDN); mpfr_mul(r22302, r22288, r22301, MPFR_RNDN); return mpfr_get_d(r22302, MPFR_RNDN); } static mpfr_t r22303, r22304, r22305, r22306, r22307, r22308, r22309, r22310, r22311, r22312, r22313, r22314, r22315, r22316, r22317, r22318, r22319, r22320, r22321, r22322, r22323, r22324, r22325, r22326, r22327, r22328, r22329, r22330; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22303); mpfr_init(r22304); mpfr_init(r22305); mpfr_init(r22306); mpfr_init(r22307); mpfr_init(r22308); mpfr_init_set_str(r22309, "27464.7644705", 10, MPFR_RNDN); mpfr_init(r22310); mpfr_init(r22311); mpfr_init_set_str(r22312, "230661.510616", 10, MPFR_RNDN); mpfr_init(r22313); mpfr_init(r22314); mpfr_init(r22315); mpfr_init(r22316); mpfr_init_set_str(r22317, "1", 10, MPFR_RNDN); mpfr_init(r22318); mpfr_init(r22319); mpfr_init(r22320); mpfr_init(r22321); mpfr_init(r22322); mpfr_init(r22323); mpfr_init(r22324); mpfr_init(r22325); mpfr_init(r22326); mpfr_init(r22327); mpfr_init(r22328); mpfr_init(r22329); mpfr_init(r22330); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r22303, x, MPFR_RNDN); mpfr_set_d(r22304, y, MPFR_RNDN); mpfr_mul(r22305, r22303, r22304, MPFR_RNDN); mpfr_set_d(r22306, z, MPFR_RNDN); mpfr_add(r22307, r22305, r22306, MPFR_RNDN); mpfr_mul(r22308, r22307, r22304, MPFR_RNDN); ; mpfr_add(r22310, r22308, r22309, MPFR_RNDN); mpfr_mul(r22311, r22310, r22304, MPFR_RNDN); ; mpfr_add(r22313, r22311, r22312, MPFR_RNDN); mpfr_mul(r22314, r22313, r22304, MPFR_RNDN); mpfr_set_d(r22315, t, MPFR_RNDN); mpfr_add(r22316, r22314, r22315, MPFR_RNDN); ; mpfr_set_d(r22318, a, MPFR_RNDN); mpfr_add(r22319, r22304, r22318, MPFR_RNDN); mpfr_mul(r22320, r22319, r22304, MPFR_RNDN); mpfr_set_d(r22321, b, MPFR_RNDN); mpfr_add(r22322, r22320, r22321, MPFR_RNDN); mpfr_mul(r22323, r22322, r22304, MPFR_RNDN); mpfr_set_d(r22324, c, MPFR_RNDN); mpfr_add(r22325, r22323, r22324, MPFR_RNDN); mpfr_mul(r22326, r22325, r22304, MPFR_RNDN); mpfr_set_d(r22327, i, MPFR_RNDN); mpfr_add(r22328, r22326, r22327, MPFR_RNDN); mpfr_div(r22329, r22317, r22328, MPFR_RNDN); mpfr_mul(r22330, r22316, r22329, MPFR_RNDN); return mpfr_get_d(r22330, MPFR_RNDN); }