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

char *name = "expq3 (problem 3.4.2)";

double f_if(float a, float b, float eps) {
        float r10133 = eps;
        float r10134 = a;
        float r10135 = b;
        float r10136 = r10134 + r10135;
        float r10137 = r10136 * r10133;
        float r10138 = exp(r10137);
        float r10139 = 1.0f;
        float r10140 = r10138 - r10139;
        float r10141 = r10133 * r10140;
        float r10142 = r10134 * r10133;
        float r10143 = exp(r10142);
        float r10144 = r10143 - r10139;
        float r10145 = r10135 * r10133;
        float r10146 = exp(r10145);
        float r10147 = r10146 - r10139;
        float r10148 = r10144 * r10147;
        float r10149 = r10141 / r10148;
        return r10149;
}

double f_id(double a, double b, double eps) {
        double r10150 = eps;
        double r10151 = a;
        double r10152 = b;
        double r10153 = r10151 + r10152;
        double r10154 = r10153 * r10150;
        double r10155 = exp(r10154);
        double r10156 = 1.0;
        double r10157 = r10155 - r10156;
        double r10158 = r10150 * r10157;
        double r10159 = r10151 * r10150;
        double r10160 = exp(r10159);
        double r10161 = r10160 - r10156;
        double r10162 = r10152 * r10150;
        double r10163 = exp(r10162);
        double r10164 = r10163 - r10156;
        double r10165 = r10161 * r10164;
        double r10166 = r10158 / r10165;
        return r10166;
}


double f_of(float a, float b, float eps) {
        float r10167 = a;
        float r10168 = b;
        float r10169 = r10167 + r10168;
        float r10170 = eps;
        float r10171 = r10169 * r10170;
        float r10172 = exp(r10171);
        float r10173 = 1.0f;
        float r10174 = r10172 - r10173;
        float r10175 = r10174 * r10170;
        float r10176 = r10170 * r10168;
        float r10177 = exp(r10176);
        float r10178 = r10177 - r10173;
        float r10179 = r10170 * r10167;
        float r10180 = exp(r10179);
        float r10181 = r10180 - r10173;
        float r10182 = r10178 * r10181;
        float r10183 = r10175 / r10182;
        float r10184 = -3.23912328517888e-14f;
        bool r10185 = r10183 <= r10184;
        float r10186 = r10173 / r10168;
        float r10187 = r10173 / r10167;
        float r10188 = r10186 + r10187;
        float r10189 = 4.3017334472689873e+132f;
        bool r10190 = r10183 <= r10189;
        float r10191 = expm1(r10179);
        float r10192 = r10170 / r10191;
        float r10193 = expm1(r10171);
        float r10194 = expm1(r10176);
        float r10195 = r10193 / r10194;
        float r10196 = r10192 * r10195;
        float r10197 = r10190 ? r10196 : r10188;
        float r10198 = r10185 ? r10188 : r10197;
        return r10198;
}

double f_od(double a, double b, double eps) {
        double r10199 = a;
        double r10200 = b;
        double r10201 = r10199 + r10200;
        double r10202 = eps;
        double r10203 = r10201 * r10202;
        double r10204 = exp(r10203);
        double r10205 = 1.0;
        double r10206 = r10204 - r10205;
        double r10207 = r10206 * r10202;
        double r10208 = r10202 * r10200;
        double r10209 = exp(r10208);
        double r10210 = r10209 - r10205;
        double r10211 = r10202 * r10199;
        double r10212 = exp(r10211);
        double r10213 = r10212 - r10205;
        double r10214 = r10210 * r10213;
        double r10215 = r10207 / r10214;
        double r10216 = -3.23912328517888e-14;
        bool r10217 = r10215 <= r10216;
        double r10218 = r10205 / r10200;
        double r10219 = r10205 / r10199;
        double r10220 = r10218 + r10219;
        double r10221 = 4.3017334472689873e+132;
        bool r10222 = r10215 <= r10221;
        double r10223 = expm1(r10211);
        double r10224 = r10202 / r10223;
        double r10225 = expm1(r10203);
        double r10226 = expm1(r10208);
        double r10227 = r10225 / r10226;
        double r10228 = r10224 * r10227;
        double r10229 = r10222 ? r10228 : r10220;
        double r10230 = r10217 ? r10220 : r10229;
        return r10230;
}

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 r10231, r10232, r10233, r10234, r10235, r10236, r10237, r10238, r10239, r10240, r10241, r10242, r10243, r10244, r10245, r10246, r10247;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init(r10231);
        mpfr_init(r10232);
        mpfr_init(r10233);
        mpfr_init(r10234);
        mpfr_init(r10235);
        mpfr_init(r10236);
        mpfr_init_set_str(r10237, "1", 10, MPFR_RNDN);
        mpfr_init(r10238);
        mpfr_init(r10239);
        mpfr_init(r10240);
        mpfr_init(r10241);
        mpfr_init(r10242);
        mpfr_init(r10243);
        mpfr_init(r10244);
        mpfr_init(r10245);
        mpfr_init(r10246);
        mpfr_init(r10247);
}

double f_im(double a, double b, double eps) {
        mpfr_set_d(r10231, eps, MPFR_RNDN);
        mpfr_set_d(r10232, a, MPFR_RNDN);
        mpfr_set_d(r10233, b, MPFR_RNDN);
        mpfr_add(r10234, r10232, r10233, MPFR_RNDN);
        mpfr_mul(r10235, r10234, r10231, MPFR_RNDN);
        mpfr_exp(r10236, r10235, MPFR_RNDN);
        ;
        mpfr_sub(r10238, r10236, r10237, MPFR_RNDN);
        mpfr_mul(r10239, r10231, r10238, MPFR_RNDN);
        mpfr_mul(r10240, r10232, r10231, MPFR_RNDN);
        mpfr_exp(r10241, r10240, MPFR_RNDN);
        mpfr_sub(r10242, r10241, r10237, MPFR_RNDN);
        mpfr_mul(r10243, r10233, r10231, MPFR_RNDN);
        mpfr_exp(r10244, r10243, MPFR_RNDN);
        mpfr_sub(r10245, r10244, r10237, MPFR_RNDN);
        mpfr_mul(r10246, r10242, r10245, MPFR_RNDN);
        mpfr_div(r10247, r10239, r10246, MPFR_RNDN);
        return mpfr_get_d(r10247, MPFR_RNDN);
}

static mpfr_t r10248, r10249, r10250, r10251, r10252, r10253, r10254, r10255, r10256, r10257, r10258, r10259, r10260, r10261, r10262, r10263, r10264, r10265, r10266, r10267, r10268, r10269, r10270, r10271, r10272, r10273, r10274, r10275, r10276, r10277, r10278, r10279;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r10248);
        mpfr_init(r10249);
        mpfr_init(r10250);
        mpfr_init(r10251);
        mpfr_init(r10252);
        mpfr_init(r10253);
        mpfr_init_set_str(r10254, "1", 10, MPFR_RNDN);
        mpfr_init(r10255);
        mpfr_init(r10256);
        mpfr_init(r10257);
        mpfr_init(r10258);
        mpfr_init(r10259);
        mpfr_init(r10260);
        mpfr_init(r10261);
        mpfr_init(r10262);
        mpfr_init(r10263);
        mpfr_init(r10264);
        mpfr_init_set_str(r10265, "-3.23912328517888e-14", 10, MPFR_RNDN);
        mpfr_init(r10266);
        mpfr_init(r10267);
        mpfr_init(r10268);
        mpfr_init(r10269);
        mpfr_init_set_str(r10270, "4.3017334472689873e+132", 10, MPFR_RNDN);
        mpfr_init(r10271);
        mpfr_init(r10272);
        mpfr_init(r10273);
        mpfr_init(r10274);
        mpfr_init(r10275);
        mpfr_init(r10276);
        mpfr_init(r10277);
        mpfr_init(r10278);
        mpfr_init(r10279);
}

double f_fm(double a, double b, double eps) {
        mpfr_set_d(r10248, a, MPFR_RNDN);
        mpfr_set_d(r10249, b, MPFR_RNDN);
        mpfr_add(r10250, r10248, r10249, MPFR_RNDN);
        mpfr_set_d(r10251, eps, MPFR_RNDN);
        mpfr_mul(r10252, r10250, r10251, MPFR_RNDN);
        mpfr_exp(r10253, r10252, MPFR_RNDN);
        ;
        mpfr_sub(r10255, r10253, r10254, MPFR_RNDN);
        mpfr_mul(r10256, r10255, r10251, MPFR_RNDN);
        mpfr_mul(r10257, r10251, r10249, MPFR_RNDN);
        mpfr_exp(r10258, r10257, MPFR_RNDN);
        mpfr_sub(r10259, r10258, r10254, MPFR_RNDN);
        mpfr_mul(r10260, r10251, r10248, MPFR_RNDN);
        mpfr_exp(r10261, r10260, MPFR_RNDN);
        mpfr_sub(r10262, r10261, r10254, MPFR_RNDN);
        mpfr_mul(r10263, r10259, r10262, MPFR_RNDN);
        mpfr_div(r10264, r10256, r10263, MPFR_RNDN);
        ;
        mpfr_set_si(r10266, mpfr_cmp(r10264, r10265) <= 0, MPFR_RNDN);
        mpfr_div(r10267, r10254, r10249, MPFR_RNDN);
        mpfr_div(r10268, r10254, r10248, MPFR_RNDN);
        mpfr_add(r10269, r10267, r10268, MPFR_RNDN);
        ;
        mpfr_set_si(r10271, mpfr_cmp(r10264, r10270) <= 0, MPFR_RNDN);
        mpfr_expm1(r10272, r10260, MPFR_RNDN);
        mpfr_div(r10273, r10251, r10272, MPFR_RNDN);
        mpfr_expm1(r10274, r10252, MPFR_RNDN);
        mpfr_expm1(r10275, r10257, MPFR_RNDN);
        mpfr_div(r10276, r10274, r10275, MPFR_RNDN);
        mpfr_mul(r10277, r10273, r10276, MPFR_RNDN);
        if (mpfr_get_si(r10271, MPFR_RNDN)) { mpfr_set(r10278, r10277, MPFR_RNDN); } else { mpfr_set(r10278, r10269, MPFR_RNDN); };
        if (mpfr_get_si(r10266, MPFR_RNDN)) { mpfr_set(r10279, r10269, MPFR_RNDN); } else { mpfr_set(r10279, r10278, MPFR_RNDN); };
        return mpfr_get_d(r10279, MPFR_RNDN);
}

static mpfr_t r10280, r10281, r10282, r10283, r10284, r10285, r10286, r10287, r10288, r10289, r10290, r10291, r10292, r10293, r10294, r10295, r10296, r10297, r10298, r10299, r10300, r10301, r10302, r10303, r10304, r10305, r10306, r10307, r10308, r10309, r10310, r10311;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r10280);
        mpfr_init(r10281);
        mpfr_init(r10282);
        mpfr_init(r10283);
        mpfr_init(r10284);
        mpfr_init(r10285);
        mpfr_init_set_str(r10286, "1", 10, MPFR_RNDN);
        mpfr_init(r10287);
        mpfr_init(r10288);
        mpfr_init(r10289);
        mpfr_init(r10290);
        mpfr_init(r10291);
        mpfr_init(r10292);
        mpfr_init(r10293);
        mpfr_init(r10294);
        mpfr_init(r10295);
        mpfr_init(r10296);
        mpfr_init_set_str(r10297, "-3.23912328517888e-14", 10, MPFR_RNDN);
        mpfr_init(r10298);
        mpfr_init(r10299);
        mpfr_init(r10300);
        mpfr_init(r10301);
        mpfr_init_set_str(r10302, "4.3017334472689873e+132", 10, MPFR_RNDN);
        mpfr_init(r10303);
        mpfr_init(r10304);
        mpfr_init(r10305);
        mpfr_init(r10306);
        mpfr_init(r10307);
        mpfr_init(r10308);
        mpfr_init(r10309);
        mpfr_init(r10310);
        mpfr_init(r10311);
}

double f_dm(double a, double b, double eps) {
        mpfr_set_d(r10280, a, MPFR_RNDN);
        mpfr_set_d(r10281, b, MPFR_RNDN);
        mpfr_add(r10282, r10280, r10281, MPFR_RNDN);
        mpfr_set_d(r10283, eps, MPFR_RNDN);
        mpfr_mul(r10284, r10282, r10283, MPFR_RNDN);
        mpfr_exp(r10285, r10284, MPFR_RNDN);
        ;
        mpfr_sub(r10287, r10285, r10286, MPFR_RNDN);
        mpfr_mul(r10288, r10287, r10283, MPFR_RNDN);
        mpfr_mul(r10289, r10283, r10281, MPFR_RNDN);
        mpfr_exp(r10290, r10289, MPFR_RNDN);
        mpfr_sub(r10291, r10290, r10286, MPFR_RNDN);
        mpfr_mul(r10292, r10283, r10280, MPFR_RNDN);
        mpfr_exp(r10293, r10292, MPFR_RNDN);
        mpfr_sub(r10294, r10293, r10286, MPFR_RNDN);
        mpfr_mul(r10295, r10291, r10294, MPFR_RNDN);
        mpfr_div(r10296, r10288, r10295, MPFR_RNDN);
        ;
        mpfr_set_si(r10298, mpfr_cmp(r10296, r10297) <= 0, MPFR_RNDN);
        mpfr_div(r10299, r10286, r10281, MPFR_RNDN);
        mpfr_div(r10300, r10286, r10280, MPFR_RNDN);
        mpfr_add(r10301, r10299, r10300, MPFR_RNDN);
        ;
        mpfr_set_si(r10303, mpfr_cmp(r10296, r10302) <= 0, MPFR_RNDN);
        mpfr_expm1(r10304, r10292, MPFR_RNDN);
        mpfr_div(r10305, r10283, r10304, MPFR_RNDN);
        mpfr_expm1(r10306, r10284, MPFR_RNDN);
        mpfr_expm1(r10307, r10289, MPFR_RNDN);
        mpfr_div(r10308, r10306, r10307, MPFR_RNDN);
        mpfr_mul(r10309, r10305, r10308, MPFR_RNDN);
        if (mpfr_get_si(r10303, MPFR_RNDN)) { mpfr_set(r10310, r10309, MPFR_RNDN); } else { mpfr_set(r10310, r10301, MPFR_RNDN); };
        if (mpfr_get_si(r10298, MPFR_RNDN)) { mpfr_set(r10311, r10301, MPFR_RNDN); } else { mpfr_set(r10311, r10310, MPFR_RNDN); };
        return mpfr_get_d(r10311, MPFR_RNDN);
}