#include #include #include #include #include char *name = "kepler0"; double f_if(float x1, float x2, float x3, float x4, float x5, float x6) { float r22393 = x2; float r22394 = x5; float r22395 = r22393 * r22394; float r22396 = x3; float r22397 = x6; float r22398 = r22396 * r22397; float r22399 = r22395 + r22398; float r22400 = r22393 * r22396; float r22401 = r22399 - r22400; float r22402 = r22394 * r22397; float r22403 = r22401 - r22402; float r22404 = x1; float r22405 = -r22404; float r22406 = r22405 + r22393; float r22407 = r22406 + r22396; float r22408 = x4; float r22409 = r22407 - r22408; float r22410 = r22409 + r22394; float r22411 = r22410 + r22397; float r22412 = r22404 * r22411; float r22413 = r22403 + r22412; return r22413; } double f_id(double x1, double x2, double x3, double x4, double x5, double x6) { double r22414 = x2; double r22415 = x5; double r22416 = r22414 * r22415; double r22417 = x3; double r22418 = x6; double r22419 = r22417 * r22418; double r22420 = r22416 + r22419; double r22421 = r22414 * r22417; double r22422 = r22420 - r22421; double r22423 = r22415 * r22418; double r22424 = r22422 - r22423; double r22425 = x1; double r22426 = -r22425; double r22427 = r22426 + r22414; double r22428 = r22427 + r22417; double r22429 = x4; double r22430 = r22428 - r22429; double r22431 = r22430 + r22415; double r22432 = r22431 + r22418; double r22433 = r22425 * r22432; double r22434 = r22424 + r22433; return r22434; } double f_of(float x1, float x2, float x3, float x4, float x5, float x6) { float r22435 = x2; float r22436 = x4; float r22437 = x3; float r22438 = r22436 - r22437; float r22439 = r22435 - r22438; float r22440 = x1; float r22441 = -r22440; float r22442 = x5; float r22443 = r22441 + r22442; float r22444 = x6; float r22445 = r22443 + r22444; float r22446 = r22439 + r22445; float r22447 = r22446 * r22440; float r22448 = r22437 - r22442; float r22449 = r22444 - r22435; float r22450 = r22448 * r22449; float r22451 = r22447 + r22450; return r22451; } double f_od(double x1, double x2, double x3, double x4, double x5, double x6) { double r22452 = x2; double r22453 = x4; double r22454 = x3; double r22455 = r22453 - r22454; double r22456 = r22452 - r22455; double r22457 = x1; double r22458 = -r22457; double r22459 = x5; double r22460 = r22458 + r22459; double r22461 = x6; double r22462 = r22460 + r22461; double r22463 = r22456 + r22462; double r22464 = r22463 * r22457; double r22465 = r22454 - r22459; double r22466 = r22461 - r22452; double r22467 = r22465 * r22466; double r22468 = r22464 + r22467; return r22468; } 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 r22469, r22470, r22471, r22472, r22473, r22474, r22475, r22476, r22477, r22478, r22479, r22480, r22481, r22482, r22483, r22484, r22485, r22486, r22487, r22488, r22489; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22469); mpfr_init(r22470); mpfr_init(r22471); mpfr_init(r22472); mpfr_init(r22473); mpfr_init(r22474); mpfr_init(r22475); mpfr_init(r22476); mpfr_init(r22477); mpfr_init(r22478); mpfr_init(r22479); mpfr_init(r22480); mpfr_init(r22481); mpfr_init(r22482); mpfr_init(r22483); mpfr_init(r22484); mpfr_init(r22485); mpfr_init(r22486); mpfr_init(r22487); mpfr_init(r22488); mpfr_init(r22489); } double f_im(double x1, double x2, double x3, double x4, double x5, double x6) { mpfr_set_d(r22469, x2, MPFR_RNDN); mpfr_set_d(r22470, x5, MPFR_RNDN); mpfr_mul(r22471, r22469, r22470, MPFR_RNDN); mpfr_set_d(r22472, x3, MPFR_RNDN); mpfr_set_d(r22473, x6, MPFR_RNDN); mpfr_mul(r22474, r22472, r22473, MPFR_RNDN); mpfr_add(r22475, r22471, r22474, MPFR_RNDN); mpfr_mul(r22476, r22469, r22472, MPFR_RNDN); mpfr_sub(r22477, r22475, r22476, MPFR_RNDN); mpfr_mul(r22478, r22470, r22473, MPFR_RNDN); mpfr_sub(r22479, r22477, r22478, MPFR_RNDN); mpfr_set_d(r22480, x1, MPFR_RNDN); mpfr_neg(r22481, r22480, MPFR_RNDN); mpfr_add(r22482, r22481, r22469, MPFR_RNDN); mpfr_add(r22483, r22482, r22472, MPFR_RNDN); mpfr_set_d(r22484, x4, MPFR_RNDN); mpfr_sub(r22485, r22483, r22484, MPFR_RNDN); mpfr_add(r22486, r22485, r22470, MPFR_RNDN); mpfr_add(r22487, r22486, r22473, MPFR_RNDN); mpfr_mul(r22488, r22480, r22487, MPFR_RNDN); mpfr_add(r22489, r22479, r22488, MPFR_RNDN); return mpfr_get_d(r22489, MPFR_RNDN); } static mpfr_t r22490, r22491, r22492, r22493, r22494, r22495, r22496, r22497, r22498, r22499, r22500, r22501, r22502, r22503, r22504, r22505, r22506; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22490); mpfr_init(r22491); mpfr_init(r22492); mpfr_init(r22493); mpfr_init(r22494); mpfr_init(r22495); mpfr_init(r22496); mpfr_init(r22497); mpfr_init(r22498); mpfr_init(r22499); mpfr_init(r22500); mpfr_init(r22501); mpfr_init(r22502); mpfr_init(r22503); mpfr_init(r22504); mpfr_init(r22505); mpfr_init(r22506); } double f_fm(double x1, double x2, double x3, double x4, double x5, double x6) { mpfr_set_d(r22490, x2, MPFR_RNDN); mpfr_set_d(r22491, x4, MPFR_RNDN); mpfr_set_d(r22492, x3, MPFR_RNDN); mpfr_sub(r22493, r22491, r22492, MPFR_RNDN); mpfr_sub(r22494, r22490, r22493, MPFR_RNDN); mpfr_set_d(r22495, x1, MPFR_RNDN); mpfr_neg(r22496, r22495, MPFR_RNDN); mpfr_set_d(r22497, x5, MPFR_RNDN); mpfr_add(r22498, r22496, r22497, MPFR_RNDN); mpfr_set_d(r22499, x6, MPFR_RNDN); mpfr_add(r22500, r22498, r22499, MPFR_RNDN); mpfr_add(r22501, r22494, r22500, MPFR_RNDN); mpfr_mul(r22502, r22501, r22495, MPFR_RNDN); mpfr_sub(r22503, r22492, r22497, MPFR_RNDN); mpfr_sub(r22504, r22499, r22490, MPFR_RNDN); mpfr_mul(r22505, r22503, r22504, MPFR_RNDN); mpfr_add(r22506, r22502, r22505, MPFR_RNDN); return mpfr_get_d(r22506, MPFR_RNDN); } static mpfr_t r22507, r22508, r22509, r22510, r22511, r22512, r22513, r22514, r22515, r22516, r22517, r22518, r22519, r22520, r22521, r22522, r22523; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22507); mpfr_init(r22508); mpfr_init(r22509); mpfr_init(r22510); mpfr_init(r22511); mpfr_init(r22512); mpfr_init(r22513); mpfr_init(r22514); mpfr_init(r22515); mpfr_init(r22516); mpfr_init(r22517); mpfr_init(r22518); mpfr_init(r22519); mpfr_init(r22520); mpfr_init(r22521); mpfr_init(r22522); mpfr_init(r22523); } double f_dm(double x1, double x2, double x3, double x4, double x5, double x6) { mpfr_set_d(r22507, x2, MPFR_RNDN); mpfr_set_d(r22508, x4, MPFR_RNDN); mpfr_set_d(r22509, x3, MPFR_RNDN); mpfr_sub(r22510, r22508, r22509, MPFR_RNDN); mpfr_sub(r22511, r22507, r22510, MPFR_RNDN); mpfr_set_d(r22512, x1, MPFR_RNDN); mpfr_neg(r22513, r22512, MPFR_RNDN); mpfr_set_d(r22514, x5, MPFR_RNDN); mpfr_add(r22515, r22513, r22514, MPFR_RNDN); mpfr_set_d(r22516, x6, MPFR_RNDN); mpfr_add(r22517, r22515, r22516, MPFR_RNDN); mpfr_add(r22518, r22511, r22517, MPFR_RNDN); mpfr_mul(r22519, r22518, r22512, MPFR_RNDN); mpfr_sub(r22520, r22509, r22514, MPFR_RNDN); mpfr_sub(r22521, r22516, r22507, MPFR_RNDN); mpfr_mul(r22522, r22520, r22521, MPFR_RNDN); mpfr_add(r22523, r22519, r22522, MPFR_RNDN); return mpfr_get_d(r22523, MPFR_RNDN); }