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