#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 r46844 = 1;
        float r46845 = x;
        float r46846 = exp(r46845);
        float r46847 = r46844 + r46846;
        float r46848 = log(r46847);
        float r46849 = y;
        float r46850 = r46845 * r46849;
        float r46851 = r46848 - r46850;
        return r46851;
}

double f_id(double x, double y) {
        double r46852 = 1;
        double r46853 = x;
        double r46854 = exp(r46853);
        double r46855 = r46852 + r46854;
        double r46856 = log(r46855);
        double r46857 = y;
        double r46858 = r46853 * r46857;
        double r46859 = r46856 - r46858;
        return r46859;
}


double f_of(float x, float y) {
        float r46860 = 1;
        float r46861 = x;
        float r46862 = exp(r46861);
        float r46863 = 3;
        float r46864 = pow(r46862, r46863);
        float r46865 = r46860 + r46864;
        float r46866 = r46861 + r46861;
        float r46867 = exp(r46866);
        float r46868 = r46860 - r46862;
        float r46869 = r46867 + r46868;
        float r46870 = r46865 / r46869;
        float r46871 = log(r46870);
        float r46872 = y;
        float r46873 = r46861 * r46872;
        float r46874 = r46871 - r46873;
        return r46874;
}

double f_od(double x, double y) {
        double r46875 = 1;
        double r46876 = x;
        double r46877 = exp(r46876);
        double r46878 = 3;
        double r46879 = pow(r46877, r46878);
        double r46880 = r46875 + r46879;
        double r46881 = r46876 + r46876;
        double r46882 = exp(r46881);
        double r46883 = r46875 - r46877;
        double r46884 = r46882 + r46883;
        double r46885 = r46880 / r46884;
        double r46886 = log(r46885);
        double r46887 = y;
        double r46888 = r46876 * r46887;
        double r46889 = r46886 - r46888;
        return r46889;
}

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 r46890, r46891, r46892, r46893, r46894, r46895, r46896, r46897;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(656);
        mpfr_init_set_str(r46890, "1", 10, MPFR_RNDN);
        mpfr_init(r46891);
        mpfr_init(r46892);
        mpfr_init(r46893);
        mpfr_init(r46894);
        mpfr_init(r46895);
        mpfr_init(r46896);
        mpfr_init(r46897);
}

double f_im(double x, double y) {
        ;
        mpfr_set_d(r46891, x, MPFR_RNDN);
        mpfr_exp(r46892, r46891, MPFR_RNDN);
        mpfr_add(r46893, r46890, r46892, MPFR_RNDN);
        mpfr_log(r46894, r46893, MPFR_RNDN);
        mpfr_set_d(r46895, y, MPFR_RNDN);
        mpfr_mul(r46896, r46891, r46895, MPFR_RNDN);
        mpfr_sub(r46897, r46894, r46896, MPFR_RNDN);
        return mpfr_get_d(r46897, MPFR_RNDN);
}

static mpfr_t r46898, r46899, r46900, r46901, r46902, r46903, r46904, r46905, r46906, r46907, r46908, r46909, r46910, r46911, r46912;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(656);
        mpfr_init_set_str(r46898, "1", 10, MPFR_RNDN);
        mpfr_init(r46899);
        mpfr_init(r46900);
        mpfr_init_set_str(r46901, "3", 10, MPFR_RNDN);
        mpfr_init(r46902);
        mpfr_init(r46903);
        mpfr_init(r46904);
        mpfr_init(r46905);
        mpfr_init(r46906);
        mpfr_init(r46907);
        mpfr_init(r46908);
        mpfr_init(r46909);
        mpfr_init(r46910);
        mpfr_init(r46911);
        mpfr_init(r46912);
}

double f_fm(double x, double y) {
        ;
        mpfr_set_d(r46899, x, MPFR_RNDN);
        mpfr_exp(r46900, r46899, MPFR_RNDN);
        ;
        mpfr_pow(r46902, r46900, r46901, MPFR_RNDN);
        mpfr_add(r46903, r46898, r46902, MPFR_RNDN);
        mpfr_add(r46904, r46899, r46899, MPFR_RNDN);
        mpfr_exp(r46905, r46904, MPFR_RNDN);
        mpfr_sub(r46906, r46898, r46900, MPFR_RNDN);
        mpfr_add(r46907, r46905, r46906, MPFR_RNDN);
        mpfr_div(r46908, r46903, r46907, MPFR_RNDN);
        mpfr_log(r46909, r46908, MPFR_RNDN);
        mpfr_set_d(r46910, y, MPFR_RNDN);
        mpfr_mul(r46911, r46899, r46910, MPFR_RNDN);
        mpfr_sub(r46912, r46909, r46911, MPFR_RNDN);
        return mpfr_get_d(r46912, MPFR_RNDN);
}

static mpfr_t r46913, r46914, r46915, r46916, r46917, r46918, r46919, r46920, r46921, r46922, r46923, r46924, r46925, r46926, r46927;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(656);
        mpfr_init_set_str(r46913, "1", 10, MPFR_RNDN);
        mpfr_init(r46914);
        mpfr_init(r46915);
        mpfr_init_set_str(r46916, "3", 10, MPFR_RNDN);
        mpfr_init(r46917);
        mpfr_init(r46918);
        mpfr_init(r46919);
        mpfr_init(r46920);
        mpfr_init(r46921);
        mpfr_init(r46922);
        mpfr_init(r46923);
        mpfr_init(r46924);
        mpfr_init(r46925);
        mpfr_init(r46926);
        mpfr_init(r46927);
}

double f_dm(double x, double y) {
        ;
        mpfr_set_d(r46914, x, MPFR_RNDN);
        mpfr_exp(r46915, r46914, MPFR_RNDN);
        ;
        mpfr_pow(r46917, r46915, r46916, MPFR_RNDN);
        mpfr_add(r46918, r46913, r46917, MPFR_RNDN);
        mpfr_add(r46919, r46914, r46914, MPFR_RNDN);
        mpfr_exp(r46920, r46919, MPFR_RNDN);
        mpfr_sub(r46921, r46913, r46915, MPFR_RNDN);
        mpfr_add(r46922, r46920, r46921, MPFR_RNDN);
        mpfr_div(r46923, r46918, r46922, MPFR_RNDN);
        mpfr_log(r46924, r46923, MPFR_RNDN);
        mpfr_set_d(r46925, y, MPFR_RNDN);
        mpfr_mul(r46926, r46914, r46925, MPFR_RNDN);
        mpfr_sub(r46927, r46924, r46926, MPFR_RNDN);
        return mpfr_get_d(r46927, MPFR_RNDN);
}