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

char *name = "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, H";

double f_if(float x, float y, float z) {
        float r26109 = x;
        float r26110 = y;
        float r26111 = r26109 + r26110;
        float r26112 = 1.0;
        float r26113 = z;
        float r26114 = r26112 - r26113;
        float r26115 = r26111 * r26114;
        return r26115;
}

double f_id(double x, double y, double z) {
        double r26116 = x;
        double r26117 = y;
        double r26118 = r26116 + r26117;
        double r26119 = 1.0;
        double r26120 = z;
        double r26121 = r26119 - r26120;
        double r26122 = r26118 * r26121;
        return r26122;
}


double f_of(float x, float y, float z) {
        float r26123 = x;
        float r26124 = y;
        float r26125 = r26123 + r26124;
        float r26126 = 1.0;
        float r26127 = z;
        float r26128 = r26126 - r26127;
        float r26129 = r26125 * r26128;
        return r26129;
}

double f_od(double x, double y, double z) {
        double r26130 = x;
        double r26131 = y;
        double r26132 = r26130 + r26131;
        double r26133 = 1.0;
        double r26134 = z;
        double r26135 = r26133 - r26134;
        double r26136 = r26132 * r26135;
        return r26136;
}

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 r26137, r26138, r26139, r26140, r26141, r26142, r26143;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r26137);
        mpfr_init(r26138);
        mpfr_init(r26139);
        mpfr_init_set_str(r26140, "1.0", 10, MPFR_RNDN);
        mpfr_init(r26141);
        mpfr_init(r26142);
        mpfr_init(r26143);
}

double f_im(double x, double y, double z) {
        mpfr_set_d(r26137, x, MPFR_RNDN);
        mpfr_set_d(r26138, y, MPFR_RNDN);
        mpfr_add(r26139, r26137, r26138, MPFR_RNDN);
        ;
        mpfr_set_d(r26141, z, MPFR_RNDN);
        mpfr_sub(r26142, r26140, r26141, MPFR_RNDN);
        mpfr_mul(r26143, r26139, r26142, MPFR_RNDN);
        return mpfr_get_d(r26143, MPFR_RNDN);
}

static mpfr_t r26144, r26145, r26146, r26147, r26148, r26149, r26150;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r26144);
        mpfr_init(r26145);
        mpfr_init(r26146);
        mpfr_init_set_str(r26147, "1.0", 10, MPFR_RNDN);
        mpfr_init(r26148);
        mpfr_init(r26149);
        mpfr_init(r26150);
}

double f_fm(double x, double y, double z) {
        mpfr_set_d(r26144, x, MPFR_RNDN);
        mpfr_set_d(r26145, y, MPFR_RNDN);
        mpfr_add(r26146, r26144, r26145, MPFR_RNDN);
        ;
        mpfr_set_d(r26148, z, MPFR_RNDN);
        mpfr_sub(r26149, r26147, r26148, MPFR_RNDN);
        mpfr_mul(r26150, r26146, r26149, MPFR_RNDN);
        return mpfr_get_d(r26150, MPFR_RNDN);
}

static mpfr_t r26151, r26152, r26153, r26154, r26155, r26156, r26157;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r26151);
        mpfr_init(r26152);
        mpfr_init(r26153);
        mpfr_init_set_str(r26154, "1.0", 10, MPFR_RNDN);
        mpfr_init(r26155);
        mpfr_init(r26156);
        mpfr_init(r26157);
}

double f_dm(double x, double y, double z) {
        mpfr_set_d(r26151, x, MPFR_RNDN);
        mpfr_set_d(r26152, y, MPFR_RNDN);
        mpfr_add(r26153, r26151, r26152, MPFR_RNDN);
        ;
        mpfr_set_d(r26155, z, MPFR_RNDN);
        mpfr_sub(r26156, r26154, r26155, MPFR_RNDN);
        mpfr_mul(r26157, r26153, r26156, MPFR_RNDN);
        return mpfr_get_d(r26157, MPFR_RNDN);
}