#include #include #include #include #include char *name = "Jmat.Real.gamma, branch z greater than 0.5"; double f_if(float z) { float r17231 = atan2(1.0, 0.0); float r17232 = 2.0f; float r17233 = r17231 * r17232; float r17234 = sqrt(r17233); float r17235 = z; float r17236 = 1.0f; float r17237 = r17235 - r17236; float r17238 = 7.0f; float r17239 = r17237 + r17238; float r17240 = 0.5f; float r17241 = r17239 + r17240; float r17242 = r17237 + r17240; float r17243 = pow(r17241, r17242); float r17244 = r17234 * r17243; float r17245 = -r17241; float r17246 = exp(r17245); float r17247 = r17244 * r17246; float r17248 = 1.0f; float r17249 = 676.5203857421875f; float r17250 = r17237 + r17236; float r17251 = r17249 / r17250; float r17252 = r17248 + r17251; float r17253 = -1259.13916015625f; float r17254 = r17237 + r17232; float r17255 = r17253 / r17254; float r17256 = r17252 + r17255; float r17257 = 771.3234252929688f; float r17258 = 3.0f; float r17259 = r17237 + r17258; float r17260 = r17257 / r17259; float r17261 = r17256 + r17260; float r17262 = -176.6150360107422f; float r17263 = 4.0f; float r17264 = r17237 + r17263; float r17265 = r17262 / r17264; float r17266 = r17261 + r17265; float r17267 = 12.507343292236328f; float r17268 = 5.0f; float r17269 = r17237 + r17268; float r17270 = r17267 / r17269; float r17271 = r17266 + r17270; float r17272 = -0.138571098446846f; float r17273 = 6.0f; float r17274 = r17237 + r17273; float r17275 = r17272 / r17274; float r17276 = r17271 + r17275; float r17277 = 9.984369171434082e-06f; float r17278 = r17277 / r17239; float r17279 = r17276 + r17278; float r17280 = 1.5056326674312004e-07f; float r17281 = 8.0f; float r17282 = r17237 + r17281; float r17283 = r17280 / r17282; float r17284 = r17279 + r17283; float r17285 = r17247 * r17284; return r17285; } double f_id(double z) { double r17286 = atan2(1.0, 0.0); double r17287 = 2.0; double r17288 = r17286 * r17287; double r17289 = sqrt(r17288); double r17290 = z; double r17291 = 1.0; double r17292 = r17290 - r17291; double r17293 = 7.0; double r17294 = r17292 + r17293; double r17295 = 0.5; double r17296 = r17294 + r17295; double r17297 = r17292 + r17295; double r17298 = pow(r17296, r17297); double r17299 = r17289 * r17298; double r17300 = -r17296; double r17301 = exp(r17300); double r17302 = r17299 * r17301; double r17303 = 1.0; double r17304 = 676.5203857421875; double r17305 = r17292 + r17291; double r17306 = r17304 / r17305; double r17307 = r17303 + r17306; double r17308 = -1259.13916015625; double r17309 = r17292 + r17287; double r17310 = r17308 / r17309; double r17311 = r17307 + r17310; double r17312 = 771.3234252929688; double r17313 = 3.0; double r17314 = r17292 + r17313; double r17315 = r17312 / r17314; double r17316 = r17311 + r17315; double r17317 = -176.6150360107422; double r17318 = 4.0; double r17319 = r17292 + r17318; double r17320 = r17317 / r17319; double r17321 = r17316 + r17320; double r17322 = 12.507343292236328; double r17323 = 5.0; double r17324 = r17292 + r17323; double r17325 = r17322 / r17324; double r17326 = r17321 + r17325; double r17327 = -0.138571098446846; double r17328 = 6.0; double r17329 = r17292 + r17328; double r17330 = r17327 / r17329; double r17331 = r17326 + r17330; double r17332 = 9.984369171434082e-06; double r17333 = r17332 / r17294; double r17334 = r17331 + r17333; double r17335 = 1.5056326674312004e-07; double r17336 = 8.0; double r17337 = r17292 + r17336; double r17338 = r17335 / r17337; double r17339 = r17334 + r17338; double r17340 = r17302 * r17339; return r17340; } double f_of(float z) { float r17341 = 2585.19482421875f; float r17342 = z; float r17343 = atan2(1.0, 0.0); float r17344 = sqrt(r17343); float r17345 = r17342 * r17344; float r17346 = 6.5f; float r17347 = exp(r17346); float r17348 = 2.0f; float r17349 = sqrt(r17348); float r17350 = r17347 / r17349; float r17351 = r17345 / r17350; float r17352 = 1.0f; float r17353 = 1.0f; float r17354 = pow(r17346, r17353); float r17355 = r17352 / r17354; float r17356 = 0.5f; float r17357 = pow(r17355, r17356); float r17358 = r17351 * r17357; float r17359 = r17349 * r17342; float r17360 = log(r17346); float r17361 = r17360 * r17360; float r17362 = r17359 * r17361; float r17363 = r17362 / r17347; float r17364 = r17344 * r17357; float r17365 = r17363 * r17364; float r17366 = 338.26019287109375f; float r17367 = r17365 * r17366; float r17368 = fma(r17341, r17358, r17367); float r17369 = 676.5203857421875f; float r17370 = r17349 / r17347; float r17371 = r17360 * r17370; float r17372 = r17370 / r17342; float r17373 = r17371 + r17372; float r17374 = r17364 * r17373; float r17375 = r17369 * r17374; float r17376 = r17360 * r17342; float r17377 = r17376 * r17349; float r17378 = r17377 / r17347; float r17379 = r17378 * r17364; float r17380 = 1656.8104248046875f; float r17381 = r17370 * r17380; float r17382 = r17364 * r17381; float r17383 = fma(r17379, r17380, r17382); float r17384 = r17375 - r17383; float r17385 = r17368 + r17384; return r17385; } double f_od(double z) { double r17386 = 2585.19482421875; double r17387 = z; double r17388 = atan2(1.0, 0.0); double r17389 = sqrt(r17388); double r17390 = r17387 * r17389; double r17391 = 6.5; double r17392 = exp(r17391); double r17393 = 2.0; double r17394 = sqrt(r17393); double r17395 = r17392 / r17394; double r17396 = r17390 / r17395; double r17397 = 1.0; double r17398 = 1.0; double r17399 = pow(r17391, r17398); double r17400 = r17397 / r17399; double r17401 = 0.5; double r17402 = pow(r17400, r17401); double r17403 = r17396 * r17402; double r17404 = r17394 * r17387; double r17405 = log(r17391); double r17406 = r17405 * r17405; double r17407 = r17404 * r17406; double r17408 = r17407 / r17392; double r17409 = r17389 * r17402; double r17410 = r17408 * r17409; double r17411 = 338.26019287109375; double r17412 = r17410 * r17411; double r17413 = fma(r17386, r17403, r17412); double r17414 = 676.5203857421875; double r17415 = r17394 / r17392; double r17416 = r17405 * r17415; double r17417 = r17415 / r17387; double r17418 = r17416 + r17417; double r17419 = r17409 * r17418; double r17420 = r17414 * r17419; double r17421 = r17405 * r17387; double r17422 = r17421 * r17394; double r17423 = r17422 / r17392; double r17424 = r17423 * r17409; double r17425 = 1656.8104248046875; double r17426 = r17415 * r17425; double r17427 = r17409 * r17426; double r17428 = fma(r17424, r17425, r17427); double r17429 = r17420 - r17428; double r17430 = r17413 + r17429; return r17430; } 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 r17431, r17432, r17433, r17434, r17435, r17436, r17437, r17438, r17439, r17440, r17441, r17442, r17443, r17444, r17445, r17446, r17447, r17448, r17449, r17450, r17451, r17452, r17453, r17454, r17455, r17456, r17457, r17458, r17459, r17460, r17461, r17462, r17463, r17464, r17465, r17466, r17467, r17468, r17469, r17470, r17471, r17472, r17473, r17474, r17475, r17476, r17477, r17478, r17479, r17480, r17481, r17482, r17483, r17484, r17485; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17431); mpfr_init_set_str(r17432, "2", 10, MPFR_RNDN); mpfr_init(r17433); mpfr_init(r17434); mpfr_init(r17435); mpfr_init_set_str(r17436, "1", 10, MPFR_RNDN); mpfr_init(r17437); mpfr_init_set_str(r17438, "7", 10, MPFR_RNDN); mpfr_init(r17439); mpfr_init_set_str(r17440, "0.5", 10, MPFR_RNDN); mpfr_init(r17441); mpfr_init(r17442); mpfr_init(r17443); mpfr_init(r17444); mpfr_init(r17445); mpfr_init(r17446); mpfr_init(r17447); mpfr_init_set_str(r17448, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init_set_str(r17449, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r17450); mpfr_init(r17451); mpfr_init(r17452); mpfr_init_set_str(r17453, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init(r17454); mpfr_init(r17455); mpfr_init(r17456); mpfr_init_set_str(r17457, "771.3234287776531", 10, MPFR_RNDN); mpfr_init_set_str(r17458, "3", 10, MPFR_RNDN); mpfr_init(r17459); mpfr_init(r17460); mpfr_init(r17461); mpfr_init_set_str(r17462, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r17463, "4", 10, MPFR_RNDN); mpfr_init(r17464); mpfr_init(r17465); mpfr_init(r17466); mpfr_init_set_str(r17467, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r17468, "5", 10, MPFR_RNDN); mpfr_init(r17469); mpfr_init(r17470); mpfr_init(r17471); mpfr_init_set_str(r17472, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r17473, "6", 10, MPFR_RNDN); mpfr_init(r17474); mpfr_init(r17475); mpfr_init(r17476); mpfr_init_set_str(r17477, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init(r17478); mpfr_init(r17479); mpfr_init_set_str(r17480, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r17481, "8", 10, MPFR_RNDN); mpfr_init(r17482); mpfr_init(r17483); mpfr_init(r17484); mpfr_init(r17485); } double f_im(double z) { mpfr_const_pi(r17431, MPFR_RNDN); ; mpfr_mul(r17433, r17431, r17432, MPFR_RNDN); mpfr_sqrt(r17434, r17433, MPFR_RNDN); mpfr_set_d(r17435, z, MPFR_RNDN); ; mpfr_sub(r17437, r17435, r17436, MPFR_RNDN); ; mpfr_add(r17439, r17437, r17438, MPFR_RNDN); ; mpfr_add(r17441, r17439, r17440, MPFR_RNDN); mpfr_add(r17442, r17437, r17440, MPFR_RNDN); mpfr_pow(r17443, r17441, r17442, MPFR_RNDN); mpfr_mul(r17444, r17434, r17443, MPFR_RNDN); mpfr_neg(r17445, r17441, MPFR_RNDN); mpfr_exp(r17446, r17445, MPFR_RNDN); mpfr_mul(r17447, r17444, r17446, MPFR_RNDN); ; ; mpfr_add(r17450, r17437, r17436, MPFR_RNDN); mpfr_div(r17451, r17449, r17450, MPFR_RNDN); mpfr_add(r17452, r17448, r17451, MPFR_RNDN); ; mpfr_add(r17454, r17437, r17432, MPFR_RNDN); mpfr_div(r17455, r17453, r17454, MPFR_RNDN); mpfr_add(r17456, r17452, r17455, MPFR_RNDN); ; ; mpfr_add(r17459, r17437, r17458, MPFR_RNDN); mpfr_div(r17460, r17457, r17459, MPFR_RNDN); mpfr_add(r17461, r17456, r17460, MPFR_RNDN); ; ; mpfr_add(r17464, r17437, r17463, MPFR_RNDN); mpfr_div(r17465, r17462, r17464, MPFR_RNDN); mpfr_add(r17466, r17461, r17465, MPFR_RNDN); ; ; mpfr_add(r17469, r17437, r17468, MPFR_RNDN); mpfr_div(r17470, r17467, r17469, MPFR_RNDN); mpfr_add(r17471, r17466, r17470, MPFR_RNDN); ; ; mpfr_add(r17474, r17437, r17473, MPFR_RNDN); mpfr_div(r17475, r17472, r17474, MPFR_RNDN); mpfr_add(r17476, r17471, r17475, MPFR_RNDN); ; mpfr_div(r17478, r17477, r17439, MPFR_RNDN); mpfr_add(r17479, r17476, r17478, MPFR_RNDN); ; ; mpfr_add(r17482, r17437, r17481, MPFR_RNDN); mpfr_div(r17483, r17480, r17482, MPFR_RNDN); mpfr_add(r17484, r17479, r17483, MPFR_RNDN); mpfr_mul(r17485, r17447, r17484, MPFR_RNDN); return mpfr_get_d(r17485, MPFR_RNDN); } static mpfr_t r17486, r17487, r17488, r17489, r17490, r17491, r17492, r17493, r17494, r17495, r17496, r17497, r17498, r17499, r17500, r17501, r17502, r17503, r17504, r17505, r17506, r17507, r17508, r17509, r17510, r17511, r17512, r17513, r17514, r17515, r17516, r17517, r17518, r17519, r17520, r17521, r17522, r17523, r17524, r17525, r17526, r17527, r17528, r17529, r17530; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17486, "2585.1948787825354", 10, MPFR_RNDN); mpfr_init(r17487); mpfr_init(r17488); mpfr_init(r17489); mpfr_init(r17490); mpfr_init_set_str(r17491, "6.5", 10, MPFR_RNDN); mpfr_init(r17492); mpfr_init_set_str(r17493, "2", 10, MPFR_RNDN); mpfr_init(r17494); mpfr_init(r17495); mpfr_init(r17496); mpfr_init_set_str(r17497, "1", 10, MPFR_RNDN); mpfr_init_set_str(r17498, "1.0", 10, MPFR_RNDN); mpfr_init(r17499); mpfr_init(r17500); mpfr_init_set_str(r17501, "0.5", 10, MPFR_RNDN); mpfr_init(r17502); mpfr_init(r17503); mpfr_init(r17504); mpfr_init(r17505); mpfr_init(r17506); mpfr_init(r17507); mpfr_init(r17508); mpfr_init(r17509); mpfr_init(r17510); mpfr_init_set_str(r17511, "338.26018406094255", 10, MPFR_RNDN); mpfr_init(r17512); mpfr_init(r17513); mpfr_init_set_str(r17514, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r17515); mpfr_init(r17516); mpfr_init(r17517); mpfr_init(r17518); mpfr_init(r17519); mpfr_init(r17520); mpfr_init(r17521); mpfr_init(r17522); mpfr_init(r17523); mpfr_init(r17524); mpfr_init_set_str(r17525, "1656.8104518737205", 10, MPFR_RNDN); mpfr_init(r17526); mpfr_init(r17527); mpfr_init(r17528); mpfr_init(r17529); mpfr_init(r17530); } double f_fm(double z) { ; mpfr_set_d(r17487, z, MPFR_RNDN); mpfr_const_pi(r17488, MPFR_RNDN); mpfr_sqrt(r17489, r17488, MPFR_RNDN); mpfr_mul(r17490, r17487, r17489, MPFR_RNDN); ; mpfr_exp(r17492, r17491, MPFR_RNDN); ; mpfr_sqrt(r17494, r17493, MPFR_RNDN); mpfr_div(r17495, r17492, r17494, MPFR_RNDN); mpfr_div(r17496, r17490, r17495, MPFR_RNDN); ; ; mpfr_pow(r17499, r17491, r17498, MPFR_RNDN); mpfr_div(r17500, r17497, r17499, MPFR_RNDN); ; mpfr_pow(r17502, r17500, r17501, MPFR_RNDN); mpfr_mul(r17503, r17496, r17502, MPFR_RNDN); mpfr_mul(r17504, r17494, r17487, MPFR_RNDN); mpfr_log(r17505, r17491, MPFR_RNDN); mpfr_sqr(r17506, r17505, MPFR_RNDN); mpfr_mul(r17507, r17504, r17506, MPFR_RNDN); mpfr_div(r17508, r17507, r17492, MPFR_RNDN); mpfr_mul(r17509, r17489, r17502, MPFR_RNDN); mpfr_mul(r17510, r17508, r17509, MPFR_RNDN); ; mpfr_mul(r17512, r17510, r17511, MPFR_RNDN); mpfr_fma(r17513, r17486, r17503, r17512, MPFR_RNDN); ; mpfr_div(r17515, r17494, r17492, MPFR_RNDN); mpfr_mul(r17516, r17505, r17515, MPFR_RNDN); mpfr_div(r17517, r17515, r17487, MPFR_RNDN); mpfr_add(r17518, r17516, r17517, MPFR_RNDN); mpfr_mul(r17519, r17509, r17518, MPFR_RNDN); mpfr_mul(r17520, r17514, r17519, MPFR_RNDN); mpfr_mul(r17521, r17505, r17487, MPFR_RNDN); mpfr_mul(r17522, r17521, r17494, MPFR_RNDN); mpfr_div(r17523, r17522, r17492, MPFR_RNDN); mpfr_mul(r17524, r17523, r17509, MPFR_RNDN); ; mpfr_mul(r17526, r17515, r17525, MPFR_RNDN); mpfr_mul(r17527, r17509, r17526, MPFR_RNDN); mpfr_fma(r17528, r17524, r17525, r17527, MPFR_RNDN); mpfr_sub(r17529, r17520, r17528, MPFR_RNDN); mpfr_add(r17530, r17513, r17529, MPFR_RNDN); return mpfr_get_d(r17530, MPFR_RNDN); } static mpfr_t r17531, r17532, r17533, r17534, r17535, r17536, r17537, r17538, r17539, r17540, r17541, r17542, r17543, r17544, r17545, r17546, r17547, r17548, r17549, r17550, r17551, r17552, r17553, r17554, r17555, r17556, r17557, r17558, r17559, r17560, r17561, r17562, r17563, r17564, r17565, r17566, r17567, r17568, r17569, r17570, r17571, r17572, r17573, r17574, r17575; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17531, "2585.1948787825354", 10, MPFR_RNDN); mpfr_init(r17532); mpfr_init(r17533); mpfr_init(r17534); mpfr_init(r17535); mpfr_init_set_str(r17536, "6.5", 10, MPFR_RNDN); mpfr_init(r17537); mpfr_init_set_str(r17538, "2", 10, MPFR_RNDN); mpfr_init(r17539); mpfr_init(r17540); mpfr_init(r17541); mpfr_init_set_str(r17542, "1", 10, MPFR_RNDN); mpfr_init_set_str(r17543, "1.0", 10, MPFR_RNDN); mpfr_init(r17544); mpfr_init(r17545); mpfr_init_set_str(r17546, "0.5", 10, MPFR_RNDN); mpfr_init(r17547); mpfr_init(r17548); mpfr_init(r17549); mpfr_init(r17550); mpfr_init(r17551); mpfr_init(r17552); mpfr_init(r17553); mpfr_init(r17554); mpfr_init(r17555); mpfr_init_set_str(r17556, "338.26018406094255", 10, MPFR_RNDN); mpfr_init(r17557); mpfr_init(r17558); mpfr_init_set_str(r17559, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r17560); mpfr_init(r17561); mpfr_init(r17562); mpfr_init(r17563); mpfr_init(r17564); mpfr_init(r17565); mpfr_init(r17566); mpfr_init(r17567); mpfr_init(r17568); mpfr_init(r17569); mpfr_init_set_str(r17570, "1656.8104518737205", 10, MPFR_RNDN); mpfr_init(r17571); mpfr_init(r17572); mpfr_init(r17573); mpfr_init(r17574); mpfr_init(r17575); } double f_dm(double z) { ; mpfr_set_d(r17532, z, MPFR_RNDN); mpfr_const_pi(r17533, MPFR_RNDN); mpfr_sqrt(r17534, r17533, MPFR_RNDN); mpfr_mul(r17535, r17532, r17534, MPFR_RNDN); ; mpfr_exp(r17537, r17536, MPFR_RNDN); ; mpfr_sqrt(r17539, r17538, MPFR_RNDN); mpfr_div(r17540, r17537, r17539, MPFR_RNDN); mpfr_div(r17541, r17535, r17540, MPFR_RNDN); ; ; mpfr_pow(r17544, r17536, r17543, MPFR_RNDN); mpfr_div(r17545, r17542, r17544, MPFR_RNDN); ; mpfr_pow(r17547, r17545, r17546, MPFR_RNDN); mpfr_mul(r17548, r17541, r17547, MPFR_RNDN); mpfr_mul(r17549, r17539, r17532, MPFR_RNDN); mpfr_log(r17550, r17536, MPFR_RNDN); mpfr_sqr(r17551, r17550, MPFR_RNDN); mpfr_mul(r17552, r17549, r17551, MPFR_RNDN); mpfr_div(r17553, r17552, r17537, MPFR_RNDN); mpfr_mul(r17554, r17534, r17547, MPFR_RNDN); mpfr_mul(r17555, r17553, r17554, MPFR_RNDN); ; mpfr_mul(r17557, r17555, r17556, MPFR_RNDN); mpfr_fma(r17558, r17531, r17548, r17557, MPFR_RNDN); ; mpfr_div(r17560, r17539, r17537, MPFR_RNDN); mpfr_mul(r17561, r17550, r17560, MPFR_RNDN); mpfr_div(r17562, r17560, r17532, MPFR_RNDN); mpfr_add(r17563, r17561, r17562, MPFR_RNDN); mpfr_mul(r17564, r17554, r17563, MPFR_RNDN); mpfr_mul(r17565, r17559, r17564, MPFR_RNDN); mpfr_mul(r17566, r17550, r17532, MPFR_RNDN); mpfr_mul(r17567, r17566, r17539, MPFR_RNDN); mpfr_div(r17568, r17567, r17537, MPFR_RNDN); mpfr_mul(r17569, r17568, r17554, MPFR_RNDN); ; mpfr_mul(r17571, r17560, r17570, MPFR_RNDN); mpfr_mul(r17572, r17554, r17571, MPFR_RNDN); mpfr_fma(r17573, r17569, r17570, r17572, MPFR_RNDN); mpfr_sub(r17574, r17565, r17573, MPFR_RNDN); mpfr_add(r17575, r17558, r17574, MPFR_RNDN); return mpfr_get_d(r17575, MPFR_RNDN); }