#include #include #include #include #include char *name = "Octave 3.8, jcobi/1"; double f_if(float alpha, float beta) { float r9445 = beta; float r9446 = alpha; float r9447 = r9445 - r9446; float r9448 = r9446 + r9445; float r9449 = 2.0f; float r9450 = r9448 + r9449; float r9451 = r9447 / r9450; float r9452 = 1.0f; float r9453 = r9451 + r9452; float r9454 = r9453 / r9449; return r9454; } double f_id(double alpha, double beta) { double r9455 = beta; double r9456 = alpha; double r9457 = r9455 - r9456; double r9458 = r9456 + r9455; double r9459 = 2.0; double r9460 = r9458 + r9459; double r9461 = r9457 / r9460; double r9462 = 1.0; double r9463 = r9461 + r9462; double r9464 = r9463 / r9459; return r9464; } double f_of(float alpha, float beta) { float r9465 = 1.0f; float r9466 = beta; float r9467 = alpha; float r9468 = r9466 - r9467; float r9469 = r9467 + r9466; float r9470 = 2.0f; float r9471 = r9469 + r9470; float r9472 = r9468 / r9471; float r9473 = r9465 + r9472; float r9474 = r9473 / r9470; float r9475 = 1.641367264033235e-13f; bool r9476 = r9474 <= r9475; float r9477 = r9466 / r9471; float r9478 = 1.0f; float r9479 = r9467 * r9467; float r9480 = r9478 / r9479; float r9481 = 4.0f; float r9482 = 8.0f; float r9483 = r9482 / r9467; float r9484 = r9481 - r9483; float r9485 = -r9470; float r9486 = r9485 / r9467; float r9487 = fma(r9480, r9484, r9486); float r9488 = r9477 - r9487; float r9489 = r9488 / r9470; float r9490 = 0.5000001989457653f; bool r9491 = r9474 <= r9490; float r9492 = 3.0f; float r9493 = pow(r9470, r9492); float r9494 = pow(r9469, r9492); float r9495 = r9493 + r9494; float r9496 = r9468 / r9495; float r9497 = r9470 * r9470; float r9498 = r9469 * r9470; float r9499 = r9497 - r9498; float r9500 = r9469 * r9469; float r9501 = r9499 + r9500; float r9502 = fma(r9496, r9501, r9465); float r9503 = r9502 / r9470; float r9504 = log(r9477); float r9505 = exp(r9504); float r9506 = r9467 / r9471; float r9507 = r9506 - r9465; float r9508 = r9505 - r9507; float r9509 = r9508 / r9470; float r9510 = r9491 ? r9503 : r9509; float r9511 = r9476 ? r9489 : r9510; return r9511; } double f_od(double alpha, double beta) { double r9512 = 1.0; double r9513 = beta; double r9514 = alpha; double r9515 = r9513 - r9514; double r9516 = r9514 + r9513; double r9517 = 2.0; double r9518 = r9516 + r9517; double r9519 = r9515 / r9518; double r9520 = r9512 + r9519; double r9521 = r9520 / r9517; double r9522 = 1.641367264033235e-13; bool r9523 = r9521 <= r9522; double r9524 = r9513 / r9518; double r9525 = 1.0; double r9526 = r9514 * r9514; double r9527 = r9525 / r9526; double r9528 = 4.0; double r9529 = 8.0; double r9530 = r9529 / r9514; double r9531 = r9528 - r9530; double r9532 = -r9517; double r9533 = r9532 / r9514; double r9534 = fma(r9527, r9531, r9533); double r9535 = r9524 - r9534; double r9536 = r9535 / r9517; double r9537 = 0.5000001989457653; bool r9538 = r9521 <= r9537; double r9539 = 3.0; double r9540 = pow(r9517, r9539); double r9541 = pow(r9516, r9539); double r9542 = r9540 + r9541; double r9543 = r9515 / r9542; double r9544 = r9517 * r9517; double r9545 = r9516 * r9517; double r9546 = r9544 - r9545; double r9547 = r9516 * r9516; double r9548 = r9546 + r9547; double r9549 = fma(r9543, r9548, r9512); double r9550 = r9549 / r9517; double r9551 = log(r9524); double r9552 = exp(r9551); double r9553 = r9514 / r9518; double r9554 = r9553 - r9512; double r9555 = r9552 - r9554; double r9556 = r9555 / r9517; double r9557 = r9538 ? r9550 : r9556; double r9558 = r9523 ? r9536 : r9557; return r9558; } 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 r9559, r9560, r9561, r9562, r9563, r9564, r9565, r9566, r9567, r9568; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r9559); mpfr_init(r9560); mpfr_init(r9561); mpfr_init(r9562); mpfr_init_set_str(r9563, "2.0", 10, MPFR_RNDN); mpfr_init(r9564); mpfr_init(r9565); mpfr_init_set_str(r9566, "1.0", 10, MPFR_RNDN); mpfr_init(r9567); mpfr_init(r9568); } double f_im(double alpha, double beta) { mpfr_set_d(r9559, beta, MPFR_RNDN); mpfr_set_d(r9560, alpha, MPFR_RNDN); mpfr_sub(r9561, r9559, r9560, MPFR_RNDN); mpfr_add(r9562, r9560, r9559, MPFR_RNDN); ; mpfr_add(r9564, r9562, r9563, MPFR_RNDN); mpfr_div(r9565, r9561, r9564, MPFR_RNDN); ; mpfr_add(r9567, r9565, r9566, MPFR_RNDN); mpfr_div(r9568, r9567, r9563, MPFR_RNDN); return mpfr_get_d(r9568, MPFR_RNDN); } static mpfr_t r9569, r9570, r9571, r9572, r9573, r9574, r9575, r9576, r9577, r9578, r9579, r9580, r9581, r9582, r9583, r9584, r9585, r9586, r9587, r9588, r9589, r9590, r9591, r9592, r9593, r9594, r9595, r9596, r9597, r9598, r9599, r9600, r9601, r9602, r9603, r9604, r9605, r9606, r9607, r9608, r9609, r9610, r9611, r9612, r9613, r9614, r9615; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r9569, "1.0", 10, MPFR_RNDN); mpfr_init(r9570); mpfr_init(r9571); mpfr_init(r9572); mpfr_init(r9573); mpfr_init_set_str(r9574, "2.0", 10, MPFR_RNDN); mpfr_init(r9575); mpfr_init(r9576); mpfr_init(r9577); mpfr_init(r9578); mpfr_init_set_str(r9579, "1.641367264033235e-13", 10, MPFR_RNDN); mpfr_init(r9580); mpfr_init(r9581); mpfr_init_set_str(r9582, "1", 10, MPFR_RNDN); mpfr_init(r9583); mpfr_init(r9584); mpfr_init_set_str(r9585, "4.0", 10, MPFR_RNDN); mpfr_init_set_str(r9586, "8.0", 10, MPFR_RNDN); mpfr_init(r9587); mpfr_init(r9588); mpfr_init(r9589); mpfr_init(r9590); mpfr_init(r9591); mpfr_init(r9592); mpfr_init(r9593); mpfr_init_set_str(r9594, "0.5000001989457653", 10, MPFR_RNDN); mpfr_init(r9595); mpfr_init_set_str(r9596, "3", 10, MPFR_RNDN); mpfr_init(r9597); mpfr_init(r9598); mpfr_init(r9599); mpfr_init(r9600); mpfr_init(r9601); mpfr_init(r9602); mpfr_init(r9603); mpfr_init(r9604); mpfr_init(r9605); mpfr_init(r9606); mpfr_init(r9607); mpfr_init(r9608); mpfr_init(r9609); mpfr_init(r9610); mpfr_init(r9611); mpfr_init(r9612); mpfr_init(r9613); mpfr_init(r9614); mpfr_init(r9615); } double f_fm(double alpha, double beta) { ; mpfr_set_d(r9570, beta, MPFR_RNDN); mpfr_set_d(r9571, alpha, MPFR_RNDN); mpfr_sub(r9572, r9570, r9571, MPFR_RNDN); mpfr_add(r9573, r9571, r9570, MPFR_RNDN); ; mpfr_add(r9575, r9573, r9574, MPFR_RNDN); mpfr_div(r9576, r9572, r9575, MPFR_RNDN); mpfr_add(r9577, r9569, r9576, MPFR_RNDN); mpfr_div(r9578, r9577, r9574, MPFR_RNDN); ; mpfr_set_si(r9580, mpfr_cmp(r9578, r9579) <= 0, MPFR_RNDN); mpfr_div(r9581, r9570, r9575, MPFR_RNDN); ; mpfr_mul(r9583, r9571, r9571, MPFR_RNDN); mpfr_div(r9584, r9582, r9583, MPFR_RNDN); ; ; mpfr_div(r9587, r9586, r9571, MPFR_RNDN); mpfr_sub(r9588, r9585, r9587, MPFR_RNDN); mpfr_neg(r9589, r9574, MPFR_RNDN); mpfr_div(r9590, r9589, r9571, MPFR_RNDN); mpfr_fma(r9591, r9584, r9588, r9590, MPFR_RNDN); mpfr_sub(r9592, r9581, r9591, MPFR_RNDN); mpfr_div(r9593, r9592, r9574, MPFR_RNDN); ; mpfr_set_si(r9595, mpfr_cmp(r9578, r9594) <= 0, MPFR_RNDN); ; mpfr_pow(r9597, r9574, r9596, MPFR_RNDN); mpfr_pow(r9598, r9573, r9596, MPFR_RNDN); mpfr_add(r9599, r9597, r9598, MPFR_RNDN); mpfr_div(r9600, r9572, r9599, MPFR_RNDN); mpfr_mul(r9601, r9574, r9574, MPFR_RNDN); mpfr_mul(r9602, r9573, r9574, MPFR_RNDN); mpfr_sub(r9603, r9601, r9602, MPFR_RNDN); mpfr_mul(r9604, r9573, r9573, MPFR_RNDN); mpfr_add(r9605, r9603, r9604, MPFR_RNDN); mpfr_fma(r9606, r9600, r9605, r9569, MPFR_RNDN); mpfr_div(r9607, r9606, r9574, MPFR_RNDN); mpfr_log(r9608, r9581, MPFR_RNDN); mpfr_exp(r9609, r9608, MPFR_RNDN); mpfr_div(r9610, r9571, r9575, MPFR_RNDN); mpfr_sub(r9611, r9610, r9569, MPFR_RNDN); mpfr_sub(r9612, r9609, r9611, MPFR_RNDN); mpfr_div(r9613, r9612, r9574, MPFR_RNDN); if (mpfr_get_si(r9595, MPFR_RNDN)) { mpfr_set(r9614, r9607, MPFR_RNDN); } else { mpfr_set(r9614, r9613, MPFR_RNDN); }; if (mpfr_get_si(r9580, MPFR_RNDN)) { mpfr_set(r9615, r9593, MPFR_RNDN); } else { mpfr_set(r9615, r9614, MPFR_RNDN); }; return mpfr_get_d(r9615, MPFR_RNDN); } static mpfr_t r9616, r9617, r9618, r9619, r9620, r9621, r9622, r9623, r9624, r9625, r9626, r9627, r9628, r9629, r9630, r9631, r9632, r9633, r9634, r9635, r9636, r9637, r9638, r9639, r9640, r9641, r9642, r9643, r9644, r9645, r9646, r9647, r9648, r9649, r9650, r9651, r9652, r9653, r9654, r9655, r9656, r9657, r9658, r9659, r9660, r9661, r9662; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r9616, "1.0", 10, MPFR_RNDN); mpfr_init(r9617); mpfr_init(r9618); mpfr_init(r9619); mpfr_init(r9620); mpfr_init_set_str(r9621, "2.0", 10, MPFR_RNDN); mpfr_init(r9622); mpfr_init(r9623); mpfr_init(r9624); mpfr_init(r9625); mpfr_init_set_str(r9626, "1.641367264033235e-13", 10, MPFR_RNDN); mpfr_init(r9627); mpfr_init(r9628); mpfr_init_set_str(r9629, "1", 10, MPFR_RNDN); mpfr_init(r9630); mpfr_init(r9631); mpfr_init_set_str(r9632, "4.0", 10, MPFR_RNDN); mpfr_init_set_str(r9633, "8.0", 10, MPFR_RNDN); mpfr_init(r9634); mpfr_init(r9635); mpfr_init(r9636); mpfr_init(r9637); mpfr_init(r9638); mpfr_init(r9639); mpfr_init(r9640); mpfr_init_set_str(r9641, "0.5000001989457653", 10, MPFR_RNDN); mpfr_init(r9642); mpfr_init_set_str(r9643, "3", 10, MPFR_RNDN); mpfr_init(r9644); mpfr_init(r9645); mpfr_init(r9646); mpfr_init(r9647); mpfr_init(r9648); mpfr_init(r9649); mpfr_init(r9650); mpfr_init(r9651); mpfr_init(r9652); mpfr_init(r9653); mpfr_init(r9654); mpfr_init(r9655); mpfr_init(r9656); mpfr_init(r9657); mpfr_init(r9658); mpfr_init(r9659); mpfr_init(r9660); mpfr_init(r9661); mpfr_init(r9662); } double f_dm(double alpha, double beta) { ; mpfr_set_d(r9617, beta, MPFR_RNDN); mpfr_set_d(r9618, alpha, MPFR_RNDN); mpfr_sub(r9619, r9617, r9618, MPFR_RNDN); mpfr_add(r9620, r9618, r9617, MPFR_RNDN); ; mpfr_add(r9622, r9620, r9621, MPFR_RNDN); mpfr_div(r9623, r9619, r9622, MPFR_RNDN); mpfr_add(r9624, r9616, r9623, MPFR_RNDN); mpfr_div(r9625, r9624, r9621, MPFR_RNDN); ; mpfr_set_si(r9627, mpfr_cmp(r9625, r9626) <= 0, MPFR_RNDN); mpfr_div(r9628, r9617, r9622, MPFR_RNDN); ; mpfr_mul(r9630, r9618, r9618, MPFR_RNDN); mpfr_div(r9631, r9629, r9630, MPFR_RNDN); ; ; mpfr_div(r9634, r9633, r9618, MPFR_RNDN); mpfr_sub(r9635, r9632, r9634, MPFR_RNDN); mpfr_neg(r9636, r9621, MPFR_RNDN); mpfr_div(r9637, r9636, r9618, MPFR_RNDN); mpfr_fma(r9638, r9631, r9635, r9637, MPFR_RNDN); mpfr_sub(r9639, r9628, r9638, MPFR_RNDN); mpfr_div(r9640, r9639, r9621, MPFR_RNDN); ; mpfr_set_si(r9642, mpfr_cmp(r9625, r9641) <= 0, MPFR_RNDN); ; mpfr_pow(r9644, r9621, r9643, MPFR_RNDN); mpfr_pow(r9645, r9620, r9643, MPFR_RNDN); mpfr_add(r9646, r9644, r9645, MPFR_RNDN); mpfr_div(r9647, r9619, r9646, MPFR_RNDN); mpfr_mul(r9648, r9621, r9621, MPFR_RNDN); mpfr_mul(r9649, r9620, r9621, MPFR_RNDN); mpfr_sub(r9650, r9648, r9649, MPFR_RNDN); mpfr_mul(r9651, r9620, r9620, MPFR_RNDN); mpfr_add(r9652, r9650, r9651, MPFR_RNDN); mpfr_fma(r9653, r9647, r9652, r9616, MPFR_RNDN); mpfr_div(r9654, r9653, r9621, MPFR_RNDN); mpfr_log(r9655, r9628, MPFR_RNDN); mpfr_exp(r9656, r9655, MPFR_RNDN); mpfr_div(r9657, r9618, r9622, MPFR_RNDN); mpfr_sub(r9658, r9657, r9616, MPFR_RNDN); mpfr_sub(r9659, r9656, r9658, MPFR_RNDN); mpfr_div(r9660, r9659, r9621, MPFR_RNDN); if (mpfr_get_si(r9642, MPFR_RNDN)) { mpfr_set(r9661, r9654, MPFR_RNDN); } else { mpfr_set(r9661, r9660, MPFR_RNDN); }; if (mpfr_get_si(r9627, MPFR_RNDN)) { mpfr_set(r9662, r9640, MPFR_RNDN); } else { mpfr_set(r9662, r9661, MPFR_RNDN); }; return mpfr_get_d(r9662, MPFR_RNDN); }