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

char *name = "Compound Interest";

double f_if(float i, float n) {
        float r16340 = 100.0f;
        float r16341 = 1.0f;
        float r16342 = i;
        float r16343 = n;
        float r16344 = r16342 / r16343;
        float r16345 = r16341 + r16344;
        float r16346 = pow(r16345, r16343);
        float r16347 = r16346 - r16341;
        float r16348 = r16347 / r16344;
        float r16349 = r16340 * r16348;
        return r16349;
}

double f_id(double i, double n) {
        double r16350 = 100.0;
        double r16351 = 1.0;
        double r16352 = i;
        double r16353 = n;
        double r16354 = r16352 / r16353;
        double r16355 = r16351 + r16354;
        double r16356 = pow(r16355, r16353);
        double r16357 = r16356 - r16351;
        double r16358 = r16357 / r16354;
        double r16359 = r16350 * r16358;
        return r16359;
}


double f_of(float i, float n) {
        float r16360 = i;
        float r16361 = -1.180495039389048e-113f;
        bool r16362 = r16360 <= r16361;
        float r16363 = n;
        float r16364 = 100.0f;
        float r16365 = r16363 * r16364;
        float r16366 = r16365 / r16360;
        float r16367 = r16360 / r16363;
        float r16368 = log1p(r16367);
        float r16369 = r16368 * r16363;
        float r16370 = expm1(r16369);
        float r16371 = r16366 * r16370;
        float r16372 = 1.9365931721357208e-122f;
        bool r16373 = r16360 <= r16372;
        float r16374 = 0.5f;
        float r16375 = r16374 * r16360;
        float r16376 = fma(r16360, r16375, r16360);
        float r16377 = r16376 / r16360;
        float r16378 = r16364 * r16377;
        float r16379 = r16378 * r16363;
        float r16380 = 3.135592761500106e-05f;
        bool r16381 = r16360 <= r16380;
        float r16382 = log(r16363);
        float r16383 = -r16382;
        float r16384 = log(r16360);
        float r16385 = -r16384;
        float r16386 = r16383 - r16385;
        float r16387 = r16386 * r16363;
        float r16388 = expm1(r16387);
        float r16389 = r16366 * r16388;
        float r16390 = r16381 ? r16371 : r16389;
        float r16391 = r16373 ? r16379 : r16390;
        float r16392 = r16362 ? r16371 : r16391;
        return r16392;
}

double f_od(double i, double n) {
        double r16393 = i;
        double r16394 = -1.180495039389048e-113;
        bool r16395 = r16393 <= r16394;
        double r16396 = n;
        double r16397 = 100.0;
        double r16398 = r16396 * r16397;
        double r16399 = r16398 / r16393;
        double r16400 = r16393 / r16396;
        double r16401 = log1p(r16400);
        double r16402 = r16401 * r16396;
        double r16403 = expm1(r16402);
        double r16404 = r16399 * r16403;
        double r16405 = 1.9365931721357208e-122;
        bool r16406 = r16393 <= r16405;
        double r16407 = 0.5;
        double r16408 = r16407 * r16393;
        double r16409 = fma(r16393, r16408, r16393);
        double r16410 = r16409 / r16393;
        double r16411 = r16397 * r16410;
        double r16412 = r16411 * r16396;
        double r16413 = 3.135592761500106e-05;
        bool r16414 = r16393 <= r16413;
        double r16415 = log(r16396);
        double r16416 = -r16415;
        double r16417 = log(r16393);
        double r16418 = -r16417;
        double r16419 = r16416 - r16418;
        double r16420 = r16419 * r16396;
        double r16421 = expm1(r16420);
        double r16422 = r16399 * r16421;
        double r16423 = r16414 ? r16404 : r16422;
        double r16424 = r16406 ? r16412 : r16423;
        double r16425 = r16395 ? r16404 : r16424;
        return r16425;
}

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 r16426, r16427, r16428, r16429, r16430, r16431, r16432, r16433, r16434, r16435;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r16426, "100", 10, MPFR_RNDN);
        mpfr_init_set_str(r16427, "1", 10, MPFR_RNDN);
        mpfr_init(r16428);
        mpfr_init(r16429);
        mpfr_init(r16430);
        mpfr_init(r16431);
        mpfr_init(r16432);
        mpfr_init(r16433);
        mpfr_init(r16434);
        mpfr_init(r16435);
}

double f_im(double i, double n) {
        ;
        ;
        mpfr_set_d(r16428, i, MPFR_RNDN);
        mpfr_set_d(r16429, n, MPFR_RNDN);
        mpfr_div(r16430, r16428, r16429, MPFR_RNDN);
        mpfr_add(r16431, r16427, r16430, MPFR_RNDN);
        mpfr_pow(r16432, r16431, r16429, MPFR_RNDN);
        mpfr_sub(r16433, r16432, r16427, MPFR_RNDN);
        mpfr_div(r16434, r16433, r16430, MPFR_RNDN);
        mpfr_mul(r16435, r16426, r16434, MPFR_RNDN);
        return mpfr_get_d(r16435, MPFR_RNDN);
}

static mpfr_t r16436, r16437, r16438, r16439, r16440, r16441, r16442, r16443, r16444, r16445, r16446, r16447, r16448, r16449, r16450, r16451, r16452, r16453, r16454, r16455, r16456, r16457, r16458, r16459, r16460, r16461, r16462, r16463, r16464, r16465, r16466, r16467, r16468;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r16436);
        mpfr_init_set_str(r16437, "-1.180495039389048e-113", 10, MPFR_RNDN);
        mpfr_init(r16438);
        mpfr_init(r16439);
        mpfr_init_set_str(r16440, "100", 10, MPFR_RNDN);
        mpfr_init(r16441);
        mpfr_init(r16442);
        mpfr_init(r16443);
        mpfr_init(r16444);
        mpfr_init(r16445);
        mpfr_init(r16446);
        mpfr_init(r16447);
        mpfr_init_set_str(r16448, "1.9365931721357208e-122", 10, MPFR_RNDN);
        mpfr_init(r16449);
        mpfr_init_set_str(r16450, "1/2", 10, MPFR_RNDN);
        mpfr_init(r16451);
        mpfr_init(r16452);
        mpfr_init(r16453);
        mpfr_init(r16454);
        mpfr_init(r16455);
        mpfr_init_set_str(r16456, "3.135592761500106e-05", 10, MPFR_RNDN);
        mpfr_init(r16457);
        mpfr_init(r16458);
        mpfr_init(r16459);
        mpfr_init(r16460);
        mpfr_init(r16461);
        mpfr_init(r16462);
        mpfr_init(r16463);
        mpfr_init(r16464);
        mpfr_init(r16465);
        mpfr_init(r16466);
        mpfr_init(r16467);
        mpfr_init(r16468);
}

double f_fm(double i, double n) {
        mpfr_set_d(r16436, i, MPFR_RNDN);
        ;
        mpfr_set_si(r16438, mpfr_cmp(r16436, r16437) <= 0, MPFR_RNDN);
        mpfr_set_d(r16439, n, MPFR_RNDN);
        ;
        mpfr_mul(r16441, r16439, r16440, MPFR_RNDN);
        mpfr_div(r16442, r16441, r16436, MPFR_RNDN);
        mpfr_div(r16443, r16436, r16439, MPFR_RNDN);
        mpfr_log1p(r16444, r16443, MPFR_RNDN);
        mpfr_mul(r16445, r16444, r16439, MPFR_RNDN);
        mpfr_expm1(r16446, r16445, MPFR_RNDN);
        mpfr_mul(r16447, r16442, r16446, MPFR_RNDN);
        ;
        mpfr_set_si(r16449, mpfr_cmp(r16436, r16448) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r16451, r16450, r16436, MPFR_RNDN);
        mpfr_fma(r16452, r16436, r16451, r16436, MPFR_RNDN);
        mpfr_div(r16453, r16452, r16436, MPFR_RNDN);
        mpfr_mul(r16454, r16440, r16453, MPFR_RNDN);
        mpfr_mul(r16455, r16454, r16439, MPFR_RNDN);
        ;
        mpfr_set_si(r16457, mpfr_cmp(r16436, r16456) <= 0, MPFR_RNDN);
        mpfr_log(r16458, r16439, MPFR_RNDN);
        mpfr_neg(r16459, r16458, MPFR_RNDN);
        mpfr_log(r16460, r16436, MPFR_RNDN);
        mpfr_neg(r16461, r16460, MPFR_RNDN);
        mpfr_sub(r16462, r16459, r16461, MPFR_RNDN);
        mpfr_mul(r16463, r16462, r16439, MPFR_RNDN);
        mpfr_expm1(r16464, r16463, MPFR_RNDN);
        mpfr_mul(r16465, r16442, r16464, MPFR_RNDN);
        if (mpfr_get_si(r16457, MPFR_RNDN)) { mpfr_set(r16466, r16447, MPFR_RNDN); } else { mpfr_set(r16466, r16465, MPFR_RNDN); };
        if (mpfr_get_si(r16449, MPFR_RNDN)) { mpfr_set(r16467, r16455, MPFR_RNDN); } else { mpfr_set(r16467, r16466, MPFR_RNDN); };
        if (mpfr_get_si(r16438, MPFR_RNDN)) { mpfr_set(r16468, r16447, MPFR_RNDN); } else { mpfr_set(r16468, r16467, MPFR_RNDN); };
        return mpfr_get_d(r16468, MPFR_RNDN);
}

static mpfr_t r16469, r16470, r16471, r16472, r16473, r16474, r16475, r16476, r16477, r16478, r16479, r16480, r16481, r16482, r16483, r16484, r16485, r16486, r16487, r16488, r16489, r16490, r16491, r16492, r16493, r16494, r16495, r16496, r16497, r16498, r16499, r16500, r16501;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r16469);
        mpfr_init_set_str(r16470, "-1.180495039389048e-113", 10, MPFR_RNDN);
        mpfr_init(r16471);
        mpfr_init(r16472);
        mpfr_init_set_str(r16473, "100", 10, MPFR_RNDN);
        mpfr_init(r16474);
        mpfr_init(r16475);
        mpfr_init(r16476);
        mpfr_init(r16477);
        mpfr_init(r16478);
        mpfr_init(r16479);
        mpfr_init(r16480);
        mpfr_init_set_str(r16481, "1.9365931721357208e-122", 10, MPFR_RNDN);
        mpfr_init(r16482);
        mpfr_init_set_str(r16483, "1/2", 10, MPFR_RNDN);
        mpfr_init(r16484);
        mpfr_init(r16485);
        mpfr_init(r16486);
        mpfr_init(r16487);
        mpfr_init(r16488);
        mpfr_init_set_str(r16489, "3.135592761500106e-05", 10, MPFR_RNDN);
        mpfr_init(r16490);
        mpfr_init(r16491);
        mpfr_init(r16492);
        mpfr_init(r16493);
        mpfr_init(r16494);
        mpfr_init(r16495);
        mpfr_init(r16496);
        mpfr_init(r16497);
        mpfr_init(r16498);
        mpfr_init(r16499);
        mpfr_init(r16500);
        mpfr_init(r16501);
}

double f_dm(double i, double n) {
        mpfr_set_d(r16469, i, MPFR_RNDN);
        ;
        mpfr_set_si(r16471, mpfr_cmp(r16469, r16470) <= 0, MPFR_RNDN);
        mpfr_set_d(r16472, n, MPFR_RNDN);
        ;
        mpfr_mul(r16474, r16472, r16473, MPFR_RNDN);
        mpfr_div(r16475, r16474, r16469, MPFR_RNDN);
        mpfr_div(r16476, r16469, r16472, MPFR_RNDN);
        mpfr_log1p(r16477, r16476, MPFR_RNDN);
        mpfr_mul(r16478, r16477, r16472, MPFR_RNDN);
        mpfr_expm1(r16479, r16478, MPFR_RNDN);
        mpfr_mul(r16480, r16475, r16479, MPFR_RNDN);
        ;
        mpfr_set_si(r16482, mpfr_cmp(r16469, r16481) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r16484, r16483, r16469, MPFR_RNDN);
        mpfr_fma(r16485, r16469, r16484, r16469, MPFR_RNDN);
        mpfr_div(r16486, r16485, r16469, MPFR_RNDN);
        mpfr_mul(r16487, r16473, r16486, MPFR_RNDN);
        mpfr_mul(r16488, r16487, r16472, MPFR_RNDN);
        ;
        mpfr_set_si(r16490, mpfr_cmp(r16469, r16489) <= 0, MPFR_RNDN);
        mpfr_log(r16491, r16472, MPFR_RNDN);
        mpfr_neg(r16492, r16491, MPFR_RNDN);
        mpfr_log(r16493, r16469, MPFR_RNDN);
        mpfr_neg(r16494, r16493, MPFR_RNDN);
        mpfr_sub(r16495, r16492, r16494, MPFR_RNDN);
        mpfr_mul(r16496, r16495, r16472, MPFR_RNDN);
        mpfr_expm1(r16497, r16496, MPFR_RNDN);
        mpfr_mul(r16498, r16475, r16497, MPFR_RNDN);
        if (mpfr_get_si(r16490, MPFR_RNDN)) { mpfr_set(r16499, r16480, MPFR_RNDN); } else { mpfr_set(r16499, r16498, MPFR_RNDN); };
        if (mpfr_get_si(r16482, MPFR_RNDN)) { mpfr_set(r16500, r16488, MPFR_RNDN); } else { mpfr_set(r16500, r16499, MPFR_RNDN); };
        if (mpfr_get_si(r16471, MPFR_RNDN)) { mpfr_set(r16501, r16480, MPFR_RNDN); } else { mpfr_set(r16501, r16500, MPFR_RNDN); };
        return mpfr_get_d(r16501, MPFR_RNDN);
}

