400   CeedInt  dim, num_comp, num_q_comp, num_nodes, P_1d = 1, Q_1d = 1, total_num_points = 0;  in CeedBasisApplyAtPointsCheckDims()  local
404   CeedCall(CeedBasisGetNumNodes1D(basis, &P_1d));  in CeedBasisApplyAtPointsCheckDims()
478   CeedInt dim, num_comp, P_1d = 1, Q_1d = 1, total_num_points = num_points[0];  in CeedBasisApplyAtPoints_Core()  local
483   CeedCall(CeedBasisGetNumNodes1D(basis, &P_1d));  in CeedBasisApplyAtPoints_Core()
507     CeedCall(CeedCalloc(P_1d * Q_1d, &chebyshev_interp_1d));  in CeedBasisApplyAtPoints_Core()
508     CeedCall(CeedCalloc(P_1d * Q_1d, &chebyshev_grad_1d));  in CeedBasisApplyAtPoints_Core()
515 …CeedCall(CeedBasisCreateTensorH1(ceed, dim, num_comp, P_1d, Q_1d, chebyshev_interp_1d, chebyshev_g…  in CeedBasisApplyAtPoints_Core()
532 …CeedCall(CeedBasisCreateTensorH1Lagrange(ceed_ref, dim, num_comp, P_1d, Q_1d, CEED_GAUSS, &basis_r…  in CeedBasisApplyAtPoints_Core()
735   CeedInt           P_1d, Q_1d;  in CeedBasisGetCollocatedGrad()  local
741   CeedCall(CeedBasisGetNumNodes1D(basis, &P_1d));  in CeedBasisGetCollocatedGrad()
745   CeedCall(CeedCalloc(P_1d * Q_1d, &interp_1d_pinv));  in CeedBasisGetCollocatedGrad()
747   CeedCall(CeedMatrixPseudoinverse(ceed, interp_1d, Q_1d, P_1d, interp_1d_pinv));  in CeedBasisGetCollocatedGrad()
749 …atrixMultiply(ceed, grad_1d, (const CeedScalar *)interp_1d_pinv, collo_grad_1d, Q_1d, Q_1d, P_1d));  in CeedBasisGetCollocatedGrad()
767   CeedInt           P_1d, Q_1d;  in CeedBasisGetChebyshevInterp1D()  local
773   CeedCall(CeedBasisGetNumNodes1D(basis, &P_1d));  in CeedBasisGetChebyshevInterp1D()
778 …CeedCheck((P_1d > 0) && (Q_1d > 0), ceed, CEED_ERROR_INCOMPATIBLE, "CeedBasis dimensions are malfo…  in CeedBasisGetChebyshevInterp1D()
789 …ixMatrixMultiply(ceed, chebyshev_coeffs_1d_inv, interp_1d, chebyshev_interp_1d, Q_1d, P_1d, Q_1d));  in CeedBasisGetChebyshevInterp1D()
824   if (basis->is_tensor_basis && (basis->Q_1d == basis->P_1d)) {  in CeedBasisIsCollocated()
827     for (CeedInt i = 0; i < basis->P_1d; i++) {  in CeedBasisIsCollocated()
828 …*is_collocated = *is_collocated && (fabs(basis->interp_1d[i + basis->P_1d * i] - 1.0) < 10 * CEED_…  in CeedBasisIsCollocated()
830 …if (j != i) *is_collocated = *is_collocated && (fabs(basis->interp_1d[j + basis->P_1d * i]) < 10 *…  in CeedBasisIsCollocated()
944     CeedInt dim, num_comp, P_1d, Q_1d;  in CeedBasisGetFlopsEstimate()  local
948     CeedCall(CeedBasisGetNumNodes1D(basis, &P_1d));  in CeedBasisGetFlopsEstimate()
951       P_1d = Q_1d;  in CeedBasisGetFlopsEstimate()
952       Q_1d = P_1d;  in CeedBasisGetFlopsEstimate()
954     CeedInt tensor_flops = 0, pre = num_comp * CeedIntPow(P_1d, dim - 1), post = 1;  in CeedBasisGetFlopsEstimate()
957       tensor_flops += 2 * pre * P_1d * post * Q_1d;  in CeedBasisGetFlopsEstimate()
958       pre /= P_1d;  in CeedBasisGetFlopsEstimate()
976         pre /= P_1d;  in CeedBasisGetFlopsEstimate()
1536 int CeedBasisCreateTensorH1(Ceed ceed, CeedInt dim, CeedInt num_comp, CeedInt P_1d, CeedInt Q_1d, c…  in CeedBasisCreateTensorH1()  argument
1543 …CeedCall(CeedBasisCreateTensorH1(delegate, dim, num_comp, P_1d, Q_1d, interp_1d, grad_1d, q_ref_1d…  in CeedBasisCreateTensorH1()
1550   CeedCheck(P_1d > 0, ceed, CEED_ERROR_DIMENSION, "CeedBasis must have at least 1 node");  in CeedBasisCreateTensorH1()
1561   (*basis)->P_1d            = P_1d;  in CeedBasisCreateTensorH1()
1563   (*basis)->P               = CeedIntPow(P_1d, dim);  in CeedBasisCreateTensorH1()
1570   CeedCall(CeedCalloc(Q_1d * P_1d, &(*basis)->interp_1d));  in CeedBasisCreateTensorH1()
1571   CeedCall(CeedCalloc(Q_1d * P_1d, &(*basis)->grad_1d));  in CeedBasisCreateTensorH1()
1572   if (interp_1d) memcpy((*basis)->interp_1d, interp_1d, Q_1d * P_1d * sizeof(interp_1d[0]));  in CeedBasisCreateTensorH1()
1573   if (grad_1d) memcpy((*basis)->grad_1d, grad_1d, Q_1d * P_1d * sizeof(grad_1d[0]));  in CeedBasisCreateTensorH1()
1574 …CeedCall(ceed->BasisCreateTensorH1(dim, P_1d, Q_1d, interp_1d, grad_1d, q_ref_1d, q_weight_1d, *ba…  in CeedBasisCreateTensorH1()
1983 …fprintf(stream, "%s  P: %" CeedInt_FMT "\n%s  Q: %" CeedInt_FMT "\n", tabs, basis->P_1d, tabs, bas…  in CeedBasisView()
1990     CeedInt           P_1d, Q_1d;  in CeedBasisView()  local
1993     CeedCall(CeedBasisGetNumNodes1D(basis, &P_1d));  in CeedBasisView()
2002     CeedCall(CeedScalarView("interp1d", "\t% 12.8f", Q_1d, P_1d, interp_1d, tabs, stream));  in CeedBasisView()
2003     CeedCall(CeedScalarView("grad1d", "\t% 12.8f", Q_1d, P_1d, grad_1d, tabs, stream));  in CeedBasisView()
2255 int CeedBasisGetNumNodes1D(CeedBasis basis, CeedInt *P_1d) {  in CeedBasisGetNumNodes1D()  argument
2257   *P_1d = basis->P_1d;  in CeedBasisGetNumNodes1D()
2344           CeedInt p = (node / CeedIntPow(basis->P_1d, d)) % basis->P_1d;  in CeedBasisGetInterp()
2347           basis->interp[qpt * (basis->P) + node] *= basis->interp_1d[q * basis->P_1d + p];  in CeedBasisGetInterp()
2398             CeedInt p = (node / CeedIntPow(basis->P_1d, d)) % basis->P_1d;  in CeedBasisGetGrad()
2401 …== d) basis->grad[(i * basis->Q + qpt) * (basis->P) + node] *= basis->grad_1d[q * basis->P_1d + p];  in CeedBasisGetGrad()
2402 …lse basis->grad[(i * basis->Q + qpt) * (basis->P) + node] *= basis->interp_1d[q * basis->P_1d + p];  in CeedBasisGetGrad()