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

char *name = "Asymptote C";

double f_if(float x) {
        float r26431 = x;
        float r26432 = 1;
        float r26433 = r26431 + r26432;
        float r26434 = r26431 / r26433;
        float r26435 = r26431 - r26432;
        float r26436 = r26433 / r26435;
        float r26437 = r26434 - r26436;
        return r26437;
}

double f_id(double x) {
        double r26438 = x;
        double r26439 = 1;
        double r26440 = r26438 + r26439;
        double r26441 = r26438 / r26440;
        double r26442 = r26438 - r26439;
        double r26443 = r26440 / r26442;
        double r26444 = r26441 - r26443;
        return r26444;
}


double f_of(float x) {
        float r26445 = x;
        float r26446 = 1;
        float r26447 = r26446 + r26445;
        float r26448 = r26445 / r26447;
        float r26449 = r26445 - r26446;
        float r26450 = r26447 / r26449;
        float r26451 = r26448 - r26450;
        float r26452 = 4.178112254527817e-10;
        bool r26453 = r26451 <= r26452;
        float r26454 = 3;
        float r26455 = r26454 / r26445;
        float r26456 = r26455 + r26446;
        float r26457 = -r26446;
        float r26458 = r26445 * r26445;
        float r26459 = r26457 / r26458;
        float r26460 = -r26455;
        float r26461 = fma(r26456, r26459, r26460);
        float r26462 = r26446 - r26445;
        float r26463 = fma(r26445, r26445, r26462);
        float r26464 = fma(r26458, r26445, r26446);
        float r26465 = r26445 / r26464;
        float r26466 = -r26447;
        float r26467 = r26466 / r26449;
        float r26468 = fma(r26463, r26465, r26467);
        float r26469 = r26453 ? r26461 : r26468;
        return r26469;
}

double f_od(double x) {
        double r26470 = x;
        double r26471 = 1;
        double r26472 = r26471 + r26470;
        double r26473 = r26470 / r26472;
        double r26474 = r26470 - r26471;
        double r26475 = r26472 / r26474;
        double r26476 = r26473 - r26475;
        double r26477 = 4.178112254527817e-10;
        bool r26478 = r26476 <= r26477;
        double r26479 = 3;
        double r26480 = r26479 / r26470;
        double r26481 = r26480 + r26471;
        double r26482 = -r26471;
        double r26483 = r26470 * r26470;
        double r26484 = r26482 / r26483;
        double r26485 = -r26480;
        double r26486 = fma(r26481, r26484, r26485);
        double r26487 = r26471 - r26470;
        double r26488 = fma(r26470, r26470, r26487);
        double r26489 = fma(r26483, r26470, r26471);
        double r26490 = r26470 / r26489;
        double r26491 = -r26472;
        double r26492 = r26491 / r26474;
        double r26493 = fma(r26488, r26490, r26492);
        double r26494 = r26478 ? r26486 : r26493;
        return r26494;
}

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 r26495, r26496, r26497, r26498, r26499, r26500, r26501;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26495);
        mpfr_init_set_str(r26496, "1", 10, MPFR_RNDN);
        mpfr_init(r26497);
        mpfr_init(r26498);
        mpfr_init(r26499);
        mpfr_init(r26500);
        mpfr_init(r26501);
}

double f_im(double x) {
        mpfr_set_d(r26495, x, MPFR_RNDN);
        ;
        mpfr_add(r26497, r26495, r26496, MPFR_RNDN);
        mpfr_div(r26498, r26495, r26497, MPFR_RNDN);
        mpfr_sub(r26499, r26495, r26496, MPFR_RNDN);
        mpfr_div(r26500, r26497, r26499, MPFR_RNDN);
        mpfr_sub(r26501, r26498, r26500, MPFR_RNDN);
        return mpfr_get_d(r26501, MPFR_RNDN);
}

static mpfr_t r26502, r26503, r26504, r26505, r26506, r26507, r26508, r26509, r26510, r26511, r26512, r26513, r26514, r26515, r26516, r26517, r26518, r26519, r26520, r26521, r26522, r26523, r26524, r26525, r26526;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26502);
        mpfr_init_set_str(r26503, "1", 10, MPFR_RNDN);
        mpfr_init(r26504);
        mpfr_init(r26505);
        mpfr_init(r26506);
        mpfr_init(r26507);
        mpfr_init(r26508);
        mpfr_init_set_str(r26509, "4.178112254527817e-10", 10, MPFR_RNDN);
        mpfr_init(r26510);
        mpfr_init_set_str(r26511, "3", 10, MPFR_RNDN);
        mpfr_init(r26512);
        mpfr_init(r26513);
        mpfr_init(r26514);
        mpfr_init(r26515);
        mpfr_init(r26516);
        mpfr_init(r26517);
        mpfr_init(r26518);
        mpfr_init(r26519);
        mpfr_init(r26520);
        mpfr_init(r26521);
        mpfr_init(r26522);
        mpfr_init(r26523);
        mpfr_init(r26524);
        mpfr_init(r26525);
        mpfr_init(r26526);
}

double f_fm(double x) {
        mpfr_set_d(r26502, x, MPFR_RNDN);
        ;
        mpfr_add(r26504, r26503, r26502, MPFR_RNDN);
        mpfr_div(r26505, r26502, r26504, MPFR_RNDN);
        mpfr_sub(r26506, r26502, r26503, MPFR_RNDN);
        mpfr_div(r26507, r26504, r26506, MPFR_RNDN);
        mpfr_sub(r26508, r26505, r26507, MPFR_RNDN);
        ;
        mpfr_set_si(r26510, mpfr_cmp(r26508, r26509) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r26512, r26511, r26502, MPFR_RNDN);
        mpfr_add(r26513, r26512, r26503, MPFR_RNDN);
        mpfr_neg(r26514, r26503, MPFR_RNDN);
        mpfr_mul(r26515, r26502, r26502, MPFR_RNDN);
        mpfr_div(r26516, r26514, r26515, MPFR_RNDN);
        mpfr_neg(r26517, r26512, MPFR_RNDN);
        mpfr_fma(r26518, r26513, r26516, r26517, MPFR_RNDN);
        mpfr_sub(r26519, r26503, r26502, MPFR_RNDN);
        mpfr_fma(r26520, r26502, r26502, r26519, MPFR_RNDN);
        mpfr_fma(r26521, r26515, r26502, r26503, MPFR_RNDN);
        mpfr_div(r26522, r26502, r26521, MPFR_RNDN);
        mpfr_neg(r26523, r26504, MPFR_RNDN);
        mpfr_div(r26524, r26523, r26506, MPFR_RNDN);
        mpfr_fma(r26525, r26520, r26522, r26524, MPFR_RNDN);
        if (mpfr_get_si(r26510, MPFR_RNDN)) { mpfr_set(r26526, r26518, MPFR_RNDN); } else { mpfr_set(r26526, r26525, MPFR_RNDN); };
        return mpfr_get_d(r26526, MPFR_RNDN);
}

static mpfr_t r26527, r26528, r26529, r26530, r26531, r26532, r26533, r26534, r26535, r26536, r26537, r26538, r26539, r26540, r26541, r26542, r26543, r26544, r26545, r26546, r26547, r26548, r26549, r26550, r26551;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26527);
        mpfr_init_set_str(r26528, "1", 10, MPFR_RNDN);
        mpfr_init(r26529);
        mpfr_init(r26530);
        mpfr_init(r26531);
        mpfr_init(r26532);
        mpfr_init(r26533);
        mpfr_init_set_str(r26534, "4.178112254527817e-10", 10, MPFR_RNDN);
        mpfr_init(r26535);
        mpfr_init_set_str(r26536, "3", 10, MPFR_RNDN);
        mpfr_init(r26537);
        mpfr_init(r26538);
        mpfr_init(r26539);
        mpfr_init(r26540);
        mpfr_init(r26541);
        mpfr_init(r26542);
        mpfr_init(r26543);
        mpfr_init(r26544);
        mpfr_init(r26545);
        mpfr_init(r26546);
        mpfr_init(r26547);
        mpfr_init(r26548);
        mpfr_init(r26549);
        mpfr_init(r26550);
        mpfr_init(r26551);
}

double f_dm(double x) {
        mpfr_set_d(r26527, x, MPFR_RNDN);
        ;
        mpfr_add(r26529, r26528, r26527, MPFR_RNDN);
        mpfr_div(r26530, r26527, r26529, MPFR_RNDN);
        mpfr_sub(r26531, r26527, r26528, MPFR_RNDN);
        mpfr_div(r26532, r26529, r26531, MPFR_RNDN);
        mpfr_sub(r26533, r26530, r26532, MPFR_RNDN);
        ;
        mpfr_set_si(r26535, mpfr_cmp(r26533, r26534) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r26537, r26536, r26527, MPFR_RNDN);
        mpfr_add(r26538, r26537, r26528, MPFR_RNDN);
        mpfr_neg(r26539, r26528, MPFR_RNDN);
        mpfr_mul(r26540, r26527, r26527, MPFR_RNDN);
        mpfr_div(r26541, r26539, r26540, MPFR_RNDN);
        mpfr_neg(r26542, r26537, MPFR_RNDN);
        mpfr_fma(r26543, r26538, r26541, r26542, MPFR_RNDN);
        mpfr_sub(r26544, r26528, r26527, MPFR_RNDN);
        mpfr_fma(r26545, r26527, r26527, r26544, MPFR_RNDN);
        mpfr_fma(r26546, r26540, r26527, r26528, MPFR_RNDN);
        mpfr_div(r26547, r26527, r26546, MPFR_RNDN);
        mpfr_neg(r26548, r26529, MPFR_RNDN);
        mpfr_div(r26549, r26548, r26531, MPFR_RNDN);
        mpfr_fma(r26550, r26545, r26547, r26549, MPFR_RNDN);
        if (mpfr_get_si(r26535, MPFR_RNDN)) { mpfr_set(r26551, r26543, MPFR_RNDN); } else { mpfr_set(r26551, r26550, MPFR_RNDN); };
        return mpfr_get_d(r26551, MPFR_RNDN);
}