#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 r23043 = x; float r23044 = y; float r23045 = r23043 * r23044; float r23046 = z; float r23047 = r23045 + r23046; float r23048 = r23047 * r23044; float r23049 = 27464.7644705; float r23050 = r23048 + r23049; float r23051 = r23050 * r23044; float r23052 = 230661.510616; float r23053 = r23051 + r23052; float r23054 = r23053 * r23044; float r23055 = t; float r23056 = r23054 + r23055; float r23057 = a; float r23058 = r23044 + r23057; float r23059 = r23058 * r23044; float r23060 = b; float r23061 = r23059 + r23060; float r23062 = r23061 * r23044; float r23063 = c; float r23064 = r23062 + r23063; float r23065 = r23064 * r23044; float r23066 = i; float r23067 = r23065 + r23066; float r23068 = r23056 / r23067; return r23068; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i) { double r23069 = x; double r23070 = y; double r23071 = r23069 * r23070; double r23072 = z; double r23073 = r23071 + r23072; double r23074 = r23073 * r23070; double r23075 = 27464.7644705; double r23076 = r23074 + r23075; double r23077 = r23076 * r23070; double r23078 = 230661.510616; double r23079 = r23077 + r23078; double r23080 = r23079 * r23070; double r23081 = t; double r23082 = r23080 + r23081; double r23083 = a; double r23084 = r23070 + r23083; double r23085 = r23084 * r23070; double r23086 = b; double r23087 = r23085 + r23086; double r23088 = r23087 * r23070; double r23089 = c; double r23090 = r23088 + r23089; double r23091 = r23090 * r23070; double r23092 = i; double r23093 = r23091 + r23092; double r23094 = r23082 / r23093; return r23094; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i) { float r23095 = x; float r23096 = y; float r23097 = r23095 * r23096; float r23098 = z; float r23099 = r23097 + r23098; float r23100 = r23099 * r23096; float r23101 = 27464.7644705; float r23102 = r23100 + r23101; float r23103 = r23102 * r23096; float r23104 = 230661.510616; float r23105 = r23103 + r23104; float r23106 = r23105 * r23096; float r23107 = t; float r23108 = r23106 + r23107; float r23109 = 1; float r23110 = a; float r23111 = r23096 + r23110; float r23112 = r23111 * r23096; float r23113 = b; float r23114 = r23112 + r23113; float r23115 = r23114 * r23096; float r23116 = c; float r23117 = r23115 + r23116; float r23118 = r23117 * r23096; float r23119 = i; float r23120 = r23118 + r23119; float r23121 = r23109 / r23120; float r23122 = r23108 * r23121; return r23122; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i) { double r23123 = x; double r23124 = y; double r23125 = r23123 * r23124; double r23126 = z; double r23127 = r23125 + r23126; double r23128 = r23127 * r23124; double r23129 = 27464.7644705; double r23130 = r23128 + r23129; double r23131 = r23130 * r23124; double r23132 = 230661.510616; double r23133 = r23131 + r23132; double r23134 = r23133 * r23124; double r23135 = t; double r23136 = r23134 + r23135; double r23137 = 1; double r23138 = a; double r23139 = r23124 + r23138; double r23140 = r23139 * r23124; double r23141 = b; double r23142 = r23140 + r23141; double r23143 = r23142 * r23124; double r23144 = c; double r23145 = r23143 + r23144; double r23146 = r23145 * r23124; double r23147 = i; double r23148 = r23146 + r23147; double r23149 = r23137 / r23148; double r23150 = r23136 * r23149; return r23150; } 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 r23151, r23152, r23153, r23154, r23155, r23156, r23157, r23158, r23159, r23160, r23161, r23162, r23163, r23164, r23165, r23166, r23167, r23168, r23169, r23170, r23171, r23172, r23173, r23174, r23175, r23176; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r23151); mpfr_init(r23152); mpfr_init(r23153); mpfr_init(r23154); mpfr_init(r23155); mpfr_init(r23156); mpfr_init_set_str(r23157, "27464.7644705", 10, MPFR_RNDN); mpfr_init(r23158); mpfr_init(r23159); mpfr_init_set_str(r23160, "230661.510616", 10, MPFR_RNDN); mpfr_init(r23161); mpfr_init(r23162); mpfr_init(r23163); mpfr_init(r23164); mpfr_init(r23165); mpfr_init(r23166); mpfr_init(r23167); mpfr_init(r23168); mpfr_init(r23169); mpfr_init(r23170); mpfr_init(r23171); mpfr_init(r23172); mpfr_init(r23173); mpfr_init(r23174); mpfr_init(r23175); mpfr_init(r23176); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r23151, x, MPFR_RNDN); mpfr_set_d(r23152, y, MPFR_RNDN); mpfr_mul(r23153, r23151, r23152, MPFR_RNDN); mpfr_set_d(r23154, z, MPFR_RNDN); mpfr_add(r23155, r23153, r23154, MPFR_RNDN); mpfr_mul(r23156, r23155, r23152, MPFR_RNDN); ; mpfr_add(r23158, r23156, r23157, MPFR_RNDN); mpfr_mul(r23159, r23158, r23152, MPFR_RNDN); ; mpfr_add(r23161, r23159, r23160, MPFR_RNDN); mpfr_mul(r23162, r23161, r23152, MPFR_RNDN); mpfr_set_d(r23163, t, MPFR_RNDN); mpfr_add(r23164, r23162, r23163, MPFR_RNDN); mpfr_set_d(r23165, a, MPFR_RNDN); mpfr_add(r23166, r23152, r23165, MPFR_RNDN); mpfr_mul(r23167, r23166, r23152, MPFR_RNDN); mpfr_set_d(r23168, b, MPFR_RNDN); mpfr_add(r23169, r23167, r23168, MPFR_RNDN); mpfr_mul(r23170, r23169, r23152, MPFR_RNDN); mpfr_set_d(r23171, c, MPFR_RNDN); mpfr_add(r23172, r23170, r23171, MPFR_RNDN); mpfr_mul(r23173, r23172, r23152, MPFR_RNDN); mpfr_set_d(r23174, i, MPFR_RNDN); mpfr_add(r23175, r23173, r23174, MPFR_RNDN); mpfr_div(r23176, r23164, r23175, MPFR_RNDN); return mpfr_get_d(r23176, MPFR_RNDN); } static mpfr_t r23177, r23178, r23179, r23180, r23181, r23182, r23183, r23184, r23185, r23186, r23187, r23188, r23189, r23190, r23191, r23192, r23193, r23194, r23195, r23196, r23197, r23198, r23199, r23200, r23201, r23202, r23203, r23204; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r23177); mpfr_init(r23178); mpfr_init(r23179); mpfr_init(r23180); mpfr_init(r23181); mpfr_init(r23182); mpfr_init_set_str(r23183, "27464.7644705", 10, MPFR_RNDN); mpfr_init(r23184); mpfr_init(r23185); mpfr_init_set_str(r23186, "230661.510616", 10, MPFR_RNDN); mpfr_init(r23187); mpfr_init(r23188); mpfr_init(r23189); mpfr_init(r23190); mpfr_init_set_str(r23191, "1", 10, MPFR_RNDN); mpfr_init(r23192); mpfr_init(r23193); mpfr_init(r23194); mpfr_init(r23195); mpfr_init(r23196); mpfr_init(r23197); mpfr_init(r23198); mpfr_init(r23199); mpfr_init(r23200); mpfr_init(r23201); mpfr_init(r23202); mpfr_init(r23203); mpfr_init(r23204); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r23177, x, MPFR_RNDN); mpfr_set_d(r23178, y, MPFR_RNDN); mpfr_mul(r23179, r23177, r23178, MPFR_RNDN); mpfr_set_d(r23180, z, MPFR_RNDN); mpfr_add(r23181, r23179, r23180, MPFR_RNDN); mpfr_mul(r23182, r23181, r23178, MPFR_RNDN); ; mpfr_add(r23184, r23182, r23183, MPFR_RNDN); mpfr_mul(r23185, r23184, r23178, MPFR_RNDN); ; mpfr_add(r23187, r23185, r23186, MPFR_RNDN); mpfr_mul(r23188, r23187, r23178, MPFR_RNDN); mpfr_set_d(r23189, t, MPFR_RNDN); mpfr_add(r23190, r23188, r23189, MPFR_RNDN); ; mpfr_set_d(r23192, a, MPFR_RNDN); mpfr_add(r23193, r23178, r23192, MPFR_RNDN); mpfr_mul(r23194, r23193, r23178, MPFR_RNDN); mpfr_set_d(r23195, b, MPFR_RNDN); mpfr_add(r23196, r23194, r23195, MPFR_RNDN); mpfr_mul(r23197, r23196, r23178, MPFR_RNDN); mpfr_set_d(r23198, c, MPFR_RNDN); mpfr_add(r23199, r23197, r23198, MPFR_RNDN); mpfr_mul(r23200, r23199, r23178, MPFR_RNDN); mpfr_set_d(r23201, i, MPFR_RNDN); mpfr_add(r23202, r23200, r23201, MPFR_RNDN); mpfr_div(r23203, r23191, r23202, MPFR_RNDN); mpfr_mul(r23204, r23190, r23203, MPFR_RNDN); return mpfr_get_d(r23204, MPFR_RNDN); } static mpfr_t r23205, r23206, r23207, r23208, r23209, r23210, r23211, r23212, r23213, r23214, r23215, r23216, r23217, r23218, r23219, r23220, r23221, r23222, r23223, r23224, r23225, r23226, r23227, r23228, r23229, r23230, r23231, r23232; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r23205); mpfr_init(r23206); mpfr_init(r23207); mpfr_init(r23208); mpfr_init(r23209); mpfr_init(r23210); mpfr_init_set_str(r23211, "27464.7644705", 10, MPFR_RNDN); mpfr_init(r23212); mpfr_init(r23213); mpfr_init_set_str(r23214, "230661.510616", 10, MPFR_RNDN); mpfr_init(r23215); mpfr_init(r23216); mpfr_init(r23217); mpfr_init(r23218); mpfr_init_set_str(r23219, "1", 10, MPFR_RNDN); 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(r23228); mpfr_init(r23229); mpfr_init(r23230); mpfr_init(r23231); mpfr_init(r23232); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r23205, x, MPFR_RNDN); mpfr_set_d(r23206, y, MPFR_RNDN); mpfr_mul(r23207, r23205, r23206, MPFR_RNDN); mpfr_set_d(r23208, z, MPFR_RNDN); mpfr_add(r23209, r23207, r23208, MPFR_RNDN); mpfr_mul(r23210, r23209, r23206, MPFR_RNDN); ; mpfr_add(r23212, r23210, r23211, MPFR_RNDN); mpfr_mul(r23213, r23212, r23206, MPFR_RNDN); ; mpfr_add(r23215, r23213, r23214, MPFR_RNDN); mpfr_mul(r23216, r23215, r23206, MPFR_RNDN); mpfr_set_d(r23217, t, MPFR_RNDN); mpfr_add(r23218, r23216, r23217, MPFR_RNDN); ; mpfr_set_d(r23220, a, MPFR_RNDN); mpfr_add(r23221, r23206, r23220, MPFR_RNDN); mpfr_mul(r23222, r23221, r23206, MPFR_RNDN); mpfr_set_d(r23223, b, MPFR_RNDN); mpfr_add(r23224, r23222, r23223, MPFR_RNDN); mpfr_mul(r23225, r23224, r23206, MPFR_RNDN); mpfr_set_d(r23226, c, MPFR_RNDN); mpfr_add(r23227, r23225, r23226, MPFR_RNDN); mpfr_mul(r23228, r23227, r23206, MPFR_RNDN); mpfr_set_d(r23229, i, MPFR_RNDN); mpfr_add(r23230, r23228, r23229, MPFR_RNDN); mpfr_div(r23231, r23219, r23230, MPFR_RNDN); mpfr_mul(r23232, r23218, r23231, MPFR_RNDN); return mpfr_get_d(r23232, MPFR_RNDN); }