#include #include #include #include #include char *name = "Octave 3.8, jcobi/3"; double f_if(float alpha, float beta) { float r10311 = alpha; float r10312 = beta; float r10313 = r10311 + r10312; float r10314 = r10312 * r10311; float r10315 = r10313 + r10314; float r10316 = 1.0; float r10317 = r10315 + r10316; float r10318 = 2; float r10319 = 1; float r10320 = r10318 * r10319; float r10321 = r10313 + r10320; float r10322 = r10317 / r10321; float r10323 = r10322 / r10321; float r10324 = r10321 + r10316; float r10325 = r10323 / r10324; return r10325; } double f_id(double alpha, double beta) { double r10326 = alpha; double r10327 = beta; double r10328 = r10326 + r10327; double r10329 = r10327 * r10326; double r10330 = r10328 + r10329; double r10331 = 1.0; double r10332 = r10330 + r10331; double r10333 = 2; double r10334 = 1; double r10335 = r10333 * r10334; double r10336 = r10328 + r10335; double r10337 = r10332 / r10336; double r10338 = r10337 / r10336; double r10339 = r10336 + r10331; double r10340 = r10338 / r10339; return r10340; } double f_of(float alpha, float beta) { float r10341 = alpha; float r10342 = 3.2400543536458e+168; bool r10343 = r10341 <= r10342; float r10344 = 1.0; float r10345 = beta; float r10346 = r10345 * r10341; float r10347 = r10345 + r10341; float r10348 = r10346 + r10347; float r10349 = r10344 + r10348; float r10350 = 2; float r10351 = r10347 + r10350; float r10352 = r10349 / r10351; float r10353 = r10352 / r10351; float r10354 = r10344 + r10351; float r10355 = r10353 / r10354; float r10356 = 0.5; float r10357 = 0.25; float r10358 = r10347 * r10357; float r10359 = r10356 + r10358; float r10360 = r10341 + r10350; float r10361 = r10345 + r10360; float r10362 = r10345 + r10350; float r10363 = r10341 + r10344; float r10364 = r10362 + r10363; float r10365 = r10361 * r10364; float r10366 = r10359 / r10365; float r10367 = r10343 ? r10355 : r10366; return r10367; } double f_od(double alpha, double beta) { double r10368 = alpha; double r10369 = 3.2400543536458e+168; bool r10370 = r10368 <= r10369; double r10371 = 1.0; double r10372 = beta; double r10373 = r10372 * r10368; double r10374 = r10372 + r10368; double r10375 = r10373 + r10374; double r10376 = r10371 + r10375; double r10377 = 2; double r10378 = r10374 + r10377; double r10379 = r10376 / r10378; double r10380 = r10379 / r10378; double r10381 = r10371 + r10378; double r10382 = r10380 / r10381; double r10383 = 0.5; double r10384 = 0.25; double r10385 = r10374 * r10384; double r10386 = r10383 + r10385; double r10387 = r10368 + r10377; double r10388 = r10372 + r10387; double r10389 = r10372 + r10377; double r10390 = r10368 + r10371; double r10391 = r10389 + r10390; double r10392 = r10388 * r10391; double r10393 = r10386 / r10392; double r10394 = r10370 ? r10382 : r10393; return r10394; } 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 r10395, r10396, r10397, r10398, r10399, r10400, r10401, r10402, r10403, r10404, r10405, r10406, r10407, r10408, r10409; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r10395); mpfr_init(r10396); mpfr_init(r10397); mpfr_init(r10398); mpfr_init(r10399); mpfr_init_set_str(r10400, "1.0", 10, MPFR_RNDN); mpfr_init(r10401); mpfr_init_set_str(r10402, "2", 10, MPFR_RNDN); mpfr_init_set_str(r10403, "1", 10, MPFR_RNDN); mpfr_init(r10404); mpfr_init(r10405); mpfr_init(r10406); mpfr_init(r10407); mpfr_init(r10408); mpfr_init(r10409); } double f_im(double alpha, double beta) { mpfr_set_d(r10395, alpha, MPFR_RNDN); mpfr_set_d(r10396, beta, MPFR_RNDN); mpfr_add(r10397, r10395, r10396, MPFR_RNDN); mpfr_mul(r10398, r10396, r10395, MPFR_RNDN); mpfr_add(r10399, r10397, r10398, MPFR_RNDN); ; mpfr_add(r10401, r10399, r10400, MPFR_RNDN); ; ; mpfr_mul(r10404, r10402, r10403, MPFR_RNDN); mpfr_add(r10405, r10397, r10404, MPFR_RNDN); mpfr_div(r10406, r10401, r10405, MPFR_RNDN); mpfr_div(r10407, r10406, r10405, MPFR_RNDN); mpfr_add(r10408, r10405, r10400, MPFR_RNDN); mpfr_div(r10409, r10407, r10408, MPFR_RNDN); return mpfr_get_d(r10409, MPFR_RNDN); } static mpfr_t r10410, r10411, r10412, r10413, r10414, r10415, r10416, r10417, r10418, r10419, r10420, r10421, r10422, r10423, r10424, r10425, r10426, r10427, r10428, r10429, r10430, r10431, r10432, r10433, r10434, r10435, r10436; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r10410); mpfr_init_set_str(r10411, "3.2400543536458e+168", 10, MPFR_RNDN); mpfr_init(r10412); mpfr_init_set_str(r10413, "1.0", 10, MPFR_RNDN); mpfr_init(r10414); mpfr_init(r10415); mpfr_init(r10416); mpfr_init(r10417); mpfr_init(r10418); mpfr_init_set_str(r10419, "2", 10, MPFR_RNDN); mpfr_init(r10420); mpfr_init(r10421); mpfr_init(r10422); mpfr_init(r10423); mpfr_init(r10424); mpfr_init_set_str(r10425, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r10426, "0.25", 10, MPFR_RNDN); mpfr_init(r10427); mpfr_init(r10428); mpfr_init(r10429); mpfr_init(r10430); mpfr_init(r10431); mpfr_init(r10432); mpfr_init(r10433); mpfr_init(r10434); mpfr_init(r10435); mpfr_init(r10436); } double f_fm(double alpha, double beta) { mpfr_set_d(r10410, alpha, MPFR_RNDN); ; mpfr_set_si(r10412, mpfr_cmp(r10410, r10411) <= 0, MPFR_RNDN); ; mpfr_set_d(r10414, beta, MPFR_RNDN); mpfr_mul(r10415, r10414, r10410, MPFR_RNDN); mpfr_add(r10416, r10414, r10410, MPFR_RNDN); mpfr_add(r10417, r10415, r10416, MPFR_RNDN); mpfr_add(r10418, r10413, r10417, MPFR_RNDN); ; mpfr_add(r10420, r10416, r10419, MPFR_RNDN); mpfr_div(r10421, r10418, r10420, MPFR_RNDN); mpfr_div(r10422, r10421, r10420, MPFR_RNDN); mpfr_add(r10423, r10413, r10420, MPFR_RNDN); mpfr_div(r10424, r10422, r10423, MPFR_RNDN); ; ; mpfr_mul(r10427, r10416, r10426, MPFR_RNDN); mpfr_add(r10428, r10425, r10427, MPFR_RNDN); mpfr_add(r10429, r10410, r10419, MPFR_RNDN); mpfr_add(r10430, r10414, r10429, MPFR_RNDN); mpfr_add(r10431, r10414, r10419, MPFR_RNDN); mpfr_add(r10432, r10410, r10413, MPFR_RNDN); mpfr_add(r10433, r10431, r10432, MPFR_RNDN); mpfr_mul(r10434, r10430, r10433, MPFR_RNDN); mpfr_div(r10435, r10428, r10434, MPFR_RNDN); if (mpfr_get_si(r10412, MPFR_RNDN)) { mpfr_set(r10436, r10424, MPFR_RNDN); } else { mpfr_set(r10436, r10435, MPFR_RNDN); }; return mpfr_get_d(r10436, MPFR_RNDN); } static mpfr_t r10437, r10438, r10439, r10440, r10441, r10442, r10443, r10444, r10445, r10446, r10447, r10448, r10449, r10450, r10451, r10452, r10453, r10454, r10455, r10456, r10457, r10458, r10459, r10460, r10461, r10462, r10463; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r10437); mpfr_init_set_str(r10438, "3.2400543536458e+168", 10, MPFR_RNDN); mpfr_init(r10439); mpfr_init_set_str(r10440, "1.0", 10, MPFR_RNDN); mpfr_init(r10441); mpfr_init(r10442); mpfr_init(r10443); mpfr_init(r10444); mpfr_init(r10445); mpfr_init_set_str(r10446, "2", 10, MPFR_RNDN); mpfr_init(r10447); mpfr_init(r10448); mpfr_init(r10449); mpfr_init(r10450); mpfr_init(r10451); mpfr_init_set_str(r10452, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r10453, "0.25", 10, MPFR_RNDN); mpfr_init(r10454); mpfr_init(r10455); mpfr_init(r10456); mpfr_init(r10457); mpfr_init(r10458); mpfr_init(r10459); mpfr_init(r10460); mpfr_init(r10461); mpfr_init(r10462); mpfr_init(r10463); } double f_dm(double alpha, double beta) { mpfr_set_d(r10437, alpha, MPFR_RNDN); ; mpfr_set_si(r10439, mpfr_cmp(r10437, r10438) <= 0, MPFR_RNDN); ; mpfr_set_d(r10441, beta, MPFR_RNDN); mpfr_mul(r10442, r10441, r10437, MPFR_RNDN); mpfr_add(r10443, r10441, r10437, MPFR_RNDN); mpfr_add(r10444, r10442, r10443, MPFR_RNDN); mpfr_add(r10445, r10440, r10444, MPFR_RNDN); ; mpfr_add(r10447, r10443, r10446, MPFR_RNDN); mpfr_div(r10448, r10445, r10447, MPFR_RNDN); mpfr_div(r10449, r10448, r10447, MPFR_RNDN); mpfr_add(r10450, r10440, r10447, MPFR_RNDN); mpfr_div(r10451, r10449, r10450, MPFR_RNDN); ; ; mpfr_mul(r10454, r10443, r10453, MPFR_RNDN); mpfr_add(r10455, r10452, r10454, MPFR_RNDN); mpfr_add(r10456, r10437, r10446, MPFR_RNDN); mpfr_add(r10457, r10441, r10456, MPFR_RNDN); mpfr_add(r10458, r10441, r10446, MPFR_RNDN); mpfr_add(r10459, r10437, r10440, MPFR_RNDN); mpfr_add(r10460, r10458, r10459, MPFR_RNDN); mpfr_mul(r10461, r10457, r10460, MPFR_RNDN); mpfr_div(r10462, r10455, r10461, MPFR_RNDN); if (mpfr_get_si(r10439, MPFR_RNDN)) { mpfr_set(r10463, r10451, MPFR_RNDN); } else { mpfr_set(r10463, r10462, MPFR_RNDN); }; return mpfr_get_d(r10463, MPFR_RNDN); }