#include #include #include #include #include char *name = "Henrywood and Agarwal, Equation (9a)"; double f_if(float w0, float M, float D, float h, float l, float d) { float r9589 = w0; float r9590 = 1.0f; float r9591 = M; float r9592 = D; float r9593 = r9591 * r9592; float r9594 = 2.0f; float r9595 = d; float r9596 = r9594 * r9595; float r9597 = r9593 / r9596; float r9598 = pow(r9597, r9594); float r9599 = h; float r9600 = l; float r9601 = r9599 / r9600; float r9602 = r9598 * r9601; float r9603 = r9590 - r9602; float r9604 = sqrt(r9603); float r9605 = r9589 * r9604; return r9605; } double f_id(double w0, double M, double D, double h, double l, double d) { double r9606 = w0; double r9607 = 1.0; double r9608 = M; double r9609 = D; double r9610 = r9608 * r9609; double r9611 = 2.0; double r9612 = d; double r9613 = r9611 * r9612; double r9614 = r9610 / r9613; double r9615 = pow(r9614, r9611); double r9616 = h; double r9617 = l; double r9618 = r9616 / r9617; double r9619 = r9615 * r9618; double r9620 = r9607 - r9619; double r9621 = sqrt(r9620); double r9622 = r9606 * r9621; return r9622; } double f_of(float w0, float M, float D, float h, float l, float d) { float r9623 = h; float r9624 = l; float r9625 = r9623 / r9624; float r9626 = -1.775245375421158e+308f; bool r9627 = r9625 <= r9626; float r9628 = -7.24625914193812e-213f; bool r9629 = r9625 <= r9628; float r9630 = !r9629; bool r9631 = r9627 || r9630; float r9632 = w0; float r9633 = -r9623; float r9634 = r9633 / r9624; float r9635 = M; float r9636 = D; float r9637 = r9635 * r9636; float r9638 = d; float r9639 = 2.0f; float r9640 = r9638 * r9639; float r9641 = r9637 / r9640; float r9642 = cbrt(r9641); float r9643 = r9642 * r9642; float r9644 = r9642 * r9643; float r9645 = r9644 * r9641; float r9646 = 1.0f; float r9647 = fma(r9634, r9645, r9646); float r9648 = sqrt(r9647); float r9649 = r9648 * r9632; float r9650 = r9631 ? r9632 : r9649; return r9650; } double f_od(double w0, double M, double D, double h, double l, double d) { double r9651 = h; double r9652 = l; double r9653 = r9651 / r9652; double r9654 = -1.775245375421158e+308; bool r9655 = r9653 <= r9654; double r9656 = -7.24625914193812e-213; bool r9657 = r9653 <= r9656; double r9658 = !r9657; bool r9659 = r9655 || r9658; double r9660 = w0; double r9661 = -r9651; double r9662 = r9661 / r9652; double r9663 = M; double r9664 = D; double r9665 = r9663 * r9664; double r9666 = d; double r9667 = 2.0; double r9668 = r9666 * r9667; double r9669 = r9665 / r9668; double r9670 = cbrt(r9669); double r9671 = r9670 * r9670; double r9672 = r9670 * r9671; double r9673 = r9672 * r9669; double r9674 = 1.0; double r9675 = fma(r9662, r9673, r9674); double r9676 = sqrt(r9675); double r9677 = r9676 * r9660; double r9678 = r9659 ? r9660 : r9677; return r9678; } 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 r9679, r9680, r9681, r9682, r9683, r9684, r9685, r9686, r9687, r9688, r9689, r9690, r9691, r9692, r9693, r9694, r9695; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r9679); mpfr_init_set_str(r9680, "1", 10, MPFR_RNDN); mpfr_init(r9681); mpfr_init(r9682); mpfr_init(r9683); mpfr_init_set_str(r9684, "2", 10, MPFR_RNDN); mpfr_init(r9685); mpfr_init(r9686); mpfr_init(r9687); mpfr_init(r9688); mpfr_init(r9689); mpfr_init(r9690); mpfr_init(r9691); mpfr_init(r9692); mpfr_init(r9693); mpfr_init(r9694); mpfr_init(r9695); } double f_im(double w0, double M, double D, double h, double l, double d) { mpfr_set_d(r9679, w0, MPFR_RNDN); ; mpfr_set_d(r9681, M, MPFR_RNDN); mpfr_set_d(r9682, D, MPFR_RNDN); mpfr_mul(r9683, r9681, r9682, MPFR_RNDN); ; mpfr_set_d(r9685, d, MPFR_RNDN); mpfr_mul(r9686, r9684, r9685, MPFR_RNDN); mpfr_div(r9687, r9683, r9686, MPFR_RNDN); mpfr_pow(r9688, r9687, r9684, MPFR_RNDN); mpfr_set_d(r9689, h, MPFR_RNDN); mpfr_set_d(r9690, l, MPFR_RNDN); mpfr_div(r9691, r9689, r9690, MPFR_RNDN); mpfr_mul(r9692, r9688, r9691, MPFR_RNDN); mpfr_sub(r9693, r9680, r9692, MPFR_RNDN); mpfr_sqrt(r9694, r9693, MPFR_RNDN); mpfr_mul(r9695, r9679, r9694, MPFR_RNDN); return mpfr_get_d(r9695, MPFR_RNDN); } static mpfr_t r9696, r9697, r9698, r9699, r9700, r9701, r9702, r9703, r9704, r9705, r9706, r9707, r9708, r9709, r9710, r9711, r9712, r9713, r9714, r9715, r9716, r9717, r9718, r9719, r9720, r9721, r9722, r9723; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r9696); mpfr_init(r9697); mpfr_init(r9698); mpfr_init_set_str(r9699, "-1.775245375421158e+308", 10, MPFR_RNDN); mpfr_init(r9700); mpfr_init_set_str(r9701, "-7.24625914193812e-213", 10, MPFR_RNDN); mpfr_init(r9702); mpfr_init(r9703); mpfr_init(r9704); mpfr_init(r9705); mpfr_init(r9706); mpfr_init(r9707); mpfr_init(r9708); mpfr_init(r9709); mpfr_init(r9710); mpfr_init(r9711); mpfr_init_set_str(r9712, "2", 10, MPFR_RNDN); mpfr_init(r9713); mpfr_init(r9714); mpfr_init(r9715); mpfr_init(r9716); mpfr_init(r9717); mpfr_init(r9718); mpfr_init_set_str(r9719, "1", 10, MPFR_RNDN); mpfr_init(r9720); mpfr_init(r9721); mpfr_init(r9722); mpfr_init(r9723); } double f_fm(double w0, double M, double D, double h, double l, double d) { mpfr_set_d(r9696, h, MPFR_RNDN); mpfr_set_d(r9697, l, MPFR_RNDN); mpfr_div(r9698, r9696, r9697, MPFR_RNDN); ; mpfr_set_si(r9700, mpfr_cmp(r9698, r9699) <= 0, MPFR_RNDN); ; mpfr_set_si(r9702, mpfr_cmp(r9698, r9701) <= 0, MPFR_RNDN); mpfr_set_si(r9703, !mpfr_get_si(r9702, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r9704, mpfr_get_si(r9700, MPFR_RNDN) || mpfr_get_si(r9703, MPFR_RNDN), MPFR_RNDN); mpfr_set_d(r9705, w0, MPFR_RNDN); mpfr_neg(r9706, r9696, MPFR_RNDN); mpfr_div(r9707, r9706, r9697, MPFR_RNDN); mpfr_set_d(r9708, M, MPFR_RNDN); mpfr_set_d(r9709, D, MPFR_RNDN); mpfr_mul(r9710, r9708, r9709, MPFR_RNDN); mpfr_set_d(r9711, d, MPFR_RNDN); ; mpfr_mul(r9713, r9711, r9712, MPFR_RNDN); mpfr_div(r9714, r9710, r9713, MPFR_RNDN); mpfr_cbrt(r9715, r9714, MPFR_RNDN); mpfr_mul(r9716, r9715, r9715, MPFR_RNDN); mpfr_mul(r9717, r9715, r9716, MPFR_RNDN); mpfr_mul(r9718, r9717, r9714, MPFR_RNDN); ; mpfr_fma(r9720, r9707, r9718, r9719, MPFR_RNDN); mpfr_sqrt(r9721, r9720, MPFR_RNDN); mpfr_mul(r9722, r9721, r9705, MPFR_RNDN); if (mpfr_get_si(r9704, MPFR_RNDN)) { mpfr_set(r9723, r9705, MPFR_RNDN); } else { mpfr_set(r9723, r9722, MPFR_RNDN); }; return mpfr_get_d(r9723, MPFR_RNDN); } static mpfr_t r9724, r9725, r9726, r9727, r9728, r9729, r9730, r9731, r9732, r9733, r9734, r9735, r9736, r9737, r9738, r9739, r9740, r9741, r9742, r9743, r9744, r9745, r9746, r9747, r9748, r9749, r9750, r9751; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r9724); mpfr_init(r9725); mpfr_init(r9726); mpfr_init_set_str(r9727, "-1.775245375421158e+308", 10, MPFR_RNDN); mpfr_init(r9728); mpfr_init_set_str(r9729, "-7.24625914193812e-213", 10, MPFR_RNDN); mpfr_init(r9730); mpfr_init(r9731); mpfr_init(r9732); mpfr_init(r9733); mpfr_init(r9734); mpfr_init(r9735); mpfr_init(r9736); mpfr_init(r9737); mpfr_init(r9738); mpfr_init(r9739); mpfr_init_set_str(r9740, "2", 10, MPFR_RNDN); mpfr_init(r9741); mpfr_init(r9742); mpfr_init(r9743); mpfr_init(r9744); mpfr_init(r9745); mpfr_init(r9746); mpfr_init_set_str(r9747, "1", 10, MPFR_RNDN); mpfr_init(r9748); mpfr_init(r9749); mpfr_init(r9750); mpfr_init(r9751); } double f_dm(double w0, double M, double D, double h, double l, double d) { mpfr_set_d(r9724, h, MPFR_RNDN); mpfr_set_d(r9725, l, MPFR_RNDN); mpfr_div(r9726, r9724, r9725, MPFR_RNDN); ; mpfr_set_si(r9728, mpfr_cmp(r9726, r9727) <= 0, MPFR_RNDN); ; mpfr_set_si(r9730, mpfr_cmp(r9726, r9729) <= 0, MPFR_RNDN); mpfr_set_si(r9731, !mpfr_get_si(r9730, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r9732, mpfr_get_si(r9728, MPFR_RNDN) || mpfr_get_si(r9731, MPFR_RNDN), MPFR_RNDN); mpfr_set_d(r9733, w0, MPFR_RNDN); mpfr_neg(r9734, r9724, MPFR_RNDN); mpfr_div(r9735, r9734, r9725, MPFR_RNDN); mpfr_set_d(r9736, M, MPFR_RNDN); mpfr_set_d(r9737, D, MPFR_RNDN); mpfr_mul(r9738, r9736, r9737, MPFR_RNDN); mpfr_set_d(r9739, d, MPFR_RNDN); ; mpfr_mul(r9741, r9739, r9740, MPFR_RNDN); mpfr_div(r9742, r9738, r9741, MPFR_RNDN); mpfr_cbrt(r9743, r9742, MPFR_RNDN); mpfr_mul(r9744, r9743, r9743, MPFR_RNDN); mpfr_mul(r9745, r9743, r9744, MPFR_RNDN); mpfr_mul(r9746, r9745, r9742, MPFR_RNDN); ; mpfr_fma(r9748, r9735, r9746, r9747, MPFR_RNDN); mpfr_sqrt(r9749, r9748, MPFR_RNDN); mpfr_mul(r9750, r9749, r9733, MPFR_RNDN); if (mpfr_get_si(r9732, MPFR_RNDN)) { mpfr_set(r9751, r9733, MPFR_RNDN); } else { mpfr_set(r9751, r9750, MPFR_RNDN); }; return mpfr_get_d(r9751, MPFR_RNDN); }