#include #include #include #include #include char *name = "Octave 3.8, jcobi/4"; double f_if(float alpha, float beta, float i) { float r26597 = i; float r26598 = alpha; float r26599 = beta; float r26600 = r26598 + r26599; float r26601 = r26600 + r26597; float r26602 = r26597 * r26601; float r26603 = r26599 * r26598; float r26604 = r26603 + r26602; float r26605 = r26602 * r26604; float r26606 = 2; float r26607 = r26606 * r26597; float r26608 = r26600 + r26607; float r26609 = r26608 * r26608; float r26610 = r26605 / r26609; float r26611 = 1.0; float r26612 = r26609 - r26611; float r26613 = r26610 / r26612; return r26613; } double f_id(double alpha, double beta, double i) { double r26614 = i; double r26615 = alpha; double r26616 = beta; double r26617 = r26615 + r26616; double r26618 = r26617 + r26614; double r26619 = r26614 * r26618; double r26620 = r26616 * r26615; double r26621 = r26620 + r26619; double r26622 = r26619 * r26621; double r26623 = 2; double r26624 = r26623 * r26614; double r26625 = r26617 + r26624; double r26626 = r26625 * r26625; double r26627 = r26622 / r26626; double r26628 = 1.0; double r26629 = r26626 - r26628; double r26630 = r26627 / r26629; return r26630; } double f_of(float alpha, float beta, float i) { float r26631 = beta; float r26632 = 1.5126697473829053e+144; bool r26633 = r26631 <= r26632; float r26634 = i; float r26635 = alpha; float r26636 = r26634 + r26635; float r26637 = r26636 + r26631; float r26638 = r26634 + r26634; float r26639 = r26635 + r26631; float r26640 = r26638 + r26639; float r26641 = r26640 / r26634; float r26642 = r26637 / r26641; float r26643 = r26639 * r26634; float r26644 = r26634 * r26634; float r26645 = r26631 * r26635; float r26646 = r26644 + r26645; float r26647 = r26643 + r26646; float r26648 = r26647 / r26640; float r26649 = r26642 * r26648; float r26650 = r26640 * r26640; float r26651 = 1.0; float r26652 = r26650 - r26651; float r26653 = r26649 / r26652; float r26654 = 0; float r26655 = r26633 ? r26653 : r26654; return r26655; } double f_od(double alpha, double beta, double i) { double r26656 = beta; double r26657 = 1.5126697473829053e+144; bool r26658 = r26656 <= r26657; double r26659 = i; double r26660 = alpha; double r26661 = r26659 + r26660; double r26662 = r26661 + r26656; double r26663 = r26659 + r26659; double r26664 = r26660 + r26656; double r26665 = r26663 + r26664; double r26666 = r26665 / r26659; double r26667 = r26662 / r26666; double r26668 = r26664 * r26659; double r26669 = r26659 * r26659; double r26670 = r26656 * r26660; double r26671 = r26669 + r26670; double r26672 = r26668 + r26671; double r26673 = r26672 / r26665; double r26674 = r26667 * r26673; double r26675 = r26665 * r26665; double r26676 = 1.0; double r26677 = r26675 - r26676; double r26678 = r26674 / r26677; double r26679 = 0; double r26680 = r26658 ? r26678 : r26679; return r26680; } 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 r26681, r26682, r26683, r26684, r26685, r26686, r26687, r26688, r26689, r26690, r26691, r26692, r26693, r26694, r26695, r26696, r26697; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r26681); mpfr_init(r26682); mpfr_init(r26683); mpfr_init(r26684); mpfr_init(r26685); mpfr_init(r26686); mpfr_init(r26687); mpfr_init(r26688); mpfr_init(r26689); mpfr_init_set_str(r26690, "2", 10, MPFR_RNDN); mpfr_init(r26691); mpfr_init(r26692); mpfr_init(r26693); mpfr_init(r26694); mpfr_init_set_str(r26695, "1.0", 10, MPFR_RNDN); mpfr_init(r26696); mpfr_init(r26697); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r26681, i, MPFR_RNDN); mpfr_set_d(r26682, alpha, MPFR_RNDN); mpfr_set_d(r26683, beta, MPFR_RNDN); mpfr_add(r26684, r26682, r26683, MPFR_RNDN); mpfr_add(r26685, r26684, r26681, MPFR_RNDN); mpfr_mul(r26686, r26681, r26685, MPFR_RNDN); mpfr_mul(r26687, r26683, r26682, MPFR_RNDN); mpfr_add(r26688, r26687, r26686, MPFR_RNDN); mpfr_mul(r26689, r26686, r26688, MPFR_RNDN); ; mpfr_mul(r26691, r26690, r26681, MPFR_RNDN); mpfr_add(r26692, r26684, r26691, MPFR_RNDN); mpfr_mul(r26693, r26692, r26692, MPFR_RNDN); mpfr_div(r26694, r26689, r26693, MPFR_RNDN); ; mpfr_sub(r26696, r26693, r26695, MPFR_RNDN); mpfr_div(r26697, r26694, r26696, MPFR_RNDN); return mpfr_get_d(r26697, MPFR_RNDN); } static mpfr_t r26698, r26699, r26700, r26701, r26702, r26703, r26704, r26705, r26706, r26707, r26708, r26709, r26710, r26711, r26712, r26713, r26714, r26715, r26716, r26717, r26718, r26719, r26720, r26721, r26722; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r26698); mpfr_init_set_str(r26699, "1.5126697473829053e+144", 10, MPFR_RNDN); mpfr_init(r26700); mpfr_init(r26701); mpfr_init(r26702); mpfr_init(r26703); mpfr_init(r26704); mpfr_init(r26705); mpfr_init(r26706); mpfr_init(r26707); mpfr_init(r26708); mpfr_init(r26709); mpfr_init(r26710); mpfr_init(r26711); mpfr_init(r26712); mpfr_init(r26713); mpfr_init(r26714); mpfr_init(r26715); mpfr_init(r26716); mpfr_init(r26717); mpfr_init_set_str(r26718, "1.0", 10, MPFR_RNDN); mpfr_init(r26719); mpfr_init(r26720); mpfr_init_set_str(r26721, "0", 10, MPFR_RNDN); mpfr_init(r26722); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r26698, beta, MPFR_RNDN); ; mpfr_set_si(r26700, mpfr_cmp(r26698, r26699) <= 0, MPFR_RNDN); mpfr_set_d(r26701, i, MPFR_RNDN); mpfr_set_d(r26702, alpha, MPFR_RNDN); mpfr_add(r26703, r26701, r26702, MPFR_RNDN); mpfr_add(r26704, r26703, r26698, MPFR_RNDN); mpfr_add(r26705, r26701, r26701, MPFR_RNDN); mpfr_add(r26706, r26702, r26698, MPFR_RNDN); mpfr_add(r26707, r26705, r26706, MPFR_RNDN); mpfr_div(r26708, r26707, r26701, MPFR_RNDN); mpfr_div(r26709, r26704, r26708, MPFR_RNDN); mpfr_mul(r26710, r26706, r26701, MPFR_RNDN); mpfr_mul(r26711, r26701, r26701, MPFR_RNDN); mpfr_mul(r26712, r26698, r26702, MPFR_RNDN); mpfr_add(r26713, r26711, r26712, MPFR_RNDN); mpfr_add(r26714, r26710, r26713, MPFR_RNDN); mpfr_div(r26715, r26714, r26707, MPFR_RNDN); mpfr_mul(r26716, r26709, r26715, MPFR_RNDN); mpfr_mul(r26717, r26707, r26707, MPFR_RNDN); ; mpfr_sub(r26719, r26717, r26718, MPFR_RNDN); mpfr_div(r26720, r26716, r26719, MPFR_RNDN); ; if (mpfr_get_si(r26700, MPFR_RNDN)) { mpfr_set(r26722, r26720, MPFR_RNDN); } else { mpfr_set(r26722, r26721, MPFR_RNDN); }; return mpfr_get_d(r26722, MPFR_RNDN); } static mpfr_t r26723, r26724, r26725, r26726, r26727, r26728, r26729, r26730, r26731, r26732, r26733, r26734, r26735, r26736, r26737, r26738, r26739, r26740, r26741, r26742, r26743, r26744, r26745, r26746, r26747; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r26723); mpfr_init_set_str(r26724, "1.5126697473829053e+144", 10, MPFR_RNDN); mpfr_init(r26725); mpfr_init(r26726); mpfr_init(r26727); mpfr_init(r26728); mpfr_init(r26729); mpfr_init(r26730); mpfr_init(r26731); mpfr_init(r26732); mpfr_init(r26733); mpfr_init(r26734); mpfr_init(r26735); mpfr_init(r26736); mpfr_init(r26737); mpfr_init(r26738); mpfr_init(r26739); mpfr_init(r26740); mpfr_init(r26741); mpfr_init(r26742); mpfr_init_set_str(r26743, "1.0", 10, MPFR_RNDN); mpfr_init(r26744); mpfr_init(r26745); mpfr_init_set_str(r26746, "0", 10, MPFR_RNDN); mpfr_init(r26747); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r26723, beta, MPFR_RNDN); ; mpfr_set_si(r26725, mpfr_cmp(r26723, r26724) <= 0, MPFR_RNDN); mpfr_set_d(r26726, i, MPFR_RNDN); mpfr_set_d(r26727, alpha, MPFR_RNDN); mpfr_add(r26728, r26726, r26727, MPFR_RNDN); mpfr_add(r26729, r26728, r26723, MPFR_RNDN); mpfr_add(r26730, r26726, r26726, MPFR_RNDN); mpfr_add(r26731, r26727, r26723, MPFR_RNDN); mpfr_add(r26732, r26730, r26731, MPFR_RNDN); mpfr_div(r26733, r26732, r26726, MPFR_RNDN); mpfr_div(r26734, r26729, r26733, MPFR_RNDN); mpfr_mul(r26735, r26731, r26726, MPFR_RNDN); mpfr_mul(r26736, r26726, r26726, MPFR_RNDN); mpfr_mul(r26737, r26723, r26727, MPFR_RNDN); mpfr_add(r26738, r26736, r26737, MPFR_RNDN); mpfr_add(r26739, r26735, r26738, MPFR_RNDN); mpfr_div(r26740, r26739, r26732, MPFR_RNDN); mpfr_mul(r26741, r26734, r26740, MPFR_RNDN); mpfr_mul(r26742, r26732, r26732, MPFR_RNDN); ; mpfr_sub(r26744, r26742, r26743, MPFR_RNDN); mpfr_div(r26745, r26741, r26744, MPFR_RNDN); ; if (mpfr_get_si(r26725, MPFR_RNDN)) { mpfr_set(r26747, r26745, MPFR_RNDN); } else { mpfr_set(r26747, r26746, MPFR_RNDN); }; return mpfr_get_d(r26747, MPFR_RNDN); }