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