#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 r10107 = d; float r10108 = h; float r10109 = r10107 / r10108; float r10110 = 1.0f; float r10111 = 2.0f; float r10112 = r10110 / r10111; float r10113 = pow(r10109, r10112); float r10114 = l; float r10115 = r10107 / r10114; float r10116 = pow(r10115, r10112); float r10117 = r10113 * r10116; float r10118 = M; float r10119 = D; float r10120 = r10118 * r10119; float r10121 = r10111 * r10107; float r10122 = r10120 / r10121; float r10123 = pow(r10122, r10111); float r10124 = r10112 * r10123; float r10125 = r10108 / r10114; float r10126 = r10124 * r10125; float r10127 = r10110 - r10126; float r10128 = r10117 * r10127; return r10128; } double f_id(double d, double h, double l, double M, double D) { double r10129 = d; double r10130 = h; double r10131 = r10129 / r10130; double r10132 = 1.0; double r10133 = 2.0; double r10134 = r10132 / r10133; double r10135 = pow(r10131, r10134); double r10136 = l; double r10137 = r10129 / r10136; double r10138 = pow(r10137, r10134); double r10139 = r10135 * r10138; double r10140 = M; double r10141 = D; double r10142 = r10140 * r10141; double r10143 = r10133 * r10129; double r10144 = r10142 / r10143; double r10145 = pow(r10144, r10133); double r10146 = r10134 * r10145; double r10147 = r10130 / r10136; double r10148 = r10146 * r10147; double r10149 = r10132 - r10148; double r10150 = r10139 * r10149; return r10150; } double f_of(float d, float h, float l, float M, float D) { float r10151 = d; float r10152 = -2.278379854642029e-25f; bool r10153 = r10151 <= r10152; float r10154 = -1.0f; float r10155 = h; float r10156 = r10154 / r10155; float r10157 = sqrt(r10156); float r10158 = r10154 / r10151; float r10159 = 0.5f; float r10160 = pow(r10158, r10159); float r10161 = r10157 / r10160; float r10162 = l; float r10163 = r10151 / r10162; float r10164 = sqrt(r10163); float r10165 = 1.0f; float r10166 = 2.0f; float r10167 = r10165 / r10166; float r10168 = pow(r10164, r10167); float r10169 = cbrt(r10163); float r10170 = fabs(r10169); float r10171 = sqrt(r10169); float r10172 = r10170 * r10171; float r10173 = pow(r10172, r10167); float r10174 = r10168 * r10173; float r10175 = r10161 * r10174; float r10176 = M; float r10177 = D; float r10178 = r10176 * r10177; float r10179 = r10166 * r10151; float r10180 = r10178 / r10179; float r10181 = pow(r10180, r10166); float r10182 = r10167 * r10181; float r10183 = r10155 / r10162; float r10184 = r10182 * r10183; float r10185 = r10165 - r10184; float r10186 = r10175 * r10185; float r10187 = -1.09704788971517e-310f; bool r10188 = r10151 <= r10187; float r10189 = r10155 * r10181; float r10190 = r10162 * r10166; float r10191 = r10189 / r10190; float r10192 = r10165 - r10191; float r10193 = r10154 / r10162; float r10194 = sqrt(r10193); float r10195 = r10194 / r10160; float r10196 = r10151 / r10155; float r10197 = pow(r10196, r10167); float r10198 = r10195 * r10197; float r10199 = r10192 * r10198; float r10200 = 4.2030112357850945e-186f; bool r10201 = r10151 <= r10200; float r10202 = r10165 / r10162; float r10203 = pow(r10202, r10167); float r10204 = pow(r10151, r10167); float r10205 = r10203 * r10204; float r10206 = r10197 * r10205; float r10207 = r10206 * r10192; float r10208 = sqrt(r10151); float r10209 = pow(r10155, r10159); float r10210 = r10208 / r10209; float r10211 = pow(r10163, r10167); float r10212 = r10210 * r10211; float r10213 = r10192 * r10212; float r10214 = r10201 ? r10207 : r10213; float r10215 = r10188 ? r10199 : r10214; float r10216 = r10153 ? r10186 : r10215; return r10216; } double f_od(double d, double h, double l, double M, double D) { double r10217 = d; double r10218 = -2.278379854642029e-25; bool r10219 = r10217 <= r10218; double r10220 = -1.0; double r10221 = h; double r10222 = r10220 / r10221; double r10223 = sqrt(r10222); double r10224 = r10220 / r10217; double r10225 = 0.5; double r10226 = pow(r10224, r10225); double r10227 = r10223 / r10226; double r10228 = l; double r10229 = r10217 / r10228; double r10230 = sqrt(r10229); double r10231 = 1.0; double r10232 = 2.0; double r10233 = r10231 / r10232; double r10234 = pow(r10230, r10233); double r10235 = cbrt(r10229); double r10236 = fabs(r10235); double r10237 = sqrt(r10235); double r10238 = r10236 * r10237; double r10239 = pow(r10238, r10233); double r10240 = r10234 * r10239; double r10241 = r10227 * r10240; double r10242 = M; double r10243 = D; double r10244 = r10242 * r10243; double r10245 = r10232 * r10217; double r10246 = r10244 / r10245; double r10247 = pow(r10246, r10232); double r10248 = r10233 * r10247; double r10249 = r10221 / r10228; double r10250 = r10248 * r10249; double r10251 = r10231 - r10250; double r10252 = r10241 * r10251; double r10253 = -1.09704788971517e-310; bool r10254 = r10217 <= r10253; double r10255 = r10221 * r10247; double r10256 = r10228 * r10232; double r10257 = r10255 / r10256; double r10258 = r10231 - r10257; double r10259 = r10220 / r10228; double r10260 = sqrt(r10259); double r10261 = r10260 / r10226; double r10262 = r10217 / r10221; double r10263 = pow(r10262, r10233); double r10264 = r10261 * r10263; double r10265 = r10258 * r10264; double r10266 = 4.2030112357850945e-186; bool r10267 = r10217 <= r10266; double r10268 = r10231 / r10228; double r10269 = pow(r10268, r10233); double r10270 = pow(r10217, r10233); double r10271 = r10269 * r10270; double r10272 = r10263 * r10271; double r10273 = r10272 * r10258; double r10274 = sqrt(r10217); double r10275 = pow(r10221, r10225); double r10276 = r10274 / r10275; double r10277 = pow(r10229, r10233); double r10278 = r10276 * r10277; double r10279 = r10258 * r10278; double r10280 = r10267 ? r10273 : r10279; double r10281 = r10254 ? r10265 : r10280; double r10282 = r10219 ? r10252 : r10281; return r10282; } 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 r10283, r10284, r10285, r10286, r10287, r10288, r10289, r10290, r10291, r10292, r10293, r10294, r10295, r10296, r10297, r10298, r10299, r10300, r10301, r10302, r10303, r10304; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r10283); mpfr_init(r10284); mpfr_init(r10285); mpfr_init_set_str(r10286, "1", 10, MPFR_RNDN); mpfr_init_set_str(r10287, "2", 10, MPFR_RNDN); mpfr_init(r10288); mpfr_init(r10289); mpfr_init(r10290); mpfr_init(r10291); mpfr_init(r10292); mpfr_init(r10293); mpfr_init(r10294); mpfr_init(r10295); mpfr_init(r10296); mpfr_init(r10297); mpfr_init(r10298); mpfr_init(r10299); mpfr_init(r10300); mpfr_init(r10301); mpfr_init(r10302); mpfr_init(r10303); mpfr_init(r10304); } double f_im(double d, double h, double l, double M, double D) { mpfr_set_d(r10283, d, MPFR_RNDN); mpfr_set_d(r10284, h, MPFR_RNDN); mpfr_div(r10285, r10283, r10284, MPFR_RNDN); ; ; mpfr_div(r10288, r10286, r10287, MPFR_RNDN); mpfr_pow(r10289, r10285, r10288, MPFR_RNDN); mpfr_set_d(r10290, l, MPFR_RNDN); mpfr_div(r10291, r10283, r10290, MPFR_RNDN); mpfr_pow(r10292, r10291, r10288, MPFR_RNDN); mpfr_mul(r10293, r10289, r10292, MPFR_RNDN); mpfr_set_d(r10294, M, MPFR_RNDN); mpfr_set_d(r10295, D, MPFR_RNDN); mpfr_mul(r10296, r10294, r10295, MPFR_RNDN); mpfr_mul(r10297, r10287, r10283, MPFR_RNDN); mpfr_div(r10298, r10296, r10297, MPFR_RNDN); mpfr_pow(r10299, r10298, r10287, MPFR_RNDN); mpfr_mul(r10300, r10288, r10299, MPFR_RNDN); mpfr_div(r10301, r10284, r10290, MPFR_RNDN); mpfr_mul(r10302, r10300, r10301, MPFR_RNDN); mpfr_sub(r10303, r10286, r10302, MPFR_RNDN); mpfr_mul(r10304, r10293, r10303, MPFR_RNDN); return mpfr_get_d(r10304, MPFR_RNDN); } static mpfr_t r10305, r10306, r10307, r10308, r10309, r10310, r10311, r10312, r10313, r10314, r10315, r10316, r10317, r10318, r10319, r10320, r10321, r10322, r10323, r10324, r10325, r10326, r10327, r10328, r10329, r10330, r10331, r10332, r10333, r10334, r10335, r10336, r10337, r10338, r10339, r10340, r10341, r10342, r10343, r10344, r10345, r10346, r10347, r10348, r10349, r10350, r10351, r10352, r10353, r10354, r10355, r10356, r10357, r10358, r10359, r10360, r10361, r10362, r10363, r10364, r10365, r10366, r10367, r10368, r10369, r10370; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r10305); mpfr_init_set_str(r10306, "-2.278379854642029e-25", 10, MPFR_RNDN); mpfr_init(r10307); mpfr_init_set_str(r10308, "-1", 10, MPFR_RNDN); mpfr_init(r10309); mpfr_init(r10310); mpfr_init(r10311); mpfr_init(r10312); mpfr_init_set_str(r10313, "1/2", 10, MPFR_RNDN); mpfr_init(r10314); mpfr_init(r10315); mpfr_init(r10316); mpfr_init(r10317); mpfr_init(r10318); mpfr_init_set_str(r10319, "1", 10, MPFR_RNDN); mpfr_init_set_str(r10320, "2", 10, MPFR_RNDN); mpfr_init(r10321); mpfr_init(r10322); mpfr_init(r10323); mpfr_init(r10324); mpfr_init(r10325); mpfr_init(r10326); mpfr_init(r10327); mpfr_init(r10328); mpfr_init(r10329); mpfr_init(r10330); mpfr_init(r10331); mpfr_init(r10332); mpfr_init(r10333); mpfr_init(r10334); mpfr_init(r10335); mpfr_init(r10336); mpfr_init(r10337); mpfr_init(r10338); mpfr_init(r10339); mpfr_init(r10340); mpfr_init_set_str(r10341, "-1.09704788971517e-310", 10, MPFR_RNDN); mpfr_init(r10342); mpfr_init(r10343); mpfr_init(r10344); mpfr_init(r10345); mpfr_init(r10346); mpfr_init(r10347); mpfr_init(r10348); mpfr_init(r10349); mpfr_init(r10350); mpfr_init(r10351); mpfr_init(r10352); mpfr_init(r10353); mpfr_init_set_str(r10354, "4.2030112357850945e-186", 10, MPFR_RNDN); mpfr_init(r10355); mpfr_init(r10356); mpfr_init(r10357); mpfr_init(r10358); mpfr_init(r10359); mpfr_init(r10360); mpfr_init(r10361); mpfr_init(r10362); mpfr_init(r10363); mpfr_init(r10364); mpfr_init(r10365); mpfr_init(r10366); mpfr_init(r10367); mpfr_init(r10368); mpfr_init(r10369); mpfr_init(r10370); } double f_fm(double d, double h, double l, double M, double D) { mpfr_set_d(r10305, d, MPFR_RNDN); ; mpfr_set_si(r10307, mpfr_cmp(r10305, r10306) <= 0, MPFR_RNDN); ; mpfr_set_d(r10309, h, MPFR_RNDN); mpfr_div(r10310, r10308, r10309, MPFR_RNDN); mpfr_sqrt(r10311, r10310, MPFR_RNDN); mpfr_div(r10312, r10308, r10305, MPFR_RNDN); ; mpfr_pow(r10314, r10312, r10313, MPFR_RNDN); mpfr_div(r10315, r10311, r10314, MPFR_RNDN); mpfr_set_d(r10316, l, MPFR_RNDN); mpfr_div(r10317, r10305, r10316, MPFR_RNDN); mpfr_sqrt(r10318, r10317, MPFR_RNDN); ; ; mpfr_div(r10321, r10319, r10320, MPFR_RNDN); mpfr_pow(r10322, r10318, r10321, MPFR_RNDN); mpfr_cbrt(r10323, r10317, MPFR_RNDN); mpfr_abs(r10324, r10323, MPFR_RNDN); mpfr_sqrt(r10325, r10323, MPFR_RNDN); mpfr_mul(r10326, r10324, r10325, MPFR_RNDN); mpfr_pow(r10327, r10326, r10321, MPFR_RNDN); mpfr_mul(r10328, r10322, r10327, MPFR_RNDN); mpfr_mul(r10329, r10315, r10328, MPFR_RNDN); mpfr_set_d(r10330, M, MPFR_RNDN); mpfr_set_d(r10331, D, MPFR_RNDN); mpfr_mul(r10332, r10330, r10331, MPFR_RNDN); mpfr_mul(r10333, r10320, r10305, MPFR_RNDN); mpfr_div(r10334, r10332, r10333, MPFR_RNDN); mpfr_pow(r10335, r10334, r10320, MPFR_RNDN); mpfr_mul(r10336, r10321, r10335, MPFR_RNDN); mpfr_div(r10337, r10309, r10316, MPFR_RNDN); mpfr_mul(r10338, r10336, r10337, MPFR_RNDN); mpfr_sub(r10339, r10319, r10338, MPFR_RNDN); mpfr_mul(r10340, r10329, r10339, MPFR_RNDN); ; mpfr_set_si(r10342, mpfr_cmp(r10305, r10341) <= 0, MPFR_RNDN); mpfr_mul(r10343, r10309, r10335, MPFR_RNDN); mpfr_mul(r10344, r10316, r10320, MPFR_RNDN); mpfr_div(r10345, r10343, r10344, MPFR_RNDN); mpfr_sub(r10346, r10319, r10345, MPFR_RNDN); mpfr_div(r10347, r10308, r10316, MPFR_RNDN); mpfr_sqrt(r10348, r10347, MPFR_RNDN); mpfr_div(r10349, r10348, r10314, MPFR_RNDN); mpfr_div(r10350, r10305, r10309, MPFR_RNDN); mpfr_pow(r10351, r10350, r10321, MPFR_RNDN); mpfr_mul(r10352, r10349, r10351, MPFR_RNDN); mpfr_mul(r10353, r10346, r10352, MPFR_RNDN); ; mpfr_set_si(r10355, mpfr_cmp(r10305, r10354) <= 0, MPFR_RNDN); mpfr_div(r10356, r10319, r10316, MPFR_RNDN); mpfr_pow(r10357, r10356, r10321, MPFR_RNDN); mpfr_pow(r10358, r10305, r10321, MPFR_RNDN); mpfr_mul(r10359, r10357, r10358, MPFR_RNDN); mpfr_mul(r10360, r10351, r10359, MPFR_RNDN); mpfr_mul(r10361, r10360, r10346, MPFR_RNDN); mpfr_sqrt(r10362, r10305, MPFR_RNDN); mpfr_pow(r10363, r10309, r10313, MPFR_RNDN); mpfr_div(r10364, r10362, r10363, MPFR_RNDN); mpfr_pow(r10365, r10317, r10321, MPFR_RNDN); mpfr_mul(r10366, r10364, r10365, MPFR_RNDN); mpfr_mul(r10367, r10346, r10366, MPFR_RNDN); if (mpfr_get_si(r10355, MPFR_RNDN)) { mpfr_set(r10368, r10361, MPFR_RNDN); } else { mpfr_set(r10368, r10367, MPFR_RNDN); }; if (mpfr_get_si(r10342, MPFR_RNDN)) { mpfr_set(r10369, r10353, MPFR_RNDN); } else { mpfr_set(r10369, r10368, MPFR_RNDN); }; if (mpfr_get_si(r10307, MPFR_RNDN)) { mpfr_set(r10370, r10340, MPFR_RNDN); } else { mpfr_set(r10370, r10369, MPFR_RNDN); }; return mpfr_get_d(r10370, MPFR_RNDN); } static mpfr_t r10371, r10372, r10373, r10374, r10375, r10376, r10377, r10378, r10379, r10380, r10381, r10382, r10383, r10384, r10385, r10386, r10387, 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; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r10371); mpfr_init_set_str(r10372, "-2.278379854642029e-25", 10, MPFR_RNDN); mpfr_init(r10373); mpfr_init_set_str(r10374, "-1", 10, MPFR_RNDN); mpfr_init(r10375); mpfr_init(r10376); mpfr_init(r10377); mpfr_init(r10378); mpfr_init_set_str(r10379, "1/2", 10, MPFR_RNDN); mpfr_init(r10380); mpfr_init(r10381); mpfr_init(r10382); mpfr_init(r10383); mpfr_init(r10384); mpfr_init_set_str(r10385, "1", 10, MPFR_RNDN); mpfr_init_set_str(r10386, "2", 10, MPFR_RNDN); mpfr_init(r10387); mpfr_init(r10388); mpfr_init(r10389); mpfr_init(r10390); 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); mpfr_init(r10402); mpfr_init(r10403); mpfr_init(r10404); mpfr_init(r10405); mpfr_init(r10406); mpfr_init_set_str(r10407, "-1.09704788971517e-310", 10, MPFR_RNDN); mpfr_init(r10408); mpfr_init(r10409); mpfr_init(r10410); mpfr_init(r10411); mpfr_init(r10412); mpfr_init(r10413); mpfr_init(r10414); mpfr_init(r10415); mpfr_init(r10416); mpfr_init(r10417); mpfr_init(r10418); mpfr_init(r10419); mpfr_init_set_str(r10420, "4.2030112357850945e-186", 10, MPFR_RNDN); 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); mpfr_init(r10431); mpfr_init(r10432); mpfr_init(r10433); mpfr_init(r10434); mpfr_init(r10435); mpfr_init(r10436); } double f_dm(double d, double h, double l, double M, double D) { mpfr_set_d(r10371, d, MPFR_RNDN); ; mpfr_set_si(r10373, mpfr_cmp(r10371, r10372) <= 0, MPFR_RNDN); ; mpfr_set_d(r10375, h, MPFR_RNDN); mpfr_div(r10376, r10374, r10375, MPFR_RNDN); mpfr_sqrt(r10377, r10376, MPFR_RNDN); mpfr_div(r10378, r10374, r10371, MPFR_RNDN); ; mpfr_pow(r10380, r10378, r10379, MPFR_RNDN); mpfr_div(r10381, r10377, r10380, MPFR_RNDN); mpfr_set_d(r10382, l, MPFR_RNDN); mpfr_div(r10383, r10371, r10382, MPFR_RNDN); mpfr_sqrt(r10384, r10383, MPFR_RNDN); ; ; mpfr_div(r10387, r10385, r10386, MPFR_RNDN); mpfr_pow(r10388, r10384, r10387, MPFR_RNDN); mpfr_cbrt(r10389, r10383, MPFR_RNDN); mpfr_abs(r10390, r10389, MPFR_RNDN); mpfr_sqrt(r10391, r10389, MPFR_RNDN); mpfr_mul(r10392, r10390, r10391, MPFR_RNDN); mpfr_pow(r10393, r10392, r10387, MPFR_RNDN); mpfr_mul(r10394, r10388, r10393, MPFR_RNDN); mpfr_mul(r10395, r10381, r10394, MPFR_RNDN); mpfr_set_d(r10396, M, MPFR_RNDN); mpfr_set_d(r10397, D, MPFR_RNDN); mpfr_mul(r10398, r10396, r10397, MPFR_RNDN); mpfr_mul(r10399, r10386, r10371, MPFR_RNDN); mpfr_div(r10400, r10398, r10399, MPFR_RNDN); mpfr_pow(r10401, r10400, r10386, MPFR_RNDN); mpfr_mul(r10402, r10387, r10401, MPFR_RNDN); mpfr_div(r10403, r10375, r10382, MPFR_RNDN); mpfr_mul(r10404, r10402, r10403, MPFR_RNDN); mpfr_sub(r10405, r10385, r10404, MPFR_RNDN); mpfr_mul(r10406, r10395, r10405, MPFR_RNDN); ; mpfr_set_si(r10408, mpfr_cmp(r10371, r10407) <= 0, MPFR_RNDN); mpfr_mul(r10409, r10375, r10401, MPFR_RNDN); mpfr_mul(r10410, r10382, r10386, MPFR_RNDN); mpfr_div(r10411, r10409, r10410, MPFR_RNDN); mpfr_sub(r10412, r10385, r10411, MPFR_RNDN); mpfr_div(r10413, r10374, r10382, MPFR_RNDN); mpfr_sqrt(r10414, r10413, MPFR_RNDN); mpfr_div(r10415, r10414, r10380, MPFR_RNDN); mpfr_div(r10416, r10371, r10375, MPFR_RNDN); mpfr_pow(r10417, r10416, r10387, MPFR_RNDN); mpfr_mul(r10418, r10415, r10417, MPFR_RNDN); mpfr_mul(r10419, r10412, r10418, MPFR_RNDN); ; mpfr_set_si(r10421, mpfr_cmp(r10371, r10420) <= 0, MPFR_RNDN); mpfr_div(r10422, r10385, r10382, MPFR_RNDN); mpfr_pow(r10423, r10422, r10387, MPFR_RNDN); mpfr_pow(r10424, r10371, r10387, MPFR_RNDN); mpfr_mul(r10425, r10423, r10424, MPFR_RNDN); mpfr_mul(r10426, r10417, r10425, MPFR_RNDN); mpfr_mul(r10427, r10426, r10412, MPFR_RNDN); mpfr_sqrt(r10428, r10371, MPFR_RNDN); mpfr_pow(r10429, r10375, r10379, MPFR_RNDN); mpfr_div(r10430, r10428, r10429, MPFR_RNDN); mpfr_pow(r10431, r10383, r10387, MPFR_RNDN); mpfr_mul(r10432, r10430, r10431, MPFR_RNDN); mpfr_mul(r10433, r10412, r10432, MPFR_RNDN); if (mpfr_get_si(r10421, MPFR_RNDN)) { mpfr_set(r10434, r10427, MPFR_RNDN); } else { mpfr_set(r10434, r10433, MPFR_RNDN); }; if (mpfr_get_si(r10408, MPFR_RNDN)) { mpfr_set(r10435, r10419, MPFR_RNDN); } else { mpfr_set(r10435, r10434, MPFR_RNDN); }; if (mpfr_get_si(r10373, MPFR_RNDN)) { mpfr_set(r10436, r10406, MPFR_RNDN); } else { mpfr_set(r10436, r10435, MPFR_RNDN); }; return mpfr_get_d(r10436, MPFR_RNDN); }