#include #include #include #include #include char *name = "Octave 3.8, jcobi/1"; double f_if(float alpha, float beta) { float r24309 = beta; float r24310 = alpha; float r24311 = r24309 - r24310; float r24312 = r24310 + r24309; float r24313 = 2.0; float r24314 = r24312 + r24313; float r24315 = r24311 / r24314; float r24316 = 1.0; float r24317 = r24315 + r24316; float r24318 = r24317 / r24313; return r24318; } double f_id(double alpha, double beta) { double r24319 = beta; double r24320 = alpha; double r24321 = r24319 - r24320; double r24322 = r24320 + r24319; double r24323 = 2.0; double r24324 = r24322 + r24323; double r24325 = r24321 / r24324; double r24326 = 1.0; double r24327 = r24325 + r24326; double r24328 = r24327 / r24323; return r24328; } double f_of(float alpha, float beta) { float r24329 = beta; float r24330 = alpha; float r24331 = r24329 - r24330; float r24332 = r24330 + r24329; float r24333 = 2.0; float r24334 = r24332 + r24333; float r24335 = r24331 / r24334; float r24336 = 1.0; float r24337 = r24335 + r24336; float r24338 = r24337 / r24333; float r24339 = 1.0824673379872252e-09; bool r24340 = r24338 <= r24339; float r24341 = r24329 / r24334; float r24342 = 4.0; float r24343 = 1; float r24344 = 2; float r24345 = pow(r24330, r24344); float r24346 = r24343 / r24345; float r24347 = r24342 * r24346; float r24348 = 8.0; float r24349 = 3; float r24350 = pow(r24330, r24349); float r24351 = r24343 / r24350; float r24352 = r24348 * r24351; float r24353 = r24343 / r24330; float r24354 = r24333 * r24353; float r24355 = r24352 + r24354; float r24356 = r24347 - r24355; float r24357 = r24341 - r24356; float r24358 = r24357 / r24333; float r24359 = r24330 / r24334; float r24360 = r24359 - r24336; float r24361 = exp(r24360); float r24362 = log(r24361); float r24363 = r24341 - r24362; float r24364 = r24363 / r24333; float r24365 = r24340 ? r24358 : r24364; return r24365; } double f_od(double alpha, double beta) { double r24366 = beta; double r24367 = alpha; double r24368 = r24366 - r24367; double r24369 = r24367 + r24366; double r24370 = 2.0; double r24371 = r24369 + r24370; double r24372 = r24368 / r24371; double r24373 = 1.0; double r24374 = r24372 + r24373; double r24375 = r24374 / r24370; double r24376 = 1.0824673379872252e-09; bool r24377 = r24375 <= r24376; double r24378 = r24366 / r24371; double r24379 = 4.0; double r24380 = 1; double r24381 = 2; double r24382 = pow(r24367, r24381); double r24383 = r24380 / r24382; double r24384 = r24379 * r24383; double r24385 = 8.0; double r24386 = 3; double r24387 = pow(r24367, r24386); double r24388 = r24380 / r24387; double r24389 = r24385 * r24388; double r24390 = r24380 / r24367; double r24391 = r24370 * r24390; double r24392 = r24389 + r24391; double r24393 = r24384 - r24392; double r24394 = r24378 - r24393; double r24395 = r24394 / r24370; double r24396 = r24367 / r24371; double r24397 = r24396 - r24373; double r24398 = exp(r24397); double r24399 = log(r24398); double r24400 = r24378 - r24399; double r24401 = r24400 / r24370; double r24402 = r24377 ? r24395 : r24401; return r24402; } 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 r24403, r24404, r24405, r24406, r24407, r24408, r24409, r24410, r24411, r24412; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r24403); mpfr_init(r24404); mpfr_init(r24405); mpfr_init(r24406); mpfr_init_set_str(r24407, "2.0", 10, MPFR_RNDN); mpfr_init(r24408); mpfr_init(r24409); mpfr_init_set_str(r24410, "1.0", 10, MPFR_RNDN); mpfr_init(r24411); mpfr_init(r24412); } double f_im(double alpha, double beta) { mpfr_set_d(r24403, beta, MPFR_RNDN); mpfr_set_d(r24404, alpha, MPFR_RNDN); mpfr_sub(r24405, r24403, r24404, MPFR_RNDN); mpfr_add(r24406, r24404, r24403, MPFR_RNDN); ; mpfr_add(r24408, r24406, r24407, MPFR_RNDN); mpfr_div(r24409, r24405, r24408, MPFR_RNDN); ; mpfr_add(r24411, r24409, r24410, MPFR_RNDN); mpfr_div(r24412, r24411, r24407, MPFR_RNDN); return mpfr_get_d(r24412, MPFR_RNDN); } static mpfr_t r24413, r24414, r24415, r24416, r24417, r24418, r24419, r24420, r24421, r24422, r24423, r24424, r24425, r24426, r24427, r24428, r24429, r24430, r24431, r24432, r24433, r24434, r24435, r24436, r24437, r24438, r24439, r24440, r24441, r24442, r24443, r24444, r24445, r24446, r24447, r24448, r24449; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r24413); mpfr_init(r24414); mpfr_init(r24415); mpfr_init(r24416); mpfr_init_set_str(r24417, "2.0", 10, MPFR_RNDN); mpfr_init(r24418); mpfr_init(r24419); mpfr_init_set_str(r24420, "1.0", 10, MPFR_RNDN); mpfr_init(r24421); mpfr_init(r24422); mpfr_init_set_str(r24423, "1.0824673379872252e-09", 10, MPFR_RNDN); mpfr_init(r24424); mpfr_init(r24425); mpfr_init_set_str(r24426, "4.0", 10, MPFR_RNDN); mpfr_init_set_str(r24427, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24428, "2", 10, MPFR_RNDN); mpfr_init(r24429); mpfr_init(r24430); mpfr_init(r24431); mpfr_init_set_str(r24432, "8.0", 10, MPFR_RNDN); mpfr_init_set_str(r24433, "3", 10, MPFR_RNDN); mpfr_init(r24434); mpfr_init(r24435); mpfr_init(r24436); mpfr_init(r24437); mpfr_init(r24438); mpfr_init(r24439); mpfr_init(r24440); mpfr_init(r24441); mpfr_init(r24442); mpfr_init(r24443); mpfr_init(r24444); mpfr_init(r24445); mpfr_init(r24446); mpfr_init(r24447); mpfr_init(r24448); mpfr_init(r24449); } double f_fm(double alpha, double beta) { mpfr_set_d(r24413, beta, MPFR_RNDN); mpfr_set_d(r24414, alpha, MPFR_RNDN); mpfr_sub(r24415, r24413, r24414, MPFR_RNDN); mpfr_add(r24416, r24414, r24413, MPFR_RNDN); ; mpfr_add(r24418, r24416, r24417, MPFR_RNDN); mpfr_div(r24419, r24415, r24418, MPFR_RNDN); ; mpfr_add(r24421, r24419, r24420, MPFR_RNDN); mpfr_div(r24422, r24421, r24417, MPFR_RNDN); ; mpfr_set_si(r24424, mpfr_cmp(r24422, r24423) <= 0, MPFR_RNDN); mpfr_div(r24425, r24413, r24418, MPFR_RNDN); ; ; ; mpfr_pow(r24429, r24414, r24428, MPFR_RNDN); mpfr_div(r24430, r24427, r24429, MPFR_RNDN); mpfr_mul(r24431, r24426, r24430, MPFR_RNDN); ; ; mpfr_pow(r24434, r24414, r24433, MPFR_RNDN); mpfr_div(r24435, r24427, r24434, MPFR_RNDN); mpfr_mul(r24436, r24432, r24435, MPFR_RNDN); mpfr_div(r24437, r24427, r24414, MPFR_RNDN); mpfr_mul(r24438, r24417, r24437, MPFR_RNDN); mpfr_add(r24439, r24436, r24438, MPFR_RNDN); mpfr_sub(r24440, r24431, r24439, MPFR_RNDN); mpfr_sub(r24441, r24425, r24440, MPFR_RNDN); mpfr_div(r24442, r24441, r24417, MPFR_RNDN); mpfr_div(r24443, r24414, r24418, MPFR_RNDN); mpfr_sub(r24444, r24443, r24420, MPFR_RNDN); mpfr_exp(r24445, r24444, MPFR_RNDN); mpfr_log(r24446, r24445, MPFR_RNDN); mpfr_sub(r24447, r24425, r24446, MPFR_RNDN); mpfr_div(r24448, r24447, r24417, MPFR_RNDN); if (mpfr_get_si(r24424, MPFR_RNDN)) { mpfr_set(r24449, r24442, MPFR_RNDN); } else { mpfr_set(r24449, r24448, MPFR_RNDN); }; return mpfr_get_d(r24449, MPFR_RNDN); } static mpfr_t r24450, r24451, r24452, r24453, r24454, r24455, r24456, r24457, r24458, r24459, r24460, r24461, r24462, r24463, r24464, r24465, r24466, r24467, r24468, r24469, r24470, r24471, r24472, r24473, r24474, r24475, r24476, r24477, r24478, r24479, r24480, r24481, r24482, r24483, r24484, r24485, r24486; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r24450); mpfr_init(r24451); mpfr_init(r24452); mpfr_init(r24453); mpfr_init_set_str(r24454, "2.0", 10, MPFR_RNDN); mpfr_init(r24455); mpfr_init(r24456); mpfr_init_set_str(r24457, "1.0", 10, MPFR_RNDN); mpfr_init(r24458); mpfr_init(r24459); mpfr_init_set_str(r24460, "1.0824673379872252e-09", 10, MPFR_RNDN); mpfr_init(r24461); mpfr_init(r24462); mpfr_init_set_str(r24463, "4.0", 10, MPFR_RNDN); mpfr_init_set_str(r24464, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24465, "2", 10, MPFR_RNDN); mpfr_init(r24466); mpfr_init(r24467); mpfr_init(r24468); mpfr_init_set_str(r24469, "8.0", 10, MPFR_RNDN); mpfr_init_set_str(r24470, "3", 10, MPFR_RNDN); mpfr_init(r24471); mpfr_init(r24472); mpfr_init(r24473); mpfr_init(r24474); mpfr_init(r24475); mpfr_init(r24476); mpfr_init(r24477); mpfr_init(r24478); mpfr_init(r24479); mpfr_init(r24480); mpfr_init(r24481); mpfr_init(r24482); mpfr_init(r24483); mpfr_init(r24484); mpfr_init(r24485); mpfr_init(r24486); } double f_dm(double alpha, double beta) { mpfr_set_d(r24450, beta, MPFR_RNDN); mpfr_set_d(r24451, alpha, MPFR_RNDN); mpfr_sub(r24452, r24450, r24451, MPFR_RNDN); mpfr_add(r24453, r24451, r24450, MPFR_RNDN); ; mpfr_add(r24455, r24453, r24454, MPFR_RNDN); mpfr_div(r24456, r24452, r24455, MPFR_RNDN); ; mpfr_add(r24458, r24456, r24457, MPFR_RNDN); mpfr_div(r24459, r24458, r24454, MPFR_RNDN); ; mpfr_set_si(r24461, mpfr_cmp(r24459, r24460) <= 0, MPFR_RNDN); mpfr_div(r24462, r24450, r24455, MPFR_RNDN); ; ; ; mpfr_pow(r24466, r24451, r24465, MPFR_RNDN); mpfr_div(r24467, r24464, r24466, MPFR_RNDN); mpfr_mul(r24468, r24463, r24467, MPFR_RNDN); ; ; mpfr_pow(r24471, r24451, r24470, MPFR_RNDN); mpfr_div(r24472, r24464, r24471, MPFR_RNDN); mpfr_mul(r24473, r24469, r24472, MPFR_RNDN); mpfr_div(r24474, r24464, r24451, MPFR_RNDN); mpfr_mul(r24475, r24454, r24474, MPFR_RNDN); mpfr_add(r24476, r24473, r24475, MPFR_RNDN); mpfr_sub(r24477, r24468, r24476, MPFR_RNDN); mpfr_sub(r24478, r24462, r24477, MPFR_RNDN); mpfr_div(r24479, r24478, r24454, MPFR_RNDN); mpfr_div(r24480, r24451, r24455, MPFR_RNDN); mpfr_sub(r24481, r24480, r24457, MPFR_RNDN); mpfr_exp(r24482, r24481, MPFR_RNDN); mpfr_log(r24483, r24482, MPFR_RNDN); mpfr_sub(r24484, r24462, r24483, MPFR_RNDN); mpfr_div(r24485, r24484, r24454, MPFR_RNDN); if (mpfr_get_si(r24461, MPFR_RNDN)) { mpfr_set(r24486, r24479, MPFR_RNDN); } else { mpfr_set(r24486, r24485, MPFR_RNDN); }; return mpfr_get_d(r24486, MPFR_RNDN); }