44 …FUNCTION_HELPER void CopyMat3(const CeedScalar A[3][3], CeedScalar B[3][3]) { CopyN((const CeedSca…  in CopyMat3()
71 CEED_QFUNCTION_HELPER void MatVecNM(const CeedScalar *A, const CeedScalar *x, const CeedInt N, cons…  in MatVecNM()  argument
75       CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) b[i] += DotN(&A[i * M], x, M);  in MatVecNM()
78 … i = 0; i < M; i++) { CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) b[i] += A[j * M + i] * x[j]; }  in MatVecNM()
84 CEED_QFUNCTION_HELPER void MatVec3(const CeedScalar A[3][3], const CeedScalar x[3], const CeedTrans…  in MatVec3()
85   MatVecNM((const CeedScalar *)A, (const CeedScalar *)x, 3, 3, transpose_A, (CeedScalar *)b);  in MatVec3()
90 CEED_QFUNCTION_HELPER void MatDiagNM(const CeedScalar *A, const CeedScalar *D, const CeedInt N, con…  in MatDiagNM()  argument
94 …i < N; i++) { CeedPragmaSIMD for (CeedInt j = 0; j < M; j++) B[i * M + j] += D[i] * A[i * M + j]; }  in MatDiagNM()
97 …i < M; i++) { CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) B[i * N + j] += D[i] * A[j * M + i]; }  in MatDiagNM()
104 CEED_QFUNCTION_HELPER void MatDiag3(const CeedScalar A[3][3], const CeedScalar D[3], const CeedTran…  in MatDiag3()
105   MatDiagNM((const CeedScalar *)A, (const CeedScalar *)D, 3, 3, transpose_A, (CeedScalar *)B);  in MatDiag3()
108 CEED_QFUNCTION_HELPER void MatMatN(const CeedScalar *A, const CeedScalar *B, const CeedInt N, const…  in MatMatN()  argument
116 …        CeedPragmaSIMD for (CeedInt k = 0; k < N; k++) C[i * N + j] += A[i * N + k] * B[k * N + j];  in MatMatN()
123 …        CeedPragmaSIMD for (CeedInt k = 0; k < N; k++) C[i * N + j] += A[i * N + k] * B[j * N + k];  in MatMatN()
134 …        CeedPragmaSIMD for (CeedInt k = 0; k < N; k++) C[i * N + j] += A[k * N + i] * B[k * N + j];  in MatMatN()
141 …        CeedPragmaSIMD for (CeedInt k = 0; k < N; k++) C[i * N + j] += A[k * N + i] * B[j * N + k];  in MatMatN()
151 CEED_QFUNCTION_HELPER void MatMat3(const CeedScalar A[3][3], const CeedScalar B[3][3], const CeedTr…  in MatMat3()
153 …MatMatN((const CeedScalar *)A, (const CeedScalar *)B, 3, transpose_A, transpose_B, (CeedScalar *)C…  in MatMat3()
157 CEED_QFUNCTION_HELPER void KMUnpack(const CeedScalar v[6], CeedScalar A[3][3]) {  in KMUnpack()
159   A[0][0]                 = v[0];  in KMUnpack()
160   A[1][1]                 = v[1];  in KMUnpack()
161   A[2][2]                 = v[2];  in KMUnpack()
162   A[2][1] = A[1][2] = weight * v[3];  in KMUnpack()
163   A[2][0] = A[0][2] = weight * v[4];  in KMUnpack()
164   A[1][0] = A[0][1] = weight * v[5];  in KMUnpack()
168 CEED_QFUNCTION_HELPER void KMPack(const CeedScalar A[3][3], CeedScalar v[6]) {  in KMPack()
170   v[0]                    = A[0][0];  in KMPack()
171   v[1]                    = A[1][1];  in KMPack()
172   v[2]                    = A[2][2];  in KMPack()
173   v[3]                    = A[2][1] * weight;  in KMPack()
174   v[4]                    = A[2][0] * weight;  in KMPack()
175   v[5]                    = A[1][0] * weight;  in KMPack()