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