#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 r23252 = x; float r23253 = y; float r23254 = 2.0; float r23255 = z; float r23256 = t; float r23257 = a; float r23258 = r23256 + r23257; float r23259 = sqrt(r23258); float r23260 = r23255 * r23259; float r23261 = r23260 / r23256; float r23262 = b; float r23263 = c; float r23264 = r23262 - r23263; float r23265 = 5.0; float r23266 = 6.0; float r23267 = r23265 / r23266; float r23268 = r23257 + r23267; float r23269 = 3.0; float r23270 = r23256 * r23269; float r23271 = r23254 / r23270; float r23272 = r23268 - r23271; float r23273 = r23264 * r23272; float r23274 = r23261 - r23273; float r23275 = r23254 * r23274; float r23276 = exp(r23275); float r23277 = r23253 * r23276; float r23278 = r23252 + r23277; float r23279 = r23252 / r23278; return r23279; } double f_id(double x, double y, double z, double t, double a, double b, double c) { double r23280 = x; double r23281 = y; double r23282 = 2.0; double r23283 = z; double r23284 = t; double r23285 = a; double r23286 = r23284 + r23285; double r23287 = sqrt(r23286); double r23288 = r23283 * r23287; double r23289 = r23288 / r23284; double r23290 = b; double r23291 = c; double r23292 = r23290 - r23291; double r23293 = 5.0; double r23294 = 6.0; double r23295 = r23293 / r23294; double r23296 = r23285 + r23295; double r23297 = 3.0; double r23298 = r23284 * r23297; double r23299 = r23282 / r23298; double r23300 = r23296 - r23299; double r23301 = r23292 * r23300; double r23302 = r23289 - r23301; double r23303 = r23282 * r23302; double r23304 = exp(r23303); double r23305 = r23281 * r23304; double r23306 = r23280 + r23305; double r23307 = r23280 / r23306; return r23307; } double f_of(float x, float y, float z, float t, float a, float b, float c) { float r23308 = x; float r23309 = y; float r23310 = 2.0; float r23311 = z; float r23312 = t; float r23313 = a; float r23314 = r23312 + r23313; float r23315 = sqrt(r23314); float r23316 = r23312 / r23315; float r23317 = r23311 / r23316; float r23318 = b; float r23319 = c; float r23320 = r23318 - r23319; float r23321 = 5.0; float r23322 = 6.0; float r23323 = r23321 / r23322; float r23324 = r23313 + r23323; float r23325 = 3.0; float r23326 = r23312 * r23325; float r23327 = r23310 / r23326; float r23328 = r23324 - r23327; float r23329 = r23320 * r23328; float r23330 = r23317 - r23329; float r23331 = r23310 * r23330; float r23332 = exp(r23331); float r23333 = r23309 * r23332; float r23334 = r23308 + r23333; float r23335 = r23308 / r23334; return r23335; } double f_od(double x, double y, double z, double t, double a, double b, double c) { double r23336 = x; double r23337 = y; double r23338 = 2.0; double r23339 = z; double r23340 = t; double r23341 = a; double r23342 = r23340 + r23341; double r23343 = sqrt(r23342); double r23344 = r23340 / r23343; double r23345 = r23339 / r23344; double r23346 = b; double r23347 = c; double r23348 = r23346 - r23347; double r23349 = 5.0; double r23350 = 6.0; double r23351 = r23349 / r23350; double r23352 = r23341 + r23351; double r23353 = 3.0; double r23354 = r23340 * r23353; double r23355 = r23338 / r23354; double r23356 = r23352 - r23355; double r23357 = r23348 * r23356; double r23358 = r23345 - r23357; double r23359 = r23338 * r23358; double r23360 = exp(r23359); double r23361 = r23337 * r23360; double r23362 = r23336 + r23361; double r23363 = r23336 / r23362; return r23363; } 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 r23364, r23365, r23366, r23367, r23368, r23369, r23370, r23371, r23372, r23373, r23374, r23375, r23376, r23377, r23378, r23379, r23380, r23381, r23382, r23383, r23384, r23385, r23386, r23387, r23388, r23389, r23390, r23391; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r23364); mpfr_init(r23365); mpfr_init_set_str(r23366, "2.0", 10, MPFR_RNDN); mpfr_init(r23367); mpfr_init(r23368); mpfr_init(r23369); mpfr_init(r23370); mpfr_init(r23371); mpfr_init(r23372); mpfr_init(r23373); mpfr_init(r23374); mpfr_init(r23375); mpfr_init(r23376); mpfr_init_set_str(r23377, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r23378, "6.0", 10, MPFR_RNDN); mpfr_init(r23379); mpfr_init(r23380); mpfr_init_set_str(r23381, "3.0", 10, MPFR_RNDN); mpfr_init(r23382); mpfr_init(r23383); mpfr_init(r23384); mpfr_init(r23385); mpfr_init(r23386); mpfr_init(r23387); mpfr_init(r23388); mpfr_init(r23389); mpfr_init(r23390); mpfr_init(r23391); } double f_im(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r23364, x, MPFR_RNDN); mpfr_set_d(r23365, y, MPFR_RNDN); ; mpfr_set_d(r23367, z, MPFR_RNDN); mpfr_set_d(r23368, t, MPFR_RNDN); mpfr_set_d(r23369, a, MPFR_RNDN); mpfr_add(r23370, r23368, r23369, MPFR_RNDN); mpfr_sqrt(r23371, r23370, MPFR_RNDN); mpfr_mul(r23372, r23367, r23371, MPFR_RNDN); mpfr_div(r23373, r23372, r23368, MPFR_RNDN); mpfr_set_d(r23374, b, MPFR_RNDN); mpfr_set_d(r23375, c, MPFR_RNDN); mpfr_sub(r23376, r23374, r23375, MPFR_RNDN); ; ; mpfr_div(r23379, r23377, r23378, MPFR_RNDN); mpfr_add(r23380, r23369, r23379, MPFR_RNDN); ; mpfr_mul(r23382, r23368, r23381, MPFR_RNDN); mpfr_div(r23383, r23366, r23382, MPFR_RNDN); mpfr_sub(r23384, r23380, r23383, MPFR_RNDN); mpfr_mul(r23385, r23376, r23384, MPFR_RNDN); mpfr_sub(r23386, r23373, r23385, MPFR_RNDN); mpfr_mul(r23387, r23366, r23386, MPFR_RNDN); mpfr_exp(r23388, r23387, MPFR_RNDN); mpfr_mul(r23389, r23365, r23388, MPFR_RNDN); mpfr_add(r23390, r23364, r23389, MPFR_RNDN); mpfr_div(r23391, r23364, r23390, MPFR_RNDN); return mpfr_get_d(r23391, MPFR_RNDN); } static mpfr_t r23392, r23393, r23394, r23395, r23396, r23397, r23398, r23399, r23400, r23401, r23402, r23403, r23404, r23405, r23406, r23407, r23408, r23409, r23410, r23411, r23412, r23413, r23414, r23415, r23416, r23417, r23418, r23419; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r23392); mpfr_init(r23393); mpfr_init_set_str(r23394, "2.0", 10, MPFR_RNDN); mpfr_init(r23395); mpfr_init(r23396); mpfr_init(r23397); mpfr_init(r23398); mpfr_init(r23399); mpfr_init(r23400); mpfr_init(r23401); mpfr_init(r23402); mpfr_init(r23403); mpfr_init(r23404); mpfr_init_set_str(r23405, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r23406, "6.0", 10, MPFR_RNDN); mpfr_init(r23407); mpfr_init(r23408); mpfr_init_set_str(r23409, "3.0", 10, MPFR_RNDN); mpfr_init(r23410); mpfr_init(r23411); mpfr_init(r23412); mpfr_init(r23413); mpfr_init(r23414); mpfr_init(r23415); mpfr_init(r23416); mpfr_init(r23417); mpfr_init(r23418); mpfr_init(r23419); } double f_fm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r23392, x, MPFR_RNDN); mpfr_set_d(r23393, y, MPFR_RNDN); ; mpfr_set_d(r23395, z, MPFR_RNDN); mpfr_set_d(r23396, t, MPFR_RNDN); mpfr_set_d(r23397, a, MPFR_RNDN); mpfr_add(r23398, r23396, r23397, MPFR_RNDN); mpfr_sqrt(r23399, r23398, MPFR_RNDN); mpfr_div(r23400, r23396, r23399, MPFR_RNDN); mpfr_div(r23401, r23395, r23400, MPFR_RNDN); mpfr_set_d(r23402, b, MPFR_RNDN); mpfr_set_d(r23403, c, MPFR_RNDN); mpfr_sub(r23404, r23402, r23403, MPFR_RNDN); ; ; mpfr_div(r23407, r23405, r23406, MPFR_RNDN); mpfr_add(r23408, r23397, r23407, MPFR_RNDN); ; mpfr_mul(r23410, r23396, r23409, MPFR_RNDN); mpfr_div(r23411, r23394, r23410, MPFR_RNDN); mpfr_sub(r23412, r23408, r23411, MPFR_RNDN); mpfr_mul(r23413, r23404, r23412, MPFR_RNDN); mpfr_sub(r23414, r23401, r23413, MPFR_RNDN); mpfr_mul(r23415, r23394, r23414, MPFR_RNDN); mpfr_exp(r23416, r23415, MPFR_RNDN); mpfr_mul(r23417, r23393, r23416, MPFR_RNDN); mpfr_add(r23418, r23392, r23417, MPFR_RNDN); mpfr_div(r23419, r23392, r23418, MPFR_RNDN); return mpfr_get_d(r23419, MPFR_RNDN); } static mpfr_t r23420, r23421, r23422, r23423, r23424, r23425, r23426, r23427, r23428, r23429, r23430, r23431, r23432, r23433, r23434, r23435, r23436, r23437, r23438, r23439, r23440, r23441, r23442, r23443, r23444, r23445, r23446, r23447; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r23420); mpfr_init(r23421); mpfr_init_set_str(r23422, "2.0", 10, MPFR_RNDN); mpfr_init(r23423); mpfr_init(r23424); mpfr_init(r23425); mpfr_init(r23426); mpfr_init(r23427); mpfr_init(r23428); mpfr_init(r23429); mpfr_init(r23430); mpfr_init(r23431); mpfr_init(r23432); mpfr_init_set_str(r23433, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r23434, "6.0", 10, MPFR_RNDN); mpfr_init(r23435); mpfr_init(r23436); mpfr_init_set_str(r23437, "3.0", 10, MPFR_RNDN); mpfr_init(r23438); mpfr_init(r23439); mpfr_init(r23440); mpfr_init(r23441); mpfr_init(r23442); mpfr_init(r23443); mpfr_init(r23444); mpfr_init(r23445); mpfr_init(r23446); mpfr_init(r23447); } double f_dm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r23420, x, MPFR_RNDN); mpfr_set_d(r23421, y, MPFR_RNDN); ; mpfr_set_d(r23423, z, MPFR_RNDN); mpfr_set_d(r23424, t, MPFR_RNDN); mpfr_set_d(r23425, a, MPFR_RNDN); mpfr_add(r23426, r23424, r23425, MPFR_RNDN); mpfr_sqrt(r23427, r23426, MPFR_RNDN); mpfr_div(r23428, r23424, r23427, MPFR_RNDN); mpfr_div(r23429, r23423, r23428, MPFR_RNDN); mpfr_set_d(r23430, b, MPFR_RNDN); mpfr_set_d(r23431, c, MPFR_RNDN); mpfr_sub(r23432, r23430, r23431, MPFR_RNDN); ; ; mpfr_div(r23435, r23433, r23434, MPFR_RNDN); mpfr_add(r23436, r23425, r23435, MPFR_RNDN); ; mpfr_mul(r23438, r23424, r23437, MPFR_RNDN); mpfr_div(r23439, r23422, r23438, MPFR_RNDN); mpfr_sub(r23440, r23436, r23439, MPFR_RNDN); mpfr_mul(r23441, r23432, r23440, MPFR_RNDN); mpfr_sub(r23442, r23429, r23441, MPFR_RNDN); mpfr_mul(r23443, r23422, r23442, MPFR_RNDN); mpfr_exp(r23444, r23443, MPFR_RNDN); mpfr_mul(r23445, r23421, r23444, MPFR_RNDN); mpfr_add(r23446, r23420, r23445, MPFR_RNDN); mpfr_div(r23447, r23420, r23446, MPFR_RNDN); return mpfr_get_d(r23447, MPFR_RNDN); }