261 amrex::Real fact, half, one, res, sum, twelve, two, xbig, xden, xinf, xminin;
281 p[0] = -1.71618513886549492533811e+0;
282 p[1] = 2.47656508055759199108314e+1;
283 p[2] = -3.79804256470945635097577e+2;
284 p[3] = 6.29331155312818442661052e+2;
285 p[4] = 8.66966202790413211295064e+2;
286 p[5] = -3.14512729688483675254357e+4;
287 p[6] = -3.61444134186911729807069e+4;
288 p[7] = 6.64561438202405440627855e+4;
290 q[0] = -3.08402300119738975254353e+1;
291 q[1] = 3.15350626979604161529144e+2;
292 q[2] = -1.01515636749021914166146e+3;
293 q[3] = -3.10777167157231109440444e+3;
294 q[4] = 2.25381184209801510330112e+4;
295 q[5] = 4.75584627752788110767815e+3;
296 q[6] = -1.34659959864969306392456e+5;
297 q[7] = -1.15132259675553483497211e+5;
300 c[0] = -1.910444077728e-03;
301 c[1] = 8.4171387781295e-04;
302 c[2] = -5.952379913043012e-04;
303 c[3] = 7.93650793500350248e-04;
304 c[4] = -2.777777777777681622553e-03;
305 c[5] = 8.333333333333333331554247e-02;
306 c[6] = 5.7083835261e-03;
323 if (y1 != std::floor(y1 * half) * two)
325 Real pi=amrex::Math::pi<Real>();
326 fact = -
pi / std::sin(
pi * res);
352 else if (y < twelve) {
364 n =
static_cast<int>(
y) - 1;
374 for (i = 0; i < 8; i++) {
375 xnum = (xnum + p[i]) * z;
376 xden = xden *
z + q[i];
378 res = xnum / xden + one;
389 for (i = 0; i < n; i++) {
402 for (i = 0; i < 6; i++) {
403 sum = sum / ysq + c[i];
406 sum = sum /
y -
y +
xxx;
407 sum = sum + (
y - half) * std::log(y);
constexpr Real xxx
Definition: ERF_AdvanceMorrison.cpp:157
real(c_double), parameter, private pi
Definition: ERF_module_mp_morr_two_moment.F90:100