#include #include #include #include #include char *name = "Jmat.Real.gamma, branch z less than 0.5"; double f_if(float z) { float r26245 = atan2(1.0, 0.0); float r26246 = z; float r26247 = r26245 * r26246; float r26248 = sin(r26247); float r26249 = r26245 / r26248; float r26250 = 2; float r26251 = r26245 * r26250; float r26252 = sqrt(r26251); float r26253 = 1; float r26254 = r26253 - r26246; float r26255 = r26254 - r26253; float r26256 = 7; float r26257 = r26255 + r26256; float r26258 = 0.5; float r26259 = r26257 + r26258; float r26260 = r26255 + r26258; float r26261 = pow(r26259, r26260); float r26262 = r26252 * r26261; float r26263 = -r26259; float r26264 = exp(r26263); float r26265 = r26262 * r26264; float r26266 = 0.9999999999998099; float r26267 = 676.5203681218851; float r26268 = r26255 + r26253; float r26269 = r26267 / r26268; float r26270 = r26266 + r26269; float r26271 = -1259.1392167224028; float r26272 = r26255 + r26250; float r26273 = r26271 / r26272; float r26274 = r26270 + r26273; float r26275 = 771.3234287776531; float r26276 = 3; float r26277 = r26255 + r26276; float r26278 = r26275 / r26277; float r26279 = r26274 + r26278; float r26280 = -176.6150291621406; float r26281 = 4; float r26282 = r26255 + r26281; float r26283 = r26280 / r26282; float r26284 = r26279 + r26283; float r26285 = 12.507343278686905; float r26286 = 5; float r26287 = r26255 + r26286; float r26288 = r26285 / r26287; float r26289 = r26284 + r26288; float r26290 = -0.13857109526572012; float r26291 = 6; float r26292 = r26255 + r26291; float r26293 = r26290 / r26292; float r26294 = r26289 + r26293; float r26295 = 9.984369578019572e-06; float r26296 = r26295 / r26257; float r26297 = r26294 + r26296; float r26298 = 1.5056327351493116e-07; float r26299 = 8; float r26300 = r26255 + r26299; float r26301 = r26298 / r26300; float r26302 = r26297 + r26301; float r26303 = r26265 * r26302; float r26304 = r26249 * r26303; return r26304; } double f_id(double z) { double r26305 = atan2(1.0, 0.0); double r26306 = z; double r26307 = r26305 * r26306; double r26308 = sin(r26307); double r26309 = r26305 / r26308; double r26310 = 2; double r26311 = r26305 * r26310; double r26312 = sqrt(r26311); double r26313 = 1; double r26314 = r26313 - r26306; double r26315 = r26314 - r26313; double r26316 = 7; double r26317 = r26315 + r26316; double r26318 = 0.5; double r26319 = r26317 + r26318; double r26320 = r26315 + r26318; double r26321 = pow(r26319, r26320); double r26322 = r26312 * r26321; double r26323 = -r26319; double r26324 = exp(r26323); double r26325 = r26322 * r26324; double r26326 = 0.9999999999998099; double r26327 = 676.5203681218851; double r26328 = r26315 + r26313; double r26329 = r26327 / r26328; double r26330 = r26326 + r26329; double r26331 = -1259.1392167224028; double r26332 = r26315 + r26310; double r26333 = r26331 / r26332; double r26334 = r26330 + r26333; double r26335 = 771.3234287776531; double r26336 = 3; double r26337 = r26315 + r26336; double r26338 = r26335 / r26337; double r26339 = r26334 + r26338; double r26340 = -176.6150291621406; double r26341 = 4; double r26342 = r26315 + r26341; double r26343 = r26340 / r26342; double r26344 = r26339 + r26343; double r26345 = 12.507343278686905; double r26346 = 5; double r26347 = r26315 + r26346; double r26348 = r26345 / r26347; double r26349 = r26344 + r26348; double r26350 = -0.13857109526572012; double r26351 = 6; double r26352 = r26315 + r26351; double r26353 = r26350 / r26352; double r26354 = r26349 + r26353; double r26355 = 9.984369578019572e-06; double r26356 = r26355 / r26317; double r26357 = r26354 + r26356; double r26358 = 1.5056327351493116e-07; double r26359 = 8; double r26360 = r26315 + r26359; double r26361 = r26358 / r26360; double r26362 = r26357 + r26361; double r26363 = r26325 * r26362; double r26364 = r26309 * r26363; return r26364; } double f_of(float z) { float r26365 = -176.6150291621406; float r26366 = 1; float r26367 = z; float r26368 = r26366 - r26367; float r26369 = 4; float r26370 = r26366 - r26369; float r26371 = r26368 - r26370; float r26372 = r26365 / r26371; float r26373 = 676.5203681218851; float r26374 = r26373 / r26368; float r26375 = 0.9999999999998099; float r26376 = r26374 + r26375; float r26377 = r26372 + r26376; float r26378 = -1259.1392167224028; float r26379 = 2; float r26380 = r26366 + r26379; float r26381 = r26366 + r26367; float r26382 = r26380 - r26381; float r26383 = r26378 / r26382; float r26384 = 771.3234287776531; float r26385 = 3; float r26386 = r26385 + r26366; float r26387 = r26386 - r26381; float r26388 = r26384 / r26387; float r26389 = r26383 + r26388; float r26390 = r26377 + r26389; float r26391 = -0.13857109526572012; float r26392 = 6; float r26393 = r26392 + r26366; float r26394 = r26393 - r26381; float r26395 = r26391 / r26394; float r26396 = 12.507343278686905; float r26397 = 5; float r26398 = r26366 - r26397; float r26399 = r26368 - r26398; float r26400 = r26396 / r26399; float r26401 = r26395 + r26400; float r26402 = 1.5056327351493116e-07; float r26403 = 8; float r26404 = r26366 + r26403; float r26405 = r26404 - r26381; float r26406 = r26402 / r26405; float r26407 = 9.984369578019572e-06; float r26408 = 7; float r26409 = r26366 - r26408; float r26410 = r26368 - r26409; float r26411 = r26407 / r26410; float r26412 = r26406 + r26411; float r26413 = r26401 + r26412; float r26414 = r26390 + r26413; float r26415 = atan2(1.0, 0.0); float r26416 = r26367 * r26415; float r26417 = sin(r26416); float r26418 = r26415 / r26417; float r26419 = r26414 * r26418; float r26420 = 0.5; float r26421 = r26420 + r26410; float r26422 = r26366 - r26420; float r26423 = r26368 - r26422; float r26424 = pow(r26421, r26423); float r26425 = r26379 * r26415; float r26426 = sqrt(r26425); float r26427 = r26426 * r26425; float r26428 = cbrt(r26427); float r26429 = r26424 * r26428; float r26430 = -r26421; float r26431 = exp(r26430); float r26432 = r26429 * r26431; float r26433 = r26419 * r26432; return r26433; } double f_od(double z) { double r26434 = -176.6150291621406; double r26435 = 1; double r26436 = z; double r26437 = r26435 - r26436; double r26438 = 4; double r26439 = r26435 - r26438; double r26440 = r26437 - r26439; double r26441 = r26434 / r26440; double r26442 = 676.5203681218851; double r26443 = r26442 / r26437; double r26444 = 0.9999999999998099; double r26445 = r26443 + r26444; double r26446 = r26441 + r26445; double r26447 = -1259.1392167224028; double r26448 = 2; double r26449 = r26435 + r26448; double r26450 = r26435 + r26436; double r26451 = r26449 - r26450; double r26452 = r26447 / r26451; double r26453 = 771.3234287776531; double r26454 = 3; double r26455 = r26454 + r26435; double r26456 = r26455 - r26450; double r26457 = r26453 / r26456; double r26458 = r26452 + r26457; double r26459 = r26446 + r26458; double r26460 = -0.13857109526572012; double r26461 = 6; double r26462 = r26461 + r26435; double r26463 = r26462 - r26450; double r26464 = r26460 / r26463; double r26465 = 12.507343278686905; double r26466 = 5; double r26467 = r26435 - r26466; double r26468 = r26437 - r26467; double r26469 = r26465 / r26468; double r26470 = r26464 + r26469; double r26471 = 1.5056327351493116e-07; double r26472 = 8; double r26473 = r26435 + r26472; double r26474 = r26473 - r26450; double r26475 = r26471 / r26474; double r26476 = 9.984369578019572e-06; double r26477 = 7; double r26478 = r26435 - r26477; double r26479 = r26437 - r26478; double r26480 = r26476 / r26479; double r26481 = r26475 + r26480; double r26482 = r26470 + r26481; double r26483 = r26459 + r26482; double r26484 = atan2(1.0, 0.0); double r26485 = r26436 * r26484; double r26486 = sin(r26485); double r26487 = r26484 / r26486; double r26488 = r26483 * r26487; double r26489 = 0.5; double r26490 = r26489 + r26479; double r26491 = r26435 - r26489; double r26492 = r26437 - r26491; double r26493 = pow(r26490, r26492); double r26494 = r26448 * r26484; double r26495 = sqrt(r26494); double r26496 = r26495 * r26494; double r26497 = cbrt(r26496); double r26498 = r26493 * r26497; double r26499 = -r26490; double r26500 = exp(r26499); double r26501 = r26498 * r26500; double r26502 = r26488 * r26501; return r26502; } 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 r26503, r26504, r26505, r26506, r26507, r26508, r26509, r26510, r26511, r26512, r26513, r26514, r26515, r26516, r26517, r26518, r26519, r26520, r26521, r26522, r26523, r26524, r26525, r26526, r26527, r26528, r26529, r26530, r26531, r26532, r26533, r26534, r26535, r26536, r26537, r26538, r26539, r26540, r26541, r26542, r26543, r26544, r26545, r26546, r26547, r26548, r26549, r26550, r26551, r26552, r26553, r26554, r26555, r26556, r26557, r26558, r26559, r26560, r26561, r26562; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r26503); mpfr_init(r26504); mpfr_init(r26505); mpfr_init(r26506); mpfr_init(r26507); mpfr_init_set_str(r26508, "2", 10, MPFR_RNDN); mpfr_init(r26509); mpfr_init(r26510); mpfr_init_set_str(r26511, "1", 10, MPFR_RNDN); mpfr_init(r26512); mpfr_init(r26513); mpfr_init_set_str(r26514, "7", 10, MPFR_RNDN); mpfr_init(r26515); mpfr_init_set_str(r26516, "0.5", 10, MPFR_RNDN); mpfr_init(r26517); mpfr_init(r26518); mpfr_init(r26519); mpfr_init(r26520); mpfr_init(r26521); mpfr_init(r26522); mpfr_init(r26523); mpfr_init_set_str(r26524, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init_set_str(r26525, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r26526); mpfr_init(r26527); mpfr_init(r26528); mpfr_init_set_str(r26529, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init(r26530); mpfr_init(r26531); mpfr_init(r26532); mpfr_init_set_str(r26533, "771.3234287776531", 10, MPFR_RNDN); mpfr_init_set_str(r26534, "3", 10, MPFR_RNDN); mpfr_init(r26535); mpfr_init(r26536); mpfr_init(r26537); mpfr_init_set_str(r26538, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r26539, "4", 10, MPFR_RNDN); mpfr_init(r26540); mpfr_init(r26541); mpfr_init(r26542); mpfr_init_set_str(r26543, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r26544, "5", 10, MPFR_RNDN); mpfr_init(r26545); mpfr_init(r26546); mpfr_init(r26547); mpfr_init_set_str(r26548, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r26549, "6", 10, MPFR_RNDN); mpfr_init(r26550); mpfr_init(r26551); mpfr_init(r26552); mpfr_init_set_str(r26553, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init(r26554); mpfr_init(r26555); mpfr_init_set_str(r26556, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r26557, "8", 10, MPFR_RNDN); mpfr_init(r26558); mpfr_init(r26559); mpfr_init(r26560); mpfr_init(r26561); mpfr_init(r26562); } double f_im(double z) { mpfr_const_pi(r26503, MPFR_RNDN); mpfr_set_d(r26504, z, MPFR_RNDN); mpfr_mul(r26505, r26503, r26504, MPFR_RNDN); mpfr_sin(r26506, r26505, MPFR_RNDN); mpfr_div(r26507, r26503, r26506, MPFR_RNDN); ; mpfr_mul(r26509, r26503, r26508, MPFR_RNDN); mpfr_sqrt(r26510, r26509, MPFR_RNDN); ; mpfr_sub(r26512, r26511, r26504, MPFR_RNDN); mpfr_sub(r26513, r26512, r26511, MPFR_RNDN); ; mpfr_add(r26515, r26513, r26514, MPFR_RNDN); ; mpfr_add(r26517, r26515, r26516, MPFR_RNDN); mpfr_add(r26518, r26513, r26516, MPFR_RNDN); mpfr_pow(r26519, r26517, r26518, MPFR_RNDN); mpfr_mul(r26520, r26510, r26519, MPFR_RNDN); mpfr_neg(r26521, r26517, MPFR_RNDN); mpfr_exp(r26522, r26521, MPFR_RNDN); mpfr_mul(r26523, r26520, r26522, MPFR_RNDN); ; ; mpfr_add(r26526, r26513, r26511, MPFR_RNDN); mpfr_div(r26527, r26525, r26526, MPFR_RNDN); mpfr_add(r26528, r26524, r26527, MPFR_RNDN); ; mpfr_add(r26530, r26513, r26508, MPFR_RNDN); mpfr_div(r26531, r26529, r26530, MPFR_RNDN); mpfr_add(r26532, r26528, r26531, MPFR_RNDN); ; ; mpfr_add(r26535, r26513, r26534, MPFR_RNDN); mpfr_div(r26536, r26533, r26535, MPFR_RNDN); mpfr_add(r26537, r26532, r26536, MPFR_RNDN); ; ; mpfr_add(r26540, r26513, r26539, MPFR_RNDN); mpfr_div(r26541, r26538, r26540, MPFR_RNDN); mpfr_add(r26542, r26537, r26541, MPFR_RNDN); ; ; mpfr_add(r26545, r26513, r26544, MPFR_RNDN); mpfr_div(r26546, r26543, r26545, MPFR_RNDN); mpfr_add(r26547, r26542, r26546, MPFR_RNDN); ; ; mpfr_add(r26550, r26513, r26549, MPFR_RNDN); mpfr_div(r26551, r26548, r26550, MPFR_RNDN); mpfr_add(r26552, r26547, r26551, MPFR_RNDN); ; mpfr_div(r26554, r26553, r26515, MPFR_RNDN); mpfr_add(r26555, r26552, r26554, MPFR_RNDN); ; ; mpfr_add(r26558, r26513, r26557, MPFR_RNDN); mpfr_div(r26559, r26556, r26558, MPFR_RNDN); mpfr_add(r26560, r26555, r26559, MPFR_RNDN); mpfr_mul(r26561, r26523, r26560, MPFR_RNDN); mpfr_mul(r26562, r26507, r26561, MPFR_RNDN); return mpfr_get_d(r26562, MPFR_RNDN); } static mpfr_t r26563, r26564, r26565, r26566, r26567, r26568, r26569, r26570, r26571, r26572, r26573, r26574, r26575, r26576, r26577, r26578, r26579, r26580, r26581, r26582, r26583, r26584, r26585, r26586, r26587, r26588, r26589, r26590, r26591, r26592, r26593, r26594, r26595, r26596, r26597, r26598, r26599, r26600, r26601, r26602, r26603, r26604, r26605, r26606, r26607, r26608, r26609, r26610, r26611, r26612, r26613, r26614, r26615, r26616, r26617, r26618, r26619, r26620, r26621, r26622, r26623, r26624, r26625, r26626, r26627, r26628, r26629, r26630, r26631; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r26563, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r26564, "1", 10, MPFR_RNDN); mpfr_init(r26565); mpfr_init(r26566); mpfr_init_set_str(r26567, "4", 10, MPFR_RNDN); mpfr_init(r26568); mpfr_init(r26569); mpfr_init(r26570); mpfr_init_set_str(r26571, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r26572); mpfr_init_set_str(r26573, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init(r26574); mpfr_init(r26575); mpfr_init_set_str(r26576, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init_set_str(r26577, "2", 10, MPFR_RNDN); mpfr_init(r26578); mpfr_init(r26579); mpfr_init(r26580); mpfr_init(r26581); mpfr_init_set_str(r26582, "771.3234287776531", 10, MPFR_RNDN); mpfr_init_set_str(r26583, "3", 10, MPFR_RNDN); mpfr_init(r26584); mpfr_init(r26585); mpfr_init(r26586); mpfr_init(r26587); mpfr_init(r26588); mpfr_init_set_str(r26589, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r26590, "6", 10, MPFR_RNDN); mpfr_init(r26591); mpfr_init(r26592); mpfr_init(r26593); mpfr_init_set_str(r26594, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r26595, "5", 10, MPFR_RNDN); mpfr_init(r26596); mpfr_init(r26597); mpfr_init(r26598); mpfr_init(r26599); mpfr_init_set_str(r26600, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r26601, "8", 10, MPFR_RNDN); mpfr_init(r26602); mpfr_init(r26603); mpfr_init(r26604); mpfr_init_set_str(r26605, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init_set_str(r26606, "7", 10, MPFR_RNDN); mpfr_init(r26607); mpfr_init(r26608); mpfr_init(r26609); mpfr_init(r26610); mpfr_init(r26611); mpfr_init(r26612); mpfr_init(r26613); mpfr_init(r26614); mpfr_init(r26615); mpfr_init(r26616); mpfr_init(r26617); mpfr_init_set_str(r26618, "0.5", 10, MPFR_RNDN); mpfr_init(r26619); mpfr_init(r26620); mpfr_init(r26621); mpfr_init(r26622); mpfr_init(r26623); mpfr_init(r26624); mpfr_init(r26625); mpfr_init(r26626); mpfr_init(r26627); mpfr_init(r26628); mpfr_init(r26629); mpfr_init(r26630); mpfr_init(r26631); } double f_fm(double z) { ; ; mpfr_set_d(r26565, z, MPFR_RNDN); mpfr_sub(r26566, r26564, r26565, MPFR_RNDN); ; mpfr_sub(r26568, r26564, r26567, MPFR_RNDN); mpfr_sub(r26569, r26566, r26568, MPFR_RNDN); mpfr_div(r26570, r26563, r26569, MPFR_RNDN); ; mpfr_div(r26572, r26571, r26566, MPFR_RNDN); ; mpfr_add(r26574, r26572, r26573, MPFR_RNDN); mpfr_add(r26575, r26570, r26574, MPFR_RNDN); ; ; mpfr_add(r26578, r26564, r26577, MPFR_RNDN); mpfr_add(r26579, r26564, r26565, MPFR_RNDN); mpfr_sub(r26580, r26578, r26579, MPFR_RNDN); mpfr_div(r26581, r26576, r26580, MPFR_RNDN); ; ; mpfr_add(r26584, r26583, r26564, MPFR_RNDN); mpfr_sub(r26585, r26584, r26579, MPFR_RNDN); mpfr_div(r26586, r26582, r26585, MPFR_RNDN); mpfr_add(r26587, r26581, r26586, MPFR_RNDN); mpfr_add(r26588, r26575, r26587, MPFR_RNDN); ; ; mpfr_add(r26591, r26590, r26564, MPFR_RNDN); mpfr_sub(r26592, r26591, r26579, MPFR_RNDN); mpfr_div(r26593, r26589, r26592, MPFR_RNDN); ; ; mpfr_sub(r26596, r26564, r26595, MPFR_RNDN); mpfr_sub(r26597, r26566, r26596, MPFR_RNDN); mpfr_div(r26598, r26594, r26597, MPFR_RNDN); mpfr_add(r26599, r26593, r26598, MPFR_RNDN); ; ; mpfr_add(r26602, r26564, r26601, MPFR_RNDN); mpfr_sub(r26603, r26602, r26579, MPFR_RNDN); mpfr_div(r26604, r26600, r26603, MPFR_RNDN); ; ; mpfr_sub(r26607, r26564, r26606, MPFR_RNDN); mpfr_sub(r26608, r26566, r26607, MPFR_RNDN); mpfr_div(r26609, r26605, r26608, MPFR_RNDN); mpfr_add(r26610, r26604, r26609, MPFR_RNDN); mpfr_add(r26611, r26599, r26610, MPFR_RNDN); mpfr_add(r26612, r26588, r26611, MPFR_RNDN); mpfr_const_pi(r26613, MPFR_RNDN); mpfr_mul(r26614, r26565, r26613, MPFR_RNDN); mpfr_sin(r26615, r26614, MPFR_RNDN); mpfr_div(r26616, r26613, r26615, MPFR_RNDN); mpfr_mul(r26617, r26612, r26616, MPFR_RNDN); ; mpfr_add(r26619, r26618, r26608, MPFR_RNDN); mpfr_sub(r26620, r26564, r26618, MPFR_RNDN); mpfr_sub(r26621, r26566, r26620, MPFR_RNDN); mpfr_pow(r26622, r26619, r26621, MPFR_RNDN); mpfr_mul(r26623, r26577, r26613, MPFR_RNDN); mpfr_sqrt(r26624, r26623, MPFR_RNDN); mpfr_mul(r26625, r26624, r26623, MPFR_RNDN); mpfr_cbrt(r26626, r26625, MPFR_RNDN); mpfr_mul(r26627, r26622, r26626, MPFR_RNDN); mpfr_neg(r26628, r26619, MPFR_RNDN); mpfr_exp(r26629, r26628, MPFR_RNDN); mpfr_mul(r26630, r26627, r26629, MPFR_RNDN); mpfr_mul(r26631, r26617, r26630, MPFR_RNDN); return mpfr_get_d(r26631, MPFR_RNDN); } static mpfr_t r26632, r26633, r26634, r26635, r26636, r26637, r26638, r26639, r26640, r26641, r26642, r26643, r26644, r26645, r26646, r26647, r26648, r26649, r26650, r26651, r26652, r26653, r26654, r26655, r26656, r26657, r26658, r26659, r26660, r26661, r26662, r26663, r26664, r26665, r26666, r26667, r26668, r26669, r26670, r26671, r26672, r26673, r26674, r26675, r26676, r26677, r26678, r26679, r26680, r26681, r26682, r26683, r26684, r26685, r26686, r26687, r26688, r26689, r26690, r26691, r26692, r26693, r26694, r26695, r26696, r26697, r26698, r26699, r26700; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r26632, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r26633, "1", 10, MPFR_RNDN); mpfr_init(r26634); mpfr_init(r26635); mpfr_init_set_str(r26636, "4", 10, MPFR_RNDN); mpfr_init(r26637); mpfr_init(r26638); mpfr_init(r26639); mpfr_init_set_str(r26640, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r26641); mpfr_init_set_str(r26642, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init(r26643); mpfr_init(r26644); mpfr_init_set_str(r26645, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init_set_str(r26646, "2", 10, MPFR_RNDN); mpfr_init(r26647); mpfr_init(r26648); mpfr_init(r26649); mpfr_init(r26650); mpfr_init_set_str(r26651, "771.3234287776531", 10, MPFR_RNDN); mpfr_init_set_str(r26652, "3", 10, MPFR_RNDN); mpfr_init(r26653); mpfr_init(r26654); mpfr_init(r26655); mpfr_init(r26656); mpfr_init(r26657); mpfr_init_set_str(r26658, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r26659, "6", 10, MPFR_RNDN); mpfr_init(r26660); mpfr_init(r26661); mpfr_init(r26662); mpfr_init_set_str(r26663, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r26664, "5", 10, MPFR_RNDN); mpfr_init(r26665); mpfr_init(r26666); mpfr_init(r26667); mpfr_init(r26668); mpfr_init_set_str(r26669, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r26670, "8", 10, MPFR_RNDN); mpfr_init(r26671); mpfr_init(r26672); mpfr_init(r26673); mpfr_init_set_str(r26674, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init_set_str(r26675, "7", 10, MPFR_RNDN); mpfr_init(r26676); mpfr_init(r26677); mpfr_init(r26678); mpfr_init(r26679); mpfr_init(r26680); mpfr_init(r26681); mpfr_init(r26682); mpfr_init(r26683); mpfr_init(r26684); mpfr_init(r26685); mpfr_init(r26686); mpfr_init_set_str(r26687, "0.5", 10, MPFR_RNDN); mpfr_init(r26688); mpfr_init(r26689); mpfr_init(r26690); mpfr_init(r26691); mpfr_init(r26692); mpfr_init(r26693); mpfr_init(r26694); mpfr_init(r26695); mpfr_init(r26696); mpfr_init(r26697); mpfr_init(r26698); mpfr_init(r26699); mpfr_init(r26700); } double f_dm(double z) { ; ; mpfr_set_d(r26634, z, MPFR_RNDN); mpfr_sub(r26635, r26633, r26634, MPFR_RNDN); ; mpfr_sub(r26637, r26633, r26636, MPFR_RNDN); mpfr_sub(r26638, r26635, r26637, MPFR_RNDN); mpfr_div(r26639, r26632, r26638, MPFR_RNDN); ; mpfr_div(r26641, r26640, r26635, MPFR_RNDN); ; mpfr_add(r26643, r26641, r26642, MPFR_RNDN); mpfr_add(r26644, r26639, r26643, MPFR_RNDN); ; ; mpfr_add(r26647, r26633, r26646, MPFR_RNDN); mpfr_add(r26648, r26633, r26634, MPFR_RNDN); mpfr_sub(r26649, r26647, r26648, MPFR_RNDN); mpfr_div(r26650, r26645, r26649, MPFR_RNDN); ; ; mpfr_add(r26653, r26652, r26633, MPFR_RNDN); mpfr_sub(r26654, r26653, r26648, MPFR_RNDN); mpfr_div(r26655, r26651, r26654, MPFR_RNDN); mpfr_add(r26656, r26650, r26655, MPFR_RNDN); mpfr_add(r26657, r26644, r26656, MPFR_RNDN); ; ; mpfr_add(r26660, r26659, r26633, MPFR_RNDN); mpfr_sub(r26661, r26660, r26648, MPFR_RNDN); mpfr_div(r26662, r26658, r26661, MPFR_RNDN); ; ; mpfr_sub(r26665, r26633, r26664, MPFR_RNDN); mpfr_sub(r26666, r26635, r26665, MPFR_RNDN); mpfr_div(r26667, r26663, r26666, MPFR_RNDN); mpfr_add(r26668, r26662, r26667, MPFR_RNDN); ; ; mpfr_add(r26671, r26633, r26670, MPFR_RNDN); mpfr_sub(r26672, r26671, r26648, MPFR_RNDN); mpfr_div(r26673, r26669, r26672, MPFR_RNDN); ; ; mpfr_sub(r26676, r26633, r26675, MPFR_RNDN); mpfr_sub(r26677, r26635, r26676, MPFR_RNDN); mpfr_div(r26678, r26674, r26677, MPFR_RNDN); mpfr_add(r26679, r26673, r26678, MPFR_RNDN); mpfr_add(r26680, r26668, r26679, MPFR_RNDN); mpfr_add(r26681, r26657, r26680, MPFR_RNDN); mpfr_const_pi(r26682, MPFR_RNDN); mpfr_mul(r26683, r26634, r26682, MPFR_RNDN); mpfr_sin(r26684, r26683, MPFR_RNDN); mpfr_div(r26685, r26682, r26684, MPFR_RNDN); mpfr_mul(r26686, r26681, r26685, MPFR_RNDN); ; mpfr_add(r26688, r26687, r26677, MPFR_RNDN); mpfr_sub(r26689, r26633, r26687, MPFR_RNDN); mpfr_sub(r26690, r26635, r26689, MPFR_RNDN); mpfr_pow(r26691, r26688, r26690, MPFR_RNDN); mpfr_mul(r26692, r26646, r26682, MPFR_RNDN); mpfr_sqrt(r26693, r26692, MPFR_RNDN); mpfr_mul(r26694, r26693, r26692, MPFR_RNDN); mpfr_cbrt(r26695, r26694, MPFR_RNDN); mpfr_mul(r26696, r26691, r26695, MPFR_RNDN); mpfr_neg(r26697, r26688, MPFR_RNDN); mpfr_exp(r26698, r26697, MPFR_RNDN); mpfr_mul(r26699, r26696, r26698, MPFR_RNDN); mpfr_mul(r26700, r26686, r26699, MPFR_RNDN); return mpfr_get_d(r26700, MPFR_RNDN); }