2504   CeedScalar P0, P1, P2, dP2, xi, wi, PI = 4.0 * atan(1.0);  in CeedGaussQuadrature()  local
2513     P2 = 0.0;  in CeedGaussQuadrature()
2515       P2 = (((CeedScalar)(2 * j - 1)) * xi * P1 - ((CeedScalar)(j - 1)) * P0) / ((CeedScalar)(j));  in CeedGaussQuadrature()
2517       P1 = P2;  in CeedGaussQuadrature()
2520     dP2 = (xi * P2 - P0) * (CeedScalar)Q / (xi * xi - 1.0);  in CeedGaussQuadrature()
2521     xi  = xi - P2 / dP2;  in CeedGaussQuadrature()
2523     for (CeedInt k = 0; k < 100 && fabs(P2) > 10 * CEED_EPSILON; k++) {  in CeedGaussQuadrature()
2527         P2 = (((CeedScalar)(2 * j - 1)) * xi * P1 - ((CeedScalar)(j - 1)) * P0) / ((CeedScalar)(j));  in CeedGaussQuadrature()
2529         P1 = P2;  in CeedGaussQuadrature()
2531       dP2 = (xi * P2 - P0) * (CeedScalar)Q / (xi * xi - 1.0);  in CeedGaussQuadrature()
2532       xi  = xi - P2 / dP2;  in CeedGaussQuadrature()
2556   CeedScalar P0, P1, P2, dP2, d2P2, xi, wi, PI = 4.0 * atan(1.0);  in CeedLobattoQuadrature()  local
2575     P2 = 0.0;  in CeedLobattoQuadrature()
2577       P2 = (((CeedScalar)(2 * j - 1)) * xi * P1 - ((CeedScalar)(j - 1)) * P0) / ((CeedScalar)(j));  in CeedLobattoQuadrature()
2579       P1 = P2;  in CeedLobattoQuadrature()
2582     dP2  = (xi * P2 - P0) * (CeedScalar)Q / (xi * xi - 1.0);  in CeedLobattoQuadrature()
2583     d2P2 = (2 * xi * dP2 - (CeedScalar)(Q * (Q - 1)) * P2) / (1.0 - xi * xi);  in CeedLobattoQuadrature()
2590         P2 = (((CeedScalar)(2 * j - 1)) * xi * P1 - ((CeedScalar)(j - 1)) * P0) / ((CeedScalar)(j));  in CeedLobattoQuadrature()
2592         P1 = P2;  in CeedLobattoQuadrature()
2594       dP2  = (xi * P2 - P0) * (CeedScalar)Q / (xi * xi - 1.0);  in CeedLobattoQuadrature()
2595       d2P2 = (2 * xi * dP2 - (CeedScalar)(Q * (Q - 1)) * P2) / (1.0 - xi * xi);  in CeedLobattoQuadrature()
2599     wi = 2.0 / (((CeedScalar)(Q * (Q - 1))) * P2 * P2);  in CeedLobattoQuadrature()