Lines Matching refs:p
13 CeedInt p = 4; in main() local
21 CeedBasisCreateTensorH1Lagrange(ceed, 1, 1, p, p, CEED_GAUSS, &basis); in main()
24 for (int i = 0; i < p; i++) { in main()
25 for (int j = 0; j < p; j++) { in main()
27 …for (int k = 0; k < p; k++) sum += interpolation[p * k + i] * quadrature_weights[k] * interpolatio… in main()
28 M[p * i + j] = sum; in main()
29 Q[p * i + j] = sum; in main()
33 CeedSymmetricSchurDecomposition(ceed, Q, lambda, p); in main()
36 for (int i = 0; i < p; i++) { in main()
37 for (int j = 0; j < p; j++) { in main()
39 for (int k = 0; k < p; k++) sum += Q[p * i + k] * lambda[k] * Q[p * j + k]; in main()
40 Q_lambda_Qt[p * i + j] = sum; in main()
43 for (int i = 0; i < p; i++) { in main()
44 for (int j = 0; j < p; j++) { in main()
45 if (fabs(M[p * i + j] - Q_lambda_Qt[p * i + j]) > 100. * CEED_EPSILON) { in main()
47 …tion [%" CeedInt_FMT ", %" CeedInt_FMT "]: %f != %f\n", i, j, M[p * i + j], Q_lambda_Qt[p * i + j]… in main()