#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r19431 = 1.0f; float r19432 = 0.10499349236488342f; float r19433 = x; float r19434 = r19433 * r19433; float r19435 = r19432 * r19434; float r19436 = r19431 + r19435; float r19437 = 0.04240605980157852f; float r19438 = r19434 * r19434; float r19439 = r19437 * r19438; float r19440 = r19436 + r19439; float r19441 = 0.007264418061822653f; float r19442 = r19438 * r19434; float r19443 = r19441 * r19442; float r19444 = r19440 + r19443; float r19445 = 0.0005064033903181553f; float r19446 = r19442 * r19434; float r19447 = r19445 * r19446; float r19448 = r19444 + r19447; float r19449 = 0.00017899709928315133f; float r19450 = r19446 * r19434; float r19451 = r19449 * r19450; float r19452 = r19448 + r19451; float r19453 = 0.7715470790863037f; float r19454 = r19453 * r19434; float r19455 = r19431 + r19454; float r19456 = 0.29097387194633484f; float r19457 = r19456 * r19438; float r19458 = r19455 + r19457; float r19459 = 0.0694555789232254f; float r19460 = r19459 * r19442; float r19461 = r19458 + r19460; float r19462 = 0.01400054432451725f; float r19463 = r19462 * r19446; float r19464 = r19461 + r19463; float r19465 = 0.0008327945251949131f; float r19466 = r19465 * r19450; float r19467 = r19464 + r19466; float r19468 = 2.0f; float r19469 = r19468 * r19449; float r19470 = r19450 * r19434; float r19471 = r19469 * r19470; float r19472 = r19467 + r19471; float r19473 = r19452 / r19472; float r19474 = r19473 * r19433; return r19474; } double f_id(double x) { double r19475 = 1.0; double r19476 = 0.10499349236488342; double r19477 = x; double r19478 = r19477 * r19477; double r19479 = r19476 * r19478; double r19480 = r19475 + r19479; double r19481 = 0.04240605980157852; double r19482 = r19478 * r19478; double r19483 = r19481 * r19482; double r19484 = r19480 + r19483; double r19485 = 0.007264418061822653; double r19486 = r19482 * r19478; double r19487 = r19485 * r19486; double r19488 = r19484 + r19487; double r19489 = 0.0005064033903181553; double r19490 = r19486 * r19478; double r19491 = r19489 * r19490; double r19492 = r19488 + r19491; double r19493 = 0.00017899709928315133; double r19494 = r19490 * r19478; double r19495 = r19493 * r19494; double r19496 = r19492 + r19495; double r19497 = 0.7715470790863037; double r19498 = r19497 * r19478; double r19499 = r19475 + r19498; double r19500 = 0.29097387194633484; double r19501 = r19500 * r19482; double r19502 = r19499 + r19501; double r19503 = 0.0694555789232254; double r19504 = r19503 * r19486; double r19505 = r19502 + r19504; double r19506 = 0.01400054432451725; double r19507 = r19506 * r19490; double r19508 = r19505 + r19507; double r19509 = 0.0008327945251949131; double r19510 = r19509 * r19494; double r19511 = r19508 + r19510; double r19512 = 2.0; double r19513 = r19512 * r19493; double r19514 = r19494 * r19478; double r19515 = r19513 * r19514; double r19516 = r19511 + r19515; double r19517 = r19496 / r19516; double r19518 = r19517 * r19477; return r19518; } double f_of(float x) { float r19519 = x; float r19520 = r19519 * r19519; float r19521 = r19520 * r19520; float r19522 = r19519 * (r19519 * r19519); float r19523 = r19522 * r19522; float r19524 = r19521 * r19523; float r19525 = 2.0f; float r19526 = 0.00017899709928315133f; float r19527 = r19525 * r19526; float r19528 = cbrt(r19527); float r19529 = cbrt(r19520); float r19530 = r19528 * r19529; float r19531 = r19530 * (r19530 * r19530); float r19532 = 0.0008327945251949131f; float r19533 = r19531 + r19532; float r19534 = r19524 * r19533; float r19535 = 1.0f; float r19536 = 0.29097387194633484f; float r19537 = r19519 * r19536; float r19538 = r19537 * r19522; float r19539 = r19535 + r19538; float r19540 = 0.7715470790863037f; float r19541 = r19540 * r19519; float r19542 = r19541 * r19519; float r19543 = r19539 + r19542; float r19544 = r19534 + r19543; float r19545 = 0.0694555789232254f; float r19546 = 0.01400054432451725f; float r19547 = r19546 * r19520; float r19548 = r19545 + r19547; float r19549 = r19521 * r19548; float r19550 = r19520 * r19549; float r19551 = r19544 + r19550; float r19552 = r19519 / r19551; float r19553 = 0.04240605980157852f; float r19554 = r19519 * r19553; float r19555 = r19554 * r19522; float r19556 = 0.10499349236488342f; float r19557 = r19520 * r19556; float r19558 = r19557 + r19535; float r19559 = r19555 + r19558; float r19560 = r19520 * r19526; float r19561 = r19521 * r19521; float r19562 = r19560 * r19561; float r19563 = r19559 + r19562; float r19564 = 0.007264418061822653f; float r19565 = 0.0005064033903181553f; float r19566 = r19565 * r19519; float r19567 = r19566 * r19519; float r19568 = r19564 + r19567; float r19569 = r19521 * r19568; float r19570 = r19520 * r19569; float r19571 = r19563 + r19570; float r19572 = r19552 * r19571; return r19572; } double f_od(double x) { double r19573 = x; double r19574 = r19573 * r19573; double r19575 = r19574 * r19574; double r19576 = r19573 * (r19573 * r19573); double r19577 = r19576 * r19576; double r19578 = r19575 * r19577; double r19579 = 2.0; double r19580 = 0.00017899709928315133; double r19581 = r19579 * r19580; double r19582 = cbrt(r19581); double r19583 = cbrt(r19574); double r19584 = r19582 * r19583; double r19585 = r19584 * (r19584 * r19584); double r19586 = 0.0008327945251949131; double r19587 = r19585 + r19586; double r19588 = r19578 * r19587; double r19589 = 1.0; double r19590 = 0.29097387194633484; double r19591 = r19573 * r19590; double r19592 = r19591 * r19576; double r19593 = r19589 + r19592; double r19594 = 0.7715470790863037; double r19595 = r19594 * r19573; double r19596 = r19595 * r19573; double r19597 = r19593 + r19596; double r19598 = r19588 + r19597; double r19599 = 0.0694555789232254; double r19600 = 0.01400054432451725; double r19601 = r19600 * r19574; double r19602 = r19599 + r19601; double r19603 = r19575 * r19602; double r19604 = r19574 * r19603; double r19605 = r19598 + r19604; double r19606 = r19573 / r19605; double r19607 = 0.04240605980157852; double r19608 = r19573 * r19607; double r19609 = r19608 * r19576; double r19610 = 0.10499349236488342; double r19611 = r19574 * r19610; double r19612 = r19611 + r19589; double r19613 = r19609 + r19612; double r19614 = r19574 * r19580; double r19615 = r19575 * r19575; double r19616 = r19614 * r19615; double r19617 = r19613 + r19616; double r19618 = 0.007264418061822653; double r19619 = 0.0005064033903181553; double r19620 = r19619 * r19573; double r19621 = r19620 * r19573; double r19622 = r19618 + r19621; double r19623 = r19575 * r19622; double r19624 = r19574 * r19623; double r19625 = r19617 + r19624; double r19626 = r19606 * r19625; return r19626; } 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 r19627, r19628, r19629, r19630, r19631, r19632, r19633, r19634, r19635, r19636, r19637, r19638, r19639, r19640, r19641, r19642, r19643, r19644, r19645, r19646, r19647, r19648, r19649, r19650, r19651, r19652, r19653, r19654, r19655, r19656, r19657, r19658, r19659, r19660, r19661, r19662, r19663, r19664, r19665, r19666, r19667, r19668, r19669, r19670; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19627, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19628, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19629); mpfr_init(r19630); mpfr_init(r19631); mpfr_init(r19632); mpfr_init_set_str(r19633, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19634); mpfr_init(r19635); mpfr_init(r19636); mpfr_init_set_str(r19637, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19638); mpfr_init(r19639); mpfr_init(r19640); mpfr_init_set_str(r19641, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19642); mpfr_init(r19643); mpfr_init(r19644); mpfr_init_set_str(r19645, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19646); mpfr_init(r19647); mpfr_init(r19648); mpfr_init_set_str(r19649, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19650); mpfr_init(r19651); mpfr_init_set_str(r19652, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19653); mpfr_init(r19654); mpfr_init_set_str(r19655, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19656); mpfr_init(r19657); mpfr_init_set_str(r19658, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19659); mpfr_init(r19660); mpfr_init_set_str(r19661, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19662); mpfr_init(r19663); mpfr_init_set_str(r19664, "2", 10, MPFR_RNDN); mpfr_init(r19665); mpfr_init(r19666); mpfr_init(r19667); mpfr_init(r19668); mpfr_init(r19669); mpfr_init(r19670); } double f_im(double x) { ; ; mpfr_set_d(r19629, x, MPFR_RNDN); mpfr_mul(r19630, r19629, r19629, MPFR_RNDN); mpfr_mul(r19631, r19628, r19630, MPFR_RNDN); mpfr_add(r19632, r19627, r19631, MPFR_RNDN); ; mpfr_mul(r19634, r19630, r19630, MPFR_RNDN); mpfr_mul(r19635, r19633, r19634, MPFR_RNDN); mpfr_add(r19636, r19632, r19635, MPFR_RNDN); ; mpfr_mul(r19638, r19634, r19630, MPFR_RNDN); mpfr_mul(r19639, r19637, r19638, MPFR_RNDN); mpfr_add(r19640, r19636, r19639, MPFR_RNDN); ; mpfr_mul(r19642, r19638, r19630, MPFR_RNDN); mpfr_mul(r19643, r19641, r19642, MPFR_RNDN); mpfr_add(r19644, r19640, r19643, MPFR_RNDN); ; mpfr_mul(r19646, r19642, r19630, MPFR_RNDN); mpfr_mul(r19647, r19645, r19646, MPFR_RNDN); mpfr_add(r19648, r19644, r19647, MPFR_RNDN); ; mpfr_mul(r19650, r19649, r19630, MPFR_RNDN); mpfr_add(r19651, r19627, r19650, MPFR_RNDN); ; mpfr_mul(r19653, r19652, r19634, MPFR_RNDN); mpfr_add(r19654, r19651, r19653, MPFR_RNDN); ; mpfr_mul(r19656, r19655, r19638, MPFR_RNDN); mpfr_add(r19657, r19654, r19656, MPFR_RNDN); ; mpfr_mul(r19659, r19658, r19642, MPFR_RNDN); mpfr_add(r19660, r19657, r19659, MPFR_RNDN); ; mpfr_mul(r19662, r19661, r19646, MPFR_RNDN); mpfr_add(r19663, r19660, r19662, MPFR_RNDN); ; mpfr_mul(r19665, r19664, r19645, MPFR_RNDN); mpfr_mul(r19666, r19646, r19630, MPFR_RNDN); mpfr_mul(r19667, r19665, r19666, MPFR_RNDN); mpfr_add(r19668, r19663, r19667, MPFR_RNDN); mpfr_div(r19669, r19648, r19668, MPFR_RNDN); mpfr_mul(r19670, r19669, r19629, MPFR_RNDN); return mpfr_get_d(r19670, MPFR_RNDN); } static mpfr_t r19671, r19672, r19673, r19674, r19675, r19676, r19677, r19678, r19679, r19680, r19681, r19682, r19683, r19684, r19685, r19686, r19687, r19688, r19689, r19690, r19691, r19692, r19693, r19694, r19695, r19696, r19697, r19698, r19699, r19700, r19701, r19702, r19703, r19704, r19705, r19706, r19707, r19708, r19709, r19710, r19711, r19712, r19713, r19714, r19715, r19716, r19717, r19718, r19719, r19720, r19721, r19722, r19723, r19724; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19671); mpfr_init(r19672); mpfr_init(r19673); mpfr_init(r19674); mpfr_init(r19675); mpfr_init(r19676); mpfr_init_set_str(r19677, "2", 10, MPFR_RNDN); mpfr_init_set_str(r19678, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19679); mpfr_init(r19680); mpfr_init(r19681); mpfr_init(r19682); mpfr_init(r19683); mpfr_init_set_str(r19684, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19685); mpfr_init(r19686); mpfr_init_set_str(r19687, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19688, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19689); mpfr_init(r19690); mpfr_init(r19691); mpfr_init_set_str(r19692, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19693); mpfr_init(r19694); mpfr_init(r19695); mpfr_init(r19696); mpfr_init_set_str(r19697, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19698, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19699); mpfr_init(r19700); mpfr_init(r19701); mpfr_init(r19702); mpfr_init(r19703); mpfr_init(r19704); mpfr_init_set_str(r19705, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19706); mpfr_init(r19707); mpfr_init_set_str(r19708, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19709); mpfr_init(r19710); mpfr_init(r19711); mpfr_init(r19712); mpfr_init(r19713); mpfr_init(r19714); mpfr_init(r19715); mpfr_init_set_str(r19716, "0.0072644182", 10, MPFR_RNDN); mpfr_init_set_str(r19717, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19718); mpfr_init(r19719); mpfr_init(r19720); mpfr_init(r19721); mpfr_init(r19722); mpfr_init(r19723); mpfr_init(r19724); } double f_fm(double x) { mpfr_set_d(r19671, x, MPFR_RNDN); mpfr_mul(r19672, r19671, r19671, MPFR_RNDN); mpfr_mul(r19673, r19672, r19672, MPFR_RNDN); mpfr_mul(r19674, r19671, r19671, MPFR_RNDN); mpfr_mul(r19674, r19674, r19671, MPFR_RNDN); mpfr_mul(r19675, r19674, r19674, MPFR_RNDN); mpfr_mul(r19676, r19673, r19675, MPFR_RNDN); ; ; mpfr_mul(r19679, r19677, r19678, MPFR_RNDN); mpfr_cbrt(r19680, r19679, MPFR_RNDN); mpfr_cbrt(r19681, r19672, MPFR_RNDN); mpfr_mul(r19682, r19680, r19681, MPFR_RNDN); mpfr_mul(r19683, r19682, r19682, MPFR_RNDN); mpfr_mul(r19683, r19683, r19682, MPFR_RNDN); ; mpfr_add(r19685, r19683, r19684, MPFR_RNDN); mpfr_mul(r19686, r19676, r19685, MPFR_RNDN); ; ; mpfr_mul(r19689, r19671, r19688, MPFR_RNDN); mpfr_mul(r19690, r19689, r19674, MPFR_RNDN); mpfr_add(r19691, r19687, r19690, MPFR_RNDN); ; mpfr_mul(r19693, r19692, r19671, MPFR_RNDN); mpfr_mul(r19694, r19693, r19671, MPFR_RNDN); mpfr_add(r19695, r19691, r19694, MPFR_RNDN); mpfr_add(r19696, r19686, r19695, MPFR_RNDN); ; ; mpfr_mul(r19699, r19698, r19672, MPFR_RNDN); mpfr_add(r19700, r19697, r19699, MPFR_RNDN); mpfr_mul(r19701, r19673, r19700, MPFR_RNDN); mpfr_mul(r19702, r19672, r19701, MPFR_RNDN); mpfr_add(r19703, r19696, r19702, MPFR_RNDN); mpfr_div(r19704, r19671, r19703, MPFR_RNDN); ; mpfr_mul(r19706, r19671, r19705, MPFR_RNDN); mpfr_mul(r19707, r19706, r19674, MPFR_RNDN); ; mpfr_mul(r19709, r19672, r19708, MPFR_RNDN); mpfr_add(r19710, r19709, r19687, MPFR_RNDN); mpfr_add(r19711, r19707, r19710, MPFR_RNDN); mpfr_mul(r19712, r19672, r19678, MPFR_RNDN); mpfr_mul(r19713, r19673, r19673, MPFR_RNDN); mpfr_mul(r19714, r19712, r19713, MPFR_RNDN); mpfr_add(r19715, r19711, r19714, MPFR_RNDN); ; ; mpfr_mul(r19718, r19717, r19671, MPFR_RNDN); mpfr_mul(r19719, r19718, r19671, MPFR_RNDN); mpfr_add(r19720, r19716, r19719, MPFR_RNDN); mpfr_mul(r19721, r19673, r19720, MPFR_RNDN); mpfr_mul(r19722, r19672, r19721, MPFR_RNDN); mpfr_add(r19723, r19715, r19722, MPFR_RNDN); mpfr_mul(r19724, r19704, r19723, MPFR_RNDN); return mpfr_get_d(r19724, MPFR_RNDN); } static mpfr_t r19725, r19726, r19727, r19728, r19729, r19730, r19731, r19732, r19733, r19734, r19735, r19736, r19737, r19738, r19739, r19740, r19741, r19742, r19743, r19744, r19745, r19746, r19747, r19748, r19749, r19750, r19751, r19752, r19753, r19754, r19755, r19756, r19757, r19758, r19759, r19760, r19761, r19762, r19763, r19764, r19765, r19766, r19767, r19768, r19769, r19770, r19771, r19772, r19773, r19774, r19775, r19776, r19777, r19778; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19725); mpfr_init(r19726); mpfr_init(r19727); mpfr_init(r19728); mpfr_init(r19729); mpfr_init(r19730); mpfr_init_set_str(r19731, "2", 10, MPFR_RNDN); mpfr_init_set_str(r19732, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19733); mpfr_init(r19734); mpfr_init(r19735); mpfr_init(r19736); mpfr_init(r19737); mpfr_init_set_str(r19738, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19739); mpfr_init(r19740); mpfr_init_set_str(r19741, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19742, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19743); mpfr_init(r19744); mpfr_init(r19745); mpfr_init_set_str(r19746, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19747); mpfr_init(r19748); mpfr_init(r19749); mpfr_init(r19750); mpfr_init_set_str(r19751, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19752, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19753); mpfr_init(r19754); mpfr_init(r19755); mpfr_init(r19756); mpfr_init(r19757); mpfr_init(r19758); mpfr_init_set_str(r19759, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19760); mpfr_init(r19761); mpfr_init_set_str(r19762, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19763); mpfr_init(r19764); mpfr_init(r19765); mpfr_init(r19766); mpfr_init(r19767); mpfr_init(r19768); mpfr_init(r19769); mpfr_init_set_str(r19770, "0.0072644182", 10, MPFR_RNDN); mpfr_init_set_str(r19771, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19772); mpfr_init(r19773); mpfr_init(r19774); mpfr_init(r19775); mpfr_init(r19776); mpfr_init(r19777); mpfr_init(r19778); } double f_dm(double x) { mpfr_set_d(r19725, x, MPFR_RNDN); mpfr_mul(r19726, r19725, r19725, MPFR_RNDN); mpfr_mul(r19727, r19726, r19726, MPFR_RNDN); mpfr_mul(r19728, r19725, r19725, MPFR_RNDN); mpfr_mul(r19728, r19728, r19725, MPFR_RNDN); mpfr_mul(r19729, r19728, r19728, MPFR_RNDN); mpfr_mul(r19730, r19727, r19729, MPFR_RNDN); ; ; mpfr_mul(r19733, r19731, r19732, MPFR_RNDN); mpfr_cbrt(r19734, r19733, MPFR_RNDN); mpfr_cbrt(r19735, r19726, MPFR_RNDN); mpfr_mul(r19736, r19734, r19735, MPFR_RNDN); mpfr_mul(r19737, r19736, r19736, MPFR_RNDN); mpfr_mul(r19737, r19737, r19736, MPFR_RNDN); ; mpfr_add(r19739, r19737, r19738, MPFR_RNDN); mpfr_mul(r19740, r19730, r19739, MPFR_RNDN); ; ; mpfr_mul(r19743, r19725, r19742, MPFR_RNDN); mpfr_mul(r19744, r19743, r19728, MPFR_RNDN); mpfr_add(r19745, r19741, r19744, MPFR_RNDN); ; mpfr_mul(r19747, r19746, r19725, MPFR_RNDN); mpfr_mul(r19748, r19747, r19725, MPFR_RNDN); mpfr_add(r19749, r19745, r19748, MPFR_RNDN); mpfr_add(r19750, r19740, r19749, MPFR_RNDN); ; ; mpfr_mul(r19753, r19752, r19726, MPFR_RNDN); mpfr_add(r19754, r19751, r19753, MPFR_RNDN); mpfr_mul(r19755, r19727, r19754, MPFR_RNDN); mpfr_mul(r19756, r19726, r19755, MPFR_RNDN); mpfr_add(r19757, r19750, r19756, MPFR_RNDN); mpfr_div(r19758, r19725, r19757, MPFR_RNDN); ; mpfr_mul(r19760, r19725, r19759, MPFR_RNDN); mpfr_mul(r19761, r19760, r19728, MPFR_RNDN); ; mpfr_mul(r19763, r19726, r19762, MPFR_RNDN); mpfr_add(r19764, r19763, r19741, MPFR_RNDN); mpfr_add(r19765, r19761, r19764, MPFR_RNDN); mpfr_mul(r19766, r19726, r19732, MPFR_RNDN); mpfr_mul(r19767, r19727, r19727, MPFR_RNDN); mpfr_mul(r19768, r19766, r19767, MPFR_RNDN); mpfr_add(r19769, r19765, r19768, MPFR_RNDN); ; ; mpfr_mul(r19772, r19771, r19725, MPFR_RNDN); mpfr_mul(r19773, r19772, r19725, MPFR_RNDN); mpfr_add(r19774, r19770, r19773, MPFR_RNDN); mpfr_mul(r19775, r19727, r19774, MPFR_RNDN); mpfr_mul(r19776, r19726, r19775, MPFR_RNDN); mpfr_add(r19777, r19769, r19776, MPFR_RNDN); mpfr_mul(r19778, r19758, r19777, MPFR_RNDN); return mpfr_get_d(r19778, MPFR_RNDN); }