#include #include #include #include #include char *name = "(+ x (- (tan (+ y z)) (tan a)))"; double f_if(float x, float y, float z, float a) { float r24391 = x; float r24392 = y; float r24393 = z; float r24394 = r24392 + r24393; float r24395 = tan(r24394); float r24396 = a; float r24397 = tan(r24396); float r24398 = r24395 - r24397; float r24399 = r24391 + r24398; return r24399; } double f_id(double x, double y, double z, double a) { double r24400 = x; double r24401 = y; double r24402 = z; double r24403 = r24401 + r24402; double r24404 = tan(r24403); double r24405 = a; double r24406 = tan(r24405); double r24407 = r24404 - r24406; double r24408 = r24400 + r24407; return r24408; } double f_of(float x, float y, float z, float a) { float r24409 = x; float r24410 = y; float r24411 = tan(r24410); float r24412 = z; float r24413 = tan(r24412); float r24414 = r24411 + r24413; float r24415 = a; float r24416 = cos(r24415); float r24417 = r24414 * r24416; float r24418 = 1; float r24419 = r24411 * r24413; float r24420 = r24418 - r24419; float r24421 = sin(r24415); float r24422 = r24420 * r24421; float r24423 = r24417 - r24422; float r24424 = r24420 * r24416; float r24425 = r24423 / r24424; float r24426 = r24409 + r24425; return r24426; } double f_od(double x, double y, double z, double a) { double r24427 = x; double r24428 = y; double r24429 = tan(r24428); double r24430 = z; double r24431 = tan(r24430); double r24432 = r24429 + r24431; double r24433 = a; double r24434 = cos(r24433); double r24435 = r24432 * r24434; double r24436 = 1; double r24437 = r24429 * r24431; double r24438 = r24436 - r24437; double r24439 = sin(r24433); double r24440 = r24438 * r24439; double r24441 = r24435 - r24440; double r24442 = r24438 * r24434; double r24443 = r24441 / r24442; double r24444 = r24427 + r24443; return r24444; } 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 r24445, r24446, r24447, r24448, r24449, r24450, r24451, r24452, r24453; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r24445); mpfr_init(r24446); mpfr_init(r24447); mpfr_init(r24448); mpfr_init(r24449); mpfr_init(r24450); mpfr_init(r24451); mpfr_init(r24452); mpfr_init(r24453); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r24445, x, MPFR_RNDN); mpfr_set_d(r24446, y, MPFR_RNDN); mpfr_set_d(r24447, z, MPFR_RNDN); mpfr_add(r24448, r24446, r24447, MPFR_RNDN); mpfr_tan(r24449, r24448, MPFR_RNDN); mpfr_set_d(r24450, a, MPFR_RNDN); mpfr_tan(r24451, r24450, MPFR_RNDN); mpfr_sub(r24452, r24449, r24451, MPFR_RNDN); mpfr_add(r24453, r24445, r24452, MPFR_RNDN); return mpfr_get_d(r24453, MPFR_RNDN); } static mpfr_t r24454, r24455, r24456, r24457, r24458, r24459, r24460, r24461, r24462, r24463, r24464, r24465, r24466, r24467, r24468, r24469, r24470, r24471; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r24454); mpfr_init(r24455); mpfr_init(r24456); mpfr_init(r24457); mpfr_init(r24458); mpfr_init(r24459); mpfr_init(r24460); mpfr_init(r24461); mpfr_init(r24462); mpfr_init_set_str(r24463, "1", 10, MPFR_RNDN); mpfr_init(r24464); mpfr_init(r24465); mpfr_init(r24466); mpfr_init(r24467); mpfr_init(r24468); mpfr_init(r24469); mpfr_init(r24470); mpfr_init(r24471); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r24454, x, MPFR_RNDN); mpfr_set_d(r24455, y, MPFR_RNDN); mpfr_tan(r24456, r24455, MPFR_RNDN); mpfr_set_d(r24457, z, MPFR_RNDN); mpfr_tan(r24458, r24457, MPFR_RNDN); mpfr_add(r24459, r24456, r24458, MPFR_RNDN); mpfr_set_d(r24460, a, MPFR_RNDN); mpfr_cos(r24461, r24460, MPFR_RNDN); mpfr_mul(r24462, r24459, r24461, MPFR_RNDN); ; mpfr_mul(r24464, r24456, r24458, MPFR_RNDN); mpfr_sub(r24465, r24463, r24464, MPFR_RNDN); mpfr_sin(r24466, r24460, MPFR_RNDN); mpfr_mul(r24467, r24465, r24466, MPFR_RNDN); mpfr_sub(r24468, r24462, r24467, MPFR_RNDN); mpfr_mul(r24469, r24465, r24461, MPFR_RNDN); mpfr_div(r24470, r24468, r24469, MPFR_RNDN); mpfr_add(r24471, r24454, r24470, MPFR_RNDN); return mpfr_get_d(r24471, MPFR_RNDN); } static mpfr_t r24472, r24473, r24474, r24475, r24476, r24477, r24478, r24479, r24480, r24481, r24482, r24483, r24484, r24485, r24486, r24487, r24488, r24489; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r24472); mpfr_init(r24473); mpfr_init(r24474); mpfr_init(r24475); mpfr_init(r24476); mpfr_init(r24477); mpfr_init(r24478); mpfr_init(r24479); mpfr_init(r24480); mpfr_init_set_str(r24481, "1", 10, MPFR_RNDN); mpfr_init(r24482); mpfr_init(r24483); mpfr_init(r24484); mpfr_init(r24485); mpfr_init(r24486); mpfr_init(r24487); mpfr_init(r24488); mpfr_init(r24489); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r24472, x, MPFR_RNDN); mpfr_set_d(r24473, y, MPFR_RNDN); mpfr_tan(r24474, r24473, MPFR_RNDN); mpfr_set_d(r24475, z, MPFR_RNDN); mpfr_tan(r24476, r24475, MPFR_RNDN); mpfr_add(r24477, r24474, r24476, MPFR_RNDN); mpfr_set_d(r24478, a, MPFR_RNDN); mpfr_cos(r24479, r24478, MPFR_RNDN); mpfr_mul(r24480, r24477, r24479, MPFR_RNDN); ; mpfr_mul(r24482, r24474, r24476, MPFR_RNDN); mpfr_sub(r24483, r24481, r24482, MPFR_RNDN); mpfr_sin(r24484, r24478, MPFR_RNDN); mpfr_mul(r24485, r24483, r24484, MPFR_RNDN); mpfr_sub(r24486, r24480, r24485, MPFR_RNDN); mpfr_mul(r24487, r24483, r24479, MPFR_RNDN); mpfr_div(r24488, r24486, r24487, MPFR_RNDN); mpfr_add(r24489, r24472, r24488, MPFR_RNDN); return mpfr_get_d(r24489, MPFR_RNDN); }