#include #include #include #include #include char *name = "Octave 3.8, jcobi/2"; double f_if(float alpha, float beta, float i) { float r9305 = alpha; float r9306 = beta; float r9307 = r9305 + r9306; float r9308 = r9306 - r9305; float r9309 = r9307 * r9308; float r9310 = 2.0f; float r9311 = i; float r9312 = r9310 * r9311; float r9313 = r9307 + r9312; float r9314 = r9309 / r9313; float r9315 = 2.0f; float r9316 = r9313 + r9315; float r9317 = r9314 / r9316; float r9318 = 1.0f; float r9319 = r9317 + r9318; float r9320 = r9319 / r9315; return r9320; } double f_id(double alpha, double beta, double i) { double r9321 = alpha; double r9322 = beta; double r9323 = r9321 + r9322; double r9324 = r9322 - r9321; double r9325 = r9323 * r9324; double r9326 = 2.0; double r9327 = i; double r9328 = r9326 * r9327; double r9329 = r9323 + r9328; double r9330 = r9325 / r9329; double r9331 = 2.0; double r9332 = r9329 + r9331; double r9333 = r9330 / r9332; double r9334 = 1.0; double r9335 = r9333 + r9334; double r9336 = r9335 / r9331; return r9336; } double f_of(float alpha, float beta, float i) { float r9337 = alpha; float r9338 = beta; float r9339 = r9337 + r9338; float r9340 = 1.0f; float r9341 = 2.0f; float r9342 = r9338 + r9341; float r9343 = 2.0f; float r9344 = i; float r9345 = fma(r9343, r9344, r9337); float r9346 = r9342 + r9345; float r9347 = r9340 / r9346; float r9348 = r9339 * r9347; float r9349 = r9338 - r9337; float r9350 = r9338 + r9345; float r9351 = r9349 / r9350; float r9352 = 1.0f; float r9353 = fma(r9348, r9351, r9352); float r9354 = r9351 * r9348; float r9355 = r9354 + r9352; float r9356 = r9353 * r9355; float r9357 = r9353 * r9356; float r9358 = cbrt(r9357); float r9359 = r9358 / r9341; return r9359; } double f_od(double alpha, double beta, double i) { double r9360 = alpha; double r9361 = beta; double r9362 = r9360 + r9361; double r9363 = 1.0; double r9364 = 2.0; double r9365 = r9361 + r9364; double r9366 = 2.0; double r9367 = i; double r9368 = fma(r9366, r9367, r9360); double r9369 = r9365 + r9368; double r9370 = r9363 / r9369; double r9371 = r9362 * r9370; double r9372 = r9361 - r9360; double r9373 = r9361 + r9368; double r9374 = r9372 / r9373; double r9375 = 1.0; double r9376 = fma(r9371, r9374, r9375); double r9377 = r9374 * r9371; double r9378 = r9377 + r9375; double r9379 = r9376 * r9378; double r9380 = r9376 * r9379; double r9381 = cbrt(r9380); double r9382 = r9381 / r9364; return r9382; } 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 r9383, r9384, r9385, r9386, r9387, r9388, r9389, r9390, r9391, r9392, r9393, r9394, r9395, r9396, r9397, r9398; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r9383); mpfr_init(r9384); mpfr_init(r9385); mpfr_init(r9386); mpfr_init(r9387); mpfr_init_set_str(r9388, "2", 10, MPFR_RNDN); mpfr_init(r9389); mpfr_init(r9390); mpfr_init(r9391); mpfr_init(r9392); mpfr_init_set_str(r9393, "2.0", 10, MPFR_RNDN); mpfr_init(r9394); mpfr_init(r9395); mpfr_init_set_str(r9396, "1.0", 10, MPFR_RNDN); mpfr_init(r9397); mpfr_init(r9398); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r9383, alpha, MPFR_RNDN); mpfr_set_d(r9384, beta, MPFR_RNDN); mpfr_add(r9385, r9383, r9384, MPFR_RNDN); mpfr_sub(r9386, r9384, r9383, MPFR_RNDN); mpfr_mul(r9387, r9385, r9386, MPFR_RNDN); ; mpfr_set_d(r9389, i, MPFR_RNDN); mpfr_mul(r9390, r9388, r9389, MPFR_RNDN); mpfr_add(r9391, r9385, r9390, MPFR_RNDN); mpfr_div(r9392, r9387, r9391, MPFR_RNDN); ; mpfr_add(r9394, r9391, r9393, MPFR_RNDN); mpfr_div(r9395, r9392, r9394, MPFR_RNDN); ; mpfr_add(r9397, r9395, r9396, MPFR_RNDN); mpfr_div(r9398, r9397, r9393, MPFR_RNDN); return mpfr_get_d(r9398, MPFR_RNDN); } static mpfr_t r9399, r9400, r9401, r9402, r9403, r9404, r9405, r9406, r9407, r9408, r9409, r9410, r9411, r9412, r9413, r9414, r9415, r9416, r9417, r9418, r9419, r9420, r9421; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r9399); mpfr_init(r9400); mpfr_init(r9401); mpfr_init_set_str(r9402, "1", 10, MPFR_RNDN); mpfr_init_set_str(r9403, "2.0", 10, MPFR_RNDN); mpfr_init(r9404); mpfr_init_set_str(r9405, "2", 10, MPFR_RNDN); mpfr_init(r9406); mpfr_init(r9407); mpfr_init(r9408); mpfr_init(r9409); mpfr_init(r9410); mpfr_init(r9411); mpfr_init(r9412); mpfr_init(r9413); mpfr_init_set_str(r9414, "1.0", 10, MPFR_RNDN); mpfr_init(r9415); mpfr_init(r9416); mpfr_init(r9417); mpfr_init(r9418); mpfr_init(r9419); mpfr_init(r9420); mpfr_init(r9421); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r9399, alpha, MPFR_RNDN); mpfr_set_d(r9400, beta, MPFR_RNDN); mpfr_add(r9401, r9399, r9400, MPFR_RNDN); ; ; mpfr_add(r9404, r9400, r9403, MPFR_RNDN); ; mpfr_set_d(r9406, i, MPFR_RNDN); mpfr_fma(r9407, r9405, r9406, r9399, MPFR_RNDN); mpfr_add(r9408, r9404, r9407, MPFR_RNDN); mpfr_div(r9409, r9402, r9408, MPFR_RNDN); mpfr_mul(r9410, r9401, r9409, MPFR_RNDN); mpfr_sub(r9411, r9400, r9399, MPFR_RNDN); mpfr_add(r9412, r9400, r9407, MPFR_RNDN); mpfr_div(r9413, r9411, r9412, MPFR_RNDN); ; mpfr_fma(r9415, r9410, r9413, r9414, MPFR_RNDN); mpfr_mul(r9416, r9413, r9410, MPFR_RNDN); mpfr_add(r9417, r9416, r9414, MPFR_RNDN); mpfr_mul(r9418, r9415, r9417, MPFR_RNDN); mpfr_mul(r9419, r9415, r9418, MPFR_RNDN); mpfr_cbrt(r9420, r9419, MPFR_RNDN); mpfr_div(r9421, r9420, r9403, MPFR_RNDN); return mpfr_get_d(r9421, MPFR_RNDN); } static mpfr_t r9422, r9423, r9424, r9425, r9426, r9427, r9428, r9429, r9430, r9431, r9432, r9433, r9434, r9435, r9436, r9437, r9438, r9439, r9440, r9441, r9442, r9443, r9444; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r9422); mpfr_init(r9423); mpfr_init(r9424); mpfr_init_set_str(r9425, "1", 10, MPFR_RNDN); mpfr_init_set_str(r9426, "2.0", 10, MPFR_RNDN); mpfr_init(r9427); mpfr_init_set_str(r9428, "2", 10, MPFR_RNDN); mpfr_init(r9429); mpfr_init(r9430); mpfr_init(r9431); mpfr_init(r9432); mpfr_init(r9433); mpfr_init(r9434); mpfr_init(r9435); mpfr_init(r9436); mpfr_init_set_str(r9437, "1.0", 10, MPFR_RNDN); mpfr_init(r9438); mpfr_init(r9439); mpfr_init(r9440); mpfr_init(r9441); mpfr_init(r9442); mpfr_init(r9443); mpfr_init(r9444); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r9422, alpha, MPFR_RNDN); mpfr_set_d(r9423, beta, MPFR_RNDN); mpfr_add(r9424, r9422, r9423, MPFR_RNDN); ; ; mpfr_add(r9427, r9423, r9426, MPFR_RNDN); ; mpfr_set_d(r9429, i, MPFR_RNDN); mpfr_fma(r9430, r9428, r9429, r9422, MPFR_RNDN); mpfr_add(r9431, r9427, r9430, MPFR_RNDN); mpfr_div(r9432, r9425, r9431, MPFR_RNDN); mpfr_mul(r9433, r9424, r9432, MPFR_RNDN); mpfr_sub(r9434, r9423, r9422, MPFR_RNDN); mpfr_add(r9435, r9423, r9430, MPFR_RNDN); mpfr_div(r9436, r9434, r9435, MPFR_RNDN); ; mpfr_fma(r9438, r9433, r9436, r9437, MPFR_RNDN); mpfr_mul(r9439, r9436, r9433, MPFR_RNDN); mpfr_add(r9440, r9439, r9437, MPFR_RNDN); mpfr_mul(r9441, r9438, r9440, MPFR_RNDN); mpfr_mul(r9442, r9438, r9441, MPFR_RNDN); mpfr_cbrt(r9443, r9442, MPFR_RNDN); mpfr_div(r9444, r9443, r9426, MPFR_RNDN); return mpfr_get_d(r9444, MPFR_RNDN); }