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