#include #include #include #include #include char *name = "Jmat.Real.erf"; double f_if(float x) { float r19077 = 1.0f; float r19078 = 0.3275911f; float r19079 = x; float r19080 = fabs(r19079); float r19081 = r19078 * r19080; float r19082 = r19077 + r19081; float r19083 = r19077 / r19082; float r19084 = 0.254829592f; float r19085 = -0.284496736f; float r19086 = 1.421413741f; float r19087 = -1.453152027f; float r19088 = 1.061405429f; float r19089 = r19083 * r19088; float r19090 = r19087 + r19089; float r19091 = r19083 * r19090; float r19092 = r19086 + r19091; float r19093 = r19083 * r19092; float r19094 = r19085 + r19093; float r19095 = r19083 * r19094; float r19096 = r19084 + r19095; float r19097 = r19083 * r19096; float r19098 = r19080 * r19080; float r19099 = -r19098; float r19100 = exp(r19099); float r19101 = r19097 * r19100; float r19102 = r19077 - r19101; return r19102; } double f_id(double x) { double r19103 = 1.0; double r19104 = 0.3275911; double r19105 = x; double r19106 = fabs(r19105); double r19107 = r19104 * r19106; double r19108 = r19103 + r19107; double r19109 = r19103 / r19108; double r19110 = 0.254829592; double r19111 = -0.284496736; double r19112 = 1.421413741; double r19113 = -1.453152027; double r19114 = 1.061405429; double r19115 = r19109 * r19114; double r19116 = r19113 + r19115; double r19117 = r19109 * r19116; double r19118 = r19112 + r19117; double r19119 = r19109 * r19118; double r19120 = r19111 + r19119; double r19121 = r19109 * r19120; double r19122 = r19110 + r19121; double r19123 = r19109 * r19122; double r19124 = r19106 * r19106; double r19125 = -r19124; double r19126 = exp(r19125); double r19127 = r19123 * r19126; double r19128 = r19103 - r19127; return r19128; } double f_of(float x) { float r19129 = x; float r19130 = -3.878602931916025e-24f; bool r19131 = r19129 <= r19130; float r19132 = 1.0f; float r19133 = 1.453152027f; float r19134 = fabs(r19129); float r19135 = r19134 * r19134; float r19136 = -r19135; float r19137 = exp(r19136); float r19138 = 0.3275911f; float r19139 = r19138 * r19134; float r19140 = r19132 + r19139; float r19141 = 4.0f; float r19142 = pow(r19140, r19141); float r19143 = r19137 / r19142; float r19144 = r19133 * r19143; float r19145 = 0.284496736f; float r19146 = r19140 * r19140; float r19147 = r19137 / r19146; float r19148 = r19145 * r19147; float r19149 = r19144 + r19148; float r19150 = r19132 + r19149; float r19151 = 3.0f; float r19152 = pow(r19150, r19151); float r19153 = 0.254829592f; float r19154 = r19137 / r19140; float r19155 = r19153 * r19154; float r19156 = 1.061405429f; float r19157 = 5.0f; float r19158 = pow(r19140, r19157); float r19159 = cbrt(r19158); float r19160 = r19159 * (r19159 * r19159); float r19161 = r19137 / r19160; float r19162 = r19156 * r19161; float r19163 = 1.421413741f; float r19164 = pow(r19140, r19151); float r19165 = r19137 / r19164; float r19166 = r19163 * r19165; float r19167 = r19162 + r19166; float r19168 = r19155 + r19167; float r19169 = pow(r19168, r19151); float r19170 = r19152 - r19169; float r19171 = r19150 * r19150; float r19172 = r19168 * r19168; float r19173 = r19150 * r19168; float r19174 = r19172 + r19173; float r19175 = r19171 + r19174; float r19176 = r19170 / r19175; float r19177 = 1.1454841341683684e-10f; bool r19178 = r19129 <= r19177; float r19179 = r19132 * r19132; float r19180 = r19149 * r19149; float r19181 = r19179 - r19180; float r19182 = r19137 / r19158; float r19183 = r19156 * r19182; float r19184 = r19183 + r19166; float r19185 = r19155 - r19184; float r19186 = r19181 * r19185; float r19187 = r19132 - r19149; float r19188 = r19155 * r19155; float r19189 = r19184 * r19184; float r19190 = r19188 - r19189; float r19191 = r19187 * r19190; float r19192 = r19186 - r19191; float r19193 = r19187 * r19185; float r19194 = r19192 / r19193; float r19195 = r19132 / r19140; float r19196 = -0.284496736f; float r19197 = -1.453152027f; float r19198 = sqrt(r19195); float r19199 = r19198 * r19198; float r19200 = r19199 * r19156; float r19201 = r19197 + r19200; float r19202 = r19195 * r19201; float r19203 = r19163 + r19202; float r19204 = r19195 * r19203; float r19205 = r19196 + r19204; float r19206 = r19195 * r19205; float r19207 = r19153 + r19206; float r19208 = r19195 * r19207; float r19209 = r19134 * r19134; float r19210 = -r19209; float r19211 = exp(r19210); float r19212 = r19208 * r19211; float r19213 = r19132 - r19212; float r19214 = r19178 ? r19194 : r19213; float r19215 = r19131 ? r19176 : r19214; return r19215; } double f_od(double x) { double r19216 = x; double r19217 = -3.878602931916025e-24; bool r19218 = r19216 <= r19217; double r19219 = 1.0; double r19220 = 1.453152027; double r19221 = fabs(r19216); double r19222 = r19221 * r19221; double r19223 = -r19222; double r19224 = exp(r19223); double r19225 = 0.3275911; double r19226 = r19225 * r19221; double r19227 = r19219 + r19226; double r19228 = 4.0; double r19229 = pow(r19227, r19228); double r19230 = r19224 / r19229; double r19231 = r19220 * r19230; double r19232 = 0.284496736; double r19233 = r19227 * r19227; double r19234 = r19224 / r19233; double r19235 = r19232 * r19234; double r19236 = r19231 + r19235; double r19237 = r19219 + r19236; double r19238 = 3.0; double r19239 = pow(r19237, r19238); double r19240 = 0.254829592; double r19241 = r19224 / r19227; double r19242 = r19240 * r19241; double r19243 = 1.061405429; double r19244 = 5.0; double r19245 = pow(r19227, r19244); double r19246 = cbrt(r19245); double r19247 = r19246 * (r19246 * r19246); double r19248 = r19224 / r19247; double r19249 = r19243 * r19248; double r19250 = 1.421413741; double r19251 = pow(r19227, r19238); double r19252 = r19224 / r19251; double r19253 = r19250 * r19252; double r19254 = r19249 + r19253; double r19255 = r19242 + r19254; double r19256 = pow(r19255, r19238); double r19257 = r19239 - r19256; double r19258 = r19237 * r19237; double r19259 = r19255 * r19255; double r19260 = r19237 * r19255; double r19261 = r19259 + r19260; double r19262 = r19258 + r19261; double r19263 = r19257 / r19262; double r19264 = 1.1454841341683684e-10; bool r19265 = r19216 <= r19264; double r19266 = r19219 * r19219; double r19267 = r19236 * r19236; double r19268 = r19266 - r19267; double r19269 = r19224 / r19245; double r19270 = r19243 * r19269; double r19271 = r19270 + r19253; double r19272 = r19242 - r19271; double r19273 = r19268 * r19272; double r19274 = r19219 - r19236; double r19275 = r19242 * r19242; double r19276 = r19271 * r19271; double r19277 = r19275 - r19276; double r19278 = r19274 * r19277; double r19279 = r19273 - r19278; double r19280 = r19274 * r19272; double r19281 = r19279 / r19280; double r19282 = r19219 / r19227; double r19283 = -0.284496736; double r19284 = -1.453152027; double r19285 = sqrt(r19282); double r19286 = r19285 * r19285; double r19287 = r19286 * r19243; double r19288 = r19284 + r19287; double r19289 = r19282 * r19288; double r19290 = r19250 + r19289; double r19291 = r19282 * r19290; double r19292 = r19283 + r19291; double r19293 = r19282 * r19292; double r19294 = r19240 + r19293; double r19295 = r19282 * r19294; double r19296 = r19221 * r19221; double r19297 = -r19296; double r19298 = exp(r19297); double r19299 = r19295 * r19298; double r19300 = r19219 - r19299; double r19301 = r19265 ? r19281 : r19300; double r19302 = r19218 ? r19263 : r19301; return r19302; } 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 r19303, r19304, r19305, r19306, r19307, r19308, r19309, r19310, r19311, r19312, r19313, r19314, r19315, r19316, r19317, r19318, r19319, r19320, r19321, r19322, r19323, r19324, r19325, r19326, r19327, r19328; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19303, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19304, "0.3275911", 10, MPFR_RNDN); mpfr_init(r19305); mpfr_init(r19306); mpfr_init(r19307); mpfr_init(r19308); mpfr_init(r19309); mpfr_init_set_str(r19310, "0.254829592", 10, MPFR_RNDN); mpfr_init_set_str(r19311, "-0.284496736", 10, MPFR_RNDN); mpfr_init_set_str(r19312, "1.421413741", 10, MPFR_RNDN); mpfr_init_set_str(r19313, "-1.453152027", 10, MPFR_RNDN); mpfr_init_set_str(r19314, "1.061405429", 10, MPFR_RNDN); mpfr_init(r19315); mpfr_init(r19316); mpfr_init(r19317); mpfr_init(r19318); mpfr_init(r19319); mpfr_init(r19320); mpfr_init(r19321); mpfr_init(r19322); mpfr_init(r19323); mpfr_init(r19324); mpfr_init(r19325); mpfr_init(r19326); mpfr_init(r19327); mpfr_init(r19328); } double f_im(double x) { ; ; mpfr_set_d(r19305, x, MPFR_RNDN); mpfr_abs(r19306, r19305, MPFR_RNDN); mpfr_mul(r19307, r19304, r19306, MPFR_RNDN); mpfr_add(r19308, r19303, r19307, MPFR_RNDN); mpfr_div(r19309, r19303, r19308, MPFR_RNDN); ; ; ; ; ; mpfr_mul(r19315, r19309, r19314, MPFR_RNDN); mpfr_add(r19316, r19313, r19315, MPFR_RNDN); mpfr_mul(r19317, r19309, r19316, MPFR_RNDN); mpfr_add(r19318, r19312, r19317, MPFR_RNDN); mpfr_mul(r19319, r19309, r19318, MPFR_RNDN); mpfr_add(r19320, r19311, r19319, MPFR_RNDN); mpfr_mul(r19321, r19309, r19320, MPFR_RNDN); mpfr_add(r19322, r19310, r19321, MPFR_RNDN); mpfr_mul(r19323, r19309, r19322, MPFR_RNDN); mpfr_mul(r19324, r19306, r19306, MPFR_RNDN); mpfr_neg(r19325, r19324, MPFR_RNDN); mpfr_exp(r19326, r19325, MPFR_RNDN); mpfr_mul(r19327, r19323, r19326, MPFR_RNDN); mpfr_sub(r19328, r19303, r19327, MPFR_RNDN); return mpfr_get_d(r19328, MPFR_RNDN); } static mpfr_t r19329, r19330, r19331, r19332, r19333, r19334, r19335, r19336, r19337, r19338, r19339, r19340, r19341, r19342, r19343, r19344, r19345, r19346, r19347, r19348, r19349, r19350, r19351, r19352, r19353, r19354, r19355, r19356, r19357, r19358, r19359, r19360, r19361, r19362, r19363, r19364, r19365, r19366, r19367, r19368, r19369, r19370, r19371, r19372, r19373, r19374, r19375, r19376, r19377, r19378, r19379, r19380, r19381, r19382, r19383, r19384, r19385, r19386, r19387, r19388, r19389, r19390, r19391, r19392, r19393, r19394, r19395, r19396, r19397, r19398, r19399, r19400, r19401, r19402, r19403, r19404, r19405, r19406, r19407, r19408, r19409, r19410, r19411, r19412, r19413, r19414, r19415; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19329); mpfr_init_set_str(r19330, "-3.878602931916025e-24", 10, MPFR_RNDN); mpfr_init(r19331); mpfr_init_set_str(r19332, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19333, "1.453152027", 10, MPFR_RNDN); mpfr_init(r19334); mpfr_init(r19335); mpfr_init(r19336); mpfr_init(r19337); mpfr_init_set_str(r19338, "0.3275911", 10, MPFR_RNDN); mpfr_init(r19339); mpfr_init(r19340); mpfr_init_set_str(r19341, "4", 10, MPFR_RNDN); mpfr_init(r19342); mpfr_init(r19343); mpfr_init(r19344); mpfr_init_set_str(r19345, "0.284496736", 10, MPFR_RNDN); mpfr_init(r19346); mpfr_init(r19347); mpfr_init(r19348); mpfr_init(r19349); mpfr_init(r19350); mpfr_init_set_str(r19351, "3", 10, MPFR_RNDN); mpfr_init(r19352); mpfr_init_set_str(r19353, "0.254829592", 10, MPFR_RNDN); mpfr_init(r19354); mpfr_init(r19355); mpfr_init_set_str(r19356, "1.061405429", 10, MPFR_RNDN); mpfr_init_set_str(r19357, "5", 10, MPFR_RNDN); mpfr_init(r19358); mpfr_init(r19359); mpfr_init(r19360); mpfr_init(r19361); mpfr_init(r19362); mpfr_init_set_str(r19363, "1.421413741", 10, MPFR_RNDN); mpfr_init(r19364); mpfr_init(r19365); mpfr_init(r19366); mpfr_init(r19367); mpfr_init(r19368); mpfr_init(r19369); mpfr_init(r19370); mpfr_init(r19371); mpfr_init(r19372); mpfr_init(r19373); mpfr_init(r19374); mpfr_init(r19375); mpfr_init(r19376); mpfr_init_set_str(r19377, "1.1454841341683684e-10", 10, MPFR_RNDN); mpfr_init(r19378); mpfr_init(r19379); mpfr_init(r19380); mpfr_init(r19381); mpfr_init(r19382); mpfr_init(r19383); mpfr_init(r19384); mpfr_init(r19385); mpfr_init(r19386); mpfr_init(r19387); mpfr_init(r19388); mpfr_init(r19389); mpfr_init(r19390); mpfr_init(r19391); mpfr_init(r19392); mpfr_init(r19393); mpfr_init(r19394); mpfr_init(r19395); mpfr_init_set_str(r19396, "-0.284496736", 10, MPFR_RNDN); mpfr_init_set_str(r19397, "-1.453152027", 10, MPFR_RNDN); mpfr_init(r19398); mpfr_init(r19399); mpfr_init(r19400); mpfr_init(r19401); mpfr_init(r19402); mpfr_init(r19403); mpfr_init(r19404); mpfr_init(r19405); mpfr_init(r19406); mpfr_init(r19407); mpfr_init(r19408); mpfr_init(r19409); mpfr_init(r19410); mpfr_init(r19411); mpfr_init(r19412); mpfr_init(r19413); mpfr_init(r19414); mpfr_init(r19415); } double f_fm(double x) { mpfr_set_d(r19329, x, MPFR_RNDN); ; mpfr_set_si(r19331, mpfr_cmp(r19329, r19330) <= 0, MPFR_RNDN); ; ; mpfr_abs(r19334, r19329, MPFR_RNDN); mpfr_sqr(r19335, r19334, MPFR_RNDN); mpfr_neg(r19336, r19335, MPFR_RNDN); mpfr_exp(r19337, r19336, MPFR_RNDN); ; mpfr_mul(r19339, r19338, r19334, MPFR_RNDN); mpfr_add(r19340, r19332, r19339, MPFR_RNDN); ; mpfr_pow(r19342, r19340, r19341, MPFR_RNDN); mpfr_div(r19343, r19337, r19342, MPFR_RNDN); mpfr_mul(r19344, r19333, r19343, MPFR_RNDN); ; mpfr_sqr(r19346, r19340, MPFR_RNDN); mpfr_div(r19347, r19337, r19346, MPFR_RNDN); mpfr_mul(r19348, r19345, r19347, MPFR_RNDN); mpfr_add(r19349, r19344, r19348, MPFR_RNDN); mpfr_add(r19350, r19332, r19349, MPFR_RNDN); ; mpfr_pow(r19352, r19350, r19351, MPFR_RNDN); ; mpfr_div(r19354, r19337, r19340, MPFR_RNDN); mpfr_mul(r19355, r19353, r19354, MPFR_RNDN); ; ; mpfr_pow(r19358, r19340, r19357, MPFR_RNDN); mpfr_cbrt(r19359, r19358, MPFR_RNDN); mpfr_mul(r19360, r19359, r19359, MPFR_RNDN); mpfr_mul(r19360, r19360, r19359, MPFR_RNDN); mpfr_div(r19361, r19337, r19360, MPFR_RNDN); mpfr_mul(r19362, r19356, r19361, MPFR_RNDN); ; mpfr_pow(r19364, r19340, r19351, MPFR_RNDN); mpfr_div(r19365, r19337, r19364, MPFR_RNDN); mpfr_mul(r19366, r19363, r19365, MPFR_RNDN); mpfr_add(r19367, r19362, r19366, MPFR_RNDN); mpfr_add(r19368, r19355, r19367, MPFR_RNDN); mpfr_pow(r19369, r19368, r19351, MPFR_RNDN); mpfr_sub(r19370, r19352, r19369, MPFR_RNDN); mpfr_sqr(r19371, r19350, MPFR_RNDN); mpfr_sqr(r19372, r19368, MPFR_RNDN); mpfr_mul(r19373, r19350, r19368, MPFR_RNDN); mpfr_add(r19374, r19372, r19373, MPFR_RNDN); mpfr_add(r19375, r19371, r19374, MPFR_RNDN); mpfr_div(r19376, r19370, r19375, MPFR_RNDN); ; mpfr_set_si(r19378, mpfr_cmp(r19329, r19377) <= 0, MPFR_RNDN); mpfr_sqr(r19379, r19332, MPFR_RNDN); mpfr_sqr(r19380, r19349, MPFR_RNDN); mpfr_sub(r19381, r19379, r19380, MPFR_RNDN); mpfr_div(r19382, r19337, r19358, MPFR_RNDN); mpfr_mul(r19383, r19356, r19382, MPFR_RNDN); mpfr_add(r19384, r19383, r19366, MPFR_RNDN); mpfr_sub(r19385, r19355, r19384, MPFR_RNDN); mpfr_mul(r19386, r19381, r19385, MPFR_RNDN); mpfr_sub(r19387, r19332, r19349, MPFR_RNDN); mpfr_sqr(r19388, r19355, MPFR_RNDN); mpfr_sqr(r19389, r19384, MPFR_RNDN); mpfr_sub(r19390, r19388, r19389, MPFR_RNDN); mpfr_mul(r19391, r19387, r19390, MPFR_RNDN); mpfr_sub(r19392, r19386, r19391, MPFR_RNDN); mpfr_mul(r19393, r19387, r19385, MPFR_RNDN); mpfr_div(r19394, r19392, r19393, MPFR_RNDN); mpfr_div(r19395, r19332, r19340, MPFR_RNDN); ; ; mpfr_sqrt(r19398, r19395, MPFR_RNDN); mpfr_sqr(r19399, r19398, MPFR_RNDN); mpfr_mul(r19400, r19399, r19356, MPFR_RNDN); mpfr_add(r19401, r19397, r19400, MPFR_RNDN); mpfr_mul(r19402, r19395, r19401, MPFR_RNDN); mpfr_add(r19403, r19363, r19402, MPFR_RNDN); mpfr_mul(r19404, r19395, r19403, MPFR_RNDN); mpfr_add(r19405, r19396, r19404, MPFR_RNDN); mpfr_mul(r19406, r19395, r19405, MPFR_RNDN); mpfr_add(r19407, r19353, r19406, MPFR_RNDN); mpfr_mul(r19408, r19395, r19407, MPFR_RNDN); mpfr_mul(r19409, r19334, r19334, MPFR_RNDN); mpfr_neg(r19410, r19409, MPFR_RNDN); mpfr_exp(r19411, r19410, MPFR_RNDN); mpfr_mul(r19412, r19408, r19411, MPFR_RNDN); mpfr_sub(r19413, r19332, r19412, MPFR_RNDN); if (mpfr_get_si(r19378, MPFR_RNDN)) { mpfr_set(r19414, r19394, MPFR_RNDN); } else { mpfr_set(r19414, r19413, MPFR_RNDN); }; if (mpfr_get_si(r19331, MPFR_RNDN)) { mpfr_set(r19415, r19376, MPFR_RNDN); } else { mpfr_set(r19415, r19414, MPFR_RNDN); }; return mpfr_get_d(r19415, MPFR_RNDN); } static mpfr_t r19416, r19417, r19418, r19419, r19420, r19421, r19422, r19423, r19424, r19425, r19426, r19427, r19428, r19429, r19430, r19431, r19432, r19433, r19434, r19435, r19436, r19437, r19438, r19439, r19440, r19441, r19442, r19443, r19444, r19445, r19446, r19447, r19448, r19449, r19450, r19451, r19452, r19453, r19454, r19455, r19456, r19457, r19458, r19459, r19460, r19461, r19462, r19463, r19464, r19465, r19466, r19467, r19468, r19469, r19470, r19471, r19472, r19473, r19474, r19475, r19476, r19477, r19478, r19479, r19480, r19481, r19482, r19483, r19484, r19485, r19486, r19487, r19488, r19489, r19490, r19491, r19492, r19493, r19494, r19495, r19496, r19497, r19498, r19499, r19500, r19501, r19502; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19416); mpfr_init_set_str(r19417, "-3.878602931916025e-24", 10, MPFR_RNDN); mpfr_init(r19418); mpfr_init_set_str(r19419, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19420, "1.453152027", 10, MPFR_RNDN); mpfr_init(r19421); mpfr_init(r19422); mpfr_init(r19423); mpfr_init(r19424); mpfr_init_set_str(r19425, "0.3275911", 10, MPFR_RNDN); mpfr_init(r19426); mpfr_init(r19427); mpfr_init_set_str(r19428, "4", 10, MPFR_RNDN); mpfr_init(r19429); mpfr_init(r19430); mpfr_init(r19431); mpfr_init_set_str(r19432, "0.284496736", 10, MPFR_RNDN); mpfr_init(r19433); mpfr_init(r19434); mpfr_init(r19435); mpfr_init(r19436); mpfr_init(r19437); mpfr_init_set_str(r19438, "3", 10, MPFR_RNDN); mpfr_init(r19439); mpfr_init_set_str(r19440, "0.254829592", 10, MPFR_RNDN); mpfr_init(r19441); mpfr_init(r19442); mpfr_init_set_str(r19443, "1.061405429", 10, MPFR_RNDN); mpfr_init_set_str(r19444, "5", 10, MPFR_RNDN); mpfr_init(r19445); mpfr_init(r19446); mpfr_init(r19447); mpfr_init(r19448); mpfr_init(r19449); mpfr_init_set_str(r19450, "1.421413741", 10, MPFR_RNDN); mpfr_init(r19451); mpfr_init(r19452); mpfr_init(r19453); mpfr_init(r19454); mpfr_init(r19455); mpfr_init(r19456); mpfr_init(r19457); mpfr_init(r19458); mpfr_init(r19459); mpfr_init(r19460); mpfr_init(r19461); mpfr_init(r19462); mpfr_init(r19463); mpfr_init_set_str(r19464, "1.1454841341683684e-10", 10, MPFR_RNDN); mpfr_init(r19465); mpfr_init(r19466); mpfr_init(r19467); mpfr_init(r19468); mpfr_init(r19469); mpfr_init(r19470); mpfr_init(r19471); mpfr_init(r19472); mpfr_init(r19473); mpfr_init(r19474); mpfr_init(r19475); mpfr_init(r19476); mpfr_init(r19477); mpfr_init(r19478); mpfr_init(r19479); mpfr_init(r19480); mpfr_init(r19481); mpfr_init(r19482); mpfr_init_set_str(r19483, "-0.284496736", 10, MPFR_RNDN); mpfr_init_set_str(r19484, "-1.453152027", 10, MPFR_RNDN); mpfr_init(r19485); mpfr_init(r19486); mpfr_init(r19487); mpfr_init(r19488); mpfr_init(r19489); mpfr_init(r19490); mpfr_init(r19491); mpfr_init(r19492); mpfr_init(r19493); mpfr_init(r19494); mpfr_init(r19495); mpfr_init(r19496); mpfr_init(r19497); mpfr_init(r19498); mpfr_init(r19499); mpfr_init(r19500); mpfr_init(r19501); mpfr_init(r19502); } double f_dm(double x) { mpfr_set_d(r19416, x, MPFR_RNDN); ; mpfr_set_si(r19418, mpfr_cmp(r19416, r19417) <= 0, MPFR_RNDN); ; ; mpfr_abs(r19421, r19416, MPFR_RNDN); mpfr_sqr(r19422, r19421, MPFR_RNDN); mpfr_neg(r19423, r19422, MPFR_RNDN); mpfr_exp(r19424, r19423, MPFR_RNDN); ; mpfr_mul(r19426, r19425, r19421, MPFR_RNDN); mpfr_add(r19427, r19419, r19426, MPFR_RNDN); ; mpfr_pow(r19429, r19427, r19428, MPFR_RNDN); mpfr_div(r19430, r19424, r19429, MPFR_RNDN); mpfr_mul(r19431, r19420, r19430, MPFR_RNDN); ; mpfr_sqr(r19433, r19427, MPFR_RNDN); mpfr_div(r19434, r19424, r19433, MPFR_RNDN); mpfr_mul(r19435, r19432, r19434, MPFR_RNDN); mpfr_add(r19436, r19431, r19435, MPFR_RNDN); mpfr_add(r19437, r19419, r19436, MPFR_RNDN); ; mpfr_pow(r19439, r19437, r19438, MPFR_RNDN); ; mpfr_div(r19441, r19424, r19427, MPFR_RNDN); mpfr_mul(r19442, r19440, r19441, MPFR_RNDN); ; ; mpfr_pow(r19445, r19427, r19444, MPFR_RNDN); mpfr_cbrt(r19446, r19445, MPFR_RNDN); mpfr_mul(r19447, r19446, r19446, MPFR_RNDN); mpfr_mul(r19447, r19447, r19446, MPFR_RNDN); mpfr_div(r19448, r19424, r19447, MPFR_RNDN); mpfr_mul(r19449, r19443, r19448, MPFR_RNDN); ; mpfr_pow(r19451, r19427, r19438, MPFR_RNDN); mpfr_div(r19452, r19424, r19451, MPFR_RNDN); mpfr_mul(r19453, r19450, r19452, MPFR_RNDN); mpfr_add(r19454, r19449, r19453, MPFR_RNDN); mpfr_add(r19455, r19442, r19454, MPFR_RNDN); mpfr_pow(r19456, r19455, r19438, MPFR_RNDN); mpfr_sub(r19457, r19439, r19456, MPFR_RNDN); mpfr_sqr(r19458, r19437, MPFR_RNDN); mpfr_sqr(r19459, r19455, MPFR_RNDN); mpfr_mul(r19460, r19437, r19455, MPFR_RNDN); mpfr_add(r19461, r19459, r19460, MPFR_RNDN); mpfr_add(r19462, r19458, r19461, MPFR_RNDN); mpfr_div(r19463, r19457, r19462, MPFR_RNDN); ; mpfr_set_si(r19465, mpfr_cmp(r19416, r19464) <= 0, MPFR_RNDN); mpfr_sqr(r19466, r19419, MPFR_RNDN); mpfr_sqr(r19467, r19436, MPFR_RNDN); mpfr_sub(r19468, r19466, r19467, MPFR_RNDN); mpfr_div(r19469, r19424, r19445, MPFR_RNDN); mpfr_mul(r19470, r19443, r19469, MPFR_RNDN); mpfr_add(r19471, r19470, r19453, MPFR_RNDN); mpfr_sub(r19472, r19442, r19471, MPFR_RNDN); mpfr_mul(r19473, r19468, r19472, MPFR_RNDN); mpfr_sub(r19474, r19419, r19436, MPFR_RNDN); mpfr_sqr(r19475, r19442, MPFR_RNDN); mpfr_sqr(r19476, r19471, MPFR_RNDN); mpfr_sub(r19477, r19475, r19476, MPFR_RNDN); mpfr_mul(r19478, r19474, r19477, MPFR_RNDN); mpfr_sub(r19479, r19473, r19478, MPFR_RNDN); mpfr_mul(r19480, r19474, r19472, MPFR_RNDN); mpfr_div(r19481, r19479, r19480, MPFR_RNDN); mpfr_div(r19482, r19419, r19427, MPFR_RNDN); ; ; mpfr_sqrt(r19485, r19482, MPFR_RNDN); mpfr_sqr(r19486, r19485, MPFR_RNDN); mpfr_mul(r19487, r19486, r19443, MPFR_RNDN); mpfr_add(r19488, r19484, r19487, MPFR_RNDN); mpfr_mul(r19489, r19482, r19488, MPFR_RNDN); mpfr_add(r19490, r19450, r19489, MPFR_RNDN); mpfr_mul(r19491, r19482, r19490, MPFR_RNDN); mpfr_add(r19492, r19483, r19491, MPFR_RNDN); mpfr_mul(r19493, r19482, r19492, MPFR_RNDN); mpfr_add(r19494, r19440, r19493, MPFR_RNDN); mpfr_mul(r19495, r19482, r19494, MPFR_RNDN); mpfr_mul(r19496, r19421, r19421, MPFR_RNDN); mpfr_neg(r19497, r19496, MPFR_RNDN); mpfr_exp(r19498, r19497, MPFR_RNDN); mpfr_mul(r19499, r19495, r19498, MPFR_RNDN); mpfr_sub(r19500, r19419, r19499, MPFR_RNDN); if (mpfr_get_si(r19465, MPFR_RNDN)) { mpfr_set(r19501, r19481, MPFR_RNDN); } else { mpfr_set(r19501, r19500, MPFR_RNDN); }; if (mpfr_get_si(r19418, MPFR_RNDN)) { mpfr_set(r19502, r19463, MPFR_RNDN); } else { mpfr_set(r19502, r19501, MPFR_RNDN); }; return mpfr_get_d(r19502, MPFR_RNDN); }