#include #include #include #include #include char *name = "Henrywood and Agarwal, Equation (9a)"; double f_if(float w0, float M, float D, float h, float l, float d) { float r26384 = w0; float r26385 = 1; float r26386 = M; float r26387 = D; float r26388 = r26386 * r26387; float r26389 = 2; float r26390 = d; float r26391 = r26389 * r26390; float r26392 = r26388 / r26391; float r26393 = pow(r26392, r26389); float r26394 = h; float r26395 = l; float r26396 = r26394 / r26395; float r26397 = r26393 * r26396; float r26398 = r26385 - r26397; float r26399 = sqrt(r26398); float r26400 = r26384 * r26399; return r26400; } double f_id(double w0, double M, double D, double h, double l, double d) { double r26401 = w0; double r26402 = 1; double r26403 = M; double r26404 = D; double r26405 = r26403 * r26404; double r26406 = 2; double r26407 = d; double r26408 = r26406 * r26407; double r26409 = r26405 / r26408; double r26410 = pow(r26409, r26406); double r26411 = h; double r26412 = l; double r26413 = r26411 / r26412; double r26414 = r26410 * r26413; double r26415 = r26402 - r26414; double r26416 = sqrt(r26415); double r26417 = r26401 * r26416; return r26417; } double f_of(float w0, float M, float D, float h, float l, float d) { float r26418 = h; float r26419 = l; float r26420 = r26418 / r26419; float r26421 = -1.7775498841371257e+308; bool r26422 = r26420 <= r26421; float r26423 = -1.958917333755248e-230; bool r26424 = r26420 <= r26423; float r26425 = !r26424; bool r26426 = r26422 || r26425; float r26427 = 1; float r26428 = D; float r26429 = M; float r26430 = r26428 * r26429; float r26431 = d; float r26432 = 2; float r26433 = r26431 * r26432; float r26434 = r26430 / r26433; float r26435 = r26434 * r26418; float r26436 = r26434 * r26435; float r26437 = r26436 / r26419; float r26438 = r26427 - r26437; float r26439 = sqrt(r26438); float r26440 = w0; float r26441 = r26439 * r26440; float r26442 = r26434 * r26420; float r26443 = r26442 * r26434; float r26444 = r26427 - r26443; float r26445 = sqrt(r26444); float r26446 = r26440 * r26445; float r26447 = r26426 ? r26441 : r26446; return r26447; } double f_od(double w0, double M, double D, double h, double l, double d) { double r26448 = h; double r26449 = l; double r26450 = r26448 / r26449; double r26451 = -1.7775498841371257e+308; bool r26452 = r26450 <= r26451; double r26453 = -1.958917333755248e-230; bool r26454 = r26450 <= r26453; double r26455 = !r26454; bool r26456 = r26452 || r26455; double r26457 = 1; double r26458 = D; double r26459 = M; double r26460 = r26458 * r26459; double r26461 = d; double r26462 = 2; double r26463 = r26461 * r26462; double r26464 = r26460 / r26463; double r26465 = r26464 * r26448; double r26466 = r26464 * r26465; double r26467 = r26466 / r26449; double r26468 = r26457 - r26467; double r26469 = sqrt(r26468); double r26470 = w0; double r26471 = r26469 * r26470; double r26472 = r26464 * r26450; double r26473 = r26472 * r26464; double r26474 = r26457 - r26473; double r26475 = sqrt(r26474); double r26476 = r26470 * r26475; double r26477 = r26456 ? r26471 : r26476; return r26477; } 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 r26478, r26479, r26480, r26481, r26482, r26483, r26484, r26485, r26486, r26487, r26488, r26489, r26490, r26491, r26492, r26493, r26494; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r26478); mpfr_init_set_str(r26479, "1", 10, MPFR_RNDN); mpfr_init(r26480); mpfr_init(r26481); mpfr_init(r26482); mpfr_init_set_str(r26483, "2", 10, MPFR_RNDN); mpfr_init(r26484); mpfr_init(r26485); mpfr_init(r26486); mpfr_init(r26487); mpfr_init(r26488); mpfr_init(r26489); mpfr_init(r26490); mpfr_init(r26491); mpfr_init(r26492); mpfr_init(r26493); mpfr_init(r26494); } double f_im(double w0, double M, double D, double h, double l, double d) { mpfr_set_d(r26478, w0, MPFR_RNDN); ; mpfr_set_d(r26480, M, MPFR_RNDN); mpfr_set_d(r26481, D, MPFR_RNDN); mpfr_mul(r26482, r26480, r26481, MPFR_RNDN); ; mpfr_set_d(r26484, d, MPFR_RNDN); mpfr_mul(r26485, r26483, r26484, MPFR_RNDN); mpfr_div(r26486, r26482, r26485, MPFR_RNDN); mpfr_pow(r26487, r26486, r26483, MPFR_RNDN); mpfr_set_d(r26488, h, MPFR_RNDN); mpfr_set_d(r26489, l, MPFR_RNDN); mpfr_div(r26490, r26488, r26489, MPFR_RNDN); mpfr_mul(r26491, r26487, r26490, MPFR_RNDN); mpfr_sub(r26492, r26479, r26491, MPFR_RNDN); mpfr_sqrt(r26493, r26492, MPFR_RNDN); mpfr_mul(r26494, r26478, r26493, MPFR_RNDN); return mpfr_get_d(r26494, MPFR_RNDN); } static mpfr_t r26495, r26496, r26497, r26498, r26499, r26500, r26501, r26502, r26503, r26504, r26505, r26506, r26507, r26508, r26509, r26510, r26511, r26512, r26513, r26514, r26515, r26516, r26517, r26518, r26519, r26520, r26521, r26522, r26523, r26524; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r26495); mpfr_init(r26496); mpfr_init(r26497); mpfr_init_set_str(r26498, "-1.7775498841371257e+308", 10, MPFR_RNDN); mpfr_init(r26499); mpfr_init_set_str(r26500, "-1.958917333755248e-230", 10, MPFR_RNDN); mpfr_init(r26501); mpfr_init(r26502); mpfr_init(r26503); mpfr_init_set_str(r26504, "1", 10, MPFR_RNDN); mpfr_init(r26505); mpfr_init(r26506); mpfr_init(r26507); mpfr_init(r26508); mpfr_init_set_str(r26509, "2", 10, MPFR_RNDN); mpfr_init(r26510); mpfr_init(r26511); 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); } double f_fm(double w0, double M, double D, double h, double l, double d) { mpfr_set_d(r26495, h, MPFR_RNDN); mpfr_set_d(r26496, l, MPFR_RNDN); mpfr_div(r26497, r26495, r26496, MPFR_RNDN); ; mpfr_set_si(r26499, mpfr_cmp(r26497, r26498) <= 0, MPFR_RNDN); ; mpfr_set_si(r26501, mpfr_cmp(r26497, r26500) <= 0, MPFR_RNDN); mpfr_set_si(r26502, !mpfr_get_si(r26501, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r26503, mpfr_get_si(r26499, MPFR_RNDN) || mpfr_get_si(r26502, MPFR_RNDN), MPFR_RNDN); ; mpfr_set_d(r26505, D, MPFR_RNDN); mpfr_set_d(r26506, M, MPFR_RNDN); mpfr_mul(r26507, r26505, r26506, MPFR_RNDN); mpfr_set_d(r26508, d, MPFR_RNDN); ; mpfr_mul(r26510, r26508, r26509, MPFR_RNDN); mpfr_div(r26511, r26507, r26510, MPFR_RNDN); mpfr_mul(r26512, r26511, r26495, MPFR_RNDN); mpfr_mul(r26513, r26511, r26512, MPFR_RNDN); mpfr_div(r26514, r26513, r26496, MPFR_RNDN); mpfr_sub(r26515, r26504, r26514, MPFR_RNDN); mpfr_sqrt(r26516, r26515, MPFR_RNDN); mpfr_set_d(r26517, w0, MPFR_RNDN); mpfr_mul(r26518, r26516, r26517, MPFR_RNDN); mpfr_mul(r26519, r26511, r26497, MPFR_RNDN); mpfr_mul(r26520, r26519, r26511, MPFR_RNDN); mpfr_sub(r26521, r26504, r26520, MPFR_RNDN); mpfr_sqrt(r26522, r26521, MPFR_RNDN); mpfr_mul(r26523, r26517, r26522, MPFR_RNDN); if (mpfr_get_si(r26503, MPFR_RNDN)) { mpfr_set(r26524, r26518, MPFR_RNDN); } else { mpfr_set(r26524, r26523, MPFR_RNDN); }; return mpfr_get_d(r26524, MPFR_RNDN); } static mpfr_t r26525, r26526, 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, r26552, r26553, r26554; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r26525); mpfr_init(r26526); mpfr_init(r26527); mpfr_init_set_str(r26528, "-1.7775498841371257e+308", 10, MPFR_RNDN); mpfr_init(r26529); mpfr_init_set_str(r26530, "-1.958917333755248e-230", 10, MPFR_RNDN); mpfr_init(r26531); mpfr_init(r26532); mpfr_init(r26533); mpfr_init_set_str(r26534, "1", 10, MPFR_RNDN); mpfr_init(r26535); mpfr_init(r26536); mpfr_init(r26537); mpfr_init(r26538); mpfr_init_set_str(r26539, "2", 10, MPFR_RNDN); 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); mpfr_init(r26552); mpfr_init(r26553); mpfr_init(r26554); } double f_dm(double w0, double M, double D, double h, double l, double d) { mpfr_set_d(r26525, h, MPFR_RNDN); mpfr_set_d(r26526, l, MPFR_RNDN); mpfr_div(r26527, r26525, r26526, MPFR_RNDN); ; mpfr_set_si(r26529, mpfr_cmp(r26527, r26528) <= 0, MPFR_RNDN); ; mpfr_set_si(r26531, mpfr_cmp(r26527, r26530) <= 0, MPFR_RNDN); mpfr_set_si(r26532, !mpfr_get_si(r26531, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r26533, mpfr_get_si(r26529, MPFR_RNDN) || mpfr_get_si(r26532, MPFR_RNDN), MPFR_RNDN); ; mpfr_set_d(r26535, D, MPFR_RNDN); mpfr_set_d(r26536, M, MPFR_RNDN); mpfr_mul(r26537, r26535, r26536, MPFR_RNDN); mpfr_set_d(r26538, d, MPFR_RNDN); ; mpfr_mul(r26540, r26538, r26539, MPFR_RNDN); mpfr_div(r26541, r26537, r26540, MPFR_RNDN); mpfr_mul(r26542, r26541, r26525, MPFR_RNDN); mpfr_mul(r26543, r26541, r26542, MPFR_RNDN); mpfr_div(r26544, r26543, r26526, MPFR_RNDN); mpfr_sub(r26545, r26534, r26544, MPFR_RNDN); mpfr_sqrt(r26546, r26545, MPFR_RNDN); mpfr_set_d(r26547, w0, MPFR_RNDN); mpfr_mul(r26548, r26546, r26547, MPFR_RNDN); mpfr_mul(r26549, r26541, r26527, MPFR_RNDN); mpfr_mul(r26550, r26549, r26541, MPFR_RNDN); mpfr_sub(r26551, r26534, r26550, MPFR_RNDN); mpfr_sqrt(r26552, r26551, MPFR_RNDN); mpfr_mul(r26553, r26547, r26552, MPFR_RNDN); if (mpfr_get_si(r26533, MPFR_RNDN)) { mpfr_set(r26554, r26548, MPFR_RNDN); } else { mpfr_set(r26554, r26553, MPFR_RNDN); }; return mpfr_get_d(r26554, MPFR_RNDN); }