32   CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) u[i] *= alpha;  in ScaleN()
37   CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) u[i] = alpha;  in SetValueN()
41 …ER void CopyN(const CeedScalar *x, CeedScalar *y, const CeedInt N) { CeedPragmaSIMD for (CeedInt i…  in CopyN()
49   CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) output += u[i] * v[i];  in DotN()
75       CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) b[i] += DotN(&A[i * M], x, M);  in MatVecNM()
78 …CeedPragmaSIMD for (CeedInt i = 0; i < M; i++) { CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) b[…  in MatVecNM()
94 …CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) { CeedPragmaSIMD for (CeedInt j = 0; j < M; j++) B[…  in MatDiagNM()
97 …CeedPragmaSIMD for (CeedInt i = 0; i < M; i++) { CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) B[…  in MatDiagNM()
114           CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) {  in MatMatN()
115             CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) {  in MatMatN()
116 …        CeedPragmaSIMD for (CeedInt k = 0; k < N; k++) C[i * N + j] += A[i * N + k] * B[k * N + j];  in MatMatN()
121           CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) {  in MatMatN()
122             CeedPragmaSIMD for (CeedInt j = 0; j < 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()
132           CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) {  in MatMatN()
133             CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) {  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()
139           CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) {  in MatMatN()
140             CeedPragmaSIMD for (CeedInt j = 0; j < 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()