Lines Matching refs:xi
2504 CeedScalar P0, P1, P2, dP2, xi, wi, PI = 4.0 * atan(1.0); in CeedGaussQuadrature() local
2509 xi = cos(PI * (CeedScalar)(2 * i + 1) / ((CeedScalar)(2 * Q))); in CeedGaussQuadrature()
2512 P1 = xi; in CeedGaussQuadrature()
2515 P2 = (((CeedScalar)(2 * j - 1)) * xi * P1 - ((CeedScalar)(j - 1)) * P0) / ((CeedScalar)(j)); in CeedGaussQuadrature()
2520 dP2 = (xi * P2 - P0) * (CeedScalar)Q / (xi * xi - 1.0); in CeedGaussQuadrature()
2521 xi = xi - P2 / dP2; in CeedGaussQuadrature()
2525 P1 = xi; in CeedGaussQuadrature()
2527 P2 = (((CeedScalar)(2 * j - 1)) * xi * P1 - ((CeedScalar)(j - 1)) * P0) / ((CeedScalar)(j)); in CeedGaussQuadrature()
2531 dP2 = (xi * P2 - P0) * (CeedScalar)Q / (xi * xi - 1.0); in CeedGaussQuadrature()
2532 xi = xi - P2 / dP2; in CeedGaussQuadrature()
2535 wi = 2.0 / ((1.0 - xi * xi) * dP2 * dP2); in CeedGaussQuadrature()
2538 q_ref_1d[i] = -xi; in CeedGaussQuadrature()
2539 q_ref_1d[Q - 1 - i] = xi; in CeedGaussQuadrature()
2556 CeedScalar P0, P1, P2, dP2, d2P2, xi, wi, PI = 4.0 * atan(1.0); in CeedLobattoQuadrature() local
2571 xi = cos(PI * (CeedScalar)(i) / (CeedScalar)(Q - 1)); in CeedLobattoQuadrature()
2574 P1 = xi; in CeedLobattoQuadrature()
2577 P2 = (((CeedScalar)(2 * j - 1)) * xi * P1 - ((CeedScalar)(j - 1)) * P0) / ((CeedScalar)(j)); 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()
2584 xi = xi - dP2 / d2P2; in CeedLobattoQuadrature()
2588 P1 = xi; in CeedLobattoQuadrature()
2590 P2 = (((CeedScalar)(2 * j - 1)) * xi * P1 - ((CeedScalar)(j - 1)) * P0) / ((CeedScalar)(j)); 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()
2596 xi = xi - dP2 / d2P2; in CeedLobattoQuadrature()
2604 q_ref_1d[i] = -xi; in CeedLobattoQuadrature()
2605 q_ref_1d[Q - 1 - i] = xi; in CeedLobattoQuadrature()