#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (23)"; double f_if(float F, float B, float x) { float r21351 = x; float r21352 = 1; float r21353 = B; float r21354 = tan(r21353); float r21355 = r21352 / r21354; float r21356 = r21351 * r21355; float r21357 = -r21356; float r21358 = F; float r21359 = sin(r21353); float r21360 = r21358 / r21359; float r21361 = r21358 * r21358; float r21362 = 2; float r21363 = r21361 + r21362; float r21364 = r21362 * r21351; float r21365 = r21363 + r21364; float r21366 = r21352 / r21362; float r21367 = -r21366; float r21368 = pow(r21365, r21367); float r21369 = r21360 * r21368; float r21370 = r21357 + r21369; return r21370; } double f_id(double F, double B, double x) { double r21371 = x; double r21372 = 1; double r21373 = B; double r21374 = tan(r21373); double r21375 = r21372 / r21374; double r21376 = r21371 * r21375; double r21377 = -r21376; double r21378 = F; double r21379 = sin(r21373); double r21380 = r21378 / r21379; double r21381 = r21378 * r21378; double r21382 = 2; double r21383 = r21381 + r21382; double r21384 = r21382 * r21371; double r21385 = r21383 + r21384; double r21386 = r21372 / r21382; double r21387 = -r21386; double r21388 = pow(r21385, r21387); double r21389 = r21380 * r21388; double r21390 = r21377 + r21389; return r21390; } double f_of(float F, float B, float x) { float r21391 = x; float r21392 = 2; float r21393 = F; float r21394 = fma(r21393, r21393, r21392); float r21395 = fma(r21391, r21392, r21394); float r21396 = 1; float r21397 = r21396 / r21392; float r21398 = -r21397; float r21399 = pow(r21395, r21398); float r21400 = B; float r21401 = sin(r21400); float r21402 = r21396 / r21401; float r21403 = r21393 * r21402; float r21404 = -r21391; float r21405 = tan(r21400); float r21406 = r21404 / r21405; float r21407 = fma(r21399, r21403, r21406); return r21407; } double f_od(double F, double B, double x) { double r21408 = x; double r21409 = 2; double r21410 = F; double r21411 = fma(r21410, r21410, r21409); double r21412 = fma(r21408, r21409, r21411); double r21413 = 1; double r21414 = r21413 / r21409; double r21415 = -r21414; double r21416 = pow(r21412, r21415); double r21417 = B; double r21418 = sin(r21417); double r21419 = r21413 / r21418; double r21420 = r21410 * r21419; double r21421 = -r21408; double r21422 = tan(r21417); double r21423 = r21421 / r21422; double r21424 = fma(r21416, r21420, r21423); return r21424; } 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 r21425, r21426, r21427, r21428, r21429, r21430, r21431, r21432, r21433, r21434, r21435, r21436, r21437, r21438, r21439, r21440, r21441, r21442, r21443, r21444; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r21425); mpfr_init_set_str(r21426, "1", 10, MPFR_RNDN); mpfr_init(r21427); mpfr_init(r21428); mpfr_init(r21429); mpfr_init(r21430); mpfr_init(r21431); mpfr_init(r21432); mpfr_init(r21433); mpfr_init(r21434); mpfr_init(r21435); mpfr_init_set_str(r21436, "2", 10, MPFR_RNDN); mpfr_init(r21437); mpfr_init(r21438); mpfr_init(r21439); mpfr_init(r21440); mpfr_init(r21441); mpfr_init(r21442); mpfr_init(r21443); mpfr_init(r21444); } double f_im(double F, double B, double x) { mpfr_set_d(r21425, x, MPFR_RNDN); ; mpfr_set_d(r21427, B, MPFR_RNDN); mpfr_tan(r21428, r21427, MPFR_RNDN); mpfr_div(r21429, r21426, r21428, MPFR_RNDN); mpfr_mul(r21430, r21425, r21429, MPFR_RNDN); mpfr_neg(r21431, r21430, MPFR_RNDN); mpfr_set_d(r21432, F, MPFR_RNDN); mpfr_sin(r21433, r21427, MPFR_RNDN); mpfr_div(r21434, r21432, r21433, MPFR_RNDN); mpfr_mul(r21435, r21432, r21432, MPFR_RNDN); ; mpfr_add(r21437, r21435, r21436, MPFR_RNDN); mpfr_mul(r21438, r21436, r21425, MPFR_RNDN); mpfr_add(r21439, r21437, r21438, MPFR_RNDN); mpfr_div(r21440, r21426, r21436, MPFR_RNDN); mpfr_neg(r21441, r21440, MPFR_RNDN); mpfr_pow(r21442, r21439, r21441, MPFR_RNDN); mpfr_mul(r21443, r21434, r21442, MPFR_RNDN); mpfr_add(r21444, r21431, r21443, MPFR_RNDN); return mpfr_get_d(r21444, MPFR_RNDN); } static mpfr_t r21445, r21446, r21447, r21448, r21449, r21450, r21451, r21452, r21453, r21454, r21455, r21456, r21457, r21458, r21459, r21460, r21461; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r21445); mpfr_init_set_str(r21446, "2", 10, MPFR_RNDN); mpfr_init(r21447); mpfr_init(r21448); mpfr_init(r21449); mpfr_init_set_str(r21450, "1", 10, MPFR_RNDN); mpfr_init(r21451); mpfr_init(r21452); mpfr_init(r21453); mpfr_init(r21454); mpfr_init(r21455); mpfr_init(r21456); mpfr_init(r21457); mpfr_init(r21458); mpfr_init(r21459); mpfr_init(r21460); mpfr_init(r21461); } double f_fm(double F, double B, double x) { mpfr_set_d(r21445, x, MPFR_RNDN); ; mpfr_set_d(r21447, F, MPFR_RNDN); mpfr_fma(r21448, r21447, r21447, r21446, MPFR_RNDN); mpfr_fma(r21449, r21445, r21446, r21448, MPFR_RNDN); ; mpfr_div(r21451, r21450, r21446, MPFR_RNDN); mpfr_neg(r21452, r21451, MPFR_RNDN); mpfr_pow(r21453, r21449, r21452, MPFR_RNDN); mpfr_set_d(r21454, B, MPFR_RNDN); mpfr_sin(r21455, r21454, MPFR_RNDN); mpfr_div(r21456, r21450, r21455, MPFR_RNDN); mpfr_mul(r21457, r21447, r21456, MPFR_RNDN); mpfr_neg(r21458, r21445, MPFR_RNDN); mpfr_tan(r21459, r21454, MPFR_RNDN); mpfr_div(r21460, r21458, r21459, MPFR_RNDN); mpfr_fma(r21461, r21453, r21457, r21460, MPFR_RNDN); return mpfr_get_d(r21461, MPFR_RNDN); } static mpfr_t r21462, r21463, r21464, r21465, r21466, r21467, r21468, r21469, r21470, r21471, r21472, r21473, r21474, r21475, r21476, r21477, r21478; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r21462); mpfr_init_set_str(r21463, "2", 10, MPFR_RNDN); mpfr_init(r21464); mpfr_init(r21465); mpfr_init(r21466); mpfr_init_set_str(r21467, "1", 10, MPFR_RNDN); mpfr_init(r21468); mpfr_init(r21469); mpfr_init(r21470); mpfr_init(r21471); mpfr_init(r21472); mpfr_init(r21473); mpfr_init(r21474); mpfr_init(r21475); mpfr_init(r21476); mpfr_init(r21477); mpfr_init(r21478); } double f_dm(double F, double B, double x) { mpfr_set_d(r21462, x, MPFR_RNDN); ; mpfr_set_d(r21464, F, MPFR_RNDN); mpfr_fma(r21465, r21464, r21464, r21463, MPFR_RNDN); mpfr_fma(r21466, r21462, r21463, r21465, MPFR_RNDN); ; mpfr_div(r21468, r21467, r21463, MPFR_RNDN); mpfr_neg(r21469, r21468, MPFR_RNDN); mpfr_pow(r21470, r21466, r21469, MPFR_RNDN); mpfr_set_d(r21471, B, MPFR_RNDN); mpfr_sin(r21472, r21471, MPFR_RNDN); mpfr_div(r21473, r21467, r21472, MPFR_RNDN); mpfr_mul(r21474, r21464, r21473, MPFR_RNDN); mpfr_neg(r21475, r21462, MPFR_RNDN); mpfr_tan(r21476, r21471, MPFR_RNDN); mpfr_div(r21477, r21475, r21476, MPFR_RNDN); mpfr_fma(r21478, r21470, r21474, r21477, MPFR_RNDN); return mpfr_get_d(r21478, MPFR_RNDN); }