#include #include #include #include #include char *name = "Octave 3.8, jcobi/4"; double f_if(float alpha, float beta, float i) { float r23370 = i; float r23371 = alpha; float r23372 = beta; float r23373 = r23371 + r23372; float r23374 = r23373 + r23370; float r23375 = r23370 * r23374; float r23376 = r23372 * r23371; float r23377 = r23376 + r23375; float r23378 = r23375 * r23377; float r23379 = 2; float r23380 = r23379 * r23370; float r23381 = r23373 + r23380; float r23382 = r23381 * r23381; float r23383 = r23378 / r23382; float r23384 = 1.0; float r23385 = r23382 - r23384; float r23386 = r23383 / r23385; return r23386; } double f_id(double alpha, double beta, double i) { double r23387 = i; double r23388 = alpha; double r23389 = beta; double r23390 = r23388 + r23389; double r23391 = r23390 + r23387; double r23392 = r23387 * r23391; double r23393 = r23389 * r23388; double r23394 = r23393 + r23392; double r23395 = r23392 * r23394; double r23396 = 2; double r23397 = r23396 * r23387; double r23398 = r23390 + r23397; double r23399 = r23398 * r23398; double r23400 = r23395 / r23399; double r23401 = 1.0; double r23402 = r23399 - r23401; double r23403 = r23400 / r23402; return r23403; } double f_of(float alpha, float beta, float i) { float r23404 = alpha; float r23405 = 6.768179654449703e+148; bool r23406 = r23404 <= r23405; float r23407 = i; float r23408 = r23404 + r23407; float r23409 = beta; float r23410 = r23408 + r23409; float r23411 = r23409 * r23404; float r23412 = fma(r23410, r23407, r23411); float r23413 = 2; float r23414 = fma(r23407, r23413, r23404); float r23415 = r23414 + r23409; float r23416 = 1.0; float r23417 = -r23416; float r23418 = fma(r23415, r23415, r23417); float r23419 = r23412 / r23418; float r23420 = r23410 / r23415; float r23421 = r23420 * r23407; float r23422 = r23419 * r23421; float r23423 = r23422 / r23415; float r23424 = 0; float r23425 = r23406 ? r23423 : r23424; return r23425; } double f_od(double alpha, double beta, double i) { double r23426 = alpha; double r23427 = 6.768179654449703e+148; bool r23428 = r23426 <= r23427; double r23429 = i; double r23430 = r23426 + r23429; double r23431 = beta; double r23432 = r23430 + r23431; double r23433 = r23431 * r23426; double r23434 = fma(r23432, r23429, r23433); double r23435 = 2; double r23436 = fma(r23429, r23435, r23426); double r23437 = r23436 + r23431; double r23438 = 1.0; double r23439 = -r23438; double r23440 = fma(r23437, r23437, r23439); double r23441 = r23434 / r23440; double r23442 = r23432 / r23437; double r23443 = r23442 * r23429; double r23444 = r23441 * r23443; double r23445 = r23444 / r23437; double r23446 = 0; double r23447 = r23428 ? r23445 : r23446; return r23447; } 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 r23448, r23449, r23450, r23451, r23452, r23453, r23454, r23455, r23456, r23457, r23458, r23459, r23460, r23461, r23462, r23463, r23464; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r23448); mpfr_init(r23449); mpfr_init(r23450); mpfr_init(r23451); mpfr_init(r23452); mpfr_init(r23453); mpfr_init(r23454); mpfr_init(r23455); mpfr_init(r23456); mpfr_init_set_str(r23457, "2", 10, MPFR_RNDN); mpfr_init(r23458); mpfr_init(r23459); mpfr_init(r23460); mpfr_init(r23461); mpfr_init_set_str(r23462, "1.0", 10, MPFR_RNDN); mpfr_init(r23463); mpfr_init(r23464); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r23448, i, MPFR_RNDN); mpfr_set_d(r23449, alpha, MPFR_RNDN); mpfr_set_d(r23450, beta, MPFR_RNDN); mpfr_add(r23451, r23449, r23450, MPFR_RNDN); mpfr_add(r23452, r23451, r23448, MPFR_RNDN); mpfr_mul(r23453, r23448, r23452, MPFR_RNDN); mpfr_mul(r23454, r23450, r23449, MPFR_RNDN); mpfr_add(r23455, r23454, r23453, MPFR_RNDN); mpfr_mul(r23456, r23453, r23455, MPFR_RNDN); ; mpfr_mul(r23458, r23457, r23448, MPFR_RNDN); mpfr_add(r23459, r23451, r23458, MPFR_RNDN); mpfr_mul(r23460, r23459, r23459, MPFR_RNDN); mpfr_div(r23461, r23456, r23460, MPFR_RNDN); ; mpfr_sub(r23463, r23460, r23462, MPFR_RNDN); mpfr_div(r23464, r23461, r23463, MPFR_RNDN); return mpfr_get_d(r23464, MPFR_RNDN); } static mpfr_t r23465, r23466, r23467, r23468, r23469, r23470, r23471, r23472, r23473, r23474, r23475, r23476, r23477, r23478, r23479, r23480, r23481, r23482, r23483, r23484, r23485, r23486; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r23465); mpfr_init_set_str(r23466, "6.768179654449703e+148", 10, MPFR_RNDN); mpfr_init(r23467); mpfr_init(r23468); mpfr_init(r23469); mpfr_init(r23470); mpfr_init(r23471); mpfr_init(r23472); mpfr_init(r23473); mpfr_init_set_str(r23474, "2", 10, MPFR_RNDN); mpfr_init(r23475); mpfr_init(r23476); mpfr_init_set_str(r23477, "1.0", 10, MPFR_RNDN); mpfr_init(r23478); mpfr_init(r23479); mpfr_init(r23480); mpfr_init(r23481); mpfr_init(r23482); mpfr_init(r23483); mpfr_init(r23484); mpfr_init_set_str(r23485, "0", 10, MPFR_RNDN); mpfr_init(r23486); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r23465, alpha, MPFR_RNDN); ; mpfr_set_si(r23467, mpfr_cmp(r23465, r23466) <= 0, MPFR_RNDN); mpfr_set_d(r23468, i, MPFR_RNDN); mpfr_add(r23469, r23465, r23468, MPFR_RNDN); mpfr_set_d(r23470, beta, MPFR_RNDN); mpfr_add(r23471, r23469, r23470, MPFR_RNDN); mpfr_mul(r23472, r23470, r23465, MPFR_RNDN); mpfr_fma(r23473, r23471, r23468, r23472, MPFR_RNDN); ; mpfr_fma(r23475, r23468, r23474, r23465, MPFR_RNDN); mpfr_add(r23476, r23475, r23470, MPFR_RNDN); ; mpfr_neg(r23478, r23477, MPFR_RNDN); mpfr_fma(r23479, r23476, r23476, r23478, MPFR_RNDN); mpfr_div(r23480, r23473, r23479, MPFR_RNDN); mpfr_div(r23481, r23471, r23476, MPFR_RNDN); mpfr_mul(r23482, r23481, r23468, MPFR_RNDN); mpfr_mul(r23483, r23480, r23482, MPFR_RNDN); mpfr_div(r23484, r23483, r23476, MPFR_RNDN); ; if (mpfr_get_si(r23467, MPFR_RNDN)) { mpfr_set(r23486, r23484, MPFR_RNDN); } else { mpfr_set(r23486, r23485, MPFR_RNDN); }; return mpfr_get_d(r23486, MPFR_RNDN); } static mpfr_t r23487, r23488, r23489, r23490, r23491, r23492, r23493, r23494, r23495, r23496, r23497, r23498, r23499, r23500, r23501, r23502, r23503, r23504, r23505, r23506, r23507, r23508; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r23487); mpfr_init_set_str(r23488, "6.768179654449703e+148", 10, MPFR_RNDN); mpfr_init(r23489); mpfr_init(r23490); mpfr_init(r23491); mpfr_init(r23492); mpfr_init(r23493); mpfr_init(r23494); mpfr_init(r23495); mpfr_init_set_str(r23496, "2", 10, MPFR_RNDN); mpfr_init(r23497); mpfr_init(r23498); mpfr_init_set_str(r23499, "1.0", 10, MPFR_RNDN); mpfr_init(r23500); mpfr_init(r23501); mpfr_init(r23502); mpfr_init(r23503); mpfr_init(r23504); mpfr_init(r23505); mpfr_init(r23506); mpfr_init_set_str(r23507, "0", 10, MPFR_RNDN); mpfr_init(r23508); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r23487, alpha, MPFR_RNDN); ; mpfr_set_si(r23489, mpfr_cmp(r23487, r23488) <= 0, MPFR_RNDN); mpfr_set_d(r23490, i, MPFR_RNDN); mpfr_add(r23491, r23487, r23490, MPFR_RNDN); mpfr_set_d(r23492, beta, MPFR_RNDN); mpfr_add(r23493, r23491, r23492, MPFR_RNDN); mpfr_mul(r23494, r23492, r23487, MPFR_RNDN); mpfr_fma(r23495, r23493, r23490, r23494, MPFR_RNDN); ; mpfr_fma(r23497, r23490, r23496, r23487, MPFR_RNDN); mpfr_add(r23498, r23497, r23492, MPFR_RNDN); ; mpfr_neg(r23500, r23499, MPFR_RNDN); mpfr_fma(r23501, r23498, r23498, r23500, MPFR_RNDN); mpfr_div(r23502, r23495, r23501, MPFR_RNDN); mpfr_div(r23503, r23493, r23498, MPFR_RNDN); mpfr_mul(r23504, r23503, r23490, MPFR_RNDN); mpfr_mul(r23505, r23502, r23504, MPFR_RNDN); mpfr_div(r23506, r23505, r23498, MPFR_RNDN); ; if (mpfr_get_si(r23489, MPFR_RNDN)) { mpfr_set(r23508, r23506, MPFR_RNDN); } else { mpfr_set(r23508, r23507, MPFR_RNDN); }; return mpfr_get_d(r23508, MPFR_RNDN); }