#include #include #include #include #include char *name = "Octave 3.8, jcobi/3"; double f_if(float alpha, float beta) { float r10271 = alpha; float r10272 = beta; float r10273 = r10271 + r10272; float r10274 = r10272 * r10271; float r10275 = r10273 + r10274; float r10276 = 1.0; float r10277 = r10275 + r10276; float r10278 = 2; float r10279 = 1; float r10280 = r10278 * r10279; float r10281 = r10273 + r10280; float r10282 = r10277 / r10281; float r10283 = r10282 / r10281; float r10284 = r10281 + r10276; float r10285 = r10283 / r10284; return r10285; } double f_id(double alpha, double beta) { double r10286 = alpha; double r10287 = beta; double r10288 = r10286 + r10287; double r10289 = r10287 * r10286; double r10290 = r10288 + r10289; double r10291 = 1.0; double r10292 = r10290 + r10291; double r10293 = 2; double r10294 = 1; double r10295 = r10293 * r10294; double r10296 = r10288 + r10295; double r10297 = r10292 / r10296; double r10298 = r10297 / r10296; double r10299 = r10296 + r10291; double r10300 = r10298 / r10299; return r10300; } double f_of(float alpha, float beta) { float r10301 = 0.5; float r10302 = 0.75; float r10303 = beta; float r10304 = alpha; float r10305 = r10303 + r10304; float r10306 = r10302 * r10305; float r10307 = r10301 + r10306; float r10308 = 2; float r10309 = r10305 + r10308; float r10310 = r10304 + r10308; float r10311 = 1.0; float r10312 = r10303 + r10311; float r10313 = r10310 + r10312; float r10314 = r10309 * r10313; float r10315 = r10307 / r10314; float r10316 = 6.9359525025468e-310; bool r10317 = r10315 <= r10316; float r10318 = r10304 + r10303; float r10319 = r10303 * r10304; float r10320 = r10319 + r10311; float r10321 = r10318 + r10320; float r10322 = r10318 + r10308; float r10323 = r10321 / r10322; float r10324 = r10308 + r10311; float r10325 = r10318 + r10324; float r10326 = r10325 * r10322; float r10327 = r10323 / r10326; float r10328 = r10317 ? r10315 : r10327; return r10328; } double f_od(double alpha, double beta) { double r10329 = 0.5; double r10330 = 0.75; double r10331 = beta; double r10332 = alpha; double r10333 = r10331 + r10332; double r10334 = r10330 * r10333; double r10335 = r10329 + r10334; double r10336 = 2; double r10337 = r10333 + r10336; double r10338 = r10332 + r10336; double r10339 = 1.0; double r10340 = r10331 + r10339; double r10341 = r10338 + r10340; double r10342 = r10337 * r10341; double r10343 = r10335 / r10342; double r10344 = 6.9359525025468e-310; bool r10345 = r10343 <= r10344; double r10346 = r10332 + r10331; double r10347 = r10331 * r10332; double r10348 = r10347 + r10339; double r10349 = r10346 + r10348; double r10350 = r10346 + r10336; double r10351 = r10349 / r10350; double r10352 = r10336 + r10339; double r10353 = r10346 + r10352; double r10354 = r10353 * r10350; double r10355 = r10351 / r10354; double r10356 = r10345 ? r10343 : r10355; return r10356; } 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 r10357, r10358, r10359, r10360, r10361, r10362, r10363, r10364, r10365, r10366, r10367, r10368, r10369, r10370, r10371; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r10357); mpfr_init(r10358); mpfr_init(r10359); mpfr_init(r10360); mpfr_init(r10361); mpfr_init_set_str(r10362, "1.0", 10, MPFR_RNDN); mpfr_init(r10363); mpfr_init_set_str(r10364, "2", 10, MPFR_RNDN); mpfr_init_set_str(r10365, "1", 10, MPFR_RNDN); mpfr_init(r10366); mpfr_init(r10367); mpfr_init(r10368); mpfr_init(r10369); mpfr_init(r10370); mpfr_init(r10371); } double f_im(double alpha, double beta) { mpfr_set_d(r10357, alpha, MPFR_RNDN); mpfr_set_d(r10358, beta, MPFR_RNDN); mpfr_add(r10359, r10357, r10358, MPFR_RNDN); mpfr_mul(r10360, r10358, r10357, MPFR_RNDN); mpfr_add(r10361, r10359, r10360, MPFR_RNDN); ; mpfr_add(r10363, r10361, r10362, MPFR_RNDN); ; ; mpfr_mul(r10366, r10364, r10365, MPFR_RNDN); mpfr_add(r10367, r10359, r10366, MPFR_RNDN); mpfr_div(r10368, r10363, r10367, MPFR_RNDN); mpfr_div(r10369, r10368, r10367, MPFR_RNDN); mpfr_add(r10370, r10367, r10362, MPFR_RNDN); mpfr_div(r10371, r10369, r10370, MPFR_RNDN); return mpfr_get_d(r10371, MPFR_RNDN); } static mpfr_t r10372, r10373, r10374, r10375, r10376, r10377, r10378, r10379, r10380, r10381, r10382, r10383, r10384, r10385, r10386, r10387, r10388, r10389, r10390, r10391, r10392, r10393, r10394, r10395, r10396, r10397, r10398, r10399; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r10372, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r10373, "0.75", 10, MPFR_RNDN); mpfr_init(r10374); mpfr_init(r10375); mpfr_init(r10376); mpfr_init(r10377); mpfr_init(r10378); mpfr_init_set_str(r10379, "2", 10, MPFR_RNDN); mpfr_init(r10380); mpfr_init(r10381); mpfr_init_set_str(r10382, "1.0", 10, MPFR_RNDN); mpfr_init(r10383); mpfr_init(r10384); mpfr_init(r10385); mpfr_init(r10386); mpfr_init_set_str(r10387, "6.9359525025468e-310", 10, MPFR_RNDN); mpfr_init(r10388); mpfr_init(r10389); mpfr_init(r10390); mpfr_init(r10391); mpfr_init(r10392); mpfr_init(r10393); mpfr_init(r10394); mpfr_init(r10395); mpfr_init(r10396); mpfr_init(r10397); mpfr_init(r10398); mpfr_init(r10399); } double f_fm(double alpha, double beta) { ; ; mpfr_set_d(r10374, beta, MPFR_RNDN); mpfr_set_d(r10375, alpha, MPFR_RNDN); mpfr_add(r10376, r10374, r10375, MPFR_RNDN); mpfr_mul(r10377, r10373, r10376, MPFR_RNDN); mpfr_add(r10378, r10372, r10377, MPFR_RNDN); ; mpfr_add(r10380, r10376, r10379, MPFR_RNDN); mpfr_add(r10381, r10375, r10379, MPFR_RNDN); ; mpfr_add(r10383, r10374, r10382, MPFR_RNDN); mpfr_add(r10384, r10381, r10383, MPFR_RNDN); mpfr_mul(r10385, r10380, r10384, MPFR_RNDN); mpfr_div(r10386, r10378, r10385, MPFR_RNDN); ; mpfr_set_si(r10388, mpfr_cmp(r10386, r10387) <= 0, MPFR_RNDN); mpfr_add(r10389, r10375, r10374, MPFR_RNDN); mpfr_mul(r10390, r10374, r10375, MPFR_RNDN); mpfr_add(r10391, r10390, r10382, MPFR_RNDN); mpfr_add(r10392, r10389, r10391, MPFR_RNDN); mpfr_add(r10393, r10389, r10379, MPFR_RNDN); mpfr_div(r10394, r10392, r10393, MPFR_RNDN); mpfr_add(r10395, r10379, r10382, MPFR_RNDN); mpfr_add(r10396, r10389, r10395, MPFR_RNDN); mpfr_mul(r10397, r10396, r10393, MPFR_RNDN); mpfr_div(r10398, r10394, r10397, MPFR_RNDN); if (mpfr_get_si(r10388, MPFR_RNDN)) { mpfr_set(r10399, r10386, MPFR_RNDN); } else { mpfr_set(r10399, r10398, MPFR_RNDN); }; return mpfr_get_d(r10399, MPFR_RNDN); } static mpfr_t r10400, r10401, r10402, r10403, r10404, r10405, r10406, r10407, r10408, r10409, r10410, r10411, r10412, r10413, r10414, r10415, r10416, r10417, r10418, r10419, r10420, r10421, r10422, r10423, r10424, r10425, r10426, r10427; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r10400, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r10401, "0.75", 10, MPFR_RNDN); mpfr_init(r10402); mpfr_init(r10403); mpfr_init(r10404); mpfr_init(r10405); mpfr_init(r10406); mpfr_init_set_str(r10407, "2", 10, MPFR_RNDN); mpfr_init(r10408); mpfr_init(r10409); mpfr_init_set_str(r10410, "1.0", 10, MPFR_RNDN); mpfr_init(r10411); mpfr_init(r10412); mpfr_init(r10413); mpfr_init(r10414); mpfr_init_set_str(r10415, "6.9359525025468e-310", 10, MPFR_RNDN); mpfr_init(r10416); mpfr_init(r10417); mpfr_init(r10418); mpfr_init(r10419); mpfr_init(r10420); mpfr_init(r10421); mpfr_init(r10422); mpfr_init(r10423); mpfr_init(r10424); mpfr_init(r10425); mpfr_init(r10426); mpfr_init(r10427); } double f_dm(double alpha, double beta) { ; ; mpfr_set_d(r10402, beta, MPFR_RNDN); mpfr_set_d(r10403, alpha, MPFR_RNDN); mpfr_add(r10404, r10402, r10403, MPFR_RNDN); mpfr_mul(r10405, r10401, r10404, MPFR_RNDN); mpfr_add(r10406, r10400, r10405, MPFR_RNDN); ; mpfr_add(r10408, r10404, r10407, MPFR_RNDN); mpfr_add(r10409, r10403, r10407, MPFR_RNDN); ; mpfr_add(r10411, r10402, r10410, MPFR_RNDN); mpfr_add(r10412, r10409, r10411, MPFR_RNDN); mpfr_mul(r10413, r10408, r10412, MPFR_RNDN); mpfr_div(r10414, r10406, r10413, MPFR_RNDN); ; mpfr_set_si(r10416, mpfr_cmp(r10414, r10415) <= 0, MPFR_RNDN); mpfr_add(r10417, r10403, r10402, MPFR_RNDN); mpfr_mul(r10418, r10402, r10403, MPFR_RNDN); mpfr_add(r10419, r10418, r10410, MPFR_RNDN); mpfr_add(r10420, r10417, r10419, MPFR_RNDN); mpfr_add(r10421, r10417, r10407, MPFR_RNDN); mpfr_div(r10422, r10420, r10421, MPFR_RNDN); mpfr_add(r10423, r10407, r10410, MPFR_RNDN); mpfr_add(r10424, r10417, r10423, MPFR_RNDN); mpfr_mul(r10425, r10424, r10421, MPFR_RNDN); mpfr_div(r10426, r10422, r10425, MPFR_RNDN); if (mpfr_get_si(r10416, MPFR_RNDN)) { mpfr_set(r10427, r10414, MPFR_RNDN); } else { mpfr_set(r10427, r10426, MPFR_RNDN); }; return mpfr_get_d(r10427, MPFR_RNDN); }