#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "Jmat.Real.lambertw, newton loop step";

double f_if(float wj, float x) {
        float r27671 = wj;
        float r27672 = exp(r27671);
        float r27673 = r27671 * r27672;
        float r27674 = x;
        float r27675 = r27673 - r27674;
        float r27676 = r27672 + r27673;
        float r27677 = r27675 / r27676;
        float r27678 = r27671 - r27677;
        return r27678;
}

double f_id(double wj, double x) {
        double r27679 = wj;
        double r27680 = exp(r27679);
        double r27681 = r27679 * r27680;
        double r27682 = x;
        double r27683 = r27681 - r27682;
        double r27684 = r27680 + r27681;
        double r27685 = r27683 / r27684;
        double r27686 = r27679 - r27685;
        return r27686;
}


double f_of(float wj, float x) {
        float r27687 = 1;
        float r27688 = wj;
        float r27689 = r27687 / r27688;
        float r27690 = r27689 + r27688;
        float r27691 = fma(r27689, r27689, r27687);
        float r27692 = r27690 - r27691;
        float r27693 = -125981247538304.98;
        bool r27694 = r27692 <= r27693;
        float r27695 = x;
        float r27696 = 2;
        float r27697 = r27695 * r27696;
        float r27698 = -r27688;
        float r27699 = fma(r27688, r27688, r27695);
        float r27700 = fma(r27697, r27698, r27699);
        float r27701 = exp(r27688);
        float r27702 = -r27695;
        float r27703 = fma(r27701, r27688, r27702);
        float r27704 = fma(r27688, r27701, r27701);
        float r27705 = r27703 / r27704;
        float r27706 = r27688 - r27705;
        float r27707 = r27694 ? r27700 : r27706;
        return r27707;
}

double f_od(double wj, double x) {
        double r27708 = 1;
        double r27709 = wj;
        double r27710 = r27708 / r27709;
        double r27711 = r27710 + r27709;
        double r27712 = fma(r27710, r27710, r27708);
        double r27713 = r27711 - r27712;
        double r27714 = -125981247538304.98;
        bool r27715 = r27713 <= r27714;
        double r27716 = x;
        double r27717 = 2;
        double r27718 = r27716 * r27717;
        double r27719 = -r27709;
        double r27720 = fma(r27709, r27709, r27716);
        double r27721 = fma(r27718, r27719, r27720);
        double r27722 = exp(r27709);
        double r27723 = -r27716;
        double r27724 = fma(r27722, r27709, r27723);
        double r27725 = fma(r27709, r27722, r27722);
        double r27726 = r27724 / r27725;
        double r27727 = r27709 - r27726;
        double r27728 = r27715 ? r27721 : r27727;
        return r27728;
}

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 r27729, r27730, r27731, r27732, r27733, r27734, r27735, r27736;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(848);
        mpfr_init(r27729);
        mpfr_init(r27730);
        mpfr_init(r27731);
        mpfr_init(r27732);
        mpfr_init(r27733);
        mpfr_init(r27734);
        mpfr_init(r27735);
        mpfr_init(r27736);
}

double f_im(double wj, double x) {
        mpfr_set_d(r27729, wj, MPFR_RNDN);
        mpfr_exp(r27730, r27729, MPFR_RNDN);
        mpfr_mul(r27731, r27729, r27730, MPFR_RNDN);
        mpfr_set_d(r27732, x, MPFR_RNDN);
        mpfr_sub(r27733, r27731, r27732, MPFR_RNDN);
        mpfr_add(r27734, r27730, r27731, MPFR_RNDN);
        mpfr_div(r27735, r27733, r27734, MPFR_RNDN);
        mpfr_sub(r27736, r27729, r27735, MPFR_RNDN);
        return mpfr_get_d(r27736, MPFR_RNDN);
}

static mpfr_t r27737, r27738, r27739, r27740, r27741, r27742, r27743, r27744, r27745, r27746, r27747, r27748, r27749, r27750, r27751, r27752, r27753, r27754, r27755, r27756, r27757;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(848);
        mpfr_init_set_str(r27737, "1", 10, MPFR_RNDN);
        mpfr_init(r27738);
        mpfr_init(r27739);
        mpfr_init(r27740);
        mpfr_init(r27741);
        mpfr_init(r27742);
        mpfr_init_set_str(r27743, "-125981247538304.98", 10, MPFR_RNDN);
        mpfr_init(r27744);
        mpfr_init(r27745);
        mpfr_init_set_str(r27746, "2", 10, MPFR_RNDN);
        mpfr_init(r27747);
        mpfr_init(r27748);
        mpfr_init(r27749);
        mpfr_init(r27750);
        mpfr_init(r27751);
        mpfr_init(r27752);
        mpfr_init(r27753);
        mpfr_init(r27754);
        mpfr_init(r27755);
        mpfr_init(r27756);
        mpfr_init(r27757);
}

double f_fm(double wj, double x) {
        ;
        mpfr_set_d(r27738, wj, MPFR_RNDN);
        mpfr_div(r27739, r27737, r27738, MPFR_RNDN);
        mpfr_add(r27740, r27739, r27738, MPFR_RNDN);
        mpfr_fma(r27741, r27739, r27739, r27737, MPFR_RNDN);
        mpfr_sub(r27742, r27740, r27741, MPFR_RNDN);
        ;
        mpfr_set_si(r27744, mpfr_cmp(r27742, r27743) <= 0, MPFR_RNDN);
        mpfr_set_d(r27745, x, MPFR_RNDN);
        ;
        mpfr_mul(r27747, r27745, r27746, MPFR_RNDN);
        mpfr_neg(r27748, r27738, MPFR_RNDN);
        mpfr_fma(r27749, r27738, r27738, r27745, MPFR_RNDN);
        mpfr_fma(r27750, r27747, r27748, r27749, MPFR_RNDN);
        mpfr_exp(r27751, r27738, MPFR_RNDN);
        mpfr_neg(r27752, r27745, MPFR_RNDN);
        mpfr_fma(r27753, r27751, r27738, r27752, MPFR_RNDN);
        mpfr_fma(r27754, r27738, r27751, r27751, MPFR_RNDN);
        mpfr_div(r27755, r27753, r27754, MPFR_RNDN);
        mpfr_sub(r27756, r27738, r27755, MPFR_RNDN);
        if (mpfr_get_si(r27744, MPFR_RNDN)) { mpfr_set(r27757, r27750, MPFR_RNDN); } else { mpfr_set(r27757, r27756, MPFR_RNDN); };
        return mpfr_get_d(r27757, MPFR_RNDN);
}

static mpfr_t r27758, r27759, r27760, r27761, r27762, r27763, r27764, r27765, r27766, r27767, r27768, r27769, r27770, r27771, r27772, r27773, r27774, r27775, r27776, r27777, r27778;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(848);
        mpfr_init_set_str(r27758, "1", 10, MPFR_RNDN);
        mpfr_init(r27759);
        mpfr_init(r27760);
        mpfr_init(r27761);
        mpfr_init(r27762);
        mpfr_init(r27763);
        mpfr_init_set_str(r27764, "-125981247538304.98", 10, MPFR_RNDN);
        mpfr_init(r27765);
        mpfr_init(r27766);
        mpfr_init_set_str(r27767, "2", 10, MPFR_RNDN);
        mpfr_init(r27768);
        mpfr_init(r27769);
        mpfr_init(r27770);
        mpfr_init(r27771);
        mpfr_init(r27772);
        mpfr_init(r27773);
        mpfr_init(r27774);
        mpfr_init(r27775);
        mpfr_init(r27776);
        mpfr_init(r27777);
        mpfr_init(r27778);
}

double f_dm(double wj, double x) {
        ;
        mpfr_set_d(r27759, wj, MPFR_RNDN);
        mpfr_div(r27760, r27758, r27759, MPFR_RNDN);
        mpfr_add(r27761, r27760, r27759, MPFR_RNDN);
        mpfr_fma(r27762, r27760, r27760, r27758, MPFR_RNDN);
        mpfr_sub(r27763, r27761, r27762, MPFR_RNDN);
        ;
        mpfr_set_si(r27765, mpfr_cmp(r27763, r27764) <= 0, MPFR_RNDN);
        mpfr_set_d(r27766, x, MPFR_RNDN);
        ;
        mpfr_mul(r27768, r27766, r27767, MPFR_RNDN);
        mpfr_neg(r27769, r27759, MPFR_RNDN);
        mpfr_fma(r27770, r27759, r27759, r27766, MPFR_RNDN);
        mpfr_fma(r27771, r27768, r27769, r27770, MPFR_RNDN);
        mpfr_exp(r27772, r27759, MPFR_RNDN);
        mpfr_neg(r27773, r27766, MPFR_RNDN);
        mpfr_fma(r27774, r27772, r27759, r27773, MPFR_RNDN);
        mpfr_fma(r27775, r27759, r27772, r27772, MPFR_RNDN);
        mpfr_div(r27776, r27774, r27775, MPFR_RNDN);
        mpfr_sub(r27777, r27759, r27776, MPFR_RNDN);
        if (mpfr_get_si(r27765, MPFR_RNDN)) { mpfr_set(r27778, r27771, MPFR_RNDN); } else { mpfr_set(r27778, r27777, MPFR_RNDN); };
        return mpfr_get_d(r27778, MPFR_RNDN);
}