#include #include #include #include #include char *name = "Henrywood and Agarwal, Equation (13)"; double f_if(float c0, float w, float h, float D, float d, float M) { float r10205 = c0; float r10206 = 2; float r10207 = w; float r10208 = r10206 * r10207; float r10209 = r10205 / r10208; float r10210 = d; float r10211 = r10210 * r10210; float r10212 = r10205 * r10211; float r10213 = h; float r10214 = r10207 * r10213; float r10215 = D; float r10216 = r10215 * r10215; float r10217 = r10214 * r10216; float r10218 = r10212 / r10217; float r10219 = r10218 * r10218; float r10220 = M; float r10221 = r10220 * r10220; float r10222 = r10219 - r10221; float r10223 = sqrt(r10222); float r10224 = r10218 + r10223; float r10225 = r10209 * r10224; return r10225; } double f_id(double c0, double w, double h, double D, double d, double M) { double r10226 = c0; double r10227 = 2; double r10228 = w; double r10229 = r10227 * r10228; double r10230 = r10226 / r10229; double r10231 = d; double r10232 = r10231 * r10231; double r10233 = r10226 * r10232; double r10234 = h; double r10235 = r10228 * r10234; double r10236 = D; double r10237 = r10236 * r10236; double r10238 = r10235 * r10237; double r10239 = r10233 / r10238; double r10240 = r10239 * r10239; double r10241 = M; double r10242 = r10241 * r10241; double r10243 = r10240 - r10242; double r10244 = sqrt(r10243); double r10245 = r10239 + r10244; double r10246 = r10230 * r10245; return r10246; } double f_of(float c0, float w, float h, float D, float d, float M) { float r10247 = d; float r10248 = D; float r10249 = r10247 / r10248; float r10250 = -2.3071785278265444e+40; bool r10251 = r10249 <= r10250; float r10252 = 0; float r10253 = -4.3953808869808786e-51; bool r10254 = r10249 <= r10253; float r10255 = c0; float r10256 = w; float r10257 = 2; float r10258 = r10256 * r10257; float r10259 = M; float r10260 = r10259 * r10259; float r10261 = r10258 / r10260; float r10262 = r10255 / r10261; float r10263 = r10249 * r10249; float r10264 = r10255 / r10256; float r10265 = h; float r10266 = r10264 / r10265; float r10267 = r10263 * r10266; float r10268 = r10267 * r10267; float r10269 = r10268 - r10260; float r10270 = sqrt(r10269); float r10271 = r10267 - r10270; float r10272 = r10262 / r10271; float r10273 = 3.0150778316863736e-12; bool r10274 = r10249 <= r10273; float r10275 = r10257 * r10256; float r10276 = r10255 / r10275; float r10277 = r10266 * r10263; float r10278 = r10259 + r10277; float r10279 = r10277 - r10259; float r10280 = r10278 * r10279; float r10281 = sqrt(r10280); float r10282 = r10277 + r10281; float r10283 = 3; float r10284 = pow(r10282, r10283); float r10285 = cbrt(r10284); float r10286 = r10276 * r10285; float r10287 = 8.43772071865126e+109; bool r10288 = r10249 <= r10287; float r10289 = r10255 * r10259; float r10290 = r10289 / r10258; float r10291 = r10259 * r10290; float r10292 = r10255 / r10265; float r10293 = r10292 / r10256; float r10294 = r10293 * r10263; float r10295 = r10294 * r10294; float r10296 = r10295 - r10260; float r10297 = sqrt(r10296); float r10298 = r10294 - r10297; float r10299 = r10291 / r10298; float r10300 = 2.6563604190423177e+110; bool r10301 = r10249 <= r10300; float r10302 = r10301 ? r10286 : r10252; float r10303 = r10288 ? r10299 : r10302; float r10304 = r10274 ? r10286 : r10303; float r10305 = r10254 ? r10272 : r10304; float r10306 = r10251 ? r10252 : r10305; return r10306; } double f_od(double c0, double w, double h, double D, double d, double M) { double r10307 = d; double r10308 = D; double r10309 = r10307 / r10308; double r10310 = -2.3071785278265444e+40; bool r10311 = r10309 <= r10310; double r10312 = 0; double r10313 = -4.3953808869808786e-51; bool r10314 = r10309 <= r10313; double r10315 = c0; double r10316 = w; double r10317 = 2; double r10318 = r10316 * r10317; double r10319 = M; double r10320 = r10319 * r10319; double r10321 = r10318 / r10320; double r10322 = r10315 / r10321; double r10323 = r10309 * r10309; double r10324 = r10315 / r10316; double r10325 = h; double r10326 = r10324 / r10325; double r10327 = r10323 * r10326; double r10328 = r10327 * r10327; double r10329 = r10328 - r10320; double r10330 = sqrt(r10329); double r10331 = r10327 - r10330; double r10332 = r10322 / r10331; double r10333 = 3.0150778316863736e-12; bool r10334 = r10309 <= r10333; double r10335 = r10317 * r10316; double r10336 = r10315 / r10335; double r10337 = r10326 * r10323; double r10338 = r10319 + r10337; double r10339 = r10337 - r10319; double r10340 = r10338 * r10339; double r10341 = sqrt(r10340); double r10342 = r10337 + r10341; double r10343 = 3; double r10344 = pow(r10342, r10343); double r10345 = cbrt(r10344); double r10346 = r10336 * r10345; double r10347 = 8.43772071865126e+109; bool r10348 = r10309 <= r10347; double r10349 = r10315 * r10319; double r10350 = r10349 / r10318; double r10351 = r10319 * r10350; double r10352 = r10315 / r10325; double r10353 = r10352 / r10316; double r10354 = r10353 * r10323; double r10355 = r10354 * r10354; double r10356 = r10355 - r10320; double r10357 = sqrt(r10356); double r10358 = r10354 - r10357; double r10359 = r10351 / r10358; double r10360 = 2.6563604190423177e+110; bool r10361 = r10309 <= r10360; double r10362 = r10361 ? r10346 : r10312; double r10363 = r10348 ? r10359 : r10362; double r10364 = r10334 ? r10346 : r10363; double r10365 = r10314 ? r10332 : r10364; double r10366 = r10311 ? r10312 : r10365; return r10366; } 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 r10367, r10368, r10369, r10370, r10371, r10372, r10373, r10374, r10375, r10376, r10377, r10378, r10379, r10380, r10381, r10382, r10383, r10384, r10385, r10386, r10387; void setup_mpfr_f_im() { mpfr_set_default_prec(6992); mpfr_init(r10367); mpfr_init_set_str(r10368, "2", 10, MPFR_RNDN); mpfr_init(r10369); mpfr_init(r10370); mpfr_init(r10371); mpfr_init(r10372); mpfr_init(r10373); mpfr_init(r10374); mpfr_init(r10375); mpfr_init(r10376); mpfr_init(r10377); mpfr_init(r10378); mpfr_init(r10379); mpfr_init(r10380); mpfr_init(r10381); mpfr_init(r10382); mpfr_init(r10383); mpfr_init(r10384); mpfr_init(r10385); mpfr_init(r10386); mpfr_init(r10387); } double f_im(double c0, double w, double h, double D, double d, double M) { mpfr_set_d(r10367, c0, MPFR_RNDN); ; mpfr_set_d(r10369, w, MPFR_RNDN); mpfr_mul(r10370, r10368, r10369, MPFR_RNDN); mpfr_div(r10371, r10367, r10370, MPFR_RNDN); mpfr_set_d(r10372, d, MPFR_RNDN); mpfr_mul(r10373, r10372, r10372, MPFR_RNDN); mpfr_mul(r10374, r10367, r10373, MPFR_RNDN); mpfr_set_d(r10375, h, MPFR_RNDN); mpfr_mul(r10376, r10369, r10375, MPFR_RNDN); mpfr_set_d(r10377, D, MPFR_RNDN); mpfr_mul(r10378, r10377, r10377, MPFR_RNDN); mpfr_mul(r10379, r10376, r10378, MPFR_RNDN); mpfr_div(r10380, r10374, r10379, MPFR_RNDN); mpfr_mul(r10381, r10380, r10380, MPFR_RNDN); mpfr_set_d(r10382, M, MPFR_RNDN); mpfr_mul(r10383, r10382, r10382, MPFR_RNDN); mpfr_sub(r10384, r10381, r10383, MPFR_RNDN); mpfr_sqrt(r10385, r10384, MPFR_RNDN); mpfr_add(r10386, r10380, r10385, MPFR_RNDN); mpfr_mul(r10387, r10371, r10386, MPFR_RNDN); return mpfr_get_d(r10387, MPFR_RNDN); } static mpfr_t r10388, r10389, r10390, r10391, r10392, r10393, r10394, r10395, r10396, r10397, r10398, r10399, r10400, r10401, r10402, r10403, r10404, r10405, r10406, r10407, r10408, r10409, r10410, r10411, r10412, r10413, r10414, r10415, r10416, r10417, r10418, r10419, r10420, r10421, r10422, r10423, r10424, r10425, r10426, r10427, r10428, r10429, r10430, r10431, r10432, r10433, r10434, r10435, r10436, r10437, r10438, r10439, r10440, r10441, r10442, r10443, r10444, r10445, r10446, r10447; void setup_mpfr_f_fm() { mpfr_set_default_prec(6992); mpfr_init(r10388); mpfr_init(r10389); mpfr_init(r10390); mpfr_init_set_str(r10391, "-2.3071785278265444e+40", 10, MPFR_RNDN); mpfr_init(r10392); mpfr_init_set_str(r10393, "0", 10, MPFR_RNDN); mpfr_init_set_str(r10394, "-4.3953808869808786e-51", 10, MPFR_RNDN); mpfr_init(r10395); mpfr_init(r10396); mpfr_init(r10397); mpfr_init_set_str(r10398, "2", 10, MPFR_RNDN); mpfr_init(r10399); mpfr_init(r10400); mpfr_init(r10401); mpfr_init(r10402); mpfr_init(r10403); mpfr_init(r10404); mpfr_init(r10405); 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, "3.0150778316863736e-12", 10, MPFR_RNDN); 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_set_str(r10424, "3", 10, MPFR_RNDN); mpfr_init(r10425); mpfr_init(r10426); mpfr_init(r10427); mpfr_init_set_str(r10428, "8.43772071865126e+109", 10, MPFR_RNDN); mpfr_init(r10429); mpfr_init(r10430); mpfr_init(r10431); mpfr_init(r10432); mpfr_init(r10433); mpfr_init(r10434); mpfr_init(r10435); mpfr_init(r10436); mpfr_init(r10437); mpfr_init(r10438); mpfr_init(r10439); mpfr_init(r10440); mpfr_init_set_str(r10441, "2.6563604190423177e+110", 10, MPFR_RNDN); mpfr_init(r10442); mpfr_init(r10443); mpfr_init(r10444); mpfr_init(r10445); mpfr_init(r10446); mpfr_init(r10447); } double f_fm(double c0, double w, double h, double D, double d, double M) { mpfr_set_d(r10388, d, MPFR_RNDN); mpfr_set_d(r10389, D, MPFR_RNDN); mpfr_div(r10390, r10388, r10389, MPFR_RNDN); ; mpfr_set_si(r10392, mpfr_cmp(r10390, r10391) <= 0, MPFR_RNDN); ; ; mpfr_set_si(r10395, mpfr_cmp(r10390, r10394) <= 0, MPFR_RNDN); mpfr_set_d(r10396, c0, MPFR_RNDN); mpfr_set_d(r10397, w, MPFR_RNDN); ; mpfr_mul(r10399, r10397, r10398, MPFR_RNDN); mpfr_set_d(r10400, M, MPFR_RNDN); mpfr_mul(r10401, r10400, r10400, MPFR_RNDN); mpfr_div(r10402, r10399, r10401, MPFR_RNDN); mpfr_div(r10403, r10396, r10402, MPFR_RNDN); mpfr_mul(r10404, r10390, r10390, MPFR_RNDN); mpfr_div(r10405, r10396, r10397, MPFR_RNDN); mpfr_set_d(r10406, h, MPFR_RNDN); mpfr_div(r10407, r10405, r10406, MPFR_RNDN); mpfr_mul(r10408, r10404, r10407, MPFR_RNDN); mpfr_mul(r10409, r10408, r10408, MPFR_RNDN); mpfr_sub(r10410, r10409, r10401, MPFR_RNDN); mpfr_sqrt(r10411, r10410, MPFR_RNDN); mpfr_sub(r10412, r10408, r10411, MPFR_RNDN); mpfr_div(r10413, r10403, r10412, MPFR_RNDN); ; mpfr_set_si(r10415, mpfr_cmp(r10390, r10414) <= 0, MPFR_RNDN); mpfr_mul(r10416, r10398, r10397, MPFR_RNDN); mpfr_div(r10417, r10396, r10416, MPFR_RNDN); mpfr_mul(r10418, r10407, r10404, MPFR_RNDN); mpfr_add(r10419, r10400, r10418, MPFR_RNDN); mpfr_sub(r10420, r10418, r10400, MPFR_RNDN); mpfr_mul(r10421, r10419, r10420, MPFR_RNDN); mpfr_sqrt(r10422, r10421, MPFR_RNDN); mpfr_add(r10423, r10418, r10422, MPFR_RNDN); ; mpfr_pow(r10425, r10423, r10424, MPFR_RNDN); mpfr_cbrt(r10426, r10425, MPFR_RNDN); mpfr_mul(r10427, r10417, r10426, MPFR_RNDN); ; mpfr_set_si(r10429, mpfr_cmp(r10390, r10428) <= 0, MPFR_RNDN); mpfr_mul(r10430, r10396, r10400, MPFR_RNDN); mpfr_div(r10431, r10430, r10399, MPFR_RNDN); mpfr_mul(r10432, r10400, r10431, MPFR_RNDN); mpfr_div(r10433, r10396, r10406, MPFR_RNDN); mpfr_div(r10434, r10433, r10397, MPFR_RNDN); mpfr_mul(r10435, r10434, r10404, MPFR_RNDN); mpfr_mul(r10436, r10435, r10435, MPFR_RNDN); mpfr_sub(r10437, r10436, r10401, MPFR_RNDN); mpfr_sqrt(r10438, r10437, MPFR_RNDN); mpfr_sub(r10439, r10435, r10438, MPFR_RNDN); mpfr_div(r10440, r10432, r10439, MPFR_RNDN); ; mpfr_set_si(r10442, mpfr_cmp(r10390, r10441) <= 0, MPFR_RNDN); if (mpfr_get_si(r10442, MPFR_RNDN)) { mpfr_set(r10443, r10427, MPFR_RNDN); } else { mpfr_set(r10443, r10393, MPFR_RNDN); }; if (mpfr_get_si(r10429, MPFR_RNDN)) { mpfr_set(r10444, r10440, MPFR_RNDN); } else { mpfr_set(r10444, r10443, MPFR_RNDN); }; if (mpfr_get_si(r10415, MPFR_RNDN)) { mpfr_set(r10445, r10427, MPFR_RNDN); } else { mpfr_set(r10445, r10444, MPFR_RNDN); }; if (mpfr_get_si(r10395, MPFR_RNDN)) { mpfr_set(r10446, r10413, MPFR_RNDN); } else { mpfr_set(r10446, r10445, MPFR_RNDN); }; if (mpfr_get_si(r10392, MPFR_RNDN)) { mpfr_set(r10447, r10393, MPFR_RNDN); } else { mpfr_set(r10447, r10446, MPFR_RNDN); }; return mpfr_get_d(r10447, MPFR_RNDN); } static mpfr_t 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, r10497, r10498, r10499, r10500, r10501, r10502, r10503, r10504, r10505, r10506, r10507; void setup_mpfr_f_dm() { mpfr_set_default_prec(6992); mpfr_init(r10448); mpfr_init(r10449); mpfr_init(r10450); mpfr_init_set_str(r10451, "-2.3071785278265444e+40", 10, MPFR_RNDN); mpfr_init(r10452); mpfr_init_set_str(r10453, "0", 10, MPFR_RNDN); mpfr_init_set_str(r10454, "-4.3953808869808786e-51", 10, MPFR_RNDN); mpfr_init(r10455); mpfr_init(r10456); mpfr_init(r10457); mpfr_init_set_str(r10458, "2", 10, MPFR_RNDN); 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(r10467); mpfr_init(r10468); mpfr_init(r10469); mpfr_init(r10470); mpfr_init(r10471); mpfr_init(r10472); mpfr_init(r10473); mpfr_init_set_str(r10474, "3.0150778316863736e-12", 10, MPFR_RNDN); mpfr_init(r10475); mpfr_init(r10476); mpfr_init(r10477); mpfr_init(r10478); mpfr_init(r10479); mpfr_init(r10480); mpfr_init(r10481); mpfr_init(r10482); mpfr_init(r10483); mpfr_init_set_str(r10484, "3", 10, MPFR_RNDN); mpfr_init(r10485); mpfr_init(r10486); mpfr_init(r10487); mpfr_init_set_str(r10488, "8.43772071865126e+109", 10, MPFR_RNDN); mpfr_init(r10489); mpfr_init(r10490); mpfr_init(r10491); mpfr_init(r10492); mpfr_init(r10493); mpfr_init(r10494); mpfr_init(r10495); mpfr_init(r10496); mpfr_init(r10497); mpfr_init(r10498); mpfr_init(r10499); mpfr_init(r10500); mpfr_init_set_str(r10501, "2.6563604190423177e+110", 10, MPFR_RNDN); mpfr_init(r10502); mpfr_init(r10503); mpfr_init(r10504); mpfr_init(r10505); mpfr_init(r10506); mpfr_init(r10507); } double f_dm(double c0, double w, double h, double D, double d, double M) { mpfr_set_d(r10448, d, MPFR_RNDN); mpfr_set_d(r10449, D, MPFR_RNDN); mpfr_div(r10450, r10448, r10449, MPFR_RNDN); ; mpfr_set_si(r10452, mpfr_cmp(r10450, r10451) <= 0, MPFR_RNDN); ; ; mpfr_set_si(r10455, mpfr_cmp(r10450, r10454) <= 0, MPFR_RNDN); mpfr_set_d(r10456, c0, MPFR_RNDN); mpfr_set_d(r10457, w, MPFR_RNDN); ; mpfr_mul(r10459, r10457, r10458, MPFR_RNDN); mpfr_set_d(r10460, M, MPFR_RNDN); mpfr_mul(r10461, r10460, r10460, MPFR_RNDN); mpfr_div(r10462, r10459, r10461, MPFR_RNDN); mpfr_div(r10463, r10456, r10462, MPFR_RNDN); mpfr_mul(r10464, r10450, r10450, MPFR_RNDN); mpfr_div(r10465, r10456, r10457, MPFR_RNDN); mpfr_set_d(r10466, h, MPFR_RNDN); mpfr_div(r10467, r10465, r10466, MPFR_RNDN); mpfr_mul(r10468, r10464, r10467, MPFR_RNDN); mpfr_mul(r10469, r10468, r10468, MPFR_RNDN); mpfr_sub(r10470, r10469, r10461, MPFR_RNDN); mpfr_sqrt(r10471, r10470, MPFR_RNDN); mpfr_sub(r10472, r10468, r10471, MPFR_RNDN); mpfr_div(r10473, r10463, r10472, MPFR_RNDN); ; mpfr_set_si(r10475, mpfr_cmp(r10450, r10474) <= 0, MPFR_RNDN); mpfr_mul(r10476, r10458, r10457, MPFR_RNDN); mpfr_div(r10477, r10456, r10476, MPFR_RNDN); mpfr_mul(r10478, r10467, r10464, MPFR_RNDN); mpfr_add(r10479, r10460, r10478, MPFR_RNDN); mpfr_sub(r10480, r10478, r10460, MPFR_RNDN); mpfr_mul(r10481, r10479, r10480, MPFR_RNDN); mpfr_sqrt(r10482, r10481, MPFR_RNDN); mpfr_add(r10483, r10478, r10482, MPFR_RNDN); ; mpfr_pow(r10485, r10483, r10484, MPFR_RNDN); mpfr_cbrt(r10486, r10485, MPFR_RNDN); mpfr_mul(r10487, r10477, r10486, MPFR_RNDN); ; mpfr_set_si(r10489, mpfr_cmp(r10450, r10488) <= 0, MPFR_RNDN); mpfr_mul(r10490, r10456, r10460, MPFR_RNDN); mpfr_div(r10491, r10490, r10459, MPFR_RNDN); mpfr_mul(r10492, r10460, r10491, MPFR_RNDN); mpfr_div(r10493, r10456, r10466, MPFR_RNDN); mpfr_div(r10494, r10493, r10457, MPFR_RNDN); mpfr_mul(r10495, r10494, r10464, MPFR_RNDN); mpfr_mul(r10496, r10495, r10495, MPFR_RNDN); mpfr_sub(r10497, r10496, r10461, MPFR_RNDN); mpfr_sqrt(r10498, r10497, MPFR_RNDN); mpfr_sub(r10499, r10495, r10498, MPFR_RNDN); mpfr_div(r10500, r10492, r10499, MPFR_RNDN); ; mpfr_set_si(r10502, mpfr_cmp(r10450, r10501) <= 0, MPFR_RNDN); if (mpfr_get_si(r10502, MPFR_RNDN)) { mpfr_set(r10503, r10487, MPFR_RNDN); } else { mpfr_set(r10503, r10453, MPFR_RNDN); }; if (mpfr_get_si(r10489, MPFR_RNDN)) { mpfr_set(r10504, r10500, MPFR_RNDN); } else { mpfr_set(r10504, r10503, MPFR_RNDN); }; if (mpfr_get_si(r10475, MPFR_RNDN)) { mpfr_set(r10505, r10487, MPFR_RNDN); } else { mpfr_set(r10505, r10504, MPFR_RNDN); }; if (mpfr_get_si(r10455, MPFR_RNDN)) { mpfr_set(r10506, r10473, MPFR_RNDN); } else { mpfr_set(r10506, r10505, MPFR_RNDN); }; if (mpfr_get_si(r10452, MPFR_RNDN)) { mpfr_set(r10507, r10453, MPFR_RNDN); } else { mpfr_set(r10507, r10506, MPFR_RNDN); }; return mpfr_get_d(r10507, MPFR_RNDN); }