#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 r9136 = x; float r9137 = y; float r9138 = r9136 * r9137; float r9139 = z; float r9140 = r9138 + r9139; float r9141 = r9140 * r9137; float r9142 = 27464.7644705; float r9143 = r9141 + r9142; float r9144 = r9143 * r9137; float r9145 = 230661.510616; float r9146 = r9144 + r9145; float r9147 = r9146 * r9137; float r9148 = t; float r9149 = r9147 + r9148; float r9150 = a; float r9151 = r9137 + r9150; float r9152 = r9151 * r9137; float r9153 = b; float r9154 = r9152 + r9153; float r9155 = r9154 * r9137; float r9156 = c; float r9157 = r9155 + r9156; float r9158 = r9157 * r9137; float r9159 = i; float r9160 = r9158 + r9159; float r9161 = r9149 / r9160; return r9161; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i) { double r9162 = x; double r9163 = y; double r9164 = r9162 * r9163; double r9165 = z; double r9166 = r9164 + r9165; double r9167 = r9166 * r9163; double r9168 = 27464.7644705; double r9169 = r9167 + r9168; double r9170 = r9169 * r9163; double r9171 = 230661.510616; double r9172 = r9170 + r9171; double r9173 = r9172 * r9163; double r9174 = t; double r9175 = r9173 + r9174; double r9176 = a; double r9177 = r9163 + r9176; double r9178 = r9177 * r9163; double r9179 = b; double r9180 = r9178 + r9179; double r9181 = r9180 * r9163; double r9182 = c; double r9183 = r9181 + r9182; double r9184 = r9183 * r9163; double r9185 = i; double r9186 = r9184 + r9185; double r9187 = r9175 / r9186; return r9187; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i) { float r9188 = y; float r9189 = r9188 * r9188; float r9190 = x; float r9191 = z; float r9192 = fma(r9188, r9190, r9191); float r9193 = 27464.7644705; float r9194 = 230661.510616; float r9195 = fma(r9188, r9193, r9194); float r9196 = fma(r9189, r9192, r9195); float r9197 = t; float r9198 = fma(r9188, r9196, r9197); float r9199 = a; float r9200 = r9188 + r9199; float r9201 = b; float r9202 = fma(r9200, r9188, r9201); float r9203 = c; float r9204 = i; float r9205 = fma(r9188, r9203, r9204); float r9206 = fma(r9189, r9202, r9205); float r9207 = r9198 / r9206; return r9207; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i) { double r9208 = y; double r9209 = r9208 * r9208; double r9210 = x; double r9211 = z; double r9212 = fma(r9208, r9210, r9211); double r9213 = 27464.7644705; double r9214 = 230661.510616; double r9215 = fma(r9208, r9213, r9214); double r9216 = fma(r9209, r9212, r9215); double r9217 = t; double r9218 = fma(r9208, r9216, r9217); double r9219 = a; double r9220 = r9208 + r9219; double r9221 = b; double r9222 = fma(r9220, r9208, r9221); double r9223 = c; double r9224 = i; double r9225 = fma(r9208, r9223, r9224); double r9226 = fma(r9209, r9222, r9225); double r9227 = r9218 / r9226; return r9227; } 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 r9228, r9229, r9230, r9231, r9232, r9233, r9234, r9235, r9236, r9237, r9238, r9239, r9240, r9241, r9242, r9243, r9244, r9245, r9246, r9247, r9248, r9249, r9250, r9251, r9252, r9253; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r9228); mpfr_init(r9229); mpfr_init(r9230); mpfr_init(r9231); mpfr_init(r9232); mpfr_init(r9233); mpfr_init_set_str(r9234, "27464.7644705", 10, MPFR_RNDN); mpfr_init(r9235); mpfr_init(r9236); mpfr_init_set_str(r9237, "230661.510616", 10, MPFR_RNDN); mpfr_init(r9238); 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(r9249); 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, double b, double c, double i) { mpfr_set_d(r9228, x, MPFR_RNDN); mpfr_set_d(r9229, y, MPFR_RNDN); mpfr_mul(r9230, r9228, r9229, MPFR_RNDN); mpfr_set_d(r9231, z, MPFR_RNDN); mpfr_add(r9232, r9230, r9231, MPFR_RNDN); mpfr_mul(r9233, r9232, r9229, MPFR_RNDN); ; mpfr_add(r9235, r9233, r9234, MPFR_RNDN); mpfr_mul(r9236, r9235, r9229, MPFR_RNDN); ; mpfr_add(r9238, r9236, r9237, MPFR_RNDN); mpfr_mul(r9239, r9238, r9229, MPFR_RNDN); mpfr_set_d(r9240, t, MPFR_RNDN); mpfr_add(r9241, r9239, r9240, MPFR_RNDN); mpfr_set_d(r9242, a, MPFR_RNDN); mpfr_add(r9243, r9229, r9242, MPFR_RNDN); mpfr_mul(r9244, r9243, r9229, MPFR_RNDN); mpfr_set_d(r9245, b, MPFR_RNDN); mpfr_add(r9246, r9244, r9245, MPFR_RNDN); mpfr_mul(r9247, r9246, r9229, MPFR_RNDN); mpfr_set_d(r9248, c, MPFR_RNDN); mpfr_add(r9249, r9247, r9248, MPFR_RNDN); mpfr_mul(r9250, r9249, r9229, MPFR_RNDN); mpfr_set_d(r9251, i, MPFR_RNDN); mpfr_add(r9252, r9250, r9251, MPFR_RNDN); mpfr_div(r9253, r9241, 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, r9273; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r9254); mpfr_init(r9255); mpfr_init(r9256); mpfr_init(r9257); mpfr_init(r9258); mpfr_init_set_str(r9259, "27464.7644705", 10, MPFR_RNDN); mpfr_init_set_str(r9260, "230661.510616", 10, MPFR_RNDN); mpfr_init(r9261); mpfr_init(r9262); mpfr_init(r9263); mpfr_init(r9264); mpfr_init(r9265); mpfr_init(r9266); mpfr_init(r9267); mpfr_init(r9268); mpfr_init(r9269); mpfr_init(r9270); mpfr_init(r9271); mpfr_init(r9272); mpfr_init(r9273); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r9254, y, MPFR_RNDN); mpfr_mul(r9255, r9254, r9254, MPFR_RNDN); mpfr_set_d(r9256, x, MPFR_RNDN); mpfr_set_d(r9257, z, MPFR_RNDN); mpfr_fma(r9258, r9254, r9256, r9257, MPFR_RNDN); ; ; mpfr_fma(r9261, r9254, r9259, r9260, MPFR_RNDN); mpfr_fma(r9262, r9255, r9258, r9261, MPFR_RNDN); mpfr_set_d(r9263, t, MPFR_RNDN); mpfr_fma(r9264, r9254, r9262, r9263, MPFR_RNDN); mpfr_set_d(r9265, a, MPFR_RNDN); mpfr_add(r9266, r9254, r9265, MPFR_RNDN); mpfr_set_d(r9267, b, MPFR_RNDN); mpfr_fma(r9268, r9266, r9254, r9267, MPFR_RNDN); mpfr_set_d(r9269, c, MPFR_RNDN); mpfr_set_d(r9270, i, MPFR_RNDN); mpfr_fma(r9271, r9254, r9269, r9270, MPFR_RNDN); mpfr_fma(r9272, r9255, r9268, r9271, MPFR_RNDN); mpfr_div(r9273, r9264, r9272, MPFR_RNDN); return mpfr_get_d(r9273, MPFR_RNDN); } static mpfr_t r9274, r9275, r9276, r9277, r9278, r9279, r9280, r9281, r9282, r9283, r9284, r9285, r9286, r9287, r9288, r9289, r9290, r9291, r9292, r9293; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r9274); mpfr_init(r9275); mpfr_init(r9276); mpfr_init(r9277); mpfr_init(r9278); mpfr_init_set_str(r9279, "27464.7644705", 10, MPFR_RNDN); mpfr_init_set_str(r9280, "230661.510616", 10, MPFR_RNDN); mpfr_init(r9281); mpfr_init(r9282); mpfr_init(r9283); mpfr_init(r9284); mpfr_init(r9285); mpfr_init(r9286); mpfr_init(r9287); mpfr_init(r9288); mpfr_init(r9289); mpfr_init(r9290); mpfr_init(r9291); mpfr_init(r9292); mpfr_init(r9293); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r9274, y, MPFR_RNDN); mpfr_mul(r9275, r9274, r9274, MPFR_RNDN); mpfr_set_d(r9276, x, MPFR_RNDN); mpfr_set_d(r9277, z, MPFR_RNDN); mpfr_fma(r9278, r9274, r9276, r9277, MPFR_RNDN); ; ; mpfr_fma(r9281, r9274, r9279, r9280, MPFR_RNDN); mpfr_fma(r9282, r9275, r9278, r9281, MPFR_RNDN); mpfr_set_d(r9283, t, MPFR_RNDN); mpfr_fma(r9284, r9274, r9282, r9283, MPFR_RNDN); mpfr_set_d(r9285, a, MPFR_RNDN); mpfr_add(r9286, r9274, r9285, MPFR_RNDN); mpfr_set_d(r9287, b, MPFR_RNDN); mpfr_fma(r9288, r9286, r9274, r9287, MPFR_RNDN); mpfr_set_d(r9289, c, MPFR_RNDN); mpfr_set_d(r9290, i, MPFR_RNDN); mpfr_fma(r9291, r9274, r9289, r9290, MPFR_RNDN); mpfr_fma(r9292, r9275, r9288, r9291, MPFR_RNDN); mpfr_div(r9293, r9284, r9292, MPFR_RNDN); return mpfr_get_d(r9293, MPFR_RNDN); }