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