#include #include #include #include #include char *name = "Octave 3.8, jcobi/4"; double f_if(float alpha, float beta, float i) { float r23038 = i; float r23039 = alpha; float r23040 = beta; float r23041 = r23039 + r23040; float r23042 = r23041 + r23038; float r23043 = r23038 * r23042; float r23044 = r23040 * r23039; float r23045 = r23044 + r23043; float r23046 = r23043 * r23045; float r23047 = 2; float r23048 = r23047 * r23038; float r23049 = r23041 + r23048; float r23050 = r23049 * r23049; float r23051 = r23046 / r23050; float r23052 = 1.0; float r23053 = r23050 - r23052; float r23054 = r23051 / r23053; return r23054; } double f_id(double alpha, double beta, double i) { double r23055 = i; double r23056 = alpha; double r23057 = beta; double r23058 = r23056 + r23057; double r23059 = r23058 + r23055; double r23060 = r23055 * r23059; double r23061 = r23057 * r23056; double r23062 = r23061 + r23060; double r23063 = r23060 * r23062; double r23064 = 2; double r23065 = r23064 * r23055; double r23066 = r23058 + r23065; double r23067 = r23066 * r23066; double r23068 = r23063 / r23067; double r23069 = 1.0; double r23070 = r23067 - r23069; double r23071 = r23068 / r23070; return r23071; } double f_of(float alpha, float beta, float i) { float r23072 = alpha; float r23073 = beta; float r23074 = r23072 + r23073; float r23075 = 1.3142276454815967e+154; bool r23076 = r23074 <= r23075; float r23077 = i; float r23078 = r23074 + r23077; float r23079 = r23077 * r23078; float r23080 = 2; float r23081 = r23080 * r23077; float r23082 = r23074 + r23081; float r23083 = r23079 / r23082; float r23084 = r23073 * r23072; float r23085 = r23084 + r23079; float r23086 = r23085 / r23082; float r23087 = r23083 * r23086; float r23088 = r23082 * r23082; float r23089 = 1.0; float r23090 = r23088 - r23089; float r23091 = r23087 / r23090; float r23092 = 0; float r23093 = r23076 ? r23091 : r23092; return r23093; } double f_od(double alpha, double beta, double i) { double r23094 = alpha; double r23095 = beta; double r23096 = r23094 + r23095; double r23097 = 1.3142276454815967e+154; bool r23098 = r23096 <= r23097; double r23099 = i; double r23100 = r23096 + r23099; double r23101 = r23099 * r23100; double r23102 = 2; double r23103 = r23102 * r23099; double r23104 = r23096 + r23103; double r23105 = r23101 / r23104; double r23106 = r23095 * r23094; double r23107 = r23106 + r23101; double r23108 = r23107 / r23104; double r23109 = r23105 * r23108; double r23110 = r23104 * r23104; double r23111 = 1.0; double r23112 = r23110 - r23111; double r23113 = r23109 / r23112; double r23114 = 0; double r23115 = r23098 ? r23113 : r23114; return r23115; } 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 r23116, r23117, r23118, r23119, r23120, r23121, r23122, r23123, r23124, r23125, r23126, r23127, r23128, r23129, r23130, r23131, r23132; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r23116); mpfr_init(r23117); mpfr_init(r23118); mpfr_init(r23119); mpfr_init(r23120); mpfr_init(r23121); mpfr_init(r23122); mpfr_init(r23123); mpfr_init(r23124); mpfr_init_set_str(r23125, "2", 10, MPFR_RNDN); mpfr_init(r23126); mpfr_init(r23127); mpfr_init(r23128); mpfr_init(r23129); mpfr_init_set_str(r23130, "1.0", 10, MPFR_RNDN); mpfr_init(r23131); mpfr_init(r23132); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r23116, i, MPFR_RNDN); mpfr_set_d(r23117, alpha, MPFR_RNDN); mpfr_set_d(r23118, beta, MPFR_RNDN); mpfr_add(r23119, r23117, r23118, MPFR_RNDN); mpfr_add(r23120, r23119, r23116, MPFR_RNDN); mpfr_mul(r23121, r23116, r23120, MPFR_RNDN); mpfr_mul(r23122, r23118, r23117, MPFR_RNDN); mpfr_add(r23123, r23122, r23121, MPFR_RNDN); mpfr_mul(r23124, r23121, r23123, MPFR_RNDN); ; mpfr_mul(r23126, r23125, r23116, MPFR_RNDN); mpfr_add(r23127, r23119, r23126, MPFR_RNDN); mpfr_mul(r23128, r23127, r23127, MPFR_RNDN); mpfr_div(r23129, r23124, r23128, MPFR_RNDN); ; mpfr_sub(r23131, r23128, r23130, MPFR_RNDN); mpfr_div(r23132, r23129, r23131, MPFR_RNDN); return mpfr_get_d(r23132, MPFR_RNDN); } static mpfr_t r23133, r23134, r23135, r23136, r23137, r23138, r23139, r23140, r23141, r23142, r23143, r23144, r23145, r23146, r23147, r23148, r23149, r23150, r23151, r23152, r23153, r23154; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r23133); mpfr_init(r23134); mpfr_init(r23135); mpfr_init_set_str(r23136, "1.3142276454815967e+154", 10, MPFR_RNDN); mpfr_init(r23137); mpfr_init(r23138); mpfr_init(r23139); mpfr_init(r23140); mpfr_init_set_str(r23141, "2", 10, MPFR_RNDN); mpfr_init(r23142); mpfr_init(r23143); mpfr_init(r23144); mpfr_init(r23145); mpfr_init(r23146); mpfr_init(r23147); mpfr_init(r23148); mpfr_init(r23149); mpfr_init_set_str(r23150, "1.0", 10, MPFR_RNDN); mpfr_init(r23151); mpfr_init(r23152); mpfr_init_set_str(r23153, "0", 10, MPFR_RNDN); mpfr_init(r23154); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r23133, alpha, MPFR_RNDN); mpfr_set_d(r23134, beta, MPFR_RNDN); mpfr_add(r23135, r23133, r23134, MPFR_RNDN); ; mpfr_set_si(r23137, mpfr_cmp(r23135, r23136) <= 0, MPFR_RNDN); mpfr_set_d(r23138, i, MPFR_RNDN); mpfr_add(r23139, r23135, r23138, MPFR_RNDN); mpfr_mul(r23140, r23138, r23139, MPFR_RNDN); ; mpfr_mul(r23142, r23141, r23138, MPFR_RNDN); mpfr_add(r23143, r23135, r23142, MPFR_RNDN); mpfr_div(r23144, r23140, r23143, MPFR_RNDN); mpfr_mul(r23145, r23134, r23133, MPFR_RNDN); mpfr_add(r23146, r23145, r23140, MPFR_RNDN); mpfr_div(r23147, r23146, r23143, MPFR_RNDN); mpfr_mul(r23148, r23144, r23147, MPFR_RNDN); mpfr_mul(r23149, r23143, r23143, MPFR_RNDN); ; mpfr_sub(r23151, r23149, r23150, MPFR_RNDN); mpfr_div(r23152, r23148, r23151, MPFR_RNDN); ; if (mpfr_get_si(r23137, MPFR_RNDN)) { mpfr_set(r23154, r23152, MPFR_RNDN); } else { mpfr_set(r23154, r23153, MPFR_RNDN); }; return mpfr_get_d(r23154, MPFR_RNDN); } static mpfr_t r23155, r23156, r23157, r23158, r23159, r23160, r23161, r23162, r23163, r23164, r23165, r23166, r23167, r23168, r23169, r23170, r23171, r23172, r23173, r23174, r23175, r23176; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r23155); mpfr_init(r23156); mpfr_init(r23157); mpfr_init_set_str(r23158, "1.3142276454815967e+154", 10, MPFR_RNDN); mpfr_init(r23159); mpfr_init(r23160); mpfr_init(r23161); mpfr_init(r23162); mpfr_init_set_str(r23163, "2", 10, MPFR_RNDN); mpfr_init(r23164); mpfr_init(r23165); mpfr_init(r23166); mpfr_init(r23167); mpfr_init(r23168); mpfr_init(r23169); mpfr_init(r23170); mpfr_init(r23171); mpfr_init_set_str(r23172, "1.0", 10, MPFR_RNDN); mpfr_init(r23173); mpfr_init(r23174); mpfr_init_set_str(r23175, "0", 10, MPFR_RNDN); mpfr_init(r23176); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r23155, alpha, MPFR_RNDN); mpfr_set_d(r23156, beta, MPFR_RNDN); mpfr_add(r23157, r23155, r23156, MPFR_RNDN); ; mpfr_set_si(r23159, mpfr_cmp(r23157, r23158) <= 0, MPFR_RNDN); mpfr_set_d(r23160, i, MPFR_RNDN); mpfr_add(r23161, r23157, r23160, MPFR_RNDN); mpfr_mul(r23162, r23160, r23161, MPFR_RNDN); ; mpfr_mul(r23164, r23163, r23160, MPFR_RNDN); mpfr_add(r23165, r23157, r23164, MPFR_RNDN); mpfr_div(r23166, r23162, r23165, MPFR_RNDN); mpfr_mul(r23167, r23156, r23155, MPFR_RNDN); mpfr_add(r23168, r23167, r23162, MPFR_RNDN); mpfr_div(r23169, r23168, r23165, MPFR_RNDN); mpfr_mul(r23170, r23166, r23169, MPFR_RNDN); mpfr_mul(r23171, r23165, r23165, MPFR_RNDN); ; mpfr_sub(r23173, r23171, r23172, MPFR_RNDN); mpfr_div(r23174, r23170, r23173, MPFR_RNDN); ; if (mpfr_get_si(r23159, MPFR_RNDN)) { mpfr_set(r23176, r23174, MPFR_RNDN); } else { mpfr_set(r23176, r23175, MPFR_RNDN); }; return mpfr_get_d(r23176, MPFR_RNDN); }