#include #include #include #include #include char *name = "Henrywood and Agarwal, Equation (12)"; double f_if(float d, float h, float l, float M, float D) { float r26226 = d; float r26227 = h; float r26228 = r26226 / r26227; float r26229 = 1; float r26230 = 2; float r26231 = r26229 / r26230; float r26232 = pow(r26228, r26231); float r26233 = l; float r26234 = r26226 / r26233; float r26235 = pow(r26234, r26231); float r26236 = r26232 * r26235; float r26237 = M; float r26238 = D; float r26239 = r26237 * r26238; float r26240 = r26230 * r26226; float r26241 = r26239 / r26240; float r26242 = pow(r26241, r26230); float r26243 = r26231 * r26242; float r26244 = r26227 / r26233; float r26245 = r26243 * r26244; float r26246 = r26229 - r26245; float r26247 = r26236 * r26246; return r26247; } double f_id(double d, double h, double l, double M, double D) { double r26248 = d; double r26249 = h; double r26250 = r26248 / r26249; double r26251 = 1; double r26252 = 2; double r26253 = r26251 / r26252; double r26254 = pow(r26250, r26253); double r26255 = l; double r26256 = r26248 / r26255; double r26257 = pow(r26256, r26253); double r26258 = r26254 * r26257; double r26259 = M; double r26260 = D; double r26261 = r26259 * r26260; double r26262 = r26252 * r26248; double r26263 = r26261 / r26262; double r26264 = pow(r26263, r26252); double r26265 = r26253 * r26264; double r26266 = r26249 / r26255; double r26267 = r26265 * r26266; double r26268 = r26251 - r26267; double r26269 = r26258 * r26268; return r26269; } double f_of(float d, float h, float l, float M, float D) { float r26270 = h; float r26271 = -6.21156791415509e-93; bool r26272 = r26270 <= r26271; float r26273 = 1; float r26274 = l; float r26275 = r26273 / r26274; float r26276 = 2; float r26277 = r26273 / r26276; float r26278 = M; float r26279 = D; float r26280 = r26278 * r26279; float r26281 = d; float r26282 = r26281 * r26276; float r26283 = r26280 / r26282; float r26284 = pow(r26283, r26276); float r26285 = r26277 * r26284; float r26286 = r26270 * r26285; float r26287 = r26275 * r26286; float r26288 = r26273 - r26287; float r26289 = r26281 / r26270; float r26290 = pow(r26289, r26277); float r26291 = 1/2; float r26292 = -1; float r26293 = r26292 / r26274; float r26294 = log(r26293); float r26295 = r26292 / r26281; float r26296 = log(r26295); float r26297 = r26294 - r26296; float r26298 = r26291 * r26297; float r26299 = exp(r26298); float r26300 = r26290 * r26299; float r26301 = r26288 * r26300; float r26302 = -5.7840309184133e-311; bool r26303 = r26270 <= r26302; float r26304 = r26292 / r26270; float r26305 = log(r26304); float r26306 = r26305 - r26296; float r26307 = r26306 * r26291; float r26308 = exp(r26307); float r26309 = r26281 / r26274; float r26310 = pow(r26309, r26277); float r26311 = r26308 * r26310; float r26312 = r26311 * r26288; float r26313 = 1.2379471543476946e-36; bool r26314 = r26270 <= r26313; float r26315 = 9.612895863884516e+116; bool r26316 = r26270 <= r26315; float r26317 = !r26316; bool r26318 = r26314 || r26317; float r26319 = r26273 / r26270; float r26320 = pow(r26319, r26277); float r26321 = pow(r26281, r26277); float r26322 = r26320 * r26321; float r26323 = r26310 * r26322; float r26324 = r26323 * r26288; float r26325 = pow(r26275, r26277); float r26326 = r26321 * r26325; float r26327 = r26326 * r26290; float r26328 = r26288 * r26327; float r26329 = r26318 ? r26324 : r26328; float r26330 = r26303 ? r26312 : r26329; float r26331 = r26272 ? r26301 : r26330; return r26331; } double f_od(double d, double h, double l, double M, double D) { double r26332 = h; double r26333 = -6.21156791415509e-93; bool r26334 = r26332 <= r26333; double r26335 = 1; double r26336 = l; double r26337 = r26335 / r26336; double r26338 = 2; double r26339 = r26335 / r26338; double r26340 = M; double r26341 = D; double r26342 = r26340 * r26341; double r26343 = d; double r26344 = r26343 * r26338; double r26345 = r26342 / r26344; double r26346 = pow(r26345, r26338); double r26347 = r26339 * r26346; double r26348 = r26332 * r26347; double r26349 = r26337 * r26348; double r26350 = r26335 - r26349; double r26351 = r26343 / r26332; double r26352 = pow(r26351, r26339); double r26353 = 1/2; double r26354 = -1; double r26355 = r26354 / r26336; double r26356 = log(r26355); double r26357 = r26354 / r26343; double r26358 = log(r26357); double r26359 = r26356 - r26358; double r26360 = r26353 * r26359; double r26361 = exp(r26360); double r26362 = r26352 * r26361; double r26363 = r26350 * r26362; double r26364 = -5.7840309184133e-311; bool r26365 = r26332 <= r26364; double r26366 = r26354 / r26332; double r26367 = log(r26366); double r26368 = r26367 - r26358; double r26369 = r26368 * r26353; double r26370 = exp(r26369); double r26371 = r26343 / r26336; double r26372 = pow(r26371, r26339); double r26373 = r26370 * r26372; double r26374 = r26373 * r26350; double r26375 = 1.2379471543476946e-36; bool r26376 = r26332 <= r26375; double r26377 = 9.612895863884516e+116; bool r26378 = r26332 <= r26377; double r26379 = !r26378; bool r26380 = r26376 || r26379; double r26381 = r26335 / r26332; double r26382 = pow(r26381, r26339); double r26383 = pow(r26343, r26339); double r26384 = r26382 * r26383; double r26385 = r26372 * r26384; double r26386 = r26385 * r26350; double r26387 = pow(r26337, r26339); double r26388 = r26383 * r26387; double r26389 = r26388 * r26352; double r26390 = r26350 * r26389; double r26391 = r26380 ? r26386 : r26390; double r26392 = r26365 ? r26374 : r26391; double r26393 = r26334 ? r26363 : r26392; return r26393; } 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 r26394, r26395, r26396, r26397, r26398, r26399, r26400, r26401, r26402, r26403, r26404, r26405, r26406, r26407, r26408, r26409, r26410, r26411, r26412, r26413, r26414, r26415; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r26394); mpfr_init(r26395); mpfr_init(r26396); mpfr_init_set_str(r26397, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26398, "2", 10, MPFR_RNDN); mpfr_init(r26399); mpfr_init(r26400); mpfr_init(r26401); mpfr_init(r26402); mpfr_init(r26403); mpfr_init(r26404); mpfr_init(r26405); mpfr_init(r26406); mpfr_init(r26407); mpfr_init(r26408); mpfr_init(r26409); mpfr_init(r26410); mpfr_init(r26411); mpfr_init(r26412); mpfr_init(r26413); mpfr_init(r26414); mpfr_init(r26415); } double f_im(double d, double h, double l, double M, double D) { mpfr_set_d(r26394, d, MPFR_RNDN); mpfr_set_d(r26395, h, MPFR_RNDN); mpfr_div(r26396, r26394, r26395, MPFR_RNDN); ; ; mpfr_div(r26399, r26397, r26398, MPFR_RNDN); mpfr_pow(r26400, r26396, r26399, MPFR_RNDN); mpfr_set_d(r26401, l, MPFR_RNDN); mpfr_div(r26402, r26394, r26401, MPFR_RNDN); mpfr_pow(r26403, r26402, r26399, MPFR_RNDN); mpfr_mul(r26404, r26400, r26403, MPFR_RNDN); mpfr_set_d(r26405, M, MPFR_RNDN); mpfr_set_d(r26406, D, MPFR_RNDN); mpfr_mul(r26407, r26405, r26406, MPFR_RNDN); mpfr_mul(r26408, r26398, r26394, MPFR_RNDN); mpfr_div(r26409, r26407, r26408, MPFR_RNDN); mpfr_pow(r26410, r26409, r26398, MPFR_RNDN); mpfr_mul(r26411, r26399, r26410, MPFR_RNDN); mpfr_div(r26412, r26395, r26401, MPFR_RNDN); mpfr_mul(r26413, r26411, r26412, MPFR_RNDN); mpfr_sub(r26414, r26397, r26413, MPFR_RNDN); mpfr_mul(r26415, r26404, r26414, MPFR_RNDN); return mpfr_get_d(r26415, MPFR_RNDN); } static mpfr_t r26416, r26417, r26418, r26419, r26420, r26421, r26422, r26423, r26424, r26425, r26426, r26427, r26428, r26429, r26430, r26431, r26432, r26433, r26434, r26435, r26436, r26437, r26438, r26439, r26440, r26441, r26442, r26443, r26444, r26445, r26446, r26447, r26448, r26449, r26450, r26451, r26452, r26453, r26454, r26455, r26456, r26457, r26458, r26459, r26460, r26461, r26462, r26463, r26464, r26465, r26466, r26467, r26468, r26469, r26470, r26471, r26472, r26473, r26474, r26475, r26476, r26477; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r26416); mpfr_init_set_str(r26417, "-6.21156791415509e-93", 10, MPFR_RNDN); mpfr_init(r26418); mpfr_init_set_str(r26419, "1", 10, MPFR_RNDN); mpfr_init(r26420); mpfr_init(r26421); mpfr_init_set_str(r26422, "2", 10, MPFR_RNDN); mpfr_init(r26423); mpfr_init(r26424); mpfr_init(r26425); mpfr_init(r26426); mpfr_init(r26427); mpfr_init(r26428); mpfr_init(r26429); mpfr_init(r26430); mpfr_init(r26431); mpfr_init(r26432); mpfr_init(r26433); mpfr_init(r26434); mpfr_init(r26435); mpfr_init(r26436); mpfr_init_set_str(r26437, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r26438, "-1", 10, MPFR_RNDN); mpfr_init(r26439); mpfr_init(r26440); mpfr_init(r26441); mpfr_init(r26442); mpfr_init(r26443); mpfr_init(r26444); mpfr_init(r26445); mpfr_init(r26446); mpfr_init(r26447); mpfr_init_set_str(r26448, "-5.7840309184133e-311", 10, MPFR_RNDN); mpfr_init(r26449); mpfr_init(r26450); mpfr_init(r26451); mpfr_init(r26452); mpfr_init(r26453); mpfr_init(r26454); mpfr_init(r26455); mpfr_init(r26456); mpfr_init(r26457); mpfr_init(r26458); mpfr_init_set_str(r26459, "1.2379471543476946e-36", 10, MPFR_RNDN); mpfr_init(r26460); mpfr_init_set_str(r26461, "9.612895863884516e+116", 10, MPFR_RNDN); mpfr_init(r26462); mpfr_init(r26463); mpfr_init(r26464); mpfr_init(r26465); mpfr_init(r26466); mpfr_init(r26467); mpfr_init(r26468); mpfr_init(r26469); mpfr_init(r26470); mpfr_init(r26471); mpfr_init(r26472); mpfr_init(r26473); mpfr_init(r26474); mpfr_init(r26475); mpfr_init(r26476); mpfr_init(r26477); } double f_fm(double d, double h, double l, double M, double D) { mpfr_set_d(r26416, h, MPFR_RNDN); ; mpfr_set_si(r26418, mpfr_cmp(r26416, r26417) <= 0, MPFR_RNDN); ; mpfr_set_d(r26420, l, MPFR_RNDN); mpfr_div(r26421, r26419, r26420, MPFR_RNDN); ; mpfr_div(r26423, r26419, r26422, MPFR_RNDN); mpfr_set_d(r26424, M, MPFR_RNDN); mpfr_set_d(r26425, D, MPFR_RNDN); mpfr_mul(r26426, r26424, r26425, MPFR_RNDN); mpfr_set_d(r26427, d, MPFR_RNDN); mpfr_mul(r26428, r26427, r26422, MPFR_RNDN); mpfr_div(r26429, r26426, r26428, MPFR_RNDN); mpfr_pow(r26430, r26429, r26422, MPFR_RNDN); mpfr_mul(r26431, r26423, r26430, MPFR_RNDN); mpfr_mul(r26432, r26416, r26431, MPFR_RNDN); mpfr_mul(r26433, r26421, r26432, MPFR_RNDN); mpfr_sub(r26434, r26419, r26433, MPFR_RNDN); mpfr_div(r26435, r26427, r26416, MPFR_RNDN); mpfr_pow(r26436, r26435, r26423, MPFR_RNDN); ; ; mpfr_div(r26439, r26438, r26420, MPFR_RNDN); mpfr_log(r26440, r26439, MPFR_RNDN); mpfr_div(r26441, r26438, r26427, MPFR_RNDN); mpfr_log(r26442, r26441, MPFR_RNDN); mpfr_sub(r26443, r26440, r26442, MPFR_RNDN); mpfr_mul(r26444, r26437, r26443, MPFR_RNDN); mpfr_exp(r26445, r26444, MPFR_RNDN); mpfr_mul(r26446, r26436, r26445, MPFR_RNDN); mpfr_mul(r26447, r26434, r26446, MPFR_RNDN); ; mpfr_set_si(r26449, mpfr_cmp(r26416, r26448) <= 0, MPFR_RNDN); mpfr_div(r26450, r26438, r26416, MPFR_RNDN); mpfr_log(r26451, r26450, MPFR_RNDN); mpfr_sub(r26452, r26451, r26442, MPFR_RNDN); mpfr_mul(r26453, r26452, r26437, MPFR_RNDN); mpfr_exp(r26454, r26453, MPFR_RNDN); mpfr_div(r26455, r26427, r26420, MPFR_RNDN); mpfr_pow(r26456, r26455, r26423, MPFR_RNDN); mpfr_mul(r26457, r26454, r26456, MPFR_RNDN); mpfr_mul(r26458, r26457, r26434, MPFR_RNDN); ; mpfr_set_si(r26460, mpfr_cmp(r26416, r26459) <= 0, MPFR_RNDN); ; mpfr_set_si(r26462, mpfr_cmp(r26416, r26461) <= 0, MPFR_RNDN); mpfr_set_si(r26463, !mpfr_get_si(r26462, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r26464, mpfr_get_si(r26460, MPFR_RNDN) || mpfr_get_si(r26463, MPFR_RNDN), MPFR_RNDN); mpfr_div(r26465, r26419, r26416, MPFR_RNDN); mpfr_pow(r26466, r26465, r26423, MPFR_RNDN); mpfr_pow(r26467, r26427, r26423, MPFR_RNDN); mpfr_mul(r26468, r26466, r26467, MPFR_RNDN); mpfr_mul(r26469, r26456, r26468, MPFR_RNDN); mpfr_mul(r26470, r26469, r26434, MPFR_RNDN); mpfr_pow(r26471, r26421, r26423, MPFR_RNDN); mpfr_mul(r26472, r26467, r26471, MPFR_RNDN); mpfr_mul(r26473, r26472, r26436, MPFR_RNDN); mpfr_mul(r26474, r26434, r26473, MPFR_RNDN); if (mpfr_get_si(r26464, MPFR_RNDN)) { mpfr_set(r26475, r26470, MPFR_RNDN); } else { mpfr_set(r26475, r26474, MPFR_RNDN); }; if (mpfr_get_si(r26449, MPFR_RNDN)) { mpfr_set(r26476, r26458, MPFR_RNDN); } else { mpfr_set(r26476, r26475, MPFR_RNDN); }; if (mpfr_get_si(r26418, MPFR_RNDN)) { mpfr_set(r26477, r26447, MPFR_RNDN); } else { mpfr_set(r26477, r26476, MPFR_RNDN); }; return mpfr_get_d(r26477, MPFR_RNDN); } static mpfr_t r26478, r26479, r26480, r26481, r26482, r26483, r26484, r26485, r26486, r26487, r26488, r26489, r26490, r26491, r26492, r26493, r26494, 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, r26525, r26526, r26527, r26528, r26529, r26530, r26531, r26532, r26533, r26534, r26535, r26536, r26537, r26538, r26539; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r26478); mpfr_init_set_str(r26479, "-6.21156791415509e-93", 10, MPFR_RNDN); mpfr_init(r26480); mpfr_init_set_str(r26481, "1", 10, MPFR_RNDN); mpfr_init(r26482); mpfr_init(r26483); mpfr_init_set_str(r26484, "2", 10, MPFR_RNDN); 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); mpfr_init(r26495); mpfr_init(r26496); mpfr_init(r26497); mpfr_init(r26498); mpfr_init_set_str(r26499, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r26500, "-1", 10, MPFR_RNDN); mpfr_init(r26501); mpfr_init(r26502); mpfr_init(r26503); mpfr_init(r26504); mpfr_init(r26505); mpfr_init(r26506); mpfr_init(r26507); mpfr_init(r26508); mpfr_init(r26509); mpfr_init_set_str(r26510, "-5.7840309184133e-311", 10, MPFR_RNDN); 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_set_str(r26521, "1.2379471543476946e-36", 10, MPFR_RNDN); mpfr_init(r26522); mpfr_init_set_str(r26523, "9.612895863884516e+116", 10, MPFR_RNDN); mpfr_init(r26524); mpfr_init(r26525); mpfr_init(r26526); mpfr_init(r26527); mpfr_init(r26528); mpfr_init(r26529); mpfr_init(r26530); mpfr_init(r26531); mpfr_init(r26532); mpfr_init(r26533); mpfr_init(r26534); mpfr_init(r26535); mpfr_init(r26536); mpfr_init(r26537); mpfr_init(r26538); mpfr_init(r26539); } double f_dm(double d, double h, double l, double M, double D) { mpfr_set_d(r26478, h, MPFR_RNDN); ; mpfr_set_si(r26480, mpfr_cmp(r26478, r26479) <= 0, MPFR_RNDN); ; mpfr_set_d(r26482, l, MPFR_RNDN); mpfr_div(r26483, r26481, r26482, MPFR_RNDN); ; mpfr_div(r26485, r26481, r26484, MPFR_RNDN); mpfr_set_d(r26486, M, MPFR_RNDN); mpfr_set_d(r26487, D, MPFR_RNDN); mpfr_mul(r26488, r26486, r26487, MPFR_RNDN); mpfr_set_d(r26489, d, MPFR_RNDN); mpfr_mul(r26490, r26489, r26484, MPFR_RNDN); mpfr_div(r26491, r26488, r26490, MPFR_RNDN); mpfr_pow(r26492, r26491, r26484, MPFR_RNDN); mpfr_mul(r26493, r26485, r26492, MPFR_RNDN); mpfr_mul(r26494, r26478, r26493, MPFR_RNDN); mpfr_mul(r26495, r26483, r26494, MPFR_RNDN); mpfr_sub(r26496, r26481, r26495, MPFR_RNDN); mpfr_div(r26497, r26489, r26478, MPFR_RNDN); mpfr_pow(r26498, r26497, r26485, MPFR_RNDN); ; ; mpfr_div(r26501, r26500, r26482, MPFR_RNDN); mpfr_log(r26502, r26501, MPFR_RNDN); mpfr_div(r26503, r26500, r26489, MPFR_RNDN); mpfr_log(r26504, r26503, MPFR_RNDN); mpfr_sub(r26505, r26502, r26504, MPFR_RNDN); mpfr_mul(r26506, r26499, r26505, MPFR_RNDN); mpfr_exp(r26507, r26506, MPFR_RNDN); mpfr_mul(r26508, r26498, r26507, MPFR_RNDN); mpfr_mul(r26509, r26496, r26508, MPFR_RNDN); ; mpfr_set_si(r26511, mpfr_cmp(r26478, r26510) <= 0, MPFR_RNDN); mpfr_div(r26512, r26500, r26478, MPFR_RNDN); mpfr_log(r26513, r26512, MPFR_RNDN); mpfr_sub(r26514, r26513, r26504, MPFR_RNDN); mpfr_mul(r26515, r26514, r26499, MPFR_RNDN); mpfr_exp(r26516, r26515, MPFR_RNDN); mpfr_div(r26517, r26489, r26482, MPFR_RNDN); mpfr_pow(r26518, r26517, r26485, MPFR_RNDN); mpfr_mul(r26519, r26516, r26518, MPFR_RNDN); mpfr_mul(r26520, r26519, r26496, MPFR_RNDN); ; mpfr_set_si(r26522, mpfr_cmp(r26478, r26521) <= 0, MPFR_RNDN); ; mpfr_set_si(r26524, mpfr_cmp(r26478, r26523) <= 0, MPFR_RNDN); mpfr_set_si(r26525, !mpfr_get_si(r26524, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r26526, mpfr_get_si(r26522, MPFR_RNDN) || mpfr_get_si(r26525, MPFR_RNDN), MPFR_RNDN); mpfr_div(r26527, r26481, r26478, MPFR_RNDN); mpfr_pow(r26528, r26527, r26485, MPFR_RNDN); mpfr_pow(r26529, r26489, r26485, MPFR_RNDN); mpfr_mul(r26530, r26528, r26529, MPFR_RNDN); mpfr_mul(r26531, r26518, r26530, MPFR_RNDN); mpfr_mul(r26532, r26531, r26496, MPFR_RNDN); mpfr_pow(r26533, r26483, r26485, MPFR_RNDN); mpfr_mul(r26534, r26529, r26533, MPFR_RNDN); mpfr_mul(r26535, r26534, r26498, MPFR_RNDN); mpfr_mul(r26536, r26496, r26535, MPFR_RNDN); if (mpfr_get_si(r26526, MPFR_RNDN)) { mpfr_set(r26537, r26532, MPFR_RNDN); } else { mpfr_set(r26537, r26536, MPFR_RNDN); }; if (mpfr_get_si(r26511, MPFR_RNDN)) { mpfr_set(r26538, r26520, MPFR_RNDN); } else { mpfr_set(r26538, r26537, MPFR_RNDN); }; if (mpfr_get_si(r26480, MPFR_RNDN)) { mpfr_set(r26539, r26509, MPFR_RNDN); } else { mpfr_set(r26539, r26538, MPFR_RNDN); }; return mpfr_get_d(r26539, MPFR_RNDN); }