#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 r10307 = w0; float r10308 = 1; float r10309 = M; float r10310 = D; float r10311 = r10309 * r10310; float r10312 = 2; float r10313 = d; float r10314 = r10312 * r10313; float r10315 = r10311 / r10314; float r10316 = pow(r10315, r10312); float r10317 = h; float r10318 = l; float r10319 = r10317 / r10318; float r10320 = r10316 * r10319; float r10321 = r10308 - r10320; float r10322 = sqrt(r10321); float r10323 = r10307 * r10322; return r10323; } double f_id(double w0, double M, double D, double h, double l, double d) { double r10324 = w0; double r10325 = 1; double r10326 = M; double r10327 = D; double r10328 = r10326 * r10327; double r10329 = 2; double r10330 = d; double r10331 = r10329 * r10330; double r10332 = r10328 / r10331; double r10333 = pow(r10332, r10329); double r10334 = h; double r10335 = l; double r10336 = r10334 / r10335; double r10337 = r10333 * r10336; double r10338 = r10325 - r10337; double r10339 = sqrt(r10338); double r10340 = r10324 * r10339; return r10340; } double f_of(float w0, float M, float D, float h, float l, float d) { float r10341 = M; float r10342 = D; float r10343 = r10341 * r10342; float r10344 = 2; float r10345 = d; float r10346 = r10344 * r10345; float r10347 = r10343 / r10346; float r10348 = pow(r10347, r10344); float r10349 = h; float r10350 = l; float r10351 = r10349 / r10350; float r10352 = r10348 * r10351; float r10353 = -6.503081529298164e-10; bool r10354 = r10352 <= r10353; float r10355 = w0; float r10356 = 1; float r10357 = r10347 * r10351; float r10358 = r10347 * r10357; float r10359 = r10356 - r10358; float r10360 = sqrt(r10359); float r10361 = r10355 * r10360; float r10362 = r10354 ? r10361 : r10355; return r10362; } double f_od(double w0, double M, double D, double h, double l, double d) { double r10363 = M; double r10364 = D; double r10365 = r10363 * r10364; double r10366 = 2; double r10367 = d; double r10368 = r10366 * r10367; double r10369 = r10365 / r10368; double r10370 = pow(r10369, r10366); double r10371 = h; double r10372 = l; double r10373 = r10371 / r10372; double r10374 = r10370 * r10373; double r10375 = -6.503081529298164e-10; bool r10376 = r10374 <= r10375; double r10377 = w0; double r10378 = 1; double r10379 = r10369 * r10373; double r10380 = r10369 * r10379; double r10381 = r10378 - r10380; double r10382 = sqrt(r10381); double r10383 = r10377 * r10382; double r10384 = r10376 ? r10383 : r10377; return r10384; } 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 r10385, r10386, r10387, r10388, r10389, r10390, r10391, r10392, r10393, r10394, r10395, r10396, r10397, r10398, r10399, r10400, r10401; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r10385); mpfr_init_set_str(r10386, "1", 10, MPFR_RNDN); mpfr_init(r10387); mpfr_init(r10388); mpfr_init(r10389); mpfr_init_set_str(r10390, "2", 10, MPFR_RNDN); mpfr_init(r10391); mpfr_init(r10392); mpfr_init(r10393); mpfr_init(r10394); mpfr_init(r10395); mpfr_init(r10396); mpfr_init(r10397); mpfr_init(r10398); mpfr_init(r10399); mpfr_init(r10400); mpfr_init(r10401); } double f_im(double w0, double M, double D, double h, double l, double d) { mpfr_set_d(r10385, w0, MPFR_RNDN); ; mpfr_set_d(r10387, M, MPFR_RNDN); mpfr_set_d(r10388, D, MPFR_RNDN); mpfr_mul(r10389, r10387, r10388, MPFR_RNDN); ; mpfr_set_d(r10391, d, MPFR_RNDN); mpfr_mul(r10392, r10390, r10391, MPFR_RNDN); mpfr_div(r10393, r10389, r10392, MPFR_RNDN); mpfr_pow(r10394, r10393, r10390, MPFR_RNDN); mpfr_set_d(r10395, h, MPFR_RNDN); mpfr_set_d(r10396, l, MPFR_RNDN); mpfr_div(r10397, r10395, r10396, MPFR_RNDN); mpfr_mul(r10398, r10394, r10397, MPFR_RNDN); mpfr_sub(r10399, r10386, r10398, MPFR_RNDN); mpfr_sqrt(r10400, r10399, MPFR_RNDN); mpfr_mul(r10401, r10385, r10400, MPFR_RNDN); return mpfr_get_d(r10401, MPFR_RNDN); } static mpfr_t r10402, r10403, r10404, r10405, r10406, r10407, r10408, r10409, r10410, r10411, r10412, r10413, r10414, r10415, r10416, r10417, r10418, r10419, r10420, r10421, r10422, r10423; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r10402); mpfr_init(r10403); mpfr_init(r10404); mpfr_init_set_str(r10405, "2", 10, MPFR_RNDN); mpfr_init(r10406); mpfr_init(r10407); mpfr_init(r10408); mpfr_init(r10409); mpfr_init(r10410); mpfr_init(r10411); mpfr_init(r10412); mpfr_init(r10413); mpfr_init_set_str(r10414, "-6.503081529298164e-10", 10, MPFR_RNDN); mpfr_init(r10415); mpfr_init(r10416); mpfr_init_set_str(r10417, "1", 10, MPFR_RNDN); mpfr_init(r10418); mpfr_init(r10419); mpfr_init(r10420); mpfr_init(r10421); mpfr_init(r10422); mpfr_init(r10423); } double f_fm(double w0, double M, double D, double h, double l, double d) { mpfr_set_d(r10402, M, MPFR_RNDN); mpfr_set_d(r10403, D, MPFR_RNDN); mpfr_mul(r10404, r10402, r10403, MPFR_RNDN); ; mpfr_set_d(r10406, d, MPFR_RNDN); mpfr_mul(r10407, r10405, r10406, MPFR_RNDN); mpfr_div(r10408, r10404, r10407, MPFR_RNDN); mpfr_pow(r10409, r10408, r10405, MPFR_RNDN); mpfr_set_d(r10410, h, MPFR_RNDN); mpfr_set_d(r10411, l, MPFR_RNDN); mpfr_div(r10412, r10410, r10411, MPFR_RNDN); mpfr_mul(r10413, r10409, r10412, MPFR_RNDN); ; mpfr_set_si(r10415, mpfr_cmp(r10413, r10414) <= 0, MPFR_RNDN); mpfr_set_d(r10416, w0, MPFR_RNDN); ; mpfr_mul(r10418, r10408, r10412, MPFR_RNDN); mpfr_mul(r10419, r10408, r10418, MPFR_RNDN); mpfr_sub(r10420, r10417, r10419, MPFR_RNDN); mpfr_sqrt(r10421, r10420, MPFR_RNDN); mpfr_mul(r10422, r10416, r10421, MPFR_RNDN); if (mpfr_get_si(r10415, MPFR_RNDN)) { mpfr_set(r10423, r10422, MPFR_RNDN); } else { mpfr_set(r10423, r10416, MPFR_RNDN); }; return mpfr_get_d(r10423, MPFR_RNDN); } static mpfr_t r10424, r10425, r10426, r10427, r10428, r10429, r10430, r10431, r10432, r10433, r10434, r10435, r10436, r10437, r10438, r10439, r10440, r10441, r10442, r10443, r10444, r10445; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r10424); mpfr_init(r10425); mpfr_init(r10426); mpfr_init_set_str(r10427, "2", 10, MPFR_RNDN); mpfr_init(r10428); mpfr_init(r10429); mpfr_init(r10430); mpfr_init(r10431); mpfr_init(r10432); mpfr_init(r10433); mpfr_init(r10434); mpfr_init(r10435); mpfr_init_set_str(r10436, "-6.503081529298164e-10", 10, MPFR_RNDN); mpfr_init(r10437); mpfr_init(r10438); mpfr_init_set_str(r10439, "1", 10, MPFR_RNDN); mpfr_init(r10440); mpfr_init(r10441); mpfr_init(r10442); mpfr_init(r10443); mpfr_init(r10444); mpfr_init(r10445); } double f_dm(double w0, double M, double D, double h, double l, double d) { mpfr_set_d(r10424, M, MPFR_RNDN); mpfr_set_d(r10425, D, MPFR_RNDN); mpfr_mul(r10426, r10424, r10425, MPFR_RNDN); ; mpfr_set_d(r10428, d, MPFR_RNDN); mpfr_mul(r10429, r10427, r10428, MPFR_RNDN); mpfr_div(r10430, r10426, r10429, MPFR_RNDN); mpfr_pow(r10431, r10430, r10427, MPFR_RNDN); mpfr_set_d(r10432, h, MPFR_RNDN); mpfr_set_d(r10433, l, MPFR_RNDN); mpfr_div(r10434, r10432, r10433, MPFR_RNDN); mpfr_mul(r10435, r10431, r10434, MPFR_RNDN); ; mpfr_set_si(r10437, mpfr_cmp(r10435, r10436) <= 0, MPFR_RNDN); mpfr_set_d(r10438, w0, MPFR_RNDN); ; mpfr_mul(r10440, r10430, r10434, MPFR_RNDN); mpfr_mul(r10441, r10430, r10440, MPFR_RNDN); mpfr_sub(r10442, r10439, r10441, MPFR_RNDN); mpfr_sqrt(r10443, r10442, MPFR_RNDN); mpfr_mul(r10444, r10438, r10443, MPFR_RNDN); if (mpfr_get_si(r10437, MPFR_RNDN)) { mpfr_set(r10445, r10444, MPFR_RNDN); } else { mpfr_set(r10445, r10438, MPFR_RNDN); }; return mpfr_get_d(r10445, MPFR_RNDN); }