#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 r10233 = d; float r10234 = h; float r10235 = r10233 / r10234; float r10236 = 1.0f; float r10237 = 2.0f; float r10238 = r10236 / r10237; float r10239 = pow(r10235, r10238); float r10240 = l; float r10241 = r10233 / r10240; float r10242 = pow(r10241, r10238); float r10243 = r10239 * r10242; float r10244 = M; float r10245 = D; float r10246 = r10244 * r10245; float r10247 = r10237 * r10233; float r10248 = r10246 / r10247; float r10249 = pow(r10248, r10237); float r10250 = r10238 * r10249; float r10251 = r10234 / r10240; float r10252 = r10250 * r10251; float r10253 = r10236 - r10252; float r10254 = r10243 * r10253; return r10254; } double f_id(double d, double h, double l, double M, double D) { double r10255 = d; double r10256 = h; double r10257 = r10255 / r10256; double r10258 = 1.0; double r10259 = 2.0; double r10260 = r10258 / r10259; double r10261 = pow(r10257, r10260); double r10262 = l; double r10263 = r10255 / r10262; double r10264 = pow(r10263, r10260); double r10265 = r10261 * r10264; double r10266 = M; double r10267 = D; double r10268 = r10266 * r10267; double r10269 = r10259 * r10255; double r10270 = r10268 / r10269; double r10271 = pow(r10270, r10259); double r10272 = r10260 * r10271; double r10273 = r10256 / r10262; double r10274 = r10272 * r10273; double r10275 = r10258 - r10274; double r10276 = r10265 * r10275; return r10276; } double f_of(float d, float h, float l, float M, float D) { float r10277 = d; float r10278 = -9.519021483566862e-25f; bool r10279 = r10277 <= r10278; float r10280 = -1.0f; float r10281 = h; float r10282 = r10280 / r10281; float r10283 = sqrt(r10282); float r10284 = r10280 / r10277; float r10285 = 0.5f; float r10286 = pow(r10284, r10285); float r10287 = r10283 / r10286; float r10288 = l; float r10289 = r10277 / r10288; float r10290 = sqrt(r10289); float r10291 = 1.0f; float r10292 = 2.0f; float r10293 = r10291 / r10292; float r10294 = pow(r10290, r10293); float r10295 = cbrt(r10289); float r10296 = fabs(r10295); float r10297 = sqrt(r10295); float r10298 = r10296 * r10297; float r10299 = pow(r10298, r10293); float r10300 = r10294 * r10299; float r10301 = r10287 * r10300; float r10302 = M; float r10303 = D; float r10304 = r10302 * r10303; float r10305 = r10292 * r10277; float r10306 = r10304 / r10305; float r10307 = pow(r10306, r10292); float r10308 = r10293 * r10307; float r10309 = r10281 / r10288; float r10310 = r10308 * r10309; float r10311 = r10291 - r10310; float r10312 = r10301 * r10311; float r10313 = -1.09704788971517e-310f; bool r10314 = r10277 <= r10313; float r10315 = r10281 * r10307; float r10316 = r10288 * r10292; float r10317 = r10315 / r10316; float r10318 = r10291 - r10317; float r10319 = r10280 / r10288; float r10320 = sqrt(r10319); float r10321 = r10320 / r10286; float r10322 = r10277 / r10281; float r10323 = pow(r10322, r10293); float r10324 = r10321 * r10323; float r10325 = r10318 * r10324; float r10326 = 2.512962136109515e-183f; bool r10327 = r10277 <= r10326; float r10328 = r10291 / r10288; float r10329 = pow(r10328, r10293); float r10330 = pow(r10277, r10293); float r10331 = r10329 * r10330; float r10332 = r10323 * r10331; float r10333 = r10332 * r10318; float r10334 = sqrt(r10277); float r10335 = pow(r10281, r10285); float r10336 = r10334 / r10335; float r10337 = pow(r10289, r10293); float r10338 = r10336 * r10337; float r10339 = r10318 * r10338; float r10340 = r10327 ? r10333 : r10339; float r10341 = r10314 ? r10325 : r10340; float r10342 = r10279 ? r10312 : r10341; return r10342; } double f_od(double d, double h, double l, double M, double D) { double r10343 = d; double r10344 = -9.519021483566862e-25; bool r10345 = r10343 <= r10344; double r10346 = -1.0; double r10347 = h; double r10348 = r10346 / r10347; double r10349 = sqrt(r10348); double r10350 = r10346 / r10343; double r10351 = 0.5; double r10352 = pow(r10350, r10351); double r10353 = r10349 / r10352; double r10354 = l; double r10355 = r10343 / r10354; double r10356 = sqrt(r10355); double r10357 = 1.0; double r10358 = 2.0; double r10359 = r10357 / r10358; double r10360 = pow(r10356, r10359); double r10361 = cbrt(r10355); double r10362 = fabs(r10361); double r10363 = sqrt(r10361); double r10364 = r10362 * r10363; double r10365 = pow(r10364, r10359); double r10366 = r10360 * r10365; double r10367 = r10353 * r10366; double r10368 = M; double r10369 = D; double r10370 = r10368 * r10369; double r10371 = r10358 * r10343; double r10372 = r10370 / r10371; double r10373 = pow(r10372, r10358); double r10374 = r10359 * r10373; double r10375 = r10347 / r10354; double r10376 = r10374 * r10375; double r10377 = r10357 - r10376; double r10378 = r10367 * r10377; double r10379 = -1.09704788971517e-310; bool r10380 = r10343 <= r10379; double r10381 = r10347 * r10373; double r10382 = r10354 * r10358; double r10383 = r10381 / r10382; double r10384 = r10357 - r10383; double r10385 = r10346 / r10354; double r10386 = sqrt(r10385); double r10387 = r10386 / r10352; double r10388 = r10343 / r10347; double r10389 = pow(r10388, r10359); double r10390 = r10387 * r10389; double r10391 = r10384 * r10390; double r10392 = 2.512962136109515e-183; bool r10393 = r10343 <= r10392; double r10394 = r10357 / r10354; double r10395 = pow(r10394, r10359); double r10396 = pow(r10343, r10359); double r10397 = r10395 * r10396; double r10398 = r10389 * r10397; double r10399 = r10398 * r10384; double r10400 = sqrt(r10343); double r10401 = pow(r10347, r10351); double r10402 = r10400 / r10401; double r10403 = pow(r10355, r10359); double r10404 = r10402 * r10403; double r10405 = r10384 * r10404; double r10406 = r10393 ? r10399 : r10405; double r10407 = r10380 ? r10391 : r10406; double r10408 = r10345 ? r10378 : r10407; return r10408; } 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 r10409, r10410, r10411, r10412, r10413, r10414, r10415, r10416, r10417, r10418, r10419, r10420, r10421, r10422, r10423, r10424, r10425, r10426, r10427, r10428, r10429, r10430; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r10409); mpfr_init(r10410); mpfr_init(r10411); mpfr_init_set_str(r10412, "1", 10, MPFR_RNDN); mpfr_init_set_str(r10413, "2", 10, MPFR_RNDN); mpfr_init(r10414); mpfr_init(r10415); mpfr_init(r10416); mpfr_init(r10417); mpfr_init(r10418); mpfr_init(r10419); mpfr_init(r10420); mpfr_init(r10421); mpfr_init(r10422); mpfr_init(r10423); mpfr_init(r10424); mpfr_init(r10425); mpfr_init(r10426); mpfr_init(r10427); mpfr_init(r10428); mpfr_init(r10429); mpfr_init(r10430); } double f_im(double d, double h, double l, double M, double D) { mpfr_set_d(r10409, d, MPFR_RNDN); mpfr_set_d(r10410, h, MPFR_RNDN); mpfr_div(r10411, r10409, r10410, MPFR_RNDN); ; ; mpfr_div(r10414, r10412, r10413, MPFR_RNDN); mpfr_pow(r10415, r10411, r10414, MPFR_RNDN); mpfr_set_d(r10416, l, MPFR_RNDN); mpfr_div(r10417, r10409, r10416, MPFR_RNDN); mpfr_pow(r10418, r10417, r10414, MPFR_RNDN); mpfr_mul(r10419, r10415, r10418, MPFR_RNDN); mpfr_set_d(r10420, M, MPFR_RNDN); mpfr_set_d(r10421, D, MPFR_RNDN); mpfr_mul(r10422, r10420, r10421, MPFR_RNDN); mpfr_mul(r10423, r10413, r10409, MPFR_RNDN); mpfr_div(r10424, r10422, r10423, MPFR_RNDN); mpfr_pow(r10425, r10424, r10413, MPFR_RNDN); mpfr_mul(r10426, r10414, r10425, MPFR_RNDN); mpfr_div(r10427, r10410, r10416, MPFR_RNDN); mpfr_mul(r10428, r10426, r10427, MPFR_RNDN); mpfr_sub(r10429, r10412, r10428, MPFR_RNDN); mpfr_mul(r10430, r10419, r10429, MPFR_RNDN); return mpfr_get_d(r10430, MPFR_RNDN); } static mpfr_t r10431, r10432, r10433, r10434, r10435, r10436, r10437, r10438, r10439, r10440, r10441, r10442, r10443, r10444, r10445, r10446, r10447, r10448, r10449, r10450, r10451, r10452, r10453, r10454, r10455, r10456, r10457, r10458, r10459, r10460, r10461, r10462, r10463, r10464, r10465, r10466, r10467, r10468, r10469, r10470, r10471, r10472, r10473, r10474, r10475, r10476, r10477, r10478, r10479, r10480, r10481, r10482, r10483, r10484, r10485, r10486, r10487, r10488, r10489, r10490, r10491, r10492, r10493, r10494, r10495, r10496; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r10431); mpfr_init_set_str(r10432, "-9.519021483566862e-25", 10, MPFR_RNDN); mpfr_init(r10433); mpfr_init_set_str(r10434, "-1", 10, MPFR_RNDN); mpfr_init(r10435); mpfr_init(r10436); mpfr_init(r10437); mpfr_init(r10438); mpfr_init_set_str(r10439, "1/2", 10, MPFR_RNDN); mpfr_init(r10440); mpfr_init(r10441); mpfr_init(r10442); mpfr_init(r10443); mpfr_init(r10444); mpfr_init_set_str(r10445, "1", 10, MPFR_RNDN); mpfr_init_set_str(r10446, "2", 10, MPFR_RNDN); mpfr_init(r10447); mpfr_init(r10448); mpfr_init(r10449); mpfr_init(r10450); mpfr_init(r10451); mpfr_init(r10452); mpfr_init(r10453); mpfr_init(r10454); mpfr_init(r10455); mpfr_init(r10456); mpfr_init(r10457); mpfr_init(r10458); mpfr_init(r10459); mpfr_init(r10460); mpfr_init(r10461); mpfr_init(r10462); mpfr_init(r10463); mpfr_init(r10464); mpfr_init(r10465); mpfr_init(r10466); mpfr_init_set_str(r10467, "-1.09704788971517e-310", 10, MPFR_RNDN); mpfr_init(r10468); mpfr_init(r10469); mpfr_init(r10470); mpfr_init(r10471); mpfr_init(r10472); mpfr_init(r10473); mpfr_init(r10474); mpfr_init(r10475); mpfr_init(r10476); mpfr_init(r10477); mpfr_init(r10478); mpfr_init(r10479); mpfr_init_set_str(r10480, "2.512962136109515e-183", 10, MPFR_RNDN); mpfr_init(r10481); mpfr_init(r10482); mpfr_init(r10483); mpfr_init(r10484); mpfr_init(r10485); mpfr_init(r10486); mpfr_init(r10487); mpfr_init(r10488); mpfr_init(r10489); mpfr_init(r10490); mpfr_init(r10491); mpfr_init(r10492); mpfr_init(r10493); mpfr_init(r10494); mpfr_init(r10495); mpfr_init(r10496); } double f_fm(double d, double h, double l, double M, double D) { mpfr_set_d(r10431, d, MPFR_RNDN); ; mpfr_set_si(r10433, mpfr_cmp(r10431, r10432) <= 0, MPFR_RNDN); ; mpfr_set_d(r10435, h, MPFR_RNDN); mpfr_div(r10436, r10434, r10435, MPFR_RNDN); mpfr_sqrt(r10437, r10436, MPFR_RNDN); mpfr_div(r10438, r10434, r10431, MPFR_RNDN); ; mpfr_pow(r10440, r10438, r10439, MPFR_RNDN); mpfr_div(r10441, r10437, r10440, MPFR_RNDN); mpfr_set_d(r10442, l, MPFR_RNDN); mpfr_div(r10443, r10431, r10442, MPFR_RNDN); mpfr_sqrt(r10444, r10443, MPFR_RNDN); ; ; mpfr_div(r10447, r10445, r10446, MPFR_RNDN); mpfr_pow(r10448, r10444, r10447, MPFR_RNDN); mpfr_cbrt(r10449, r10443, MPFR_RNDN); mpfr_abs(r10450, r10449, MPFR_RNDN); mpfr_sqrt(r10451, r10449, MPFR_RNDN); mpfr_mul(r10452, r10450, r10451, MPFR_RNDN); mpfr_pow(r10453, r10452, r10447, MPFR_RNDN); mpfr_mul(r10454, r10448, r10453, MPFR_RNDN); mpfr_mul(r10455, r10441, r10454, MPFR_RNDN); mpfr_set_d(r10456, M, MPFR_RNDN); mpfr_set_d(r10457, D, MPFR_RNDN); mpfr_mul(r10458, r10456, r10457, MPFR_RNDN); mpfr_mul(r10459, r10446, r10431, MPFR_RNDN); mpfr_div(r10460, r10458, r10459, MPFR_RNDN); mpfr_pow(r10461, r10460, r10446, MPFR_RNDN); mpfr_mul(r10462, r10447, r10461, MPFR_RNDN); mpfr_div(r10463, r10435, r10442, MPFR_RNDN); mpfr_mul(r10464, r10462, r10463, MPFR_RNDN); mpfr_sub(r10465, r10445, r10464, MPFR_RNDN); mpfr_mul(r10466, r10455, r10465, MPFR_RNDN); ; mpfr_set_si(r10468, mpfr_cmp(r10431, r10467) <= 0, MPFR_RNDN); mpfr_mul(r10469, r10435, r10461, MPFR_RNDN); mpfr_mul(r10470, r10442, r10446, MPFR_RNDN); mpfr_div(r10471, r10469, r10470, MPFR_RNDN); mpfr_sub(r10472, r10445, r10471, MPFR_RNDN); mpfr_div(r10473, r10434, r10442, MPFR_RNDN); mpfr_sqrt(r10474, r10473, MPFR_RNDN); mpfr_div(r10475, r10474, r10440, MPFR_RNDN); mpfr_div(r10476, r10431, r10435, MPFR_RNDN); mpfr_pow(r10477, r10476, r10447, MPFR_RNDN); mpfr_mul(r10478, r10475, r10477, MPFR_RNDN); mpfr_mul(r10479, r10472, r10478, MPFR_RNDN); ; mpfr_set_si(r10481, mpfr_cmp(r10431, r10480) <= 0, MPFR_RNDN); mpfr_div(r10482, r10445, r10442, MPFR_RNDN); mpfr_pow(r10483, r10482, r10447, MPFR_RNDN); mpfr_pow(r10484, r10431, r10447, MPFR_RNDN); mpfr_mul(r10485, r10483, r10484, MPFR_RNDN); mpfr_mul(r10486, r10477, r10485, MPFR_RNDN); mpfr_mul(r10487, r10486, r10472, MPFR_RNDN); mpfr_sqrt(r10488, r10431, MPFR_RNDN); mpfr_pow(r10489, r10435, r10439, MPFR_RNDN); mpfr_div(r10490, r10488, r10489, MPFR_RNDN); mpfr_pow(r10491, r10443, r10447, MPFR_RNDN); mpfr_mul(r10492, r10490, r10491, MPFR_RNDN); mpfr_mul(r10493, r10472, r10492, MPFR_RNDN); if (mpfr_get_si(r10481, MPFR_RNDN)) { mpfr_set(r10494, r10487, MPFR_RNDN); } else { mpfr_set(r10494, r10493, MPFR_RNDN); }; if (mpfr_get_si(r10468, MPFR_RNDN)) { mpfr_set(r10495, r10479, MPFR_RNDN); } else { mpfr_set(r10495, r10494, MPFR_RNDN); }; if (mpfr_get_si(r10433, MPFR_RNDN)) { mpfr_set(r10496, r10466, MPFR_RNDN); } else { mpfr_set(r10496, r10495, MPFR_RNDN); }; return mpfr_get_d(r10496, MPFR_RNDN); } static mpfr_t r10497, r10498, r10499, r10500, r10501, r10502, r10503, r10504, r10505, r10506, r10507, r10508, r10509, r10510, r10511, r10512, r10513, r10514, r10515, r10516, r10517, r10518, r10519, r10520, r10521, r10522, r10523, r10524, r10525, r10526, r10527, r10528, r10529, r10530, r10531, r10532, r10533, r10534, r10535, r10536, r10537, r10538, r10539, r10540, r10541, r10542, r10543, r10544, r10545, r10546, r10547, r10548, r10549, r10550, r10551, r10552, r10553, r10554, r10555, r10556, r10557, r10558, r10559, r10560, r10561, r10562; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r10497); mpfr_init_set_str(r10498, "-9.519021483566862e-25", 10, MPFR_RNDN); mpfr_init(r10499); mpfr_init_set_str(r10500, "-1", 10, MPFR_RNDN); mpfr_init(r10501); mpfr_init(r10502); mpfr_init(r10503); mpfr_init(r10504); mpfr_init_set_str(r10505, "1/2", 10, MPFR_RNDN); mpfr_init(r10506); mpfr_init(r10507); mpfr_init(r10508); mpfr_init(r10509); mpfr_init(r10510); mpfr_init_set_str(r10511, "1", 10, MPFR_RNDN); mpfr_init_set_str(r10512, "2", 10, MPFR_RNDN); mpfr_init(r10513); mpfr_init(r10514); mpfr_init(r10515); mpfr_init(r10516); mpfr_init(r10517); mpfr_init(r10518); mpfr_init(r10519); mpfr_init(r10520); mpfr_init(r10521); mpfr_init(r10522); mpfr_init(r10523); mpfr_init(r10524); mpfr_init(r10525); mpfr_init(r10526); mpfr_init(r10527); mpfr_init(r10528); mpfr_init(r10529); mpfr_init(r10530); mpfr_init(r10531); mpfr_init(r10532); mpfr_init_set_str(r10533, "-1.09704788971517e-310", 10, MPFR_RNDN); mpfr_init(r10534); mpfr_init(r10535); mpfr_init(r10536); mpfr_init(r10537); mpfr_init(r10538); mpfr_init(r10539); mpfr_init(r10540); mpfr_init(r10541); mpfr_init(r10542); mpfr_init(r10543); mpfr_init(r10544); mpfr_init(r10545); mpfr_init_set_str(r10546, "2.512962136109515e-183", 10, MPFR_RNDN); mpfr_init(r10547); mpfr_init(r10548); mpfr_init(r10549); mpfr_init(r10550); mpfr_init(r10551); mpfr_init(r10552); mpfr_init(r10553); mpfr_init(r10554); mpfr_init(r10555); mpfr_init(r10556); mpfr_init(r10557); mpfr_init(r10558); mpfr_init(r10559); mpfr_init(r10560); mpfr_init(r10561); mpfr_init(r10562); } double f_dm(double d, double h, double l, double M, double D) { mpfr_set_d(r10497, d, MPFR_RNDN); ; mpfr_set_si(r10499, mpfr_cmp(r10497, r10498) <= 0, MPFR_RNDN); ; mpfr_set_d(r10501, h, MPFR_RNDN); mpfr_div(r10502, r10500, r10501, MPFR_RNDN); mpfr_sqrt(r10503, r10502, MPFR_RNDN); mpfr_div(r10504, r10500, r10497, MPFR_RNDN); ; mpfr_pow(r10506, r10504, r10505, MPFR_RNDN); mpfr_div(r10507, r10503, r10506, MPFR_RNDN); mpfr_set_d(r10508, l, MPFR_RNDN); mpfr_div(r10509, r10497, r10508, MPFR_RNDN); mpfr_sqrt(r10510, r10509, MPFR_RNDN); ; ; mpfr_div(r10513, r10511, r10512, MPFR_RNDN); mpfr_pow(r10514, r10510, r10513, MPFR_RNDN); mpfr_cbrt(r10515, r10509, MPFR_RNDN); mpfr_abs(r10516, r10515, MPFR_RNDN); mpfr_sqrt(r10517, r10515, MPFR_RNDN); mpfr_mul(r10518, r10516, r10517, MPFR_RNDN); mpfr_pow(r10519, r10518, r10513, MPFR_RNDN); mpfr_mul(r10520, r10514, r10519, MPFR_RNDN); mpfr_mul(r10521, r10507, r10520, MPFR_RNDN); mpfr_set_d(r10522, M, MPFR_RNDN); mpfr_set_d(r10523, D, MPFR_RNDN); mpfr_mul(r10524, r10522, r10523, MPFR_RNDN); mpfr_mul(r10525, r10512, r10497, MPFR_RNDN); mpfr_div(r10526, r10524, r10525, MPFR_RNDN); mpfr_pow(r10527, r10526, r10512, MPFR_RNDN); mpfr_mul(r10528, r10513, r10527, MPFR_RNDN); mpfr_div(r10529, r10501, r10508, MPFR_RNDN); mpfr_mul(r10530, r10528, r10529, MPFR_RNDN); mpfr_sub(r10531, r10511, r10530, MPFR_RNDN); mpfr_mul(r10532, r10521, r10531, MPFR_RNDN); ; mpfr_set_si(r10534, mpfr_cmp(r10497, r10533) <= 0, MPFR_RNDN); mpfr_mul(r10535, r10501, r10527, MPFR_RNDN); mpfr_mul(r10536, r10508, r10512, MPFR_RNDN); mpfr_div(r10537, r10535, r10536, MPFR_RNDN); mpfr_sub(r10538, r10511, r10537, MPFR_RNDN); mpfr_div(r10539, r10500, r10508, MPFR_RNDN); mpfr_sqrt(r10540, r10539, MPFR_RNDN); mpfr_div(r10541, r10540, r10506, MPFR_RNDN); mpfr_div(r10542, r10497, r10501, MPFR_RNDN); mpfr_pow(r10543, r10542, r10513, MPFR_RNDN); mpfr_mul(r10544, r10541, r10543, MPFR_RNDN); mpfr_mul(r10545, r10538, r10544, MPFR_RNDN); ; mpfr_set_si(r10547, mpfr_cmp(r10497, r10546) <= 0, MPFR_RNDN); mpfr_div(r10548, r10511, r10508, MPFR_RNDN); mpfr_pow(r10549, r10548, r10513, MPFR_RNDN); mpfr_pow(r10550, r10497, r10513, MPFR_RNDN); mpfr_mul(r10551, r10549, r10550, MPFR_RNDN); mpfr_mul(r10552, r10543, r10551, MPFR_RNDN); mpfr_mul(r10553, r10552, r10538, MPFR_RNDN); mpfr_sqrt(r10554, r10497, MPFR_RNDN); mpfr_pow(r10555, r10501, r10505, MPFR_RNDN); mpfr_div(r10556, r10554, r10555, MPFR_RNDN); mpfr_pow(r10557, r10509, r10513, MPFR_RNDN); mpfr_mul(r10558, r10556, r10557, MPFR_RNDN); mpfr_mul(r10559, r10538, r10558, MPFR_RNDN); if (mpfr_get_si(r10547, MPFR_RNDN)) { mpfr_set(r10560, r10553, MPFR_RNDN); } else { mpfr_set(r10560, r10559, MPFR_RNDN); }; if (mpfr_get_si(r10534, MPFR_RNDN)) { mpfr_set(r10561, r10545, MPFR_RNDN); } else { mpfr_set(r10561, r10560, MPFR_RNDN); }; if (mpfr_get_si(r10499, MPFR_RNDN)) { mpfr_set(r10562, r10532, MPFR_RNDN); } else { mpfr_set(r10562, r10561, MPFR_RNDN); }; return mpfr_get_d(r10562, MPFR_RNDN); }