#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r19271 = 1.0f; float r19272 = 0.1049934947f; float r19273 = x; float r19274 = r19273 * r19273; float r19275 = r19272 * r19274; float r19276 = r19271 + r19275; float r19277 = 0.0424060604f; float r19278 = r19274 * r19274; float r19279 = r19277 * r19278; float r19280 = r19276 + r19279; float r19281 = 0.0072644182f; float r19282 = r19278 * r19274; float r19283 = r19281 * r19282; float r19284 = r19280 + r19283; float r19285 = 0.0005064034f; float r19286 = r19282 * r19274; float r19287 = r19285 * r19286; float r19288 = r19284 + r19287; float r19289 = 0.0001789971f; float r19290 = r19286 * r19274; float r19291 = r19289 * r19290; float r19292 = r19288 + r19291; float r19293 = 0.7715471019f; float r19294 = r19293 * r19274; float r19295 = r19271 + r19294; float r19296 = 0.2909738639f; float r19297 = r19296 * r19278; float r19298 = r19295 + r19297; float r19299 = 0.0694555761f; float r19300 = r19299 * r19282; float r19301 = r19298 + r19300; float r19302 = 0.0140005442f; float r19303 = r19302 * r19286; float r19304 = r19301 + r19303; float r19305 = 0.0008327945f; float r19306 = r19305 * r19290; float r19307 = r19304 + r19306; float r19308 = 2.0f; float r19309 = r19308 * r19289; float r19310 = r19290 * r19274; float r19311 = r19309 * r19310; float r19312 = r19307 + r19311; float r19313 = r19292 / r19312; float r19314 = r19313 * r19273; return r19314; } double f_id(double x) { double r19315 = 1.0; double r19316 = 0.1049934947; double r19317 = x; double r19318 = r19317 * r19317; double r19319 = r19316 * r19318; double r19320 = r19315 + r19319; double r19321 = 0.0424060604; double r19322 = r19318 * r19318; double r19323 = r19321 * r19322; double r19324 = r19320 + r19323; double r19325 = 0.0072644182; double r19326 = r19322 * r19318; double r19327 = r19325 * r19326; double r19328 = r19324 + r19327; double r19329 = 0.0005064034; double r19330 = r19326 * r19318; double r19331 = r19329 * r19330; double r19332 = r19328 + r19331; double r19333 = 0.0001789971; double r19334 = r19330 * r19318; double r19335 = r19333 * r19334; double r19336 = r19332 + r19335; double r19337 = 0.7715471019; double r19338 = r19337 * r19318; double r19339 = r19315 + r19338; double r19340 = 0.2909738639; double r19341 = r19340 * r19322; double r19342 = r19339 + r19341; double r19343 = 0.0694555761; double r19344 = r19343 * r19326; double r19345 = r19342 + r19344; double r19346 = 0.0140005442; double r19347 = r19346 * r19330; double r19348 = r19345 + r19347; double r19349 = 0.0008327945; double r19350 = r19349 * r19334; double r19351 = r19348 + r19350; double r19352 = 2.0; double r19353 = r19352 * r19333; double r19354 = r19334 * r19318; double r19355 = r19353 * r19354; double r19356 = r19351 + r19355; double r19357 = r19336 / r19356; double r19358 = r19357 * r19317; return r19358; } double f_of(float x) { float r19359 = x; float r19360 = r19359 * r19359; float r19361 = 0.0001789971f; float r19362 = r19360 * r19361; float r19363 = r19360 * r19360; float r19364 = r19363 * r19363; float r19365 = r19362 * r19364; float r19366 = 0.0424060604f; float r19367 = r19359 * r19366; float r19368 = r19359 * (r19359 * r19359); float r19369 = r19367 * r19368; float r19370 = 1.0f; float r19371 = r19369 + r19370; float r19372 = 0.1049934947f; float r19373 = r19372 * r19360; float r19374 = r19371 + r19373; float r19375 = r19365 + r19374; float r19376 = r19359 * r19375; float r19377 = 0.7715471019f; float r19378 = r19359 * r19377; float r19379 = r19378 * r19359; float r19380 = r19370 + r19379; float r19381 = 0.2909738639f; float r19382 = r19359 * r19381; float r19383 = r19368 * r19382; float r19384 = r19380 + r19383; float r19385 = r19368 * r19368; float r19386 = 0.0694555761f; float r19387 = 0.0140005442f; float r19388 = r19387 * r19359; float r19389 = r19359 * r19388; float r19390 = r19386 + r19389; float r19391 = r19385 * r19390; float r19392 = r19385 * r19363; float r19393 = 2.0f; float r19394 = r19393 * r19361; float r19395 = r19394 * r19360; float r19396 = 0.0008327945f; float r19397 = r19395 + r19396; float r19398 = r19392 * r19397; float r19399 = r19391 + r19398; float r19400 = r19384 + r19399; float r19401 = r19376 / r19400; float r19402 = r19359 * r19385; float r19403 = 0.0005064034f; float r19404 = r19403 * r19360; float r19405 = 0.0072644182f; float r19406 = r19404 + r19405; float r19407 = r19402 * r19406; float r19408 = r19360 * (r19360 * r19360); float r19409 = r19363 * r19408; float r19410 = r19360 * r19394; float r19411 = r19410 + r19396; float r19412 = r19409 * r19411; float r19413 = r19377 * r19359; float r19414 = r19359 * r19413; float r19415 = r19382 * r19368; float r19416 = r19415 + r19370; float r19417 = r19414 + r19416; float r19418 = r19412 + r19417; float r19419 = r19389 + r19386; float r19420 = r19408 * r19419; float r19421 = r19418 + r19420; float r19422 = r19407 / r19421; float r19423 = r19401 + r19422; return r19423; } double f_od(double x) { double r19424 = x; double r19425 = r19424 * r19424; double r19426 = 0.0001789971; double r19427 = r19425 * r19426; double r19428 = r19425 * r19425; double r19429 = r19428 * r19428; double r19430 = r19427 * r19429; double r19431 = 0.0424060604; double r19432 = r19424 * r19431; double r19433 = r19424 * (r19424 * r19424); double r19434 = r19432 * r19433; double r19435 = 1.0; double r19436 = r19434 + r19435; double r19437 = 0.1049934947; double r19438 = r19437 * r19425; double r19439 = r19436 + r19438; double r19440 = r19430 + r19439; double r19441 = r19424 * r19440; double r19442 = 0.7715471019; double r19443 = r19424 * r19442; double r19444 = r19443 * r19424; double r19445 = r19435 + r19444; double r19446 = 0.2909738639; double r19447 = r19424 * r19446; double r19448 = r19433 * r19447; double r19449 = r19445 + r19448; double r19450 = r19433 * r19433; double r19451 = 0.0694555761; double r19452 = 0.0140005442; double r19453 = r19452 * r19424; double r19454 = r19424 * r19453; double r19455 = r19451 + r19454; double r19456 = r19450 * r19455; double r19457 = r19450 * r19428; double r19458 = 2.0; double r19459 = r19458 * r19426; double r19460 = r19459 * r19425; double r19461 = 0.0008327945; double r19462 = r19460 + r19461; double r19463 = r19457 * r19462; double r19464 = r19456 + r19463; double r19465 = r19449 + r19464; double r19466 = r19441 / r19465; double r19467 = r19424 * r19450; double r19468 = 0.0005064034; double r19469 = r19468 * r19425; double r19470 = 0.0072644182; double r19471 = r19469 + r19470; double r19472 = r19467 * r19471; double r19473 = r19425 * (r19425 * r19425); double r19474 = r19428 * r19473; double r19475 = r19425 * r19459; double r19476 = r19475 + r19461; double r19477 = r19474 * r19476; double r19478 = r19442 * r19424; double r19479 = r19424 * r19478; double r19480 = r19447 * r19433; double r19481 = r19480 + r19435; double r19482 = r19479 + r19481; double r19483 = r19477 + r19482; double r19484 = r19454 + r19451; double r19485 = r19473 * r19484; double r19486 = r19483 + r19485; double r19487 = r19472 / r19486; double r19488 = r19466 + r19487; return r19488; } 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 r19489, r19490, r19491, r19492, r19493, r19494, r19495, r19496, r19497, r19498, r19499, r19500, r19501, r19502, r19503, r19504, r19505, r19506, r19507, r19508, r19509, r19510, r19511, r19512, r19513, r19514, r19515, r19516, r19517, r19518, r19519, r19520, r19521, r19522, r19523, r19524, r19525, r19526, r19527, r19528, r19529, r19530, r19531, r19532; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19489, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19490, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19491); mpfr_init(r19492); mpfr_init(r19493); mpfr_init(r19494); mpfr_init_set_str(r19495, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19496); mpfr_init(r19497); mpfr_init(r19498); mpfr_init_set_str(r19499, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19500); mpfr_init(r19501); mpfr_init(r19502); mpfr_init_set_str(r19503, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19504); mpfr_init(r19505); mpfr_init(r19506); mpfr_init_set_str(r19507, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19508); mpfr_init(r19509); mpfr_init(r19510); mpfr_init_set_str(r19511, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19512); mpfr_init(r19513); mpfr_init_set_str(r19514, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19515); mpfr_init(r19516); mpfr_init_set_str(r19517, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19518); mpfr_init(r19519); mpfr_init_set_str(r19520, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19521); mpfr_init(r19522); mpfr_init_set_str(r19523, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19524); mpfr_init(r19525); mpfr_init_set_str(r19526, "2", 10, MPFR_RNDN); mpfr_init(r19527); mpfr_init(r19528); mpfr_init(r19529); mpfr_init(r19530); mpfr_init(r19531); mpfr_init(r19532); } double f_im(double x) { ; ; mpfr_set_d(r19491, x, MPFR_RNDN); mpfr_mul(r19492, r19491, r19491, MPFR_RNDN); mpfr_mul(r19493, r19490, r19492, MPFR_RNDN); mpfr_add(r19494, r19489, r19493, MPFR_RNDN); ; mpfr_mul(r19496, r19492, r19492, MPFR_RNDN); mpfr_mul(r19497, r19495, r19496, MPFR_RNDN); mpfr_add(r19498, r19494, r19497, MPFR_RNDN); ; mpfr_mul(r19500, r19496, r19492, MPFR_RNDN); mpfr_mul(r19501, r19499, r19500, MPFR_RNDN); mpfr_add(r19502, r19498, r19501, MPFR_RNDN); ; mpfr_mul(r19504, r19500, r19492, MPFR_RNDN); mpfr_mul(r19505, r19503, r19504, MPFR_RNDN); mpfr_add(r19506, r19502, r19505, MPFR_RNDN); ; mpfr_mul(r19508, r19504, r19492, MPFR_RNDN); mpfr_mul(r19509, r19507, r19508, MPFR_RNDN); mpfr_add(r19510, r19506, r19509, MPFR_RNDN); ; mpfr_mul(r19512, r19511, r19492, MPFR_RNDN); mpfr_add(r19513, r19489, r19512, MPFR_RNDN); ; mpfr_mul(r19515, r19514, r19496, MPFR_RNDN); mpfr_add(r19516, r19513, r19515, MPFR_RNDN); ; mpfr_mul(r19518, r19517, r19500, MPFR_RNDN); mpfr_add(r19519, r19516, r19518, MPFR_RNDN); ; mpfr_mul(r19521, r19520, r19504, MPFR_RNDN); mpfr_add(r19522, r19519, r19521, MPFR_RNDN); ; mpfr_mul(r19524, r19523, r19508, MPFR_RNDN); mpfr_add(r19525, r19522, r19524, MPFR_RNDN); ; mpfr_mul(r19527, r19526, r19507, MPFR_RNDN); mpfr_mul(r19528, r19508, r19492, MPFR_RNDN); mpfr_mul(r19529, r19527, r19528, MPFR_RNDN); mpfr_add(r19530, r19525, r19529, MPFR_RNDN); mpfr_div(r19531, r19510, r19530, MPFR_RNDN); mpfr_mul(r19532, r19531, r19491, MPFR_RNDN); return mpfr_get_d(r19532, MPFR_RNDN); } static mpfr_t r19533, r19534, r19535, r19536, r19537, r19538, r19539, r19540, 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, r19585, r19586, r19587, r19588, r19589, r19590, r19591, r19592, r19593, r19594, r19595, r19596, r19597; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19533); mpfr_init(r19534); mpfr_init_set_str(r19535, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19536); mpfr_init(r19537); mpfr_init(r19538); mpfr_init(r19539); mpfr_init_set_str(r19540, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19541); mpfr_init(r19542); mpfr_init(r19543); mpfr_init_set_str(r19544, "1", 10, MPFR_RNDN); mpfr_init(r19545); mpfr_init_set_str(r19546, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19547); mpfr_init(r19548); mpfr_init(r19549); mpfr_init(r19550); mpfr_init_set_str(r19551, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19552); mpfr_init(r19553); mpfr_init(r19554); mpfr_init_set_str(r19555, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19556); mpfr_init(r19557); mpfr_init(r19558); mpfr_init(r19559); mpfr_init_set_str(r19560, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19561, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19562); mpfr_init(r19563); mpfr_init(r19564); mpfr_init(r19565); mpfr_init(r19566); mpfr_init_set_str(r19567, "2", 10, MPFR_RNDN); mpfr_init(r19568); mpfr_init(r19569); mpfr_init_set_str(r19570, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19571); mpfr_init(r19572); mpfr_init(r19573); mpfr_init(r19574); mpfr_init(r19575); mpfr_init(r19576); mpfr_init_set_str(r19577, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19578); mpfr_init_set_str(r19579, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19580); mpfr_init(r19581); mpfr_init(r19582); mpfr_init(r19583); mpfr_init(r19584); mpfr_init(r19585); mpfr_init(r19586); mpfr_init(r19587); mpfr_init(r19588); mpfr_init(r19589); mpfr_init(r19590); mpfr_init(r19591); mpfr_init(r19592); mpfr_init(r19593); mpfr_init(r19594); mpfr_init(r19595); mpfr_init(r19596); mpfr_init(r19597); } double f_fm(double x) { mpfr_set_d(r19533, x, MPFR_RNDN); mpfr_sqr(r19534, r19533, MPFR_RNDN); ; mpfr_mul(r19536, r19534, r19535, MPFR_RNDN); mpfr_mul(r19537, r19534, r19534, MPFR_RNDN); mpfr_mul(r19538, r19537, r19537, MPFR_RNDN); mpfr_mul(r19539, r19536, r19538, MPFR_RNDN); ; mpfr_mul(r19541, r19533, r19540, MPFR_RNDN); mpfr_mul(r19542, r19533, r19533, MPFR_RNDN); mpfr_mul(r19542, r19542, r19533, MPFR_RNDN); mpfr_mul(r19543, r19541, r19542, MPFR_RNDN); ; mpfr_add(r19545, r19543, r19544, MPFR_RNDN); ; mpfr_mul(r19547, r19546, r19534, MPFR_RNDN); mpfr_add(r19548, r19545, r19547, MPFR_RNDN); mpfr_add(r19549, r19539, r19548, MPFR_RNDN); mpfr_mul(r19550, r19533, r19549, MPFR_RNDN); ; mpfr_mul(r19552, r19533, r19551, MPFR_RNDN); mpfr_mul(r19553, r19552, r19533, MPFR_RNDN); mpfr_add(r19554, r19544, r19553, MPFR_RNDN); ; mpfr_mul(r19556, r19533, r19555, MPFR_RNDN); mpfr_mul(r19557, r19542, r19556, MPFR_RNDN); mpfr_add(r19558, r19554, r19557, MPFR_RNDN); mpfr_mul(r19559, r19542, r19542, MPFR_RNDN); ; ; mpfr_mul(r19562, r19561, r19533, MPFR_RNDN); mpfr_mul(r19563, r19533, r19562, MPFR_RNDN); mpfr_add(r19564, r19560, r19563, MPFR_RNDN); mpfr_mul(r19565, r19559, r19564, MPFR_RNDN); mpfr_mul(r19566, r19559, r19537, MPFR_RNDN); ; mpfr_mul(r19568, r19567, r19535, MPFR_RNDN); mpfr_mul(r19569, r19568, r19534, MPFR_RNDN); ; mpfr_add(r19571, r19569, r19570, MPFR_RNDN); mpfr_mul(r19572, r19566, r19571, MPFR_RNDN); mpfr_add(r19573, r19565, r19572, MPFR_RNDN); mpfr_add(r19574, r19558, r19573, MPFR_RNDN); mpfr_div(r19575, r19550, r19574, MPFR_RNDN); mpfr_mul(r19576, r19533, r19559, MPFR_RNDN); ; mpfr_mul(r19578, r19577, r19534, MPFR_RNDN); ; mpfr_add(r19580, r19578, r19579, MPFR_RNDN); mpfr_mul(r19581, r19576, r19580, MPFR_RNDN); mpfr_mul(r19582, r19534, r19534, MPFR_RNDN); mpfr_mul(r19582, r19582, r19534, MPFR_RNDN); mpfr_mul(r19583, r19537, r19582, MPFR_RNDN); mpfr_mul(r19584, r19534, r19568, MPFR_RNDN); mpfr_add(r19585, r19584, r19570, MPFR_RNDN); mpfr_mul(r19586, r19583, r19585, MPFR_RNDN); mpfr_mul(r19587, r19551, r19533, MPFR_RNDN); mpfr_mul(r19588, r19533, r19587, MPFR_RNDN); mpfr_mul(r19589, r19556, r19542, MPFR_RNDN); mpfr_add(r19590, r19589, r19544, MPFR_RNDN); mpfr_add(r19591, r19588, r19590, MPFR_RNDN); mpfr_add(r19592, r19586, r19591, MPFR_RNDN); mpfr_add(r19593, r19563, r19560, MPFR_RNDN); mpfr_mul(r19594, r19582, r19593, MPFR_RNDN); mpfr_add(r19595, r19592, r19594, MPFR_RNDN); mpfr_div(r19596, r19581, r19595, MPFR_RNDN); mpfr_add(r19597, r19575, r19596, MPFR_RNDN); return mpfr_get_d(r19597, MPFR_RNDN); } static mpfr_t 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; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19598); mpfr_init(r19599); mpfr_init_set_str(r19600, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19601); mpfr_init(r19602); mpfr_init(r19603); mpfr_init(r19604); mpfr_init_set_str(r19605, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19606); mpfr_init(r19607); mpfr_init(r19608); mpfr_init_set_str(r19609, "1", 10, MPFR_RNDN); mpfr_init(r19610); mpfr_init_set_str(r19611, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19612); mpfr_init(r19613); mpfr_init(r19614); mpfr_init(r19615); mpfr_init_set_str(r19616, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19617); mpfr_init(r19618); mpfr_init(r19619); mpfr_init_set_str(r19620, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19621); mpfr_init(r19622); mpfr_init(r19623); mpfr_init(r19624); mpfr_init_set_str(r19625, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19626, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19627); mpfr_init(r19628); mpfr_init(r19629); mpfr_init(r19630); mpfr_init(r19631); mpfr_init_set_str(r19632, "2", 10, MPFR_RNDN); mpfr_init(r19633); mpfr_init(r19634); mpfr_init_set_str(r19635, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19636); mpfr_init(r19637); mpfr_init(r19638); mpfr_init(r19639); mpfr_init(r19640); mpfr_init(r19641); mpfr_init_set_str(r19642, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19643); mpfr_init_set_str(r19644, "0.0072644182", 10, MPFR_RNDN); 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); } double f_dm(double x) { mpfr_set_d(r19598, x, MPFR_RNDN); mpfr_sqr(r19599, r19598, MPFR_RNDN); ; mpfr_mul(r19601, r19599, r19600, MPFR_RNDN); mpfr_mul(r19602, r19599, r19599, MPFR_RNDN); mpfr_mul(r19603, r19602, r19602, MPFR_RNDN); mpfr_mul(r19604, r19601, r19603, MPFR_RNDN); ; mpfr_mul(r19606, r19598, r19605, MPFR_RNDN); mpfr_mul(r19607, r19598, r19598, MPFR_RNDN); mpfr_mul(r19607, r19607, r19598, MPFR_RNDN); mpfr_mul(r19608, r19606, r19607, MPFR_RNDN); ; mpfr_add(r19610, r19608, r19609, MPFR_RNDN); ; mpfr_mul(r19612, r19611, r19599, MPFR_RNDN); mpfr_add(r19613, r19610, r19612, MPFR_RNDN); mpfr_add(r19614, r19604, r19613, MPFR_RNDN); mpfr_mul(r19615, r19598, r19614, MPFR_RNDN); ; mpfr_mul(r19617, r19598, r19616, MPFR_RNDN); mpfr_mul(r19618, r19617, r19598, MPFR_RNDN); mpfr_add(r19619, r19609, r19618, MPFR_RNDN); ; mpfr_mul(r19621, r19598, r19620, MPFR_RNDN); mpfr_mul(r19622, r19607, r19621, MPFR_RNDN); mpfr_add(r19623, r19619, r19622, MPFR_RNDN); mpfr_mul(r19624, r19607, r19607, MPFR_RNDN); ; ; mpfr_mul(r19627, r19626, r19598, MPFR_RNDN); mpfr_mul(r19628, r19598, r19627, MPFR_RNDN); mpfr_add(r19629, r19625, r19628, MPFR_RNDN); mpfr_mul(r19630, r19624, r19629, MPFR_RNDN); mpfr_mul(r19631, r19624, r19602, MPFR_RNDN); ; mpfr_mul(r19633, r19632, r19600, MPFR_RNDN); mpfr_mul(r19634, r19633, r19599, MPFR_RNDN); ; mpfr_add(r19636, r19634, r19635, MPFR_RNDN); mpfr_mul(r19637, r19631, r19636, MPFR_RNDN); mpfr_add(r19638, r19630, r19637, MPFR_RNDN); mpfr_add(r19639, r19623, r19638, MPFR_RNDN); mpfr_div(r19640, r19615, r19639, MPFR_RNDN); mpfr_mul(r19641, r19598, r19624, MPFR_RNDN); ; mpfr_mul(r19643, r19642, r19599, MPFR_RNDN); ; mpfr_add(r19645, r19643, r19644, MPFR_RNDN); mpfr_mul(r19646, r19641, r19645, MPFR_RNDN); mpfr_mul(r19647, r19599, r19599, MPFR_RNDN); mpfr_mul(r19647, r19647, r19599, MPFR_RNDN); mpfr_mul(r19648, r19602, r19647, MPFR_RNDN); mpfr_mul(r19649, r19599, r19633, MPFR_RNDN); mpfr_add(r19650, r19649, r19635, MPFR_RNDN); mpfr_mul(r19651, r19648, r19650, MPFR_RNDN); mpfr_mul(r19652, r19616, r19598, MPFR_RNDN); mpfr_mul(r19653, r19598, r19652, MPFR_RNDN); mpfr_mul(r19654, r19621, r19607, MPFR_RNDN); mpfr_add(r19655, r19654, r19609, MPFR_RNDN); mpfr_add(r19656, r19653, r19655, MPFR_RNDN); mpfr_add(r19657, r19651, r19656, MPFR_RNDN); mpfr_add(r19658, r19628, r19625, MPFR_RNDN); mpfr_mul(r19659, r19647, r19658, MPFR_RNDN); mpfr_add(r19660, r19657, r19659, MPFR_RNDN); mpfr_div(r19661, r19646, r19660, MPFR_RNDN); mpfr_add(r19662, r19640, r19661, MPFR_RNDN); return mpfr_get_d(r19662, MPFR_RNDN); }