#include #include #include #include #include char *name = "Numeric.SpecFunctions:logGammaL from math-functions-0.1.5.2"; double f_if(float x, float y, float z, float t, float a) { float r9171 = x; float r9172 = y; float r9173 = r9171 + r9172; float r9174 = log(r9173); float r9175 = z; float r9176 = log(r9175); float r9177 = r9174 + r9176; float r9178 = t; float r9179 = r9177 - r9178; float r9180 = a; float r9181 = 0.5; float r9182 = r9180 - r9181; float r9183 = log(r9178); float r9184 = r9182 * r9183; float r9185 = r9179 + r9184; return r9185; } double f_id(double x, double y, double z, double t, double a) { double r9186 = x; double r9187 = y; double r9188 = r9186 + r9187; double r9189 = log(r9188); double r9190 = z; double r9191 = log(r9190); double r9192 = r9189 + r9191; double r9193 = t; double r9194 = r9192 - r9193; double r9195 = a; double r9196 = 0.5; double r9197 = r9195 - r9196; double r9198 = log(r9193); double r9199 = r9197 * r9198; double r9200 = r9194 + r9199; return r9200; } double f_of(float x, float y, float z, float t, float a) { float r9201 = x; float r9202 = y; float r9203 = r9201 + r9202; float r9204 = log(r9203); float r9205 = r9204 * r9204; float r9206 = z; float r9207 = log(r9206); float r9208 = r9207 * r9207; float r9209 = r9205 - r9208; float r9210 = r9204 - r9207; float r9211 = r9209 / r9210; float r9212 = t; float r9213 = r9211 - r9212; float r9214 = a; float r9215 = 0.5; float r9216 = r9214 - r9215; float r9217 = log(r9212); float r9218 = r9216 * r9217; float r9219 = r9213 + r9218; return r9219; } double f_od(double x, double y, double z, double t, double a) { double r9220 = x; double r9221 = y; double r9222 = r9220 + r9221; double r9223 = log(r9222); double r9224 = r9223 * r9223; double r9225 = z; double r9226 = log(r9225); double r9227 = r9226 * r9226; double r9228 = r9224 - r9227; double r9229 = r9223 - r9226; double r9230 = r9228 / r9229; double r9231 = t; double r9232 = r9230 - r9231; double r9233 = a; double r9234 = 0.5; double r9235 = r9233 - r9234; double r9236 = log(r9231); double r9237 = r9235 * r9236; double r9238 = r9232 + r9237; return r9238; } 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 r9239, r9240, r9241, r9242, r9243, r9244, r9245, r9246, r9247, r9248, r9249, r9250, r9251, r9252, r9253; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r9239); mpfr_init(r9240); mpfr_init(r9241); mpfr_init(r9242); mpfr_init(r9243); mpfr_init(r9244); mpfr_init(r9245); mpfr_init(r9246); mpfr_init(r9247); mpfr_init(r9248); mpfr_init_set_str(r9249, "0.5", 10, MPFR_RNDN); mpfr_init(r9250); mpfr_init(r9251); mpfr_init(r9252); mpfr_init(r9253); } double f_im(double x, double y, double z, double t, double a) { mpfr_set_d(r9239, x, MPFR_RNDN); mpfr_set_d(r9240, y, MPFR_RNDN); mpfr_add(r9241, r9239, r9240, MPFR_RNDN); mpfr_log(r9242, r9241, MPFR_RNDN); mpfr_set_d(r9243, z, MPFR_RNDN); mpfr_log(r9244, r9243, MPFR_RNDN); mpfr_add(r9245, r9242, r9244, MPFR_RNDN); mpfr_set_d(r9246, t, MPFR_RNDN); mpfr_sub(r9247, r9245, r9246, MPFR_RNDN); mpfr_set_d(r9248, a, MPFR_RNDN); ; mpfr_sub(r9250, r9248, r9249, MPFR_RNDN); mpfr_log(r9251, r9246, MPFR_RNDN); mpfr_mul(r9252, r9250, r9251, MPFR_RNDN); mpfr_add(r9253, r9247, r9252, MPFR_RNDN); return mpfr_get_d(r9253, MPFR_RNDN); } static mpfr_t r9254, r9255, r9256, r9257, r9258, r9259, r9260, r9261, r9262, r9263, r9264, r9265, r9266, r9267, r9268, r9269, r9270, r9271, r9272; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r9254); mpfr_init(r9255); mpfr_init(r9256); mpfr_init(r9257); mpfr_init(r9258); mpfr_init(r9259); mpfr_init(r9260); mpfr_init(r9261); mpfr_init(r9262); mpfr_init(r9263); mpfr_init(r9264); mpfr_init(r9265); mpfr_init(r9266); mpfr_init(r9267); mpfr_init_set_str(r9268, "0.5", 10, MPFR_RNDN); mpfr_init(r9269); mpfr_init(r9270); mpfr_init(r9271); mpfr_init(r9272); } double f_fm(double x, double y, double z, double t, double a) { mpfr_set_d(r9254, x, MPFR_RNDN); mpfr_set_d(r9255, y, MPFR_RNDN); mpfr_add(r9256, r9254, r9255, MPFR_RNDN); mpfr_log(r9257, r9256, MPFR_RNDN); mpfr_mul(r9258, r9257, r9257, MPFR_RNDN); mpfr_set_d(r9259, z, MPFR_RNDN); mpfr_log(r9260, r9259, MPFR_RNDN); mpfr_mul(r9261, r9260, r9260, MPFR_RNDN); mpfr_sub(r9262, r9258, r9261, MPFR_RNDN); mpfr_sub(r9263, r9257, r9260, MPFR_RNDN); mpfr_div(r9264, r9262, r9263, MPFR_RNDN); mpfr_set_d(r9265, t, MPFR_RNDN); mpfr_sub(r9266, r9264, r9265, MPFR_RNDN); mpfr_set_d(r9267, a, MPFR_RNDN); ; mpfr_sub(r9269, r9267, r9268, MPFR_RNDN); mpfr_log(r9270, r9265, MPFR_RNDN); mpfr_mul(r9271, r9269, r9270, MPFR_RNDN); mpfr_add(r9272, r9266, r9271, MPFR_RNDN); return mpfr_get_d(r9272, MPFR_RNDN); } static mpfr_t r9273, r9274, r9275, r9276, r9277, r9278, r9279, r9280, r9281, r9282, r9283, r9284, r9285, r9286, r9287, r9288, r9289, r9290, r9291; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r9273); mpfr_init(r9274); mpfr_init(r9275); mpfr_init(r9276); mpfr_init(r9277); mpfr_init(r9278); mpfr_init(r9279); mpfr_init(r9280); mpfr_init(r9281); mpfr_init(r9282); mpfr_init(r9283); mpfr_init(r9284); mpfr_init(r9285); mpfr_init(r9286); mpfr_init_set_str(r9287, "0.5", 10, MPFR_RNDN); mpfr_init(r9288); mpfr_init(r9289); mpfr_init(r9290); mpfr_init(r9291); } double f_dm(double x, double y, double z, double t, double a) { mpfr_set_d(r9273, x, MPFR_RNDN); mpfr_set_d(r9274, y, MPFR_RNDN); mpfr_add(r9275, r9273, r9274, MPFR_RNDN); mpfr_log(r9276, r9275, MPFR_RNDN); mpfr_mul(r9277, r9276, r9276, MPFR_RNDN); mpfr_set_d(r9278, z, MPFR_RNDN); mpfr_log(r9279, r9278, MPFR_RNDN); mpfr_mul(r9280, r9279, r9279, MPFR_RNDN); mpfr_sub(r9281, r9277, r9280, MPFR_RNDN); mpfr_sub(r9282, r9276, r9279, MPFR_RNDN); mpfr_div(r9283, r9281, r9282, MPFR_RNDN); mpfr_set_d(r9284, t, MPFR_RNDN); mpfr_sub(r9285, r9283, r9284, MPFR_RNDN); mpfr_set_d(r9286, a, MPFR_RNDN); ; mpfr_sub(r9288, r9286, r9287, MPFR_RNDN); mpfr_log(r9289, r9284, MPFR_RNDN); mpfr_mul(r9290, r9288, r9289, MPFR_RNDN); mpfr_add(r9291, r9285, r9290, MPFR_RNDN); return mpfr_get_d(r9291, MPFR_RNDN); }