#include #include #include #include #include char *name = "Octave 3.8, jcobi/3"; double f_if(float alpha, float beta) { float r26322 = alpha; float r26323 = beta; float r26324 = r26322 + r26323; float r26325 = r26323 * r26322; float r26326 = r26324 + r26325; float r26327 = 1.0; float r26328 = r26326 + r26327; float r26329 = 2; float r26330 = 1; float r26331 = r26329 * r26330; float r26332 = r26324 + r26331; float r26333 = r26328 / r26332; float r26334 = r26333 / r26332; float r26335 = r26332 + r26327; float r26336 = r26334 / r26335; return r26336; } double f_id(double alpha, double beta) { double r26337 = alpha; double r26338 = beta; double r26339 = r26337 + r26338; double r26340 = r26338 * r26337; double r26341 = r26339 + r26340; double r26342 = 1.0; double r26343 = r26341 + r26342; double r26344 = 2; double r26345 = 1; double r26346 = r26344 * r26345; double r26347 = r26339 + r26346; double r26348 = r26343 / r26347; double r26349 = r26348 / r26347; double r26350 = r26347 + r26342; double r26351 = r26349 / r26350; return r26351; } double f_of(float alpha, float beta) { float r26352 = beta; float r26353 = alpha; float r26354 = r26352 * r26353; float r26355 = r26352 + r26353; float r26356 = r26354 + r26355; float r26357 = 1.0; float r26358 = r26356 + r26357; float r26359 = 2; float r26360 = r26355 + r26359; float r26361 = r26358 / r26360; float r26362 = r26361 / r26360; float r26363 = r26357 + r26360; float r26364 = r26362 / r26363; float r26365 = 1.7835983677761644e+308; bool r26366 = r26364 <= r26365; float r26367 = 0.25; float r26368 = r26367 * r26355; float r26369 = 0.5; float r26370 = r26368 + r26369; float r26371 = r26357 + r26359; float r26372 = r26355 + r26371; float r26373 = r26360 * r26372; float r26374 = r26370 / r26373; float r26375 = r26366 ? r26364 : r26374; return r26375; } double f_od(double alpha, double beta) { double r26376 = beta; double r26377 = alpha; double r26378 = r26376 * r26377; double r26379 = r26376 + r26377; double r26380 = r26378 + r26379; double r26381 = 1.0; double r26382 = r26380 + r26381; double r26383 = 2; double r26384 = r26379 + r26383; double r26385 = r26382 / r26384; double r26386 = r26385 / r26384; double r26387 = r26381 + r26384; double r26388 = r26386 / r26387; double r26389 = 1.7835983677761644e+308; bool r26390 = r26388 <= r26389; double r26391 = 0.25; double r26392 = r26391 * r26379; double r26393 = 0.5; double r26394 = r26392 + r26393; double r26395 = r26381 + r26383; double r26396 = r26379 + r26395; double r26397 = r26384 * r26396; double r26398 = r26394 / r26397; double r26399 = r26390 ? r26388 : r26398; return r26399; } 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 r26400, r26401, r26402, r26403, r26404, r26405, r26406, r26407, r26408, r26409, r26410, r26411, r26412, r26413, r26414; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r26400); mpfr_init(r26401); mpfr_init(r26402); mpfr_init(r26403); mpfr_init(r26404); mpfr_init_set_str(r26405, "1.0", 10, MPFR_RNDN); mpfr_init(r26406); mpfr_init_set_str(r26407, "2", 10, MPFR_RNDN); mpfr_init_set_str(r26408, "1", 10, MPFR_RNDN); mpfr_init(r26409); mpfr_init(r26410); mpfr_init(r26411); mpfr_init(r26412); mpfr_init(r26413); mpfr_init(r26414); } double f_im(double alpha, double beta) { mpfr_set_d(r26400, alpha, MPFR_RNDN); mpfr_set_d(r26401, beta, MPFR_RNDN); mpfr_add(r26402, r26400, r26401, MPFR_RNDN); mpfr_mul(r26403, r26401, r26400, MPFR_RNDN); mpfr_add(r26404, r26402, r26403, MPFR_RNDN); ; mpfr_add(r26406, r26404, r26405, MPFR_RNDN); ; ; mpfr_mul(r26409, r26407, r26408, MPFR_RNDN); mpfr_add(r26410, r26402, r26409, MPFR_RNDN); mpfr_div(r26411, r26406, r26410, MPFR_RNDN); mpfr_div(r26412, r26411, r26410, MPFR_RNDN); mpfr_add(r26413, r26410, r26405, MPFR_RNDN); mpfr_div(r26414, r26412, r26413, MPFR_RNDN); return mpfr_get_d(r26414, MPFR_RNDN); } static mpfr_t r26415, r26416, r26417, r26418, r26419, r26420, r26421, r26422, r26423, r26424, r26425, r26426, r26427, r26428, r26429, r26430, r26431, r26432, r26433, r26434, r26435, r26436, r26437, r26438; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r26415); mpfr_init(r26416); mpfr_init(r26417); mpfr_init(r26418); mpfr_init(r26419); mpfr_init_set_str(r26420, "1.0", 10, MPFR_RNDN); mpfr_init(r26421); mpfr_init_set_str(r26422, "2", 10, MPFR_RNDN); mpfr_init(r26423); mpfr_init(r26424); mpfr_init(r26425); mpfr_init(r26426); mpfr_init(r26427); mpfr_init_set_str(r26428, "1.7835983677761644e+308", 10, MPFR_RNDN); mpfr_init(r26429); mpfr_init_set_str(r26430, "0.25", 10, MPFR_RNDN); mpfr_init(r26431); mpfr_init_set_str(r26432, "0.5", 10, MPFR_RNDN); mpfr_init(r26433); mpfr_init(r26434); mpfr_init(r26435); mpfr_init(r26436); mpfr_init(r26437); mpfr_init(r26438); } double f_fm(double alpha, double beta) { mpfr_set_d(r26415, beta, MPFR_RNDN); mpfr_set_d(r26416, alpha, MPFR_RNDN); mpfr_mul(r26417, r26415, r26416, MPFR_RNDN); mpfr_add(r26418, r26415, r26416, MPFR_RNDN); mpfr_add(r26419, r26417, r26418, MPFR_RNDN); ; mpfr_add(r26421, r26419, r26420, MPFR_RNDN); ; mpfr_add(r26423, r26418, r26422, MPFR_RNDN); mpfr_div(r26424, r26421, r26423, MPFR_RNDN); mpfr_div(r26425, r26424, r26423, MPFR_RNDN); mpfr_add(r26426, r26420, r26423, MPFR_RNDN); mpfr_div(r26427, r26425, r26426, MPFR_RNDN); ; mpfr_set_si(r26429, mpfr_cmp(r26427, r26428) <= 0, MPFR_RNDN); ; mpfr_mul(r26431, r26430, r26418, MPFR_RNDN); ; mpfr_add(r26433, r26431, r26432, MPFR_RNDN); mpfr_add(r26434, r26420, r26422, MPFR_RNDN); mpfr_add(r26435, r26418, r26434, MPFR_RNDN); mpfr_mul(r26436, r26423, r26435, MPFR_RNDN); mpfr_div(r26437, r26433, r26436, MPFR_RNDN); if (mpfr_get_si(r26429, MPFR_RNDN)) { mpfr_set(r26438, r26427, MPFR_RNDN); } else { mpfr_set(r26438, r26437, MPFR_RNDN); }; return mpfr_get_d(r26438, MPFR_RNDN); } static mpfr_t r26439, r26440, r26441, r26442, r26443, r26444, r26445, r26446, r26447, r26448, r26449, r26450, r26451, r26452, r26453, r26454, r26455, r26456, r26457, r26458, r26459, r26460, r26461, r26462; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r26439); mpfr_init(r26440); mpfr_init(r26441); mpfr_init(r26442); mpfr_init(r26443); mpfr_init_set_str(r26444, "1.0", 10, MPFR_RNDN); mpfr_init(r26445); mpfr_init_set_str(r26446, "2", 10, MPFR_RNDN); mpfr_init(r26447); mpfr_init(r26448); mpfr_init(r26449); mpfr_init(r26450); mpfr_init(r26451); mpfr_init_set_str(r26452, "1.7835983677761644e+308", 10, MPFR_RNDN); mpfr_init(r26453); mpfr_init_set_str(r26454, "0.25", 10, MPFR_RNDN); mpfr_init(r26455); mpfr_init_set_str(r26456, "0.5", 10, MPFR_RNDN); mpfr_init(r26457); mpfr_init(r26458); mpfr_init(r26459); mpfr_init(r26460); mpfr_init(r26461); mpfr_init(r26462); } double f_dm(double alpha, double beta) { mpfr_set_d(r26439, beta, MPFR_RNDN); mpfr_set_d(r26440, alpha, MPFR_RNDN); mpfr_mul(r26441, r26439, r26440, MPFR_RNDN); mpfr_add(r26442, r26439, r26440, MPFR_RNDN); mpfr_add(r26443, r26441, r26442, MPFR_RNDN); ; mpfr_add(r26445, r26443, r26444, MPFR_RNDN); ; mpfr_add(r26447, r26442, r26446, MPFR_RNDN); mpfr_div(r26448, r26445, r26447, MPFR_RNDN); mpfr_div(r26449, r26448, r26447, MPFR_RNDN); mpfr_add(r26450, r26444, r26447, MPFR_RNDN); mpfr_div(r26451, r26449, r26450, MPFR_RNDN); ; mpfr_set_si(r26453, mpfr_cmp(r26451, r26452) <= 0, MPFR_RNDN); ; mpfr_mul(r26455, r26454, r26442, MPFR_RNDN); ; mpfr_add(r26457, r26455, r26456, MPFR_RNDN); mpfr_add(r26458, r26444, r26446, MPFR_RNDN); mpfr_add(r26459, r26442, r26458, MPFR_RNDN); mpfr_mul(r26460, r26447, r26459, MPFR_RNDN); mpfr_div(r26461, r26457, r26460, MPFR_RNDN); if (mpfr_get_si(r26453, MPFR_RNDN)) { mpfr_set(r26462, r26451, MPFR_RNDN); } else { mpfr_set(r26462, r26461, MPFR_RNDN); }; return mpfr_get_d(r26462, MPFR_RNDN); }