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

char *name = "Logistic regression 2";

double f_if(float x, float y) {
        float r27336 = 1;
        float r27337 = x;
        float r27338 = exp(r27337);
        float r27339 = r27336 + r27338;
        float r27340 = log(r27339);
        float r27341 = y;
        float r27342 = r27337 * r27341;
        float r27343 = r27340 - r27342;
        return r27343;
}

double f_id(double x, double y) {
        double r27344 = 1;
        double r27345 = x;
        double r27346 = exp(r27345);
        double r27347 = r27344 + r27346;
        double r27348 = log(r27347);
        double r27349 = y;
        double r27350 = r27345 * r27349;
        double r27351 = r27348 - r27350;
        return r27351;
}


double f_of(float x, float y) {
        float r27352 = x;
        float r27353 = exp(r27352);
        float r27354 = log1p(r27353);
        float r27355 = y;
        float r27356 = r27355 * r27352;
        float r27357 = r27354 - r27356;
        return r27357;
}

double f_od(double x, double y) {
        double r27358 = x;
        double r27359 = exp(r27358);
        double r27360 = log1p(r27359);
        double r27361 = y;
        double r27362 = r27361 * r27358;
        double r27363 = r27360 - r27362;
        return r27363;
}

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 r27364, r27365, r27366, r27367, r27368, r27369, r27370, r27371;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r27364, "1", 10, MPFR_RNDN);
        mpfr_init(r27365);
        mpfr_init(r27366);
        mpfr_init(r27367);
        mpfr_init(r27368);
        mpfr_init(r27369);
        mpfr_init(r27370);
        mpfr_init(r27371);
}

double f_im(double x, double y) {
        ;
        mpfr_set_d(r27365, x, MPFR_RNDN);
        mpfr_exp(r27366, r27365, MPFR_RNDN);
        mpfr_add(r27367, r27364, r27366, MPFR_RNDN);
        mpfr_log(r27368, r27367, MPFR_RNDN);
        mpfr_set_d(r27369, y, MPFR_RNDN);
        mpfr_mul(r27370, r27365, r27369, MPFR_RNDN);
        mpfr_sub(r27371, r27368, r27370, MPFR_RNDN);
        return mpfr_get_d(r27371, MPFR_RNDN);
}

static mpfr_t r27372, r27373, r27374, r27375, r27376, r27377;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r27372);
        mpfr_init(r27373);
        mpfr_init(r27374);
        mpfr_init(r27375);
        mpfr_init(r27376);
        mpfr_init(r27377);
}

double f_fm(double x, double y) {
        mpfr_set_d(r27372, x, MPFR_RNDN);
        mpfr_exp(r27373, r27372, MPFR_RNDN);
        mpfr_log1p(r27374, r27373, MPFR_RNDN);
        mpfr_set_d(r27375, y, MPFR_RNDN);
        mpfr_mul(r27376, r27375, r27372, MPFR_RNDN);
        mpfr_sub(r27377, r27374, r27376, MPFR_RNDN);
        return mpfr_get_d(r27377, MPFR_RNDN);
}

static mpfr_t r27378, r27379, r27380, r27381, r27382, r27383;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r27378);
        mpfr_init(r27379);
        mpfr_init(r27380);
        mpfr_init(r27381);
        mpfr_init(r27382);
        mpfr_init(r27383);
}

double f_dm(double x, double y) {
        mpfr_set_d(r27378, x, MPFR_RNDN);
        mpfr_exp(r27379, r27378, MPFR_RNDN);
        mpfr_log1p(r27380, r27379, MPFR_RNDN);
        mpfr_set_d(r27381, y, MPFR_RNDN);
        mpfr_mul(r27382, r27381, r27378, MPFR_RNDN);
        mpfr_sub(r27383, r27380, r27382, MPFR_RNDN);
        return mpfr_get_d(r27383, MPFR_RNDN);
}