#include #include #include #include #include char *name = "Octave 3.8, jcobi/4, as called"; double f_if(float i) { float r23402 = i; float r23403 = r23402 * r23402; float r23404 = r23403 * r23403; float r23405 = 2; float r23406 = r23405 * r23402; float r23407 = r23406 * r23406; float r23408 = r23404 / r23407; float r23409 = 1.0; float r23410 = r23407 - r23409; float r23411 = r23408 / r23410; return r23411; } double f_id(double i) { double r23412 = i; double r23413 = r23412 * r23412; double r23414 = r23413 * r23413; double r23415 = 2; double r23416 = r23415 * r23412; double r23417 = r23416 * r23416; double r23418 = r23414 / r23417; double r23419 = 1.0; double r23420 = r23417 - r23419; double r23421 = r23418 / r23420; return r23421; } double f_of(float i) { float r23422 = i; float r23423 = 1001.461792354028; bool r23424 = r23422 <= r23423; float r23425 = 1; float r23426 = 2; float r23427 = r23425 / r23426; float r23428 = r23427 * r23422; float r23429 = r23428 * r23428; float r23430 = r23426 * r23422; float r23431 = r23430 * r23430; float r23432 = 1.0; float r23433 = r23431 - r23432; float r23434 = r23429 / r23433; float r23435 = 0.015625; float r23436 = r23435 / r23422; float r23437 = r23436 / r23422; float r23438 = 0.00390625; float r23439 = 4; float r23440 = pow(r23422, r23439); float r23441 = r23438 / r23440; float r23442 = 1/16; float r23443 = r23441 + r23442; float r23444 = r23437 + r23443; float r23445 = r23424 ? r23434 : r23444; return r23445; } double f_od(double i) { double r23446 = i; double r23447 = 1001.461792354028; bool r23448 = r23446 <= r23447; double r23449 = 1; double r23450 = 2; double r23451 = r23449 / r23450; double r23452 = r23451 * r23446; double r23453 = r23452 * r23452; double r23454 = r23450 * r23446; double r23455 = r23454 * r23454; double r23456 = 1.0; double r23457 = r23455 - r23456; double r23458 = r23453 / r23457; double r23459 = 0.015625; double r23460 = r23459 / r23446; double r23461 = r23460 / r23446; double r23462 = 0.00390625; double r23463 = 4; double r23464 = pow(r23446, r23463); double r23465 = r23462 / r23464; double r23466 = 1/16; double r23467 = r23465 + r23466; double r23468 = r23461 + r23467; double r23469 = r23448 ? r23458 : r23468; return r23469; } 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 r23470, r23471, r23472, r23473, r23474, r23475, r23476, r23477, r23478, r23479; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r23470); mpfr_init(r23471); mpfr_init(r23472); mpfr_init_set_str(r23473, "2", 10, MPFR_RNDN); mpfr_init(r23474); mpfr_init(r23475); mpfr_init(r23476); mpfr_init_set_str(r23477, "1.0", 10, MPFR_RNDN); mpfr_init(r23478); mpfr_init(r23479); } double f_im(double i) { mpfr_set_d(r23470, i, MPFR_RNDN); mpfr_mul(r23471, r23470, r23470, MPFR_RNDN); mpfr_mul(r23472, r23471, r23471, MPFR_RNDN); ; mpfr_mul(r23474, r23473, r23470, MPFR_RNDN); mpfr_mul(r23475, r23474, r23474, MPFR_RNDN); mpfr_div(r23476, r23472, r23475, MPFR_RNDN); ; mpfr_sub(r23478, r23475, r23477, MPFR_RNDN); mpfr_div(r23479, r23476, r23478, MPFR_RNDN); return mpfr_get_d(r23479, MPFR_RNDN); } static mpfr_t r23480, r23481, r23482, r23483, r23484, r23485, r23486, r23487, r23488, r23489, r23490, r23491, r23492, r23493, r23494, r23495, r23496, r23497, r23498, r23499, r23500, r23501, r23502, r23503; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r23480); mpfr_init_set_str(r23481, "1001.461792354028", 10, MPFR_RNDN); mpfr_init(r23482); mpfr_init_set_str(r23483, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23484, "2", 10, MPFR_RNDN); mpfr_init(r23485); mpfr_init(r23486); mpfr_init(r23487); mpfr_init(r23488); mpfr_init(r23489); mpfr_init_set_str(r23490, "1.0", 10, MPFR_RNDN); mpfr_init(r23491); mpfr_init(r23492); mpfr_init_set_str(r23493, "0.015625", 10, MPFR_RNDN); mpfr_init(r23494); mpfr_init(r23495); mpfr_init_set_str(r23496, "0.00390625", 10, MPFR_RNDN); mpfr_init_set_str(r23497, "4", 10, MPFR_RNDN); mpfr_init(r23498); mpfr_init(r23499); mpfr_init_set_str(r23500, "1/16", 10, MPFR_RNDN); mpfr_init(r23501); mpfr_init(r23502); mpfr_init(r23503); } double f_fm(double i) { mpfr_set_d(r23480, i, MPFR_RNDN); ; mpfr_set_si(r23482, mpfr_cmp(r23480, r23481) <= 0, MPFR_RNDN); ; ; mpfr_div(r23485, r23483, r23484, MPFR_RNDN); mpfr_mul(r23486, r23485, r23480, MPFR_RNDN); mpfr_mul(r23487, r23486, r23486, MPFR_RNDN); mpfr_mul(r23488, r23484, r23480, MPFR_RNDN); mpfr_mul(r23489, r23488, r23488, MPFR_RNDN); ; mpfr_sub(r23491, r23489, r23490, MPFR_RNDN); mpfr_div(r23492, r23487, r23491, MPFR_RNDN); ; mpfr_div(r23494, r23493, r23480, MPFR_RNDN); mpfr_div(r23495, r23494, r23480, MPFR_RNDN); ; ; mpfr_pow(r23498, r23480, r23497, MPFR_RNDN); mpfr_div(r23499, r23496, r23498, MPFR_RNDN); ; mpfr_add(r23501, r23499, r23500, MPFR_RNDN); mpfr_add(r23502, r23495, r23501, MPFR_RNDN); if (mpfr_get_si(r23482, MPFR_RNDN)) { mpfr_set(r23503, r23492, MPFR_RNDN); } else { mpfr_set(r23503, r23502, MPFR_RNDN); }; return mpfr_get_d(r23503, MPFR_RNDN); } static mpfr_t r23504, r23505, r23506, r23507, r23508, r23509, r23510, r23511, r23512, r23513, r23514, r23515, r23516, r23517, r23518, r23519, r23520, r23521, r23522, r23523, r23524, r23525, r23526, r23527; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r23504); mpfr_init_set_str(r23505, "1001.461792354028", 10, MPFR_RNDN); mpfr_init(r23506); mpfr_init_set_str(r23507, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23508, "2", 10, MPFR_RNDN); mpfr_init(r23509); mpfr_init(r23510); mpfr_init(r23511); mpfr_init(r23512); mpfr_init(r23513); mpfr_init_set_str(r23514, "1.0", 10, MPFR_RNDN); mpfr_init(r23515); mpfr_init(r23516); mpfr_init_set_str(r23517, "0.015625", 10, MPFR_RNDN); mpfr_init(r23518); mpfr_init(r23519); mpfr_init_set_str(r23520, "0.00390625", 10, MPFR_RNDN); mpfr_init_set_str(r23521, "4", 10, MPFR_RNDN); mpfr_init(r23522); mpfr_init(r23523); mpfr_init_set_str(r23524, "1/16", 10, MPFR_RNDN); mpfr_init(r23525); mpfr_init(r23526); mpfr_init(r23527); } double f_dm(double i) { mpfr_set_d(r23504, i, MPFR_RNDN); ; mpfr_set_si(r23506, mpfr_cmp(r23504, r23505) <= 0, MPFR_RNDN); ; ; mpfr_div(r23509, r23507, r23508, MPFR_RNDN); mpfr_mul(r23510, r23509, r23504, MPFR_RNDN); mpfr_mul(r23511, r23510, r23510, MPFR_RNDN); mpfr_mul(r23512, r23508, r23504, MPFR_RNDN); mpfr_mul(r23513, r23512, r23512, MPFR_RNDN); ; mpfr_sub(r23515, r23513, r23514, MPFR_RNDN); mpfr_div(r23516, r23511, r23515, MPFR_RNDN); ; mpfr_div(r23518, r23517, r23504, MPFR_RNDN); mpfr_div(r23519, r23518, r23504, MPFR_RNDN); ; ; mpfr_pow(r23522, r23504, r23521, MPFR_RNDN); mpfr_div(r23523, r23520, r23522, MPFR_RNDN); ; mpfr_add(r23525, r23523, r23524, MPFR_RNDN); mpfr_add(r23526, r23519, r23525, MPFR_RNDN); if (mpfr_get_si(r23506, MPFR_RNDN)) { mpfr_set(r23527, r23516, MPFR_RNDN); } else { mpfr_set(r23527, r23526, MPFR_RNDN); }; return mpfr_get_d(r23527, MPFR_RNDN); }