#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r11349 = 1.0f; float r11350 = 0.1049934947f; float r11351 = x; float r11352 = r11351 * r11351; float r11353 = r11350 * r11352; float r11354 = r11349 + r11353; float r11355 = 0.0424060604f; float r11356 = r11352 * r11352; float r11357 = r11355 * r11356; float r11358 = r11354 + r11357; float r11359 = 0.0072644182f; float r11360 = r11356 * r11352; float r11361 = r11359 * r11360; float r11362 = r11358 + r11361; float r11363 = 0.0005064034f; float r11364 = r11360 * r11352; float r11365 = r11363 * r11364; float r11366 = r11362 + r11365; float r11367 = 0.0001789971f; float r11368 = r11364 * r11352; float r11369 = r11367 * r11368; float r11370 = r11366 + r11369; float r11371 = 0.7715471019f; float r11372 = r11371 * r11352; float r11373 = r11349 + r11372; float r11374 = 0.2909738639f; float r11375 = r11374 * r11356; float r11376 = r11373 + r11375; float r11377 = 0.0694555761f; float r11378 = r11377 * r11360; float r11379 = r11376 + r11378; float r11380 = 0.0140005442f; float r11381 = r11380 * r11364; float r11382 = r11379 + r11381; float r11383 = 0.0008327945f; float r11384 = r11383 * r11368; float r11385 = r11382 + r11384; float r11386 = 2.0f; float r11387 = r11386 * r11367; float r11388 = r11368 * r11352; float r11389 = r11387 * r11388; float r11390 = r11385 + r11389; float r11391 = r11370 / r11390; float r11392 = r11391 * r11351; return r11392; } double f_id(double x) { double r11393 = 1.0; double r11394 = 0.1049934947; double r11395 = x; double r11396 = r11395 * r11395; double r11397 = r11394 * r11396; double r11398 = r11393 + r11397; double r11399 = 0.0424060604; double r11400 = r11396 * r11396; double r11401 = r11399 * r11400; double r11402 = r11398 + r11401; double r11403 = 0.0072644182; double r11404 = r11400 * r11396; double r11405 = r11403 * r11404; double r11406 = r11402 + r11405; double r11407 = 0.0005064034; double r11408 = r11404 * r11396; double r11409 = r11407 * r11408; double r11410 = r11406 + r11409; double r11411 = 0.0001789971; double r11412 = r11408 * r11396; double r11413 = r11411 * r11412; double r11414 = r11410 + r11413; double r11415 = 0.7715471019; double r11416 = r11415 * r11396; double r11417 = r11393 + r11416; double r11418 = 0.2909738639; double r11419 = r11418 * r11400; double r11420 = r11417 + r11419; double r11421 = 0.0694555761; double r11422 = r11421 * r11404; double r11423 = r11420 + r11422; double r11424 = 0.0140005442; double r11425 = r11424 * r11408; double r11426 = r11423 + r11425; double r11427 = 0.0008327945; double r11428 = r11427 * r11412; double r11429 = r11426 + r11428; double r11430 = 2.0; double r11431 = r11430 * r11411; double r11432 = r11412 * r11396; double r11433 = r11431 * r11432; double r11434 = r11429 + r11433; double r11435 = r11414 / r11434; double r11436 = r11435 * r11395; return r11436; } double f_of(float x) { float r11437 = 0.2514179000665375f; float r11438 = x; float r11439 = r11437 / r11438; float r11440 = r11438 * r11438; float r11441 = r11439 / r11440; float r11442 = 0.5f; float r11443 = r11442 / r11438; float r11444 = 0.15298196345929327f; float r11445 = 5.0f; float r11446 = pow(r11438, r11445); float r11447 = r11444 / r11446; float r11448 = r11443 + r11447; float r11449 = r11441 + r11448; float r11450 = -0.00040364450155768605f; bool r11451 = r11449 <= r11450; float r11452 = 1.0f; float r11453 = 0.265709700396151f; float r11454 = 4.0f; float r11455 = pow(r11438, r11454); float r11456 = r11453 * r11455; float r11457 = r11452 + r11456; float r11458 = 2.0f; float r11459 = pow(r11438, r11458); float r11460 = 0.6665536072f; float r11461 = r11459 * r11460; float r11462 = r11457 - r11461; float r11463 = r11438 * r11462; float r11464 = 8.03349330015686e-09f; bool r11465 = r11449 <= r11464; float r11466 = r11440 * r11440; float r11467 = r11440 * r11466; float r11468 = r11440 * r11467; float r11469 = r11468 * r11440; float r11470 = 0.0001789971f; float r11471 = r11469 * r11470; float r11472 = 0.0072644182f; float r11473 = r11472 * r11467; float r11474 = 0.1049934947f; float r11475 = r11440 * r11474; float r11476 = r11475 + r11452; float r11477 = 0.0424060604f; float r11478 = r11477 * r11466; float r11479 = r11476 + r11478; float r11480 = r11473 + r11479; float r11481 = 0.0005064034f; float r11482 = r11481 * r11468; float r11483 = r11480 + r11482; float r11484 = r11471 + r11483; float r11485 = r11469 * r11440; float r11486 = r11470 * r11458; float r11487 = r11485 * r11486; float r11488 = 0.0008327945f; float r11489 = r11469 * r11488; float r11490 = 0.0694555761f; float r11491 = r11490 * r11467; float r11492 = 0.7715471019f; float r11493 = r11440 * r11492; float r11494 = r11493 + r11452; float r11495 = 0.2909738639f; float r11496 = r11495 * r11466; float r11497 = r11494 + r11496; float r11498 = r11491 + r11497; float r11499 = 0.0140005442f; float r11500 = r11468 * r11499; float r11501 = r11498 + r11500; float r11502 = r11489 + r11501; float r11503 = r11487 + r11502; float r11504 = r11484 / r11503; float r11505 = log(r11504); float r11506 = exp(r11505); float r11507 = r11438 * r11506; float r11508 = r11465 ? r11449 : r11507; float r11509 = r11451 ? r11463 : r11508; return r11509; } double f_od(double x) { double r11510 = 0.2514179000665375; double r11511 = x; double r11512 = r11510 / r11511; double r11513 = r11511 * r11511; double r11514 = r11512 / r11513; double r11515 = 0.5; double r11516 = r11515 / r11511; double r11517 = 0.15298196345929327; double r11518 = 5.0; double r11519 = pow(r11511, r11518); double r11520 = r11517 / r11519; double r11521 = r11516 + r11520; double r11522 = r11514 + r11521; double r11523 = -0.00040364450155768605; bool r11524 = r11522 <= r11523; double r11525 = 1.0; double r11526 = 0.265709700396151; double r11527 = 4.0; double r11528 = pow(r11511, r11527); double r11529 = r11526 * r11528; double r11530 = r11525 + r11529; double r11531 = 2.0; double r11532 = pow(r11511, r11531); double r11533 = 0.6665536072; double r11534 = r11532 * r11533; double r11535 = r11530 - r11534; double r11536 = r11511 * r11535; double r11537 = 8.03349330015686e-09; bool r11538 = r11522 <= r11537; double r11539 = r11513 * r11513; double r11540 = r11513 * r11539; double r11541 = r11513 * r11540; double r11542 = r11541 * r11513; double r11543 = 0.0001789971; double r11544 = r11542 * r11543; double r11545 = 0.0072644182; double r11546 = r11545 * r11540; double r11547 = 0.1049934947; double r11548 = r11513 * r11547; double r11549 = r11548 + r11525; double r11550 = 0.0424060604; double r11551 = r11550 * r11539; double r11552 = r11549 + r11551; double r11553 = r11546 + r11552; double r11554 = 0.0005064034; double r11555 = r11554 * r11541; double r11556 = r11553 + r11555; double r11557 = r11544 + r11556; double r11558 = r11542 * r11513; double r11559 = r11543 * r11531; double r11560 = r11558 * r11559; double r11561 = 0.0008327945; double r11562 = r11542 * r11561; double r11563 = 0.0694555761; double r11564 = r11563 * r11540; double r11565 = 0.7715471019; double r11566 = r11513 * r11565; double r11567 = r11566 + r11525; double r11568 = 0.2909738639; double r11569 = r11568 * r11539; double r11570 = r11567 + r11569; double r11571 = r11564 + r11570; double r11572 = 0.0140005442; double r11573 = r11541 * r11572; double r11574 = r11571 + r11573; double r11575 = r11562 + r11574; double r11576 = r11560 + r11575; double r11577 = r11557 / r11576; double r11578 = log(r11577); double r11579 = exp(r11578); double r11580 = r11511 * r11579; double r11581 = r11538 ? r11522 : r11580; double r11582 = r11524 ? r11536 : r11581; return r11582; } 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 r11583, r11584, r11585, r11586, r11587, r11588, r11589, r11590, r11591, r11592, r11593, r11594, r11595, r11596, r11597, r11598, r11599, r11600, r11601, r11602, r11603, r11604, r11605, r11606, r11607, r11608, r11609, r11610, r11611, r11612, r11613, r11614, r11615, r11616, r11617, r11618, r11619, r11620, r11621, r11622, r11623, r11624, r11625, r11626; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r11583, "1", 10, MPFR_RNDN); mpfr_init_set_str(r11584, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r11585); mpfr_init(r11586); mpfr_init(r11587); mpfr_init(r11588); mpfr_init_set_str(r11589, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r11590); mpfr_init(r11591); mpfr_init(r11592); mpfr_init_set_str(r11593, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r11594); mpfr_init(r11595); mpfr_init(r11596); mpfr_init_set_str(r11597, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r11598); mpfr_init(r11599); mpfr_init(r11600); mpfr_init_set_str(r11601, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r11602); mpfr_init(r11603); mpfr_init(r11604); mpfr_init_set_str(r11605, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r11606); mpfr_init(r11607); mpfr_init_set_str(r11608, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r11609); mpfr_init(r11610); mpfr_init_set_str(r11611, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r11612); mpfr_init(r11613); mpfr_init_set_str(r11614, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r11615); mpfr_init(r11616); mpfr_init_set_str(r11617, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r11618); mpfr_init(r11619); mpfr_init_set_str(r11620, "2", 10, MPFR_RNDN); mpfr_init(r11621); mpfr_init(r11622); mpfr_init(r11623); mpfr_init(r11624); mpfr_init(r11625); mpfr_init(r11626); } double f_im(double x) { ; ; mpfr_set_d(r11585, x, MPFR_RNDN); mpfr_mul(r11586, r11585, r11585, MPFR_RNDN); mpfr_mul(r11587, r11584, r11586, MPFR_RNDN); mpfr_add(r11588, r11583, r11587, MPFR_RNDN); ; mpfr_mul(r11590, r11586, r11586, MPFR_RNDN); mpfr_mul(r11591, r11589, r11590, MPFR_RNDN); mpfr_add(r11592, r11588, r11591, MPFR_RNDN); ; mpfr_mul(r11594, r11590, r11586, MPFR_RNDN); mpfr_mul(r11595, r11593, r11594, MPFR_RNDN); mpfr_add(r11596, r11592, r11595, MPFR_RNDN); ; mpfr_mul(r11598, r11594, r11586, MPFR_RNDN); mpfr_mul(r11599, r11597, r11598, MPFR_RNDN); mpfr_add(r11600, r11596, r11599, MPFR_RNDN); ; mpfr_mul(r11602, r11598, r11586, MPFR_RNDN); mpfr_mul(r11603, r11601, r11602, MPFR_RNDN); mpfr_add(r11604, r11600, r11603, MPFR_RNDN); ; mpfr_mul(r11606, r11605, r11586, MPFR_RNDN); mpfr_add(r11607, r11583, r11606, MPFR_RNDN); ; mpfr_mul(r11609, r11608, r11590, MPFR_RNDN); mpfr_add(r11610, r11607, r11609, MPFR_RNDN); ; mpfr_mul(r11612, r11611, r11594, MPFR_RNDN); mpfr_add(r11613, r11610, r11612, MPFR_RNDN); ; mpfr_mul(r11615, r11614, r11598, MPFR_RNDN); mpfr_add(r11616, r11613, r11615, MPFR_RNDN); ; mpfr_mul(r11618, r11617, r11602, MPFR_RNDN); mpfr_add(r11619, r11616, r11618, MPFR_RNDN); ; mpfr_mul(r11621, r11620, r11601, MPFR_RNDN); mpfr_mul(r11622, r11602, r11586, MPFR_RNDN); mpfr_mul(r11623, r11621, r11622, MPFR_RNDN); mpfr_add(r11624, r11619, r11623, MPFR_RNDN); mpfr_div(r11625, r11604, r11624, MPFR_RNDN); mpfr_mul(r11626, r11625, r11585, MPFR_RNDN); return mpfr_get_d(r11626, MPFR_RNDN); } static mpfr_t r11627, r11628, r11629, r11630, r11631, r11632, r11633, r11634, r11635, r11636, r11637, r11638, r11639, r11640, r11641, r11642, r11643, r11644, r11645, r11646, r11647, r11648, r11649, r11650, r11651, r11652, r11653, r11654, r11655, r11656, r11657, r11658, r11659, r11660, r11661, r11662, r11663, r11664, r11665, r11666, r11667, r11668, r11669, r11670, r11671, r11672, r11673, r11674, r11675, r11676, r11677, r11678, r11679, r11680, r11681, r11682, r11683, r11684, r11685, r11686, r11687, r11688, r11689, r11690, r11691, r11692, r11693, r11694, r11695, r11696, r11697, r11698, r11699; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r11627, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r11628); mpfr_init(r11629); mpfr_init(r11630); mpfr_init(r11631); mpfr_init_set_str(r11632, "0.5", 10, MPFR_RNDN); mpfr_init(r11633); mpfr_init_set_str(r11634, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init_set_str(r11635, "5", 10, MPFR_RNDN); mpfr_init(r11636); mpfr_init(r11637); mpfr_init(r11638); mpfr_init(r11639); mpfr_init_set_str(r11640, "-0.00040364450155768605", 10, MPFR_RNDN); mpfr_init(r11641); mpfr_init_set_str(r11642, "1", 10, MPFR_RNDN); mpfr_init_set_str(r11643, "0.265709700396151", 10, MPFR_RNDN); mpfr_init_set_str(r11644, "4", 10, MPFR_RNDN); mpfr_init(r11645); mpfr_init(r11646); mpfr_init(r11647); mpfr_init_set_str(r11648, "2", 10, MPFR_RNDN); mpfr_init(r11649); mpfr_init_set_str(r11650, "0.6665536072", 10, MPFR_RNDN); mpfr_init(r11651); mpfr_init(r11652); mpfr_init(r11653); mpfr_init_set_str(r11654, "8.03349330015686e-09", 10, MPFR_RNDN); mpfr_init(r11655); mpfr_init(r11656); mpfr_init(r11657); mpfr_init(r11658); mpfr_init(r11659); mpfr_init_set_str(r11660, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r11661); mpfr_init_set_str(r11662, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r11663); mpfr_init_set_str(r11664, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r11665); mpfr_init(r11666); mpfr_init_set_str(r11667, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r11668); mpfr_init(r11669); mpfr_init(r11670); mpfr_init_set_str(r11671, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r11672); mpfr_init(r11673); mpfr_init(r11674); mpfr_init(r11675); mpfr_init(r11676); mpfr_init(r11677); mpfr_init_set_str(r11678, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r11679); mpfr_init_set_str(r11680, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r11681); mpfr_init_set_str(r11682, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r11683); mpfr_init(r11684); mpfr_init_set_str(r11685, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r11686); mpfr_init(r11687); mpfr_init(r11688); mpfr_init_set_str(r11689, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r11690); mpfr_init(r11691); mpfr_init(r11692); mpfr_init(r11693); mpfr_init(r11694); mpfr_init(r11695); mpfr_init(r11696); mpfr_init(r11697); mpfr_init(r11698); mpfr_init(r11699); } double f_fm(double x) { ; mpfr_set_d(r11628, x, MPFR_RNDN); mpfr_div(r11629, r11627, r11628, MPFR_RNDN); mpfr_mul(r11630, r11628, r11628, MPFR_RNDN); mpfr_div(r11631, r11629, r11630, MPFR_RNDN); ; mpfr_div(r11633, r11632, r11628, MPFR_RNDN); ; ; mpfr_pow(r11636, r11628, r11635, MPFR_RNDN); mpfr_div(r11637, r11634, r11636, MPFR_RNDN); mpfr_add(r11638, r11633, r11637, MPFR_RNDN); mpfr_add(r11639, r11631, r11638, MPFR_RNDN); ; mpfr_set_si(r11641, mpfr_cmp(r11639, r11640) <= 0, MPFR_RNDN); ; ; ; mpfr_pow(r11645, r11628, r11644, MPFR_RNDN); mpfr_mul(r11646, r11643, r11645, MPFR_RNDN); mpfr_add(r11647, r11642, r11646, MPFR_RNDN); ; mpfr_pow(r11649, r11628, r11648, MPFR_RNDN); ; mpfr_mul(r11651, r11649, r11650, MPFR_RNDN); mpfr_sub(r11652, r11647, r11651, MPFR_RNDN); mpfr_mul(r11653, r11628, r11652, MPFR_RNDN); ; mpfr_set_si(r11655, mpfr_cmp(r11639, r11654) <= 0, MPFR_RNDN); mpfr_mul(r11656, r11630, r11630, MPFR_RNDN); mpfr_mul(r11657, r11630, r11656, MPFR_RNDN); mpfr_mul(r11658, r11630, r11657, MPFR_RNDN); mpfr_mul(r11659, r11658, r11630, MPFR_RNDN); ; mpfr_mul(r11661, r11659, r11660, MPFR_RNDN); ; mpfr_mul(r11663, r11662, r11657, MPFR_RNDN); ; mpfr_mul(r11665, r11630, r11664, MPFR_RNDN); mpfr_add(r11666, r11665, r11642, MPFR_RNDN); ; mpfr_mul(r11668, r11667, r11656, MPFR_RNDN); mpfr_add(r11669, r11666, r11668, MPFR_RNDN); mpfr_add(r11670, r11663, r11669, MPFR_RNDN); ; mpfr_mul(r11672, r11671, r11658, MPFR_RNDN); mpfr_add(r11673, r11670, r11672, MPFR_RNDN); mpfr_add(r11674, r11661, r11673, MPFR_RNDN); mpfr_mul(r11675, r11659, r11630, MPFR_RNDN); mpfr_mul(r11676, r11660, r11648, MPFR_RNDN); mpfr_mul(r11677, r11675, r11676, MPFR_RNDN); ; mpfr_mul(r11679, r11659, r11678, MPFR_RNDN); ; mpfr_mul(r11681, r11680, r11657, MPFR_RNDN); ; mpfr_mul(r11683, r11630, r11682, MPFR_RNDN); mpfr_add(r11684, r11683, r11642, MPFR_RNDN); ; mpfr_mul(r11686, r11685, r11656, MPFR_RNDN); mpfr_add(r11687, r11684, r11686, MPFR_RNDN); mpfr_add(r11688, r11681, r11687, MPFR_RNDN); ; mpfr_mul(r11690, r11658, r11689, MPFR_RNDN); mpfr_add(r11691, r11688, r11690, MPFR_RNDN); mpfr_add(r11692, r11679, r11691, MPFR_RNDN); mpfr_add(r11693, r11677, r11692, MPFR_RNDN); mpfr_div(r11694, r11674, r11693, MPFR_RNDN); mpfr_log(r11695, r11694, MPFR_RNDN); mpfr_exp(r11696, r11695, MPFR_RNDN); mpfr_mul(r11697, r11628, r11696, MPFR_RNDN); if (mpfr_get_si(r11655, MPFR_RNDN)) { mpfr_set(r11698, r11639, MPFR_RNDN); } else { mpfr_set(r11698, r11697, MPFR_RNDN); }; if (mpfr_get_si(r11641, MPFR_RNDN)) { mpfr_set(r11699, r11653, MPFR_RNDN); } else { mpfr_set(r11699, r11698, MPFR_RNDN); }; return mpfr_get_d(r11699, MPFR_RNDN); } static mpfr_t r11700, r11701, r11702, r11703, r11704, r11705, r11706, r11707, r11708, r11709, r11710, r11711, r11712, r11713, r11714, r11715, r11716, r11717, r11718, r11719, r11720, r11721, r11722, r11723, r11724, r11725, r11726, r11727, r11728, r11729, r11730, r11731, r11732, r11733, r11734, r11735, r11736, r11737, r11738, r11739, r11740, r11741, r11742, r11743, r11744, r11745, r11746, r11747, r11748, r11749, r11750, r11751, r11752, r11753, r11754, r11755, r11756, r11757, r11758, r11759, r11760, r11761, r11762, r11763, r11764, r11765, r11766, r11767, r11768, r11769, r11770, r11771, r11772; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r11700, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r11701); mpfr_init(r11702); mpfr_init(r11703); mpfr_init(r11704); mpfr_init_set_str(r11705, "0.5", 10, MPFR_RNDN); mpfr_init(r11706); mpfr_init_set_str(r11707, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init_set_str(r11708, "5", 10, MPFR_RNDN); mpfr_init(r11709); mpfr_init(r11710); mpfr_init(r11711); mpfr_init(r11712); mpfr_init_set_str(r11713, "-0.00040364450155768605", 10, MPFR_RNDN); mpfr_init(r11714); mpfr_init_set_str(r11715, "1", 10, MPFR_RNDN); mpfr_init_set_str(r11716, "0.265709700396151", 10, MPFR_RNDN); mpfr_init_set_str(r11717, "4", 10, MPFR_RNDN); mpfr_init(r11718); mpfr_init(r11719); mpfr_init(r11720); mpfr_init_set_str(r11721, "2", 10, MPFR_RNDN); mpfr_init(r11722); mpfr_init_set_str(r11723, "0.6665536072", 10, MPFR_RNDN); mpfr_init(r11724); mpfr_init(r11725); mpfr_init(r11726); mpfr_init_set_str(r11727, "8.03349330015686e-09", 10, MPFR_RNDN); mpfr_init(r11728); mpfr_init(r11729); mpfr_init(r11730); mpfr_init(r11731); mpfr_init(r11732); mpfr_init_set_str(r11733, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r11734); mpfr_init_set_str(r11735, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r11736); mpfr_init_set_str(r11737, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r11738); mpfr_init(r11739); mpfr_init_set_str(r11740, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r11741); mpfr_init(r11742); mpfr_init(r11743); mpfr_init_set_str(r11744, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r11745); mpfr_init(r11746); mpfr_init(r11747); mpfr_init(r11748); mpfr_init(r11749); mpfr_init(r11750); mpfr_init_set_str(r11751, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r11752); mpfr_init_set_str(r11753, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r11754); mpfr_init_set_str(r11755, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r11756); mpfr_init(r11757); mpfr_init_set_str(r11758, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r11759); mpfr_init(r11760); mpfr_init(r11761); mpfr_init_set_str(r11762, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r11763); mpfr_init(r11764); mpfr_init(r11765); mpfr_init(r11766); mpfr_init(r11767); mpfr_init(r11768); mpfr_init(r11769); mpfr_init(r11770); mpfr_init(r11771); mpfr_init(r11772); } double f_dm(double x) { ; mpfr_set_d(r11701, x, MPFR_RNDN); mpfr_div(r11702, r11700, r11701, MPFR_RNDN); mpfr_mul(r11703, r11701, r11701, MPFR_RNDN); mpfr_div(r11704, r11702, r11703, MPFR_RNDN); ; mpfr_div(r11706, r11705, r11701, MPFR_RNDN); ; ; mpfr_pow(r11709, r11701, r11708, MPFR_RNDN); mpfr_div(r11710, r11707, r11709, MPFR_RNDN); mpfr_add(r11711, r11706, r11710, MPFR_RNDN); mpfr_add(r11712, r11704, r11711, MPFR_RNDN); ; mpfr_set_si(r11714, mpfr_cmp(r11712, r11713) <= 0, MPFR_RNDN); ; ; ; mpfr_pow(r11718, r11701, r11717, MPFR_RNDN); mpfr_mul(r11719, r11716, r11718, MPFR_RNDN); mpfr_add(r11720, r11715, r11719, MPFR_RNDN); ; mpfr_pow(r11722, r11701, r11721, MPFR_RNDN); ; mpfr_mul(r11724, r11722, r11723, MPFR_RNDN); mpfr_sub(r11725, r11720, r11724, MPFR_RNDN); mpfr_mul(r11726, r11701, r11725, MPFR_RNDN); ; mpfr_set_si(r11728, mpfr_cmp(r11712, r11727) <= 0, MPFR_RNDN); mpfr_mul(r11729, r11703, r11703, MPFR_RNDN); mpfr_mul(r11730, r11703, r11729, MPFR_RNDN); mpfr_mul(r11731, r11703, r11730, MPFR_RNDN); mpfr_mul(r11732, r11731, r11703, MPFR_RNDN); ; mpfr_mul(r11734, r11732, r11733, MPFR_RNDN); ; mpfr_mul(r11736, r11735, r11730, MPFR_RNDN); ; mpfr_mul(r11738, r11703, r11737, MPFR_RNDN); mpfr_add(r11739, r11738, r11715, MPFR_RNDN); ; mpfr_mul(r11741, r11740, r11729, MPFR_RNDN); mpfr_add(r11742, r11739, r11741, MPFR_RNDN); mpfr_add(r11743, r11736, r11742, MPFR_RNDN); ; mpfr_mul(r11745, r11744, r11731, MPFR_RNDN); mpfr_add(r11746, r11743, r11745, MPFR_RNDN); mpfr_add(r11747, r11734, r11746, MPFR_RNDN); mpfr_mul(r11748, r11732, r11703, MPFR_RNDN); mpfr_mul(r11749, r11733, r11721, MPFR_RNDN); mpfr_mul(r11750, r11748, r11749, MPFR_RNDN); ; mpfr_mul(r11752, r11732, r11751, MPFR_RNDN); ; mpfr_mul(r11754, r11753, r11730, MPFR_RNDN); ; mpfr_mul(r11756, r11703, r11755, MPFR_RNDN); mpfr_add(r11757, r11756, r11715, MPFR_RNDN); ; mpfr_mul(r11759, r11758, r11729, MPFR_RNDN); mpfr_add(r11760, r11757, r11759, MPFR_RNDN); mpfr_add(r11761, r11754, r11760, MPFR_RNDN); ; mpfr_mul(r11763, r11731, r11762, MPFR_RNDN); mpfr_add(r11764, r11761, r11763, MPFR_RNDN); mpfr_add(r11765, r11752, r11764, MPFR_RNDN); mpfr_add(r11766, r11750, r11765, MPFR_RNDN); mpfr_div(r11767, r11747, r11766, MPFR_RNDN); mpfr_log(r11768, r11767, MPFR_RNDN); mpfr_exp(r11769, r11768, MPFR_RNDN); mpfr_mul(r11770, r11701, r11769, MPFR_RNDN); if (mpfr_get_si(r11728, MPFR_RNDN)) { mpfr_set(r11771, r11712, MPFR_RNDN); } else { mpfr_set(r11771, r11770, MPFR_RNDN); }; if (mpfr_get_si(r11714, MPFR_RNDN)) { mpfr_set(r11772, r11726, MPFR_RNDN); } else { mpfr_set(r11772, r11771, MPFR_RNDN); }; return mpfr_get_d(r11772, MPFR_RNDN); }