#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r19300 = 1.0f; float r19301 = 0.1049934947f; float r19302 = x; float r19303 = r19302 * r19302; float r19304 = r19301 * r19303; float r19305 = r19300 + r19304; float r19306 = 0.0424060604f; float r19307 = r19303 * r19303; float r19308 = r19306 * r19307; float r19309 = r19305 + r19308; float r19310 = 0.0072644182f; float r19311 = r19307 * r19303; float r19312 = r19310 * r19311; float r19313 = r19309 + r19312; float r19314 = 0.0005064034f; float r19315 = r19311 * r19303; float r19316 = r19314 * r19315; float r19317 = r19313 + r19316; float r19318 = 0.0001789971f; float r19319 = r19315 * r19303; float r19320 = r19318 * r19319; float r19321 = r19317 + r19320; float r19322 = 0.7715471019f; float r19323 = r19322 * r19303; float r19324 = r19300 + r19323; float r19325 = 0.2909738639f; float r19326 = r19325 * r19307; float r19327 = r19324 + r19326; float r19328 = 0.0694555761f; float r19329 = r19328 * r19311; float r19330 = r19327 + r19329; float r19331 = 0.0140005442f; float r19332 = r19331 * r19315; float r19333 = r19330 + r19332; float r19334 = 0.0008327945f; float r19335 = r19334 * r19319; float r19336 = r19333 + r19335; float r19337 = 2.0f; float r19338 = r19337 * r19318; float r19339 = r19319 * r19303; float r19340 = r19338 * r19339; float r19341 = r19336 + r19340; float r19342 = r19321 / r19341; float r19343 = r19342 * r19302; return r19343; } double f_id(double x) { double r19344 = 1.0; double r19345 = 0.1049934947; double r19346 = x; double r19347 = r19346 * r19346; double r19348 = r19345 * r19347; double r19349 = r19344 + r19348; double r19350 = 0.0424060604; double r19351 = r19347 * r19347; double r19352 = r19350 * r19351; double r19353 = r19349 + r19352; double r19354 = 0.0072644182; double r19355 = r19351 * r19347; double r19356 = r19354 * r19355; double r19357 = r19353 + r19356; double r19358 = 0.0005064034; double r19359 = r19355 * r19347; double r19360 = r19358 * r19359; double r19361 = r19357 + r19360; double r19362 = 0.0001789971; double r19363 = r19359 * r19347; double r19364 = r19362 * r19363; double r19365 = r19361 + r19364; double r19366 = 0.7715471019; double r19367 = r19366 * r19347; double r19368 = r19344 + r19367; double r19369 = 0.2909738639; double r19370 = r19369 * r19351; double r19371 = r19368 + r19370; double r19372 = 0.0694555761; double r19373 = r19372 * r19355; double r19374 = r19371 + r19373; double r19375 = 0.0140005442; double r19376 = r19375 * r19359; double r19377 = r19374 + r19376; double r19378 = 0.0008327945; double r19379 = r19378 * r19363; double r19380 = r19377 + r19379; double r19381 = 2.0; double r19382 = r19381 * r19362; double r19383 = r19363 * r19347; double r19384 = r19382 * r19383; double r19385 = r19380 + r19384; double r19386 = r19365 / r19385; double r19387 = r19386 * r19346; return r19387; } double f_of(float x) { float r19388 = x; float r19389 = -1.1392497874101738e+24f; bool r19390 = r19388 <= r19389; float r19391 = 1.0f; float r19392 = 0.1049934947f; float r19393 = r19388 * r19388; float r19394 = r19392 / r19393; float r19395 = r19391 + r19394; float r19396 = 0.0072644182f; float r19397 = r19396 * r19391; float r19398 = r19388 * (r19388 * r19388); float r19399 = r19398 * r19398; float r19400 = r19397 / r19399; float r19401 = r19395 + r19400; float r19402 = 0.0424060604f; float r19403 = r19402 * r19391; float r19404 = r19403 / r19398; float r19405 = r19404 / r19388; float r19406 = 0.0001789971f; float r19407 = 10.0f; float r19408 = pow(r19388, r19407); float r19409 = r19406 / r19408; float r19410 = r19405 + r19409; float r19411 = 0.0005064034f; float r19412 = r19391 * r19411; float r19413 = r19412 / r19399; float r19414 = r19413 / r19393; float r19415 = r19410 + r19414; float r19416 = r19401 + r19415; float r19417 = r19391 / r19388; float r19418 = r19417 / r19398; float r19419 = 0.2909738639f; float r19420 = r19418 * r19419; float r19421 = r19391 + r19420; float r19422 = 0.0694555761f; float r19423 = r19422 * r19391; float r19424 = r19423 / r19399; float r19425 = r19421 + r19424; float r19426 = 0.0140005442f; float r19427 = r19391 * r19426; float r19428 = r19427 / r19399; float r19429 = r19428 / r19393; float r19430 = r19391 / r19398; float r19431 = r19430 * r19430; float r19432 = 6.0f; float r19433 = pow(r19388, r19432); float r19434 = 0.0003579942f; float r19435 = r19433 / r19434; float r19436 = r19431 / r19435; float r19437 = r19429 + r19436; float r19438 = r19425 + r19437; float r19439 = 0.7715471019f; float r19440 = r19439 / r19393; float r19441 = 4.0f; float r19442 = pow(r19388, r19441); float r19443 = 0.0008327945f; float r19444 = r19442 / r19443; float r19445 = r19431 / r19444; float r19446 = r19440 + r19445; float r19447 = r19438 + r19446; float r19448 = r19447 * r19388; float r19449 = r19416 / r19448; float r19450 = 1.653196896254347e+29f; bool r19451 = r19388 <= r19450; float r19452 = r19388 * r19388; float r19453 = r19452 * (r19452 * r19452); float r19454 = r19426 * r19388; float r19455 = r19388 * r19454; float r19456 = r19455 + r19422; float r19457 = r19453 * r19456; float r19458 = r19452 * r19452; float r19459 = r19458 * r19453; float r19460 = 2.0f; float r19461 = r19460 * r19406; float r19462 = r19452 * r19461; float r19463 = r19462 + r19443; float r19464 = r19459 * r19463; float r19465 = r19457 + r19464; float r19466 = r19388 * r19419; float r19467 = r19466 * r19398; float r19468 = r19452 * r19439; float r19469 = r19467 + r19468; float r19470 = r19469 + r19391; float r19471 = r19465 + r19470; float r19472 = r19458 * r19458; float r19473 = r19452 * r19406; float r19474 = r19472 * r19473; float r19475 = r19388 * r19402; float r19476 = r19475 * r19398; float r19477 = r19392 * r19388; float r19478 = r19388 * r19477; float r19479 = r19391 + r19478; float r19480 = r19476 + r19479; float r19481 = r19474 + r19480; float r19482 = r19388 * r19411; float r19483 = r19482 * r19388; float r19484 = r19396 + r19483; float r19485 = r19484 * r19453; float r19486 = r19481 + r19485; float r19487 = r19471 / r19486; float r19488 = r19388 / r19487; float r19489 = pow(r19488, r19391); float r19490 = r19451 ? r19489 : r19449; float r19491 = r19390 ? r19449 : r19490; return r19491; } double f_od(double x) { double r19492 = x; double r19493 = -1.1392497874101738e+24; bool r19494 = r19492 <= r19493; double r19495 = 1.0; double r19496 = 0.1049934947; double r19497 = r19492 * r19492; double r19498 = r19496 / r19497; double r19499 = r19495 + r19498; double r19500 = 0.0072644182; double r19501 = r19500 * r19495; double r19502 = r19492 * (r19492 * r19492); double r19503 = r19502 * r19502; double r19504 = r19501 / r19503; double r19505 = r19499 + r19504; double r19506 = 0.0424060604; double r19507 = r19506 * r19495; double r19508 = r19507 / r19502; double r19509 = r19508 / r19492; double r19510 = 0.0001789971; double r19511 = 10.0; double r19512 = pow(r19492, r19511); double r19513 = r19510 / r19512; double r19514 = r19509 + r19513; double r19515 = 0.0005064034; double r19516 = r19495 * r19515; double r19517 = r19516 / r19503; double r19518 = r19517 / r19497; double r19519 = r19514 + r19518; double r19520 = r19505 + r19519; double r19521 = r19495 / r19492; double r19522 = r19521 / r19502; double r19523 = 0.2909738639; double r19524 = r19522 * r19523; double r19525 = r19495 + r19524; double r19526 = 0.0694555761; double r19527 = r19526 * r19495; double r19528 = r19527 / r19503; double r19529 = r19525 + r19528; double r19530 = 0.0140005442; double r19531 = r19495 * r19530; double r19532 = r19531 / r19503; double r19533 = r19532 / r19497; double r19534 = r19495 / r19502; double r19535 = r19534 * r19534; double r19536 = 6.0; double r19537 = pow(r19492, r19536); double r19538 = 0.0003579942; double r19539 = r19537 / r19538; double r19540 = r19535 / r19539; double r19541 = r19533 + r19540; double r19542 = r19529 + r19541; double r19543 = 0.7715471019; double r19544 = r19543 / r19497; double r19545 = 4.0; double r19546 = pow(r19492, r19545); double r19547 = 0.0008327945; double r19548 = r19546 / r19547; double r19549 = r19535 / r19548; double r19550 = r19544 + r19549; double r19551 = r19542 + r19550; double r19552 = r19551 * r19492; double r19553 = r19520 / r19552; double r19554 = 1.653196896254347e+29; bool r19555 = r19492 <= r19554; double r19556 = r19492 * r19492; double r19557 = r19556 * (r19556 * r19556); double r19558 = r19530 * r19492; double r19559 = r19492 * r19558; double r19560 = r19559 + r19526; double r19561 = r19557 * r19560; double r19562 = r19556 * r19556; double r19563 = r19562 * r19557; double r19564 = 2.0; double r19565 = r19564 * r19510; double r19566 = r19556 * r19565; double r19567 = r19566 + r19547; double r19568 = r19563 * r19567; double r19569 = r19561 + r19568; double r19570 = r19492 * r19523; double r19571 = r19570 * r19502; double r19572 = r19556 * r19543; double r19573 = r19571 + r19572; double r19574 = r19573 + r19495; double r19575 = r19569 + r19574; double r19576 = r19562 * r19562; double r19577 = r19556 * r19510; double r19578 = r19576 * r19577; double r19579 = r19492 * r19506; double r19580 = r19579 * r19502; double r19581 = r19496 * r19492; double r19582 = r19492 * r19581; double r19583 = r19495 + r19582; double r19584 = r19580 + r19583; double r19585 = r19578 + r19584; double r19586 = r19492 * r19515; double r19587 = r19586 * r19492; double r19588 = r19500 + r19587; double r19589 = r19588 * r19557; double r19590 = r19585 + r19589; double r19591 = r19575 / r19590; double r19592 = r19492 / r19591; double r19593 = pow(r19592, r19495); double r19594 = r19555 ? r19593 : r19553; double r19595 = r19494 ? r19553 : r19594; return r19595; } 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 r19596, r19597, r19598, r19599, r19600, r19601, r19602, r19603, r19604, r19605, r19606, r19607, r19608, r19609, r19610, r19611, r19612, r19613, r19614, r19615, r19616, r19617, r19618, r19619, r19620, r19621, r19622, r19623, r19624, r19625, r19626, r19627, r19628, r19629, r19630, r19631, r19632, r19633, r19634, r19635, r19636, r19637, r19638, r19639; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19596, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19597, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19598); mpfr_init(r19599); mpfr_init(r19600); mpfr_init(r19601); mpfr_init_set_str(r19602, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19603); mpfr_init(r19604); mpfr_init(r19605); mpfr_init_set_str(r19606, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19607); mpfr_init(r19608); mpfr_init(r19609); mpfr_init_set_str(r19610, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19611); mpfr_init(r19612); mpfr_init(r19613); mpfr_init_set_str(r19614, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19615); mpfr_init(r19616); mpfr_init(r19617); mpfr_init_set_str(r19618, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19619); mpfr_init(r19620); mpfr_init_set_str(r19621, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19622); mpfr_init(r19623); mpfr_init_set_str(r19624, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19625); mpfr_init(r19626); mpfr_init_set_str(r19627, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19628); mpfr_init(r19629); mpfr_init_set_str(r19630, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19631); mpfr_init(r19632); mpfr_init_set_str(r19633, "2", 10, MPFR_RNDN); mpfr_init(r19634); mpfr_init(r19635); mpfr_init(r19636); mpfr_init(r19637); mpfr_init(r19638); mpfr_init(r19639); } double f_im(double x) { ; ; mpfr_set_d(r19598, x, MPFR_RNDN); mpfr_mul(r19599, r19598, r19598, MPFR_RNDN); mpfr_mul(r19600, r19597, r19599, MPFR_RNDN); mpfr_add(r19601, r19596, r19600, MPFR_RNDN); ; mpfr_mul(r19603, r19599, r19599, MPFR_RNDN); mpfr_mul(r19604, r19602, r19603, MPFR_RNDN); mpfr_add(r19605, r19601, r19604, MPFR_RNDN); ; mpfr_mul(r19607, r19603, r19599, MPFR_RNDN); mpfr_mul(r19608, r19606, r19607, MPFR_RNDN); mpfr_add(r19609, r19605, r19608, MPFR_RNDN); ; mpfr_mul(r19611, r19607, r19599, MPFR_RNDN); mpfr_mul(r19612, r19610, r19611, MPFR_RNDN); mpfr_add(r19613, r19609, r19612, MPFR_RNDN); ; mpfr_mul(r19615, r19611, r19599, MPFR_RNDN); mpfr_mul(r19616, r19614, r19615, MPFR_RNDN); mpfr_add(r19617, r19613, r19616, MPFR_RNDN); ; mpfr_mul(r19619, r19618, r19599, MPFR_RNDN); mpfr_add(r19620, r19596, r19619, MPFR_RNDN); ; mpfr_mul(r19622, r19621, r19603, MPFR_RNDN); mpfr_add(r19623, r19620, r19622, MPFR_RNDN); ; mpfr_mul(r19625, r19624, r19607, MPFR_RNDN); mpfr_add(r19626, r19623, r19625, MPFR_RNDN); ; mpfr_mul(r19628, r19627, r19611, MPFR_RNDN); mpfr_add(r19629, r19626, r19628, MPFR_RNDN); ; mpfr_mul(r19631, r19630, r19615, MPFR_RNDN); mpfr_add(r19632, r19629, r19631, MPFR_RNDN); ; mpfr_mul(r19634, r19633, r19614, MPFR_RNDN); mpfr_mul(r19635, r19615, r19599, MPFR_RNDN); mpfr_mul(r19636, r19634, r19635, MPFR_RNDN); mpfr_add(r19637, r19632, r19636, MPFR_RNDN); mpfr_div(r19638, r19617, r19637, MPFR_RNDN); mpfr_mul(r19639, r19638, r19598, MPFR_RNDN); return mpfr_get_d(r19639, MPFR_RNDN); } static mpfr_t 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, 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, r19725, r19726, r19727, r19728, r19729, r19730, r19731, r19732, r19733, r19734, r19735, r19736, r19737, r19738, r19739, r19740, r19741, r19742, r19743; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19640); mpfr_init_set_str(r19641, "-1.1392497874101738e+24", 10, MPFR_RNDN); mpfr_init(r19642); mpfr_init_set_str(r19643, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19644, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19645); mpfr_init(r19646); mpfr_init(r19647); mpfr_init_set_str(r19648, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19649); mpfr_init(r19650); mpfr_init(r19651); mpfr_init(r19652); mpfr_init(r19653); mpfr_init_set_str(r19654, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19655); mpfr_init(r19656); mpfr_init(r19657); mpfr_init_set_str(r19658, "0.0001789971", 10, MPFR_RNDN); mpfr_init_set_str(r19659, "10", 10, MPFR_RNDN); mpfr_init(r19660); mpfr_init(r19661); mpfr_init(r19662); mpfr_init_set_str(r19663, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19664); mpfr_init(r19665); mpfr_init(r19666); mpfr_init(r19667); mpfr_init(r19668); mpfr_init(r19669); mpfr_init(r19670); mpfr_init_set_str(r19671, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19672); mpfr_init(r19673); mpfr_init_set_str(r19674, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19675); mpfr_init(r19676); mpfr_init(r19677); mpfr_init_set_str(r19678, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19679); mpfr_init(r19680); mpfr_init(r19681); mpfr_init(r19682); mpfr_init(r19683); mpfr_init_set_str(r19684, "6", 10, MPFR_RNDN); mpfr_init(r19685); mpfr_init_set_str(r19686, "0.0003579942", 10, MPFR_RNDN); mpfr_init(r19687); mpfr_init(r19688); mpfr_init(r19689); mpfr_init(r19690); mpfr_init_set_str(r19691, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19692); mpfr_init_set_str(r19693, "4", 10, MPFR_RNDN); mpfr_init(r19694); mpfr_init_set_str(r19695, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19696); mpfr_init(r19697); mpfr_init(r19698); mpfr_init(r19699); mpfr_init(r19700); mpfr_init(r19701); mpfr_init_set_str(r19702, "1.653196896254347e+29", 10, MPFR_RNDN); mpfr_init(r19703); mpfr_init(r19704); mpfr_init(r19705); mpfr_init(r19706); mpfr_init(r19707); mpfr_init(r19708); mpfr_init(r19709); mpfr_init(r19710); mpfr_init(r19711); mpfr_init_set_str(r19712, "2", 10, MPFR_RNDN); mpfr_init(r19713); mpfr_init(r19714); mpfr_init(r19715); mpfr_init(r19716); mpfr_init(r19717); mpfr_init(r19718); mpfr_init(r19719); mpfr_init(r19720); mpfr_init(r19721); mpfr_init(r19722); mpfr_init(r19723); mpfr_init(r19724); mpfr_init(r19725); mpfr_init(r19726); mpfr_init(r19727); mpfr_init(r19728); mpfr_init(r19729); mpfr_init(r19730); mpfr_init(r19731); mpfr_init(r19732); mpfr_init(r19733); mpfr_init(r19734); mpfr_init(r19735); mpfr_init(r19736); mpfr_init(r19737); mpfr_init(r19738); mpfr_init(r19739); mpfr_init(r19740); mpfr_init(r19741); mpfr_init(r19742); mpfr_init(r19743); } double f_fm(double x) { mpfr_set_d(r19640, x, MPFR_RNDN); ; mpfr_set_si(r19642, mpfr_cmp(r19640, r19641) <= 0, MPFR_RNDN); ; ; mpfr_mul(r19645, r19640, r19640, MPFR_RNDN); mpfr_div(r19646, r19644, r19645, MPFR_RNDN); mpfr_add(r19647, r19643, r19646, MPFR_RNDN); ; mpfr_mul(r19649, r19648, r19643, MPFR_RNDN); mpfr_mul(r19650, r19640, r19640, MPFR_RNDN); mpfr_mul(r19650, r19650, r19640, MPFR_RNDN); mpfr_sqr(r19651, r19650, MPFR_RNDN); mpfr_div(r19652, r19649, r19651, MPFR_RNDN); mpfr_add(r19653, r19647, r19652, MPFR_RNDN); ; mpfr_mul(r19655, r19654, r19643, MPFR_RNDN); mpfr_div(r19656, r19655, r19650, MPFR_RNDN); mpfr_div(r19657, r19656, r19640, MPFR_RNDN); ; ; mpfr_pow(r19660, r19640, r19659, MPFR_RNDN); mpfr_div(r19661, r19658, r19660, MPFR_RNDN); mpfr_add(r19662, r19657, r19661, MPFR_RNDN); ; mpfr_mul(r19664, r19643, r19663, MPFR_RNDN); mpfr_div(r19665, r19664, r19651, MPFR_RNDN); mpfr_div(r19666, r19665, r19645, MPFR_RNDN); mpfr_add(r19667, r19662, r19666, MPFR_RNDN); mpfr_add(r19668, r19653, r19667, MPFR_RNDN); mpfr_div(r19669, r19643, r19640, MPFR_RNDN); mpfr_div(r19670, r19669, r19650, MPFR_RNDN); ; mpfr_mul(r19672, r19670, r19671, MPFR_RNDN); mpfr_add(r19673, r19643, r19672, MPFR_RNDN); ; mpfr_mul(r19675, r19674, r19643, MPFR_RNDN); mpfr_div(r19676, r19675, r19651, MPFR_RNDN); mpfr_add(r19677, r19673, r19676, MPFR_RNDN); ; mpfr_mul(r19679, r19643, r19678, MPFR_RNDN); mpfr_div(r19680, r19679, r19651, MPFR_RNDN); mpfr_div(r19681, r19680, r19645, MPFR_RNDN); mpfr_div(r19682, r19643, r19650, MPFR_RNDN); mpfr_mul(r19683, r19682, r19682, MPFR_RNDN); ; mpfr_pow(r19685, r19640, r19684, MPFR_RNDN); ; mpfr_div(r19687, r19685, r19686, MPFR_RNDN); mpfr_div(r19688, r19683, r19687, MPFR_RNDN); mpfr_add(r19689, r19681, r19688, MPFR_RNDN); mpfr_add(r19690, r19677, r19689, MPFR_RNDN); ; mpfr_div(r19692, r19691, r19645, MPFR_RNDN); ; mpfr_pow(r19694, r19640, r19693, MPFR_RNDN); ; mpfr_div(r19696, r19694, r19695, MPFR_RNDN); mpfr_div(r19697, r19683, r19696, MPFR_RNDN); mpfr_add(r19698, r19692, r19697, MPFR_RNDN); mpfr_add(r19699, r19690, r19698, MPFR_RNDN); mpfr_mul(r19700, r19699, r19640, MPFR_RNDN); mpfr_div(r19701, r19668, r19700, MPFR_RNDN); ; mpfr_set_si(r19703, mpfr_cmp(r19640, r19702) <= 0, MPFR_RNDN); mpfr_sqr(r19704, r19640, MPFR_RNDN); mpfr_mul(r19705, r19704, r19704, MPFR_RNDN); mpfr_mul(r19705, r19705, r19704, MPFR_RNDN); mpfr_mul(r19706, r19678, r19640, MPFR_RNDN); mpfr_mul(r19707, r19640, r19706, MPFR_RNDN); mpfr_add(r19708, r19707, r19674, MPFR_RNDN); mpfr_mul(r19709, r19705, r19708, MPFR_RNDN); mpfr_mul(r19710, r19704, r19704, MPFR_RNDN); mpfr_mul(r19711, r19710, r19705, MPFR_RNDN); ; mpfr_mul(r19713, r19712, r19658, MPFR_RNDN); mpfr_mul(r19714, r19704, r19713, MPFR_RNDN); mpfr_add(r19715, r19714, r19695, MPFR_RNDN); mpfr_mul(r19716, r19711, r19715, MPFR_RNDN); mpfr_add(r19717, r19709, r19716, MPFR_RNDN); mpfr_mul(r19718, r19640, r19671, MPFR_RNDN); mpfr_mul(r19719, r19718, r19650, MPFR_RNDN); mpfr_mul(r19720, r19704, r19691, MPFR_RNDN); mpfr_add(r19721, r19719, r19720, MPFR_RNDN); mpfr_add(r19722, r19721, r19643, MPFR_RNDN); mpfr_add(r19723, r19717, r19722, MPFR_RNDN); mpfr_mul(r19724, r19710, r19710, MPFR_RNDN); mpfr_mul(r19725, r19704, r19658, MPFR_RNDN); mpfr_mul(r19726, r19724, r19725, MPFR_RNDN); mpfr_mul(r19727, r19640, r19654, MPFR_RNDN); mpfr_mul(r19728, r19727, r19650, MPFR_RNDN); mpfr_mul(r19729, r19644, r19640, MPFR_RNDN); mpfr_mul(r19730, r19640, r19729, MPFR_RNDN); mpfr_add(r19731, r19643, r19730, MPFR_RNDN); mpfr_add(r19732, r19728, r19731, MPFR_RNDN); mpfr_add(r19733, r19726, r19732, MPFR_RNDN); mpfr_mul(r19734, r19640, r19663, MPFR_RNDN); mpfr_mul(r19735, r19734, r19640, MPFR_RNDN); mpfr_add(r19736, r19648, r19735, MPFR_RNDN); mpfr_mul(r19737, r19736, r19705, MPFR_RNDN); mpfr_add(r19738, r19733, r19737, MPFR_RNDN); mpfr_div(r19739, r19723, r19738, MPFR_RNDN); mpfr_div(r19740, r19640, r19739, MPFR_RNDN); mpfr_pow(r19741, r19740, r19643, MPFR_RNDN); if (mpfr_get_si(r19703, MPFR_RNDN)) { mpfr_set(r19742, r19741, MPFR_RNDN); } else { mpfr_set(r19742, r19701, MPFR_RNDN); }; if (mpfr_get_si(r19642, MPFR_RNDN)) { mpfr_set(r19743, r19701, MPFR_RNDN); } else { mpfr_set(r19743, r19742, MPFR_RNDN); }; return mpfr_get_d(r19743, MPFR_RNDN); } static mpfr_t 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, r19779, r19780, r19781, r19782, r19783, r19784, r19785, r19786, r19787, r19788, r19789, r19790, r19791, r19792, r19793, r19794, r19795, r19796, r19797, r19798, r19799, r19800, r19801, r19802, r19803, r19804, r19805, r19806, r19807, r19808, r19809, r19810, r19811, r19812, r19813, r19814, r19815, r19816, r19817, r19818, r19819, r19820, r19821, r19822, r19823, r19824, r19825, r19826, r19827, r19828, r19829, r19830, r19831, r19832, r19833, r19834, r19835, r19836, r19837, r19838, r19839, r19840, r19841, r19842, r19843, r19844, r19845, r19846, r19847; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19744); mpfr_init_set_str(r19745, "-1.1392497874101738e+24", 10, MPFR_RNDN); mpfr_init(r19746); mpfr_init_set_str(r19747, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19748, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19749); mpfr_init(r19750); mpfr_init(r19751); mpfr_init_set_str(r19752, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19753); mpfr_init(r19754); mpfr_init(r19755); mpfr_init(r19756); mpfr_init(r19757); mpfr_init_set_str(r19758, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19759); mpfr_init(r19760); mpfr_init(r19761); mpfr_init_set_str(r19762, "0.0001789971", 10, MPFR_RNDN); mpfr_init_set_str(r19763, "10", 10, MPFR_RNDN); mpfr_init(r19764); mpfr_init(r19765); mpfr_init(r19766); mpfr_init_set_str(r19767, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19768); mpfr_init(r19769); mpfr_init(r19770); mpfr_init(r19771); mpfr_init(r19772); mpfr_init(r19773); mpfr_init(r19774); mpfr_init_set_str(r19775, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19776); mpfr_init(r19777); mpfr_init_set_str(r19778, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19779); mpfr_init(r19780); mpfr_init(r19781); mpfr_init_set_str(r19782, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19783); mpfr_init(r19784); mpfr_init(r19785); mpfr_init(r19786); mpfr_init(r19787); mpfr_init_set_str(r19788, "6", 10, MPFR_RNDN); mpfr_init(r19789); mpfr_init_set_str(r19790, "0.0003579942", 10, MPFR_RNDN); mpfr_init(r19791); mpfr_init(r19792); mpfr_init(r19793); mpfr_init(r19794); mpfr_init_set_str(r19795, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19796); mpfr_init_set_str(r19797, "4", 10, MPFR_RNDN); mpfr_init(r19798); mpfr_init_set_str(r19799, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19800); mpfr_init(r19801); mpfr_init(r19802); mpfr_init(r19803); mpfr_init(r19804); mpfr_init(r19805); mpfr_init_set_str(r19806, "1.653196896254347e+29", 10, MPFR_RNDN); mpfr_init(r19807); mpfr_init(r19808); mpfr_init(r19809); mpfr_init(r19810); mpfr_init(r19811); mpfr_init(r19812); mpfr_init(r19813); mpfr_init(r19814); mpfr_init(r19815); mpfr_init_set_str(r19816, "2", 10, MPFR_RNDN); mpfr_init(r19817); mpfr_init(r19818); mpfr_init(r19819); mpfr_init(r19820); mpfr_init(r19821); mpfr_init(r19822); mpfr_init(r19823); mpfr_init(r19824); mpfr_init(r19825); mpfr_init(r19826); mpfr_init(r19827); mpfr_init(r19828); mpfr_init(r19829); mpfr_init(r19830); mpfr_init(r19831); mpfr_init(r19832); mpfr_init(r19833); mpfr_init(r19834); mpfr_init(r19835); mpfr_init(r19836); mpfr_init(r19837); mpfr_init(r19838); mpfr_init(r19839); mpfr_init(r19840); mpfr_init(r19841); mpfr_init(r19842); mpfr_init(r19843); mpfr_init(r19844); mpfr_init(r19845); mpfr_init(r19846); mpfr_init(r19847); } double f_dm(double x) { mpfr_set_d(r19744, x, MPFR_RNDN); ; mpfr_set_si(r19746, mpfr_cmp(r19744, r19745) <= 0, MPFR_RNDN); ; ; mpfr_mul(r19749, r19744, r19744, MPFR_RNDN); mpfr_div(r19750, r19748, r19749, MPFR_RNDN); mpfr_add(r19751, r19747, r19750, MPFR_RNDN); ; mpfr_mul(r19753, r19752, r19747, MPFR_RNDN); mpfr_mul(r19754, r19744, r19744, MPFR_RNDN); mpfr_mul(r19754, r19754, r19744, MPFR_RNDN); mpfr_sqr(r19755, r19754, MPFR_RNDN); mpfr_div(r19756, r19753, r19755, MPFR_RNDN); mpfr_add(r19757, r19751, r19756, MPFR_RNDN); ; mpfr_mul(r19759, r19758, r19747, MPFR_RNDN); mpfr_div(r19760, r19759, r19754, MPFR_RNDN); mpfr_div(r19761, r19760, r19744, MPFR_RNDN); ; ; mpfr_pow(r19764, r19744, r19763, MPFR_RNDN); mpfr_div(r19765, r19762, r19764, MPFR_RNDN); mpfr_add(r19766, r19761, r19765, MPFR_RNDN); ; mpfr_mul(r19768, r19747, r19767, MPFR_RNDN); mpfr_div(r19769, r19768, r19755, MPFR_RNDN); mpfr_div(r19770, r19769, r19749, MPFR_RNDN); mpfr_add(r19771, r19766, r19770, MPFR_RNDN); mpfr_add(r19772, r19757, r19771, MPFR_RNDN); mpfr_div(r19773, r19747, r19744, MPFR_RNDN); mpfr_div(r19774, r19773, r19754, MPFR_RNDN); ; mpfr_mul(r19776, r19774, r19775, MPFR_RNDN); mpfr_add(r19777, r19747, r19776, MPFR_RNDN); ; mpfr_mul(r19779, r19778, r19747, MPFR_RNDN); mpfr_div(r19780, r19779, r19755, MPFR_RNDN); mpfr_add(r19781, r19777, r19780, MPFR_RNDN); ; mpfr_mul(r19783, r19747, r19782, MPFR_RNDN); mpfr_div(r19784, r19783, r19755, MPFR_RNDN); mpfr_div(r19785, r19784, r19749, MPFR_RNDN); mpfr_div(r19786, r19747, r19754, MPFR_RNDN); mpfr_mul(r19787, r19786, r19786, MPFR_RNDN); ; mpfr_pow(r19789, r19744, r19788, MPFR_RNDN); ; mpfr_div(r19791, r19789, r19790, MPFR_RNDN); mpfr_div(r19792, r19787, r19791, MPFR_RNDN); mpfr_add(r19793, r19785, r19792, MPFR_RNDN); mpfr_add(r19794, r19781, r19793, MPFR_RNDN); ; mpfr_div(r19796, r19795, r19749, MPFR_RNDN); ; mpfr_pow(r19798, r19744, r19797, MPFR_RNDN); ; mpfr_div(r19800, r19798, r19799, MPFR_RNDN); mpfr_div(r19801, r19787, r19800, MPFR_RNDN); mpfr_add(r19802, r19796, r19801, MPFR_RNDN); mpfr_add(r19803, r19794, r19802, MPFR_RNDN); mpfr_mul(r19804, r19803, r19744, MPFR_RNDN); mpfr_div(r19805, r19772, r19804, MPFR_RNDN); ; mpfr_set_si(r19807, mpfr_cmp(r19744, r19806) <= 0, MPFR_RNDN); mpfr_sqr(r19808, r19744, MPFR_RNDN); mpfr_mul(r19809, r19808, r19808, MPFR_RNDN); mpfr_mul(r19809, r19809, r19808, MPFR_RNDN); mpfr_mul(r19810, r19782, r19744, MPFR_RNDN); mpfr_mul(r19811, r19744, r19810, MPFR_RNDN); mpfr_add(r19812, r19811, r19778, MPFR_RNDN); mpfr_mul(r19813, r19809, r19812, MPFR_RNDN); mpfr_mul(r19814, r19808, r19808, MPFR_RNDN); mpfr_mul(r19815, r19814, r19809, MPFR_RNDN); ; mpfr_mul(r19817, r19816, r19762, MPFR_RNDN); mpfr_mul(r19818, r19808, r19817, MPFR_RNDN); mpfr_add(r19819, r19818, r19799, MPFR_RNDN); mpfr_mul(r19820, r19815, r19819, MPFR_RNDN); mpfr_add(r19821, r19813, r19820, MPFR_RNDN); mpfr_mul(r19822, r19744, r19775, MPFR_RNDN); mpfr_mul(r19823, r19822, r19754, MPFR_RNDN); mpfr_mul(r19824, r19808, r19795, MPFR_RNDN); mpfr_add(r19825, r19823, r19824, MPFR_RNDN); mpfr_add(r19826, r19825, r19747, MPFR_RNDN); mpfr_add(r19827, r19821, r19826, MPFR_RNDN); mpfr_mul(r19828, r19814, r19814, MPFR_RNDN); mpfr_mul(r19829, r19808, r19762, MPFR_RNDN); mpfr_mul(r19830, r19828, r19829, MPFR_RNDN); mpfr_mul(r19831, r19744, r19758, MPFR_RNDN); mpfr_mul(r19832, r19831, r19754, MPFR_RNDN); mpfr_mul(r19833, r19748, r19744, MPFR_RNDN); mpfr_mul(r19834, r19744, r19833, MPFR_RNDN); mpfr_add(r19835, r19747, r19834, MPFR_RNDN); mpfr_add(r19836, r19832, r19835, MPFR_RNDN); mpfr_add(r19837, r19830, r19836, MPFR_RNDN); mpfr_mul(r19838, r19744, r19767, MPFR_RNDN); mpfr_mul(r19839, r19838, r19744, MPFR_RNDN); mpfr_add(r19840, r19752, r19839, MPFR_RNDN); mpfr_mul(r19841, r19840, r19809, MPFR_RNDN); mpfr_add(r19842, r19837, r19841, MPFR_RNDN); mpfr_div(r19843, r19827, r19842, MPFR_RNDN); mpfr_div(r19844, r19744, r19843, MPFR_RNDN); mpfr_pow(r19845, r19844, r19747, MPFR_RNDN); if (mpfr_get_si(r19807, MPFR_RNDN)) { mpfr_set(r19846, r19845, MPFR_RNDN); } else { mpfr_set(r19846, r19805, MPFR_RNDN); }; if (mpfr_get_si(r19746, MPFR_RNDN)) { mpfr_set(r19847, r19805, MPFR_RNDN); } else { mpfr_set(r19847, r19846, MPFR_RNDN); }; return mpfr_get_d(r19847, MPFR_RNDN); }