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

char *name = "Graphics.Rasterific.Svg.PathConverter:arcToSegments from rasterific-svg-0.2.3.1";

double f_if(float x, float y, float z, float t) {
        float r50698 = x;
        float r50699 = r50698 * r50698;
        float r50700 = y;
        float r50701 = r50700 * r50700;
        float r50702 = r50699 / r50701;
        float r50703 = z;
        float r50704 = r50703 * r50703;
        float r50705 = t;
        float r50706 = r50705 * r50705;
        float r50707 = r50704 / r50706;
        float r50708 = r50702 + r50707;
        return r50708;
}

double f_id(double x, double y, double z, double t) {
        double r50709 = x;
        double r50710 = r50709 * r50709;
        double r50711 = y;
        double r50712 = r50711 * r50711;
        double r50713 = r50710 / r50712;
        double r50714 = z;
        double r50715 = r50714 * r50714;
        double r50716 = t;
        double r50717 = r50716 * r50716;
        double r50718 = r50715 / r50717;
        double r50719 = r50713 + r50718;
        return r50719;
}


double f_of(float x, float y, float z, float t) {
        float r50720 = x;
        float r50721 = y;
        float r50722 = r50720 / r50721;
        float r50723 = fabs(r50722);
        float r50724 = sqrt(r50723);
        float r50725 = r50723 * r50724;
        float r50726 = r50725 * r50724;
        float r50727 = z;
        float r50728 = t;
        float r50729 = r50727 / r50728;
        float r50730 = fabs(r50729);
        float r50731 = r50730 * r50730;
        float r50732 = r50726 + r50731;
        return r50732;
}

double f_od(double x, double y, double z, double t) {
        double r50733 = x;
        double r50734 = y;
        double r50735 = r50733 / r50734;
        double r50736 = fabs(r50735);
        double r50737 = sqrt(r50736);
        double r50738 = r50736 * r50737;
        double r50739 = r50738 * r50737;
        double r50740 = z;
        double r50741 = t;
        double r50742 = r50740 / r50741;
        double r50743 = fabs(r50742);
        double r50744 = r50743 * r50743;
        double r50745 = r50739 + r50744;
        return r50745;
}

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 r50746, r50747, r50748, r50749, r50750, r50751, r50752, r50753, r50754, r50755, r50756;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r50746);
        mpfr_init(r50747);
        mpfr_init(r50748);
        mpfr_init(r50749);
        mpfr_init(r50750);
        mpfr_init(r50751);
        mpfr_init(r50752);
        mpfr_init(r50753);
        mpfr_init(r50754);
        mpfr_init(r50755);
        mpfr_init(r50756);
}

double f_im(double x, double y, double z, double t) {
        mpfr_set_d(r50746, x, MPFR_RNDN);
        mpfr_mul(r50747, r50746, r50746, MPFR_RNDN);
        mpfr_set_d(r50748, y, MPFR_RNDN);
        mpfr_mul(r50749, r50748, r50748, MPFR_RNDN);
        mpfr_div(r50750, r50747, r50749, MPFR_RNDN);
        mpfr_set_d(r50751, z, MPFR_RNDN);
        mpfr_mul(r50752, r50751, r50751, MPFR_RNDN);
        mpfr_set_d(r50753, t, MPFR_RNDN);
        mpfr_mul(r50754, r50753, r50753, MPFR_RNDN);
        mpfr_div(r50755, r50752, r50754, MPFR_RNDN);
        mpfr_add(r50756, r50750, r50755, MPFR_RNDN);
        return mpfr_get_d(r50756, MPFR_RNDN);
}

static mpfr_t r50757, r50758, r50759, r50760, r50761, r50762, r50763, r50764, r50765, r50766, r50767, r50768, r50769;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r50757);
        mpfr_init(r50758);
        mpfr_init(r50759);
        mpfr_init(r50760);
        mpfr_init(r50761);
        mpfr_init(r50762);
        mpfr_init(r50763);
        mpfr_init(r50764);
        mpfr_init(r50765);
        mpfr_init(r50766);
        mpfr_init(r50767);
        mpfr_init(r50768);
        mpfr_init(r50769);
}

double f_fm(double x, double y, double z, double t) {
        mpfr_set_d(r50757, x, MPFR_RNDN);
        mpfr_set_d(r50758, y, MPFR_RNDN);
        mpfr_div(r50759, r50757, r50758, MPFR_RNDN);
        mpfr_abs(r50760, r50759, MPFR_RNDN);
        mpfr_sqrt(r50761, r50760, MPFR_RNDN);
        mpfr_mul(r50762, r50760, r50761, MPFR_RNDN);
        mpfr_mul(r50763, r50762, r50761, MPFR_RNDN);
        mpfr_set_d(r50764, z, MPFR_RNDN);
        mpfr_set_d(r50765, t, MPFR_RNDN);
        mpfr_div(r50766, r50764, r50765, MPFR_RNDN);
        mpfr_abs(r50767, r50766, MPFR_RNDN);
        mpfr_mul(r50768, r50767, r50767, MPFR_RNDN);
        mpfr_add(r50769, r50763, r50768, MPFR_RNDN);
        return mpfr_get_d(r50769, MPFR_RNDN);
}

static mpfr_t r50770, r50771, r50772, r50773, r50774, r50775, r50776, r50777, r50778, r50779, r50780, r50781, r50782;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r50770);
        mpfr_init(r50771);
        mpfr_init(r50772);
        mpfr_init(r50773);
        mpfr_init(r50774);
        mpfr_init(r50775);
        mpfr_init(r50776);
        mpfr_init(r50777);
        mpfr_init(r50778);
        mpfr_init(r50779);
        mpfr_init(r50780);
        mpfr_init(r50781);
        mpfr_init(r50782);
}

double f_dm(double x, double y, double z, double t) {
        mpfr_set_d(r50770, x, MPFR_RNDN);
        mpfr_set_d(r50771, y, MPFR_RNDN);
        mpfr_div(r50772, r50770, r50771, MPFR_RNDN);
        mpfr_abs(r50773, r50772, MPFR_RNDN);
        mpfr_sqrt(r50774, r50773, MPFR_RNDN);
        mpfr_mul(r50775, r50773, r50774, MPFR_RNDN);
        mpfr_mul(r50776, r50775, r50774, MPFR_RNDN);
        mpfr_set_d(r50777, z, MPFR_RNDN);
        mpfr_set_d(r50778, t, MPFR_RNDN);
        mpfr_div(r50779, r50777, r50778, MPFR_RNDN);
        mpfr_abs(r50780, r50779, MPFR_RNDN);
        mpfr_mul(r50781, r50780, r50780, MPFR_RNDN);
        mpfr_add(r50782, r50776, r50781, MPFR_RNDN);
        return mpfr_get_d(r50782, MPFR_RNDN);
}