#include #include #include #include #include char *name = "Octave 3.8, jcobi/2"; double f_if(float alpha, float beta, float i) { float r26538 = alpha; float r26539 = beta; float r26540 = r26538 + r26539; float r26541 = r26539 - r26538; float r26542 = r26540 * r26541; float r26543 = 2; float r26544 = i; float r26545 = r26543 * r26544; float r26546 = r26540 + r26545; float r26547 = r26542 / r26546; float r26548 = 2.0; float r26549 = r26546 + r26548; float r26550 = r26547 / r26549; float r26551 = 1.0; float r26552 = r26550 + r26551; float r26553 = r26552 / r26548; return r26553; } double f_id(double alpha, double beta, double i) { double r26554 = alpha; double r26555 = beta; double r26556 = r26554 + r26555; double r26557 = r26555 - r26554; double r26558 = r26556 * r26557; double r26559 = 2; double r26560 = i; double r26561 = r26559 * r26560; double r26562 = r26556 + r26561; double r26563 = r26558 / r26562; double r26564 = 2.0; double r26565 = r26562 + r26564; double r26566 = r26563 / r26565; double r26567 = 1.0; double r26568 = r26566 + r26567; double r26569 = r26568 / r26564; return r26569; } double f_of(float alpha, float beta, float i) { float r26570 = beta; float r26571 = alpha; float r26572 = r26570 + r26571; float r26573 = i; float r26574 = 2; float r26575 = fma(r26573, r26574, r26570); float r26576 = 2.0; float r26577 = r26576 + r26571; float r26578 = r26575 + r26577; float r26579 = r26572 / r26578; float r26580 = fma(r26574, r26573, r26572); float r26581 = r26570 / r26580; float r26582 = r26571 / r26580; float r26583 = r26581 - r26582; float r26584 = 1.0; float r26585 = fma(r26579, r26583, r26584); float r26586 = exp(r26585); float r26587 = log(r26586); float r26588 = r26587 / r26576; return r26588; } double f_od(double alpha, double beta, double i) { double r26589 = beta; double r26590 = alpha; double r26591 = r26589 + r26590; double r26592 = i; double r26593 = 2; double r26594 = fma(r26592, r26593, r26589); double r26595 = 2.0; double r26596 = r26595 + r26590; double r26597 = r26594 + r26596; double r26598 = r26591 / r26597; double r26599 = fma(r26593, r26592, r26591); double r26600 = r26589 / r26599; double r26601 = r26590 / r26599; double r26602 = r26600 - r26601; double r26603 = 1.0; double r26604 = fma(r26598, r26602, r26603); double r26605 = exp(r26604); double r26606 = log(r26605); double r26607 = r26606 / r26595; return r26607; } 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 r26608, r26609, r26610, r26611, r26612, r26613, r26614, r26615, r26616, r26617, r26618, r26619, r26620, r26621, r26622, r26623; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r26608); mpfr_init(r26609); mpfr_init(r26610); mpfr_init(r26611); mpfr_init(r26612); mpfr_init_set_str(r26613, "2", 10, MPFR_RNDN); mpfr_init(r26614); mpfr_init(r26615); mpfr_init(r26616); mpfr_init(r26617); mpfr_init_set_str(r26618, "2.0", 10, MPFR_RNDN); mpfr_init(r26619); mpfr_init(r26620); mpfr_init_set_str(r26621, "1.0", 10, MPFR_RNDN); mpfr_init(r26622); mpfr_init(r26623); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r26608, alpha, MPFR_RNDN); mpfr_set_d(r26609, beta, MPFR_RNDN); mpfr_add(r26610, r26608, r26609, MPFR_RNDN); mpfr_sub(r26611, r26609, r26608, MPFR_RNDN); mpfr_mul(r26612, r26610, r26611, MPFR_RNDN); ; mpfr_set_d(r26614, i, MPFR_RNDN); mpfr_mul(r26615, r26613, r26614, MPFR_RNDN); mpfr_add(r26616, r26610, r26615, MPFR_RNDN); mpfr_div(r26617, r26612, r26616, MPFR_RNDN); ; mpfr_add(r26619, r26616, r26618, MPFR_RNDN); mpfr_div(r26620, r26617, r26619, MPFR_RNDN); ; mpfr_add(r26622, r26620, r26621, MPFR_RNDN); mpfr_div(r26623, r26622, r26618, MPFR_RNDN); return mpfr_get_d(r26623, MPFR_RNDN); } static mpfr_t r26624, r26625, r26626, r26627, r26628, r26629, r26630, r26631, r26632, r26633, r26634, r26635, r26636, r26637, r26638, r26639, r26640, r26641, r26642; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r26624); mpfr_init(r26625); mpfr_init(r26626); mpfr_init(r26627); mpfr_init_set_str(r26628, "2", 10, MPFR_RNDN); mpfr_init(r26629); mpfr_init_set_str(r26630, "2.0", 10, MPFR_RNDN); mpfr_init(r26631); mpfr_init(r26632); mpfr_init(r26633); mpfr_init(r26634); mpfr_init(r26635); mpfr_init(r26636); mpfr_init(r26637); mpfr_init_set_str(r26638, "1.0", 10, MPFR_RNDN); mpfr_init(r26639); mpfr_init(r26640); mpfr_init(r26641); mpfr_init(r26642); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r26624, beta, MPFR_RNDN); mpfr_set_d(r26625, alpha, MPFR_RNDN); mpfr_add(r26626, r26624, r26625, MPFR_RNDN); mpfr_set_d(r26627, i, MPFR_RNDN); ; mpfr_fma(r26629, r26627, r26628, r26624, MPFR_RNDN); ; mpfr_add(r26631, r26630, r26625, MPFR_RNDN); mpfr_add(r26632, r26629, r26631, MPFR_RNDN); mpfr_div(r26633, r26626, r26632, MPFR_RNDN); mpfr_fma(r26634, r26628, r26627, r26626, MPFR_RNDN); mpfr_div(r26635, r26624, r26634, MPFR_RNDN); mpfr_div(r26636, r26625, r26634, MPFR_RNDN); mpfr_sub(r26637, r26635, r26636, MPFR_RNDN); ; mpfr_fma(r26639, r26633, r26637, r26638, MPFR_RNDN); mpfr_exp(r26640, r26639, MPFR_RNDN); mpfr_log(r26641, r26640, MPFR_RNDN); mpfr_div(r26642, r26641, r26630, MPFR_RNDN); return mpfr_get_d(r26642, MPFR_RNDN); } static mpfr_t r26643, r26644, r26645, r26646, r26647, r26648, r26649, r26650, r26651, r26652, r26653, r26654, r26655, r26656, r26657, r26658, r26659, r26660, r26661; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r26643); mpfr_init(r26644); mpfr_init(r26645); mpfr_init(r26646); mpfr_init_set_str(r26647, "2", 10, MPFR_RNDN); mpfr_init(r26648); mpfr_init_set_str(r26649, "2.0", 10, MPFR_RNDN); mpfr_init(r26650); mpfr_init(r26651); mpfr_init(r26652); mpfr_init(r26653); mpfr_init(r26654); mpfr_init(r26655); mpfr_init(r26656); mpfr_init_set_str(r26657, "1.0", 10, MPFR_RNDN); mpfr_init(r26658); mpfr_init(r26659); mpfr_init(r26660); mpfr_init(r26661); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r26643, beta, MPFR_RNDN); mpfr_set_d(r26644, alpha, MPFR_RNDN); mpfr_add(r26645, r26643, r26644, MPFR_RNDN); mpfr_set_d(r26646, i, MPFR_RNDN); ; mpfr_fma(r26648, r26646, r26647, r26643, MPFR_RNDN); ; mpfr_add(r26650, r26649, r26644, MPFR_RNDN); mpfr_add(r26651, r26648, r26650, MPFR_RNDN); mpfr_div(r26652, r26645, r26651, MPFR_RNDN); mpfr_fma(r26653, r26647, r26646, r26645, MPFR_RNDN); mpfr_div(r26654, r26643, r26653, MPFR_RNDN); mpfr_div(r26655, r26644, r26653, MPFR_RNDN); mpfr_sub(r26656, r26654, r26655, MPFR_RNDN); ; mpfr_fma(r26658, r26652, r26656, r26657, MPFR_RNDN); mpfr_exp(r26659, r26658, MPFR_RNDN); mpfr_log(r26660, r26659, MPFR_RNDN); mpfr_div(r26661, r26660, r26649, MPFR_RNDN); return mpfr_get_d(r26661, MPFR_RNDN); }