#include #include #include #include #include char *name = "Octave 3.8, jcobi/4"; double f_if(float alpha, float beta, float i) { float r26629 = i; float r26630 = alpha; float r26631 = beta; float r26632 = r26630 + r26631; float r26633 = r26632 + r26629; float r26634 = r26629 * r26633; float r26635 = r26631 * r26630; float r26636 = r26635 + r26634; float r26637 = r26634 * r26636; float r26638 = 2; float r26639 = r26638 * r26629; float r26640 = r26632 + r26639; float r26641 = r26640 * r26640; float r26642 = r26637 / r26641; float r26643 = 1.0; float r26644 = r26641 - r26643; float r26645 = r26642 / r26644; return r26645; } double f_id(double alpha, double beta, double i) { double r26646 = i; double r26647 = alpha; double r26648 = beta; double r26649 = r26647 + r26648; double r26650 = r26649 + r26646; double r26651 = r26646 * r26650; double r26652 = r26648 * r26647; double r26653 = r26652 + r26651; double r26654 = r26651 * r26653; double r26655 = 2; double r26656 = r26655 * r26646; double r26657 = r26649 + r26656; double r26658 = r26657 * r26657; double r26659 = r26654 / r26658; double r26660 = 1.0; double r26661 = r26658 - r26660; double r26662 = r26659 / r26661; return r26662; } double f_of(float alpha, float beta, float i) { float r26663 = beta; float r26664 = 1.5126697473829053e+144; bool r26665 = r26663 <= r26664; float r26666 = i; float r26667 = alpha; float r26668 = r26666 + r26667; float r26669 = r26668 + r26663; float r26670 = r26666 + r26666; float r26671 = r26667 + r26663; float r26672 = r26670 + r26671; float r26673 = r26672 / r26666; float r26674 = r26669 / r26673; float r26675 = r26671 * r26666; float r26676 = r26666 * r26666; float r26677 = r26663 * r26667; float r26678 = r26676 + r26677; float r26679 = r26675 + r26678; float r26680 = r26679 / r26672; float r26681 = r26674 * r26680; float r26682 = r26672 * r26672; float r26683 = 1.0; float r26684 = r26682 - r26683; float r26685 = r26681 / r26684; float r26686 = 0; float r26687 = r26665 ? r26685 : r26686; return r26687; } double f_od(double alpha, double beta, double i) { double r26688 = beta; double r26689 = 1.5126697473829053e+144; bool r26690 = r26688 <= r26689; double r26691 = i; double r26692 = alpha; double r26693 = r26691 + r26692; double r26694 = r26693 + r26688; double r26695 = r26691 + r26691; double r26696 = r26692 + r26688; double r26697 = r26695 + r26696; double r26698 = r26697 / r26691; double r26699 = r26694 / r26698; double r26700 = r26696 * r26691; double r26701 = r26691 * r26691; double r26702 = r26688 * r26692; double r26703 = r26701 + r26702; double r26704 = r26700 + r26703; double r26705 = r26704 / r26697; double r26706 = r26699 * r26705; double r26707 = r26697 * r26697; double r26708 = 1.0; double r26709 = r26707 - r26708; double r26710 = r26706 / r26709; double r26711 = 0; double r26712 = r26690 ? r26710 : r26711; return r26712; } 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 r26713, r26714, r26715, r26716, r26717, r26718, r26719, r26720, r26721, r26722, r26723, r26724, r26725, r26726, r26727, r26728, r26729; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r26713); mpfr_init(r26714); mpfr_init(r26715); mpfr_init(r26716); mpfr_init(r26717); mpfr_init(r26718); mpfr_init(r26719); mpfr_init(r26720); mpfr_init(r26721); mpfr_init_set_str(r26722, "2", 10, MPFR_RNDN); mpfr_init(r26723); mpfr_init(r26724); mpfr_init(r26725); mpfr_init(r26726); mpfr_init_set_str(r26727, "1.0", 10, MPFR_RNDN); mpfr_init(r26728); mpfr_init(r26729); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r26713, i, MPFR_RNDN); mpfr_set_d(r26714, alpha, MPFR_RNDN); mpfr_set_d(r26715, beta, MPFR_RNDN); mpfr_add(r26716, r26714, r26715, MPFR_RNDN); mpfr_add(r26717, r26716, r26713, MPFR_RNDN); mpfr_mul(r26718, r26713, r26717, MPFR_RNDN); mpfr_mul(r26719, r26715, r26714, MPFR_RNDN); mpfr_add(r26720, r26719, r26718, MPFR_RNDN); mpfr_mul(r26721, r26718, r26720, MPFR_RNDN); ; mpfr_mul(r26723, r26722, r26713, MPFR_RNDN); mpfr_add(r26724, r26716, r26723, MPFR_RNDN); mpfr_mul(r26725, r26724, r26724, MPFR_RNDN); mpfr_div(r26726, r26721, r26725, MPFR_RNDN); ; mpfr_sub(r26728, r26725, r26727, MPFR_RNDN); mpfr_div(r26729, r26726, r26728, MPFR_RNDN); return mpfr_get_d(r26729, MPFR_RNDN); } static mpfr_t r26730, r26731, r26732, r26733, r26734, r26735, r26736, r26737, r26738, r26739, r26740, r26741, r26742, r26743, r26744, r26745, r26746, r26747, r26748, r26749, r26750, r26751, r26752, r26753, r26754; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r26730); mpfr_init_set_str(r26731, "1.5126697473829053e+144", 10, MPFR_RNDN); 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(r26743); mpfr_init(r26744); mpfr_init(r26745); mpfr_init(r26746); mpfr_init(r26747); mpfr_init(r26748); mpfr_init(r26749); mpfr_init_set_str(r26750, "1.0", 10, MPFR_RNDN); mpfr_init(r26751); mpfr_init(r26752); mpfr_init_set_str(r26753, "0", 10, MPFR_RNDN); mpfr_init(r26754); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r26730, beta, MPFR_RNDN); ; mpfr_set_si(r26732, mpfr_cmp(r26730, r26731) <= 0, MPFR_RNDN); mpfr_set_d(r26733, i, MPFR_RNDN); mpfr_set_d(r26734, alpha, MPFR_RNDN); mpfr_add(r26735, r26733, r26734, MPFR_RNDN); mpfr_add(r26736, r26735, r26730, MPFR_RNDN); mpfr_add(r26737, r26733, r26733, MPFR_RNDN); mpfr_add(r26738, r26734, r26730, MPFR_RNDN); mpfr_add(r26739, r26737, r26738, MPFR_RNDN); mpfr_div(r26740, r26739, r26733, MPFR_RNDN); mpfr_div(r26741, r26736, r26740, MPFR_RNDN); mpfr_mul(r26742, r26738, r26733, MPFR_RNDN); mpfr_mul(r26743, r26733, r26733, MPFR_RNDN); mpfr_mul(r26744, r26730, r26734, MPFR_RNDN); mpfr_add(r26745, r26743, r26744, MPFR_RNDN); mpfr_add(r26746, r26742, r26745, MPFR_RNDN); mpfr_div(r26747, r26746, r26739, MPFR_RNDN); mpfr_mul(r26748, r26741, r26747, MPFR_RNDN); mpfr_mul(r26749, r26739, r26739, MPFR_RNDN); ; mpfr_sub(r26751, r26749, r26750, MPFR_RNDN); mpfr_div(r26752, r26748, r26751, MPFR_RNDN); ; if (mpfr_get_si(r26732, MPFR_RNDN)) { mpfr_set(r26754, r26752, MPFR_RNDN); } else { mpfr_set(r26754, r26753, MPFR_RNDN); }; return mpfr_get_d(r26754, MPFR_RNDN); } static mpfr_t r26755, r26756, r26757, r26758, r26759, r26760, r26761, r26762, r26763, r26764, r26765, r26766, r26767, r26768, r26769, r26770, r26771, r26772, r26773, r26774, r26775, r26776, r26777, r26778, r26779; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r26755); mpfr_init_set_str(r26756, "1.5126697473829053e+144", 10, MPFR_RNDN); mpfr_init(r26757); mpfr_init(r26758); mpfr_init(r26759); mpfr_init(r26760); mpfr_init(r26761); mpfr_init(r26762); mpfr_init(r26763); mpfr_init(r26764); mpfr_init(r26765); mpfr_init(r26766); mpfr_init(r26767); mpfr_init(r26768); mpfr_init(r26769); mpfr_init(r26770); mpfr_init(r26771); mpfr_init(r26772); mpfr_init(r26773); mpfr_init(r26774); mpfr_init_set_str(r26775, "1.0", 10, MPFR_RNDN); mpfr_init(r26776); mpfr_init(r26777); mpfr_init_set_str(r26778, "0", 10, MPFR_RNDN); mpfr_init(r26779); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r26755, beta, MPFR_RNDN); ; mpfr_set_si(r26757, mpfr_cmp(r26755, r26756) <= 0, MPFR_RNDN); mpfr_set_d(r26758, i, MPFR_RNDN); mpfr_set_d(r26759, alpha, MPFR_RNDN); mpfr_add(r26760, r26758, r26759, MPFR_RNDN); mpfr_add(r26761, r26760, r26755, MPFR_RNDN); mpfr_add(r26762, r26758, r26758, MPFR_RNDN); mpfr_add(r26763, r26759, r26755, MPFR_RNDN); mpfr_add(r26764, r26762, r26763, MPFR_RNDN); mpfr_div(r26765, r26764, r26758, MPFR_RNDN); mpfr_div(r26766, r26761, r26765, MPFR_RNDN); mpfr_mul(r26767, r26763, r26758, MPFR_RNDN); mpfr_mul(r26768, r26758, r26758, MPFR_RNDN); mpfr_mul(r26769, r26755, r26759, MPFR_RNDN); mpfr_add(r26770, r26768, r26769, MPFR_RNDN); mpfr_add(r26771, r26767, r26770, MPFR_RNDN); mpfr_div(r26772, r26771, r26764, MPFR_RNDN); mpfr_mul(r26773, r26766, r26772, MPFR_RNDN); mpfr_mul(r26774, r26764, r26764, MPFR_RNDN); ; mpfr_sub(r26776, r26774, r26775, MPFR_RNDN); mpfr_div(r26777, r26773, r26776, MPFR_RNDN); ; if (mpfr_get_si(r26757, MPFR_RNDN)) { mpfr_set(r26779, r26777, MPFR_RNDN); } else { mpfr_set(r26779, r26778, MPFR_RNDN); }; return mpfr_get_d(r26779, MPFR_RNDN); }