#include #include #include #include #include char *name = "Jmat.Real.gamma, branch z greater than 0.5"; double f_if(float z) { float r30319 = atan2(1.0, 0.0); float r30320 = 2; float r30321 = r30319 * r30320; float r30322 = sqrt(r30321); float r30323 = z; float r30324 = 1; float r30325 = r30323 - r30324; float r30326 = 7; float r30327 = r30325 + r30326; float r30328 = 0.5; float r30329 = r30327 + r30328; float r30330 = r30325 + r30328; float r30331 = pow(r30329, r30330); float r30332 = r30322 * r30331; float r30333 = -r30329; float r30334 = exp(r30333); float r30335 = r30332 * r30334; float r30336 = 0.9999999999998099; float r30337 = 676.5203681218851; float r30338 = r30325 + r30324; float r30339 = r30337 / r30338; float r30340 = r30336 + r30339; float r30341 = -1259.1392167224028; float r30342 = r30325 + r30320; float r30343 = r30341 / r30342; float r30344 = r30340 + r30343; float r30345 = 771.3234287776531; float r30346 = 3; float r30347 = r30325 + r30346; float r30348 = r30345 / r30347; float r30349 = r30344 + r30348; float r30350 = -176.6150291621406; float r30351 = 4; float r30352 = r30325 + r30351; float r30353 = r30350 / r30352; float r30354 = r30349 + r30353; float r30355 = 12.507343278686905; float r30356 = 5; float r30357 = r30325 + r30356; float r30358 = r30355 / r30357; float r30359 = r30354 + r30358; float r30360 = -0.13857109526572012; float r30361 = 6; float r30362 = r30325 + r30361; float r30363 = r30360 / r30362; float r30364 = r30359 + r30363; float r30365 = 9.984369578019572e-06; float r30366 = r30365 / r30327; float r30367 = r30364 + r30366; float r30368 = 1.5056327351493116e-07; float r30369 = 8; float r30370 = r30325 + r30369; float r30371 = r30368 / r30370; float r30372 = r30367 + r30371; float r30373 = r30335 * r30372; return r30373; } double f_id(double z) { double r30374 = atan2(1.0, 0.0); double r30375 = 2; double r30376 = r30374 * r30375; double r30377 = sqrt(r30376); double r30378 = z; double r30379 = 1; double r30380 = r30378 - r30379; double r30381 = 7; double r30382 = r30380 + r30381; double r30383 = 0.5; double r30384 = r30382 + r30383; double r30385 = r30380 + r30383; double r30386 = pow(r30384, r30385); double r30387 = r30377 * r30386; double r30388 = -r30384; double r30389 = exp(r30388); double r30390 = r30387 * r30389; double r30391 = 0.9999999999998099; double r30392 = 676.5203681218851; double r30393 = r30380 + r30379; double r30394 = r30392 / r30393; double r30395 = r30391 + r30394; double r30396 = -1259.1392167224028; double r30397 = r30380 + r30375; double r30398 = r30396 / r30397; double r30399 = r30395 + r30398; double r30400 = 771.3234287776531; double r30401 = 3; double r30402 = r30380 + r30401; double r30403 = r30400 / r30402; double r30404 = r30399 + r30403; double r30405 = -176.6150291621406; double r30406 = 4; double r30407 = r30380 + r30406; double r30408 = r30405 / r30407; double r30409 = r30404 + r30408; double r30410 = 12.507343278686905; double r30411 = 5; double r30412 = r30380 + r30411; double r30413 = r30410 / r30412; double r30414 = r30409 + r30413; double r30415 = -0.13857109526572012; double r30416 = 6; double r30417 = r30380 + r30416; double r30418 = r30415 / r30417; double r30419 = r30414 + r30418; double r30420 = 9.984369578019572e-06; double r30421 = r30420 / r30382; double r30422 = r30419 + r30421; double r30423 = 1.5056327351493116e-07; double r30424 = 8; double r30425 = r30380 + r30424; double r30426 = r30423 / r30425; double r30427 = r30422 + r30426; double r30428 = r30390 * r30427; return r30428; } double f_of(float z) { float r30429 = atan2(1.0, 0.0); float r30430 = r30429 + r30429; float r30431 = sqrt(r30430); float r30432 = 0.5; float r30433 = 7; float r30434 = r30432 + r30433; float r30435 = z; float r30436 = 1; float r30437 = r30435 - r30436; float r30438 = r30434 + r30437; float r30439 = r30432 + r30437; float r30440 = pow(r30438, r30439); float r30441 = exp(r30438); float r30442 = r30440 / r30441; float r30443 = r30431 * r30442; float r30444 = -0.13857109526572012; float r30445 = 6; float r30446 = r30436 - r30445; float r30447 = r30435 - r30446; float r30448 = r30444 / r30447; float r30449 = 9.984369578019572e-06; float r30450 = r30437 + r30433; float r30451 = r30449 / r30450; float r30452 = r30448 + r30451; float r30453 = 12.507343278686905; float r30454 = 5; float r30455 = r30435 + r30454; float r30456 = r30455 - r30436; float r30457 = r30453 / r30456; float r30458 = r30452 + r30457; float r30459 = -176.6150291621406; float r30460 = 4; float r30461 = r30435 + r30460; float r30462 = r30461 - r30436; float r30463 = r30459 / r30462; float r30464 = 676.5203681218851; float r30465 = 0; float r30466 = r30435 - r30465; float r30467 = r30464 / r30466; float r30468 = 0.9999999999998099; float r30469 = r30467 + r30468; float r30470 = r30463 + r30469; float r30471 = -1259.1392167224028; float r30472 = 2; float r30473 = r30436 - r30472; float r30474 = r30435 - r30473; float r30475 = r30471 / r30474; float r30476 = 771.3234287776531; float r30477 = 3; float r30478 = r30477 + r30435; float r30479 = r30478 - r30436; float r30480 = r30476 / r30479; float r30481 = r30475 + r30480; float r30482 = r30470 + r30481; float r30483 = r30458 + r30482; float r30484 = 1.5056327351493116e-07; float r30485 = 8; float r30486 = r30435 + r30485; float r30487 = r30486 - r30436; float r30488 = r30484 / r30487; float r30489 = r30483 + r30488; float r30490 = r30443 * r30489; return r30490; } double f_od(double z) { double r30491 = atan2(1.0, 0.0); double r30492 = r30491 + r30491; double r30493 = sqrt(r30492); double r30494 = 0.5; double r30495 = 7; double r30496 = r30494 + r30495; double r30497 = z; double r30498 = 1; double r30499 = r30497 - r30498; double r30500 = r30496 + r30499; double r30501 = r30494 + r30499; double r30502 = pow(r30500, r30501); double r30503 = exp(r30500); double r30504 = r30502 / r30503; double r30505 = r30493 * r30504; double r30506 = -0.13857109526572012; double r30507 = 6; double r30508 = r30498 - r30507; double r30509 = r30497 - r30508; double r30510 = r30506 / r30509; double r30511 = 9.984369578019572e-06; double r30512 = r30499 + r30495; double r30513 = r30511 / r30512; double r30514 = r30510 + r30513; double r30515 = 12.507343278686905; double r30516 = 5; double r30517 = r30497 + r30516; double r30518 = r30517 - r30498; double r30519 = r30515 / r30518; double r30520 = r30514 + r30519; double r30521 = -176.6150291621406; double r30522 = 4; double r30523 = r30497 + r30522; double r30524 = r30523 - r30498; double r30525 = r30521 / r30524; double r30526 = 676.5203681218851; double r30527 = 0; double r30528 = r30497 - r30527; double r30529 = r30526 / r30528; double r30530 = 0.9999999999998099; double r30531 = r30529 + r30530; double r30532 = r30525 + r30531; double r30533 = -1259.1392167224028; double r30534 = 2; double r30535 = r30498 - r30534; double r30536 = r30497 - r30535; double r30537 = r30533 / r30536; double r30538 = 771.3234287776531; double r30539 = 3; double r30540 = r30539 + r30497; double r30541 = r30540 - r30498; double r30542 = r30538 / r30541; double r30543 = r30537 + r30542; double r30544 = r30532 + r30543; double r30545 = r30520 + r30544; double r30546 = 1.5056327351493116e-07; double r30547 = 8; double r30548 = r30497 + r30547; double r30549 = r30548 - r30498; double r30550 = r30546 / r30549; double r30551 = r30545 + r30550; double r30552 = r30505 * r30551; return r30552; } 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 r30553, r30554, r30555, r30556, r30557, r30558, r30559, r30560, r30561, r30562, r30563, r30564, r30565, r30566, r30567, r30568, r30569, r30570, r30571, r30572, r30573, r30574, r30575, r30576, r30577, r30578, r30579, r30580, r30581, r30582, r30583, r30584, r30585, r30586, r30587, r30588, r30589, r30590, r30591, r30592, r30593, r30594, r30595, r30596, r30597, r30598, r30599, r30600, r30601, r30602, r30603, r30604, r30605, r30606, r30607; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r30553); mpfr_init_set_str(r30554, "2", 10, MPFR_RNDN); mpfr_init(r30555); mpfr_init(r30556); mpfr_init(r30557); mpfr_init_set_str(r30558, "1", 10, MPFR_RNDN); mpfr_init(r30559); mpfr_init_set_str(r30560, "7", 10, MPFR_RNDN); mpfr_init(r30561); mpfr_init_set_str(r30562, "0.5", 10, MPFR_RNDN); mpfr_init(r30563); mpfr_init(r30564); mpfr_init(r30565); mpfr_init(r30566); mpfr_init(r30567); mpfr_init(r30568); mpfr_init(r30569); mpfr_init_set_str(r30570, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init_set_str(r30571, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r30572); mpfr_init(r30573); mpfr_init(r30574); mpfr_init_set_str(r30575, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init(r30576); mpfr_init(r30577); mpfr_init(r30578); mpfr_init_set_str(r30579, "771.3234287776531", 10, MPFR_RNDN); mpfr_init_set_str(r30580, "3", 10, MPFR_RNDN); mpfr_init(r30581); mpfr_init(r30582); mpfr_init(r30583); mpfr_init_set_str(r30584, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r30585, "4", 10, MPFR_RNDN); mpfr_init(r30586); mpfr_init(r30587); mpfr_init(r30588); mpfr_init_set_str(r30589, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r30590, "5", 10, MPFR_RNDN); mpfr_init(r30591); mpfr_init(r30592); mpfr_init(r30593); mpfr_init_set_str(r30594, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r30595, "6", 10, MPFR_RNDN); mpfr_init(r30596); mpfr_init(r30597); mpfr_init(r30598); mpfr_init_set_str(r30599, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init(r30600); mpfr_init(r30601); mpfr_init_set_str(r30602, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r30603, "8", 10, MPFR_RNDN); mpfr_init(r30604); mpfr_init(r30605); mpfr_init(r30606); mpfr_init(r30607); } double f_im(double z) { mpfr_const_pi(r30553, MPFR_RNDN); ; mpfr_mul(r30555, r30553, r30554, MPFR_RNDN); mpfr_sqrt(r30556, r30555, MPFR_RNDN); mpfr_set_d(r30557, z, MPFR_RNDN); ; mpfr_sub(r30559, r30557, r30558, MPFR_RNDN); ; mpfr_add(r30561, r30559, r30560, MPFR_RNDN); ; mpfr_add(r30563, r30561, r30562, MPFR_RNDN); mpfr_add(r30564, r30559, r30562, MPFR_RNDN); mpfr_pow(r30565, r30563, r30564, MPFR_RNDN); mpfr_mul(r30566, r30556, r30565, MPFR_RNDN); mpfr_neg(r30567, r30563, MPFR_RNDN); mpfr_exp(r30568, r30567, MPFR_RNDN); mpfr_mul(r30569, r30566, r30568, MPFR_RNDN); ; ; mpfr_add(r30572, r30559, r30558, MPFR_RNDN); mpfr_div(r30573, r30571, r30572, MPFR_RNDN); mpfr_add(r30574, r30570, r30573, MPFR_RNDN); ; mpfr_add(r30576, r30559, r30554, MPFR_RNDN); mpfr_div(r30577, r30575, r30576, MPFR_RNDN); mpfr_add(r30578, r30574, r30577, MPFR_RNDN); ; ; mpfr_add(r30581, r30559, r30580, MPFR_RNDN); mpfr_div(r30582, r30579, r30581, MPFR_RNDN); mpfr_add(r30583, r30578, r30582, MPFR_RNDN); ; ; mpfr_add(r30586, r30559, r30585, MPFR_RNDN); mpfr_div(r30587, r30584, r30586, MPFR_RNDN); mpfr_add(r30588, r30583, r30587, MPFR_RNDN); ; ; mpfr_add(r30591, r30559, r30590, MPFR_RNDN); mpfr_div(r30592, r30589, r30591, MPFR_RNDN); mpfr_add(r30593, r30588, r30592, MPFR_RNDN); ; ; mpfr_add(r30596, r30559, r30595, MPFR_RNDN); mpfr_div(r30597, r30594, r30596, MPFR_RNDN); mpfr_add(r30598, r30593, r30597, MPFR_RNDN); ; mpfr_div(r30600, r30599, r30561, MPFR_RNDN); mpfr_add(r30601, r30598, r30600, MPFR_RNDN); ; ; mpfr_add(r30604, r30559, r30603, MPFR_RNDN); mpfr_div(r30605, r30602, r30604, MPFR_RNDN); mpfr_add(r30606, r30601, r30605, MPFR_RNDN); mpfr_mul(r30607, r30569, r30606, MPFR_RNDN); return mpfr_get_d(r30607, MPFR_RNDN); } static mpfr_t r30608, r30609, r30610, r30611, r30612, r30613, r30614, r30615, r30616, r30617, r30618, r30619, r30620, r30621, r30622, r30623, r30624, r30625, r30626, r30627, r30628, r30629, r30630, r30631, r30632, r30633, r30634, r30635, r30636, r30637, r30638, r30639, r30640, r30641, r30642, r30643, r30644, r30645, r30646, r30647, r30648, r30649, r30650, r30651, r30652, r30653, r30654, r30655, r30656, r30657, r30658, r30659, r30660, r30661, r30662, r30663, r30664, r30665, r30666, r30667, r30668, r30669; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r30608); mpfr_init(r30609); mpfr_init(r30610); mpfr_init_set_str(r30611, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r30612, "7", 10, MPFR_RNDN); mpfr_init(r30613); mpfr_init(r30614); mpfr_init_set_str(r30615, "1", 10, MPFR_RNDN); mpfr_init(r30616); mpfr_init(r30617); mpfr_init(r30618); mpfr_init(r30619); mpfr_init(r30620); mpfr_init(r30621); mpfr_init(r30622); mpfr_init_set_str(r30623, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r30624, "6", 10, MPFR_RNDN); mpfr_init(r30625); mpfr_init(r30626); mpfr_init(r30627); mpfr_init_set_str(r30628, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init(r30629); mpfr_init(r30630); mpfr_init(r30631); mpfr_init_set_str(r30632, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r30633, "5", 10, MPFR_RNDN); mpfr_init(r30634); mpfr_init(r30635); mpfr_init(r30636); mpfr_init(r30637); mpfr_init_set_str(r30638, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r30639, "4", 10, MPFR_RNDN); mpfr_init(r30640); mpfr_init(r30641); mpfr_init(r30642); mpfr_init_set_str(r30643, "676.5203681218851", 10, MPFR_RNDN); mpfr_init_set_str(r30644, "0", 10, MPFR_RNDN); mpfr_init(r30645); mpfr_init(r30646); mpfr_init_set_str(r30647, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init(r30648); mpfr_init(r30649); mpfr_init_set_str(r30650, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init_set_str(r30651, "2", 10, MPFR_RNDN); mpfr_init(r30652); mpfr_init(r30653); mpfr_init(r30654); mpfr_init_set_str(r30655, "771.3234287776531", 10, MPFR_RNDN); mpfr_init_set_str(r30656, "3", 10, MPFR_RNDN); mpfr_init(r30657); mpfr_init(r30658); mpfr_init(r30659); mpfr_init(r30660); mpfr_init(r30661); mpfr_init(r30662); mpfr_init_set_str(r30663, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r30664, "8", 10, MPFR_RNDN); mpfr_init(r30665); mpfr_init(r30666); mpfr_init(r30667); mpfr_init(r30668); mpfr_init(r30669); } double f_fm(double z) { mpfr_const_pi(r30608, MPFR_RNDN); mpfr_add(r30609, r30608, r30608, MPFR_RNDN); mpfr_sqrt(r30610, r30609, MPFR_RNDN); ; ; mpfr_add(r30613, r30611, r30612, MPFR_RNDN); mpfr_set_d(r30614, z, MPFR_RNDN); ; mpfr_sub(r30616, r30614, r30615, MPFR_RNDN); mpfr_add(r30617, r30613, r30616, MPFR_RNDN); mpfr_add(r30618, r30611, r30616, MPFR_RNDN); mpfr_pow(r30619, r30617, r30618, MPFR_RNDN); mpfr_exp(r30620, r30617, MPFR_RNDN); mpfr_div(r30621, r30619, r30620, MPFR_RNDN); mpfr_mul(r30622, r30610, r30621, MPFR_RNDN); ; ; mpfr_sub(r30625, r30615, r30624, MPFR_RNDN); mpfr_sub(r30626, r30614, r30625, MPFR_RNDN); mpfr_div(r30627, r30623, r30626, MPFR_RNDN); ; mpfr_add(r30629, r30616, r30612, MPFR_RNDN); mpfr_div(r30630, r30628, r30629, MPFR_RNDN); mpfr_add(r30631, r30627, r30630, MPFR_RNDN); ; ; mpfr_add(r30634, r30614, r30633, MPFR_RNDN); mpfr_sub(r30635, r30634, r30615, MPFR_RNDN); mpfr_div(r30636, r30632, r30635, MPFR_RNDN); mpfr_add(r30637, r30631, r30636, MPFR_RNDN); ; ; mpfr_add(r30640, r30614, r30639, MPFR_RNDN); mpfr_sub(r30641, r30640, r30615, MPFR_RNDN); mpfr_div(r30642, r30638, r30641, MPFR_RNDN); ; ; mpfr_sub(r30645, r30614, r30644, MPFR_RNDN); mpfr_div(r30646, r30643, r30645, MPFR_RNDN); ; mpfr_add(r30648, r30646, r30647, MPFR_RNDN); mpfr_add(r30649, r30642, r30648, MPFR_RNDN); ; ; mpfr_sub(r30652, r30615, r30651, MPFR_RNDN); mpfr_sub(r30653, r30614, r30652, MPFR_RNDN); mpfr_div(r30654, r30650, r30653, MPFR_RNDN); ; ; mpfr_add(r30657, r30656, r30614, MPFR_RNDN); mpfr_sub(r30658, r30657, r30615, MPFR_RNDN); mpfr_div(r30659, r30655, r30658, MPFR_RNDN); mpfr_add(r30660, r30654, r30659, MPFR_RNDN); mpfr_add(r30661, r30649, r30660, MPFR_RNDN); mpfr_add(r30662, r30637, r30661, MPFR_RNDN); ; ; mpfr_add(r30665, r30614, r30664, MPFR_RNDN); mpfr_sub(r30666, r30665, r30615, MPFR_RNDN); mpfr_div(r30667, r30663, r30666, MPFR_RNDN); mpfr_add(r30668, r30662, r30667, MPFR_RNDN); mpfr_mul(r30669, r30622, r30668, MPFR_RNDN); return mpfr_get_d(r30669, MPFR_RNDN); } static mpfr_t r30670, r30671, r30672, r30673, r30674, r30675, r30676, r30677, r30678, r30679, r30680, r30681, r30682, r30683, r30684, r30685, r30686, r30687, r30688, r30689, r30690, r30691, r30692, r30693, r30694, r30695, r30696, r30697, r30698, r30699, r30700, r30701, r30702, r30703, r30704, r30705, r30706, r30707, r30708, r30709, r30710, r30711, r30712, r30713, r30714, r30715, r30716, r30717, r30718, r30719, r30720, r30721, r30722, r30723, r30724, r30725, r30726, r30727, r30728, r30729, r30730, r30731; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r30670); mpfr_init(r30671); mpfr_init(r30672); mpfr_init_set_str(r30673, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r30674, "7", 10, MPFR_RNDN); mpfr_init(r30675); mpfr_init(r30676); mpfr_init_set_str(r30677, "1", 10, MPFR_RNDN); mpfr_init(r30678); mpfr_init(r30679); mpfr_init(r30680); mpfr_init(r30681); mpfr_init(r30682); mpfr_init(r30683); mpfr_init(r30684); mpfr_init_set_str(r30685, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r30686, "6", 10, MPFR_RNDN); mpfr_init(r30687); mpfr_init(r30688); mpfr_init(r30689); mpfr_init_set_str(r30690, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init(r30691); mpfr_init(r30692); mpfr_init(r30693); mpfr_init_set_str(r30694, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r30695, "5", 10, MPFR_RNDN); mpfr_init(r30696); mpfr_init(r30697); mpfr_init(r30698); mpfr_init(r30699); mpfr_init_set_str(r30700, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r30701, "4", 10, MPFR_RNDN); mpfr_init(r30702); mpfr_init(r30703); mpfr_init(r30704); mpfr_init_set_str(r30705, "676.5203681218851", 10, MPFR_RNDN); mpfr_init_set_str(r30706, "0", 10, MPFR_RNDN); mpfr_init(r30707); mpfr_init(r30708); mpfr_init_set_str(r30709, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init(r30710); mpfr_init(r30711); mpfr_init_set_str(r30712, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init_set_str(r30713, "2", 10, MPFR_RNDN); mpfr_init(r30714); mpfr_init(r30715); mpfr_init(r30716); mpfr_init_set_str(r30717, "771.3234287776531", 10, MPFR_RNDN); mpfr_init_set_str(r30718, "3", 10, MPFR_RNDN); mpfr_init(r30719); mpfr_init(r30720); mpfr_init(r30721); mpfr_init(r30722); mpfr_init(r30723); mpfr_init(r30724); mpfr_init_set_str(r30725, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r30726, "8", 10, MPFR_RNDN); mpfr_init(r30727); mpfr_init(r30728); mpfr_init(r30729); mpfr_init(r30730); mpfr_init(r30731); } double f_dm(double z) { mpfr_const_pi(r30670, MPFR_RNDN); mpfr_add(r30671, r30670, r30670, MPFR_RNDN); mpfr_sqrt(r30672, r30671, MPFR_RNDN); ; ; mpfr_add(r30675, r30673, r30674, MPFR_RNDN); mpfr_set_d(r30676, z, MPFR_RNDN); ; mpfr_sub(r30678, r30676, r30677, MPFR_RNDN); mpfr_add(r30679, r30675, r30678, MPFR_RNDN); mpfr_add(r30680, r30673, r30678, MPFR_RNDN); mpfr_pow(r30681, r30679, r30680, MPFR_RNDN); mpfr_exp(r30682, r30679, MPFR_RNDN); mpfr_div(r30683, r30681, r30682, MPFR_RNDN); mpfr_mul(r30684, r30672, r30683, MPFR_RNDN); ; ; mpfr_sub(r30687, r30677, r30686, MPFR_RNDN); mpfr_sub(r30688, r30676, r30687, MPFR_RNDN); mpfr_div(r30689, r30685, r30688, MPFR_RNDN); ; mpfr_add(r30691, r30678, r30674, MPFR_RNDN); mpfr_div(r30692, r30690, r30691, MPFR_RNDN); mpfr_add(r30693, r30689, r30692, MPFR_RNDN); ; ; mpfr_add(r30696, r30676, r30695, MPFR_RNDN); mpfr_sub(r30697, r30696, r30677, MPFR_RNDN); mpfr_div(r30698, r30694, r30697, MPFR_RNDN); mpfr_add(r30699, r30693, r30698, MPFR_RNDN); ; ; mpfr_add(r30702, r30676, r30701, MPFR_RNDN); mpfr_sub(r30703, r30702, r30677, MPFR_RNDN); mpfr_div(r30704, r30700, r30703, MPFR_RNDN); ; ; mpfr_sub(r30707, r30676, r30706, MPFR_RNDN); mpfr_div(r30708, r30705, r30707, MPFR_RNDN); ; mpfr_add(r30710, r30708, r30709, MPFR_RNDN); mpfr_add(r30711, r30704, r30710, MPFR_RNDN); ; ; mpfr_sub(r30714, r30677, r30713, MPFR_RNDN); mpfr_sub(r30715, r30676, r30714, MPFR_RNDN); mpfr_div(r30716, r30712, r30715, MPFR_RNDN); ; ; mpfr_add(r30719, r30718, r30676, MPFR_RNDN); mpfr_sub(r30720, r30719, r30677, MPFR_RNDN); mpfr_div(r30721, r30717, r30720, MPFR_RNDN); mpfr_add(r30722, r30716, r30721, MPFR_RNDN); mpfr_add(r30723, r30711, r30722, MPFR_RNDN); mpfr_add(r30724, r30699, r30723, MPFR_RNDN); ; ; mpfr_add(r30727, r30676, r30726, MPFR_RNDN); mpfr_sub(r30728, r30727, r30677, MPFR_RNDN); mpfr_div(r30729, r30725, r30728, MPFR_RNDN); mpfr_add(r30730, r30724, r30729, MPFR_RNDN); mpfr_mul(r30731, r30684, r30730, MPFR_RNDN); return mpfr_get_d(r30731, MPFR_RNDN); }