Searched refs:P2 (Results 1 – 1 of 1) sorted by relevance
2504 CeedScalar P0, P1, P2, dP2, xi, wi, PI = 4.0 * atan(1.0); in CeedGaussQuadrature() local2513 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()[all …]