\frac{\left(\left(\left(\left(1 + 0.1049934947 \cdot \left(x \cdot x\right)\right) + 0.0424060604 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + 0.0072644182 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + 0.0005064034 \cdot \left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + 0.0001789971 \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)}{\left(\left(\left(\left(\left(1 + 0.7715471019 \cdot \left(x \cdot x\right)\right) + 0.2909738639 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + 0.0694555761 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + 0.0140005442 \cdot \left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + 0.0008327945 \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + \left(2 \cdot 0.0001789971\right) \cdot \left(\left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)} \cdot x\begin{array}{l}
\mathbf{if}\;x \le -610.1882030804575:\\
\;\;\;\;\frac{\frac{0.2514179000665375}{x}}{x \cdot x} + \left(\frac{0.15298196345929327}{\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot x\right)} + \frac{0.5}{x}\right)\\
\mathbf{elif}\;x \le 975.8488061330117:\\
\;\;\;\;x \cdot \sqrt[3]{\frac{\mathsf{fma}\left(x \cdot x, 0.0001789971 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right), \mathsf{fma}\left(0.0005064034, \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), \mathsf{fma}\left(0.0072644182 \cdot \left(x \cdot x\right), \left(x \cdot x\right) \cdot \left(x \cdot x\right), \mathsf{fma}\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right), 0.0424060604, \mathsf{fma}\left(0.1049934947, x \cdot x, 1\right)\right)\right)\right)\right)}{\mathsf{fma}\left(0.0140005442, \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), \mathsf{fma}\left(x \cdot x, 0.0008327945 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right), \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right) \cdot \left(0.0001789971 + 0.0001789971\right)\right) + \mathsf{fma}\left(x \cdot x, \mathsf{fma}\left(x \cdot x, 0.2909738639, \left(0.0694555761 \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right), \mathsf{fma}\left(x \cdot x, 0.7715471019, 1\right)\right)\right)} \cdot \left(\frac{\mathsf{fma}\left(x \cdot x, 0.0001789971 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right), \mathsf{fma}\left(0.0005064034, \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), \mathsf{fma}\left(0.0072644182 \cdot \left(x \cdot x\right), \left(x \cdot x\right) \cdot \left(x \cdot x\right), \mathsf{fma}\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right), 0.0424060604, \mathsf{fma}\left(0.1049934947, x \cdot x, 1\right)\right)\right)\right)\right)}{\mathsf{fma}\left(0.0140005442, \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), \mathsf{fma}\left(x \cdot x, 0.0008327945 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right), \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right) \cdot \left(0.0001789971 + 0.0001789971\right)\right) + \mathsf{fma}\left(x \cdot x, \mathsf{fma}\left(x \cdot x, 0.2909738639, \left(0.0694555761 \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right), \mathsf{fma}\left(x \cdot x, 0.7715471019, 1\right)\right)\right)} \cdot \frac{\mathsf{fma}\left(x \cdot x, 0.0001789971 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right), \mathsf{fma}\left(0.0005064034, \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), \mathsf{fma}\left(0.0072644182 \cdot \left(x \cdot x\right), \left(x \cdot x\right) \cdot \left(x \cdot x\right), \mathsf{fma}\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right), 0.0424060604, \mathsf{fma}\left(0.1049934947, x \cdot x, 1\right)\right)\right)\right)\right)}{\mathsf{fma}\left(0.0140005442, \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), \mathsf{fma}\left(x \cdot x, 0.0008327945 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right), \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right) \cdot \left(0.0001789971 + 0.0001789971\right)\right) + \mathsf{fma}\left(x \cdot x, \mathsf{fma}\left(x \cdot x, 0.2909738639, \left(0.0694555761 \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right), \mathsf{fma}\left(x \cdot x, 0.7715471019, 1\right)\right)\right)}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{0.2514179000665375}{x}}{x \cdot x} + \left(\frac{0.15298196345929327}{\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot x\right)} + \frac{0.5}{x}\right)\\
\end{array}double f(double x) {
double r5822055 = 1.0;
double r5822056 = 0.1049934947;
double r5822057 = x;
double r5822058 = r5822057 * r5822057;
double r5822059 = r5822056 * r5822058;
double r5822060 = r5822055 + r5822059;
double r5822061 = 0.0424060604;
double r5822062 = r5822058 * r5822058;
double r5822063 = r5822061 * r5822062;
double r5822064 = r5822060 + r5822063;
double r5822065 = 0.0072644182;
double r5822066 = r5822062 * r5822058;
double r5822067 = r5822065 * r5822066;
double r5822068 = r5822064 + r5822067;
double r5822069 = 0.0005064034;
double r5822070 = r5822066 * r5822058;
double r5822071 = r5822069 * r5822070;
double r5822072 = r5822068 + r5822071;
double r5822073 = 0.0001789971;
double r5822074 = r5822070 * r5822058;
double r5822075 = r5822073 * r5822074;
double r5822076 = r5822072 + r5822075;
double r5822077 = 0.7715471019;
double r5822078 = r5822077 * r5822058;
double r5822079 = r5822055 + r5822078;
double r5822080 = 0.2909738639;
double r5822081 = r5822080 * r5822062;
double r5822082 = r5822079 + r5822081;
double r5822083 = 0.0694555761;
double r5822084 = r5822083 * r5822066;
double r5822085 = r5822082 + r5822084;
double r5822086 = 0.0140005442;
double r5822087 = r5822086 * r5822070;
double r5822088 = r5822085 + r5822087;
double r5822089 = 0.0008327945;
double r5822090 = r5822089 * r5822074;
double r5822091 = r5822088 + r5822090;
double r5822092 = 2.0;
double r5822093 = r5822092 * r5822073;
double r5822094 = r5822074 * r5822058;
double r5822095 = r5822093 * r5822094;
double r5822096 = r5822091 + r5822095;
double r5822097 = r5822076 / r5822096;
double r5822098 = r5822097 * r5822057;
return r5822098;
}
double f(double x) {
double r5822099 = x;
double r5822100 = -610.1882030804575;
bool r5822101 = r5822099 <= r5822100;
double r5822102 = 0.2514179000665375;
double r5822103 = r5822102 / r5822099;
double r5822104 = r5822099 * r5822099;
double r5822105 = r5822103 / r5822104;
double r5822106 = 0.15298196345929327;
double r5822107 = r5822104 * r5822099;
double r5822108 = r5822104 * r5822107;
double r5822109 = r5822106 / r5822108;
double r5822110 = 0.5;
double r5822111 = r5822110 / r5822099;
double r5822112 = r5822109 + r5822111;
double r5822113 = r5822105 + r5822112;
double r5822114 = 975.8488061330117;
bool r5822115 = r5822099 <= r5822114;
double r5822116 = 0.0001789971;
double r5822117 = r5822104 * r5822104;
double r5822118 = r5822117 * r5822117;
double r5822119 = r5822116 * r5822118;
double r5822120 = 0.0005064034;
double r5822121 = 0.0072644182;
double r5822122 = r5822121 * r5822104;
double r5822123 = 0.0424060604;
double r5822124 = 0.1049934947;
double r5822125 = 1.0;
double r5822126 = fma(r5822124, r5822104, r5822125);
double r5822127 = fma(r5822117, r5822123, r5822126);
double r5822128 = fma(r5822122, r5822117, r5822127);
double r5822129 = fma(r5822120, r5822118, r5822128);
double r5822130 = fma(r5822104, r5822119, r5822129);
double r5822131 = 0.0140005442;
double r5822132 = 0.0008327945;
double r5822133 = r5822132 * r5822118;
double r5822134 = r5822117 * r5822118;
double r5822135 = r5822116 + r5822116;
double r5822136 = r5822134 * r5822135;
double r5822137 = fma(r5822104, r5822133, r5822136);
double r5822138 = 0.2909738639;
double r5822139 = 0.0694555761;
double r5822140 = r5822139 * r5822104;
double r5822141 = r5822140 * r5822104;
double r5822142 = fma(r5822104, r5822138, r5822141);
double r5822143 = 0.7715471019;
double r5822144 = fma(r5822104, r5822143, r5822125);
double r5822145 = fma(r5822104, r5822142, r5822144);
double r5822146 = r5822137 + r5822145;
double r5822147 = fma(r5822131, r5822118, r5822146);
double r5822148 = r5822130 / r5822147;
double r5822149 = r5822148 * r5822148;
double r5822150 = r5822148 * r5822149;
double r5822151 = cbrt(r5822150);
double r5822152 = r5822099 * r5822151;
double r5822153 = r5822115 ? r5822152 : r5822113;
double r5822154 = r5822101 ? r5822113 : r5822153;
return r5822154;
}



Bits error versus x
if x < -610.1882030804575 or 975.8488061330117 < x Initial program 57.8
Taylor expanded around inf 0.0
Simplified0.0
if -610.1882030804575 < x < 975.8488061330117Initial program 0.0
rmApplied add-cbrt-cube0.0
Applied add-cbrt-cube0.0
Applied cbrt-undiv0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2019158 +o rules:numerics
(FPCore (x)
:name "Jmat.Real.dawson"
(* (/ (+ (+ (+ (+ (+ 1 (* 0.1049934947 (* x x))) (* 0.0424060604 (* (* x x) (* x x)))) (* 0.0072644182 (* (* (* x x) (* x x)) (* x x)))) (* 0.0005064034 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))) (* 0.0001789971 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)))) (+ (+ (+ (+ (+ (+ 1 (* 0.7715471019 (* x x))) (* 0.2909738639 (* (* x x) (* x x)))) (* 0.0694555761 (* (* (* x x) (* x x)) (* x x)))) (* 0.0140005442 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))) (* 0.0008327945 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)))) (* (* 2 0.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x))))) x))