1146 …(Ceed ceed, const CeedScalar *mat_A, const CeedScalar *mat_B, CeedScalar *mat_C, CeedInt m, CeedIn…  in CeedMatrixMatrixMultiply()  argument
1152       mat_C[j + i * n] = sum;  in CeedMatrixMatrixMultiply()
1448   CeedScalar *mat_C, *mat_G, *vec_D;  in CeedSimultaneousDiagonalization()  local
1450   CeedCall(CeedCalloc(n * n, &mat_C));  in CeedSimultaneousDiagonalization()
1477       mat_C[j * n + i] = mat_G[i * n + j];  in CeedSimultaneousDiagonalization()
1481 …CeedCall(CeedMatrixMatrixMultiply(ceed, (const CeedScalar *)mat_C, (const CeedScalar *)mat_A, mat_…  in CeedSimultaneousDiagonalization()
1483 …dMatrixMatrixMultiply(ceed, (const CeedScalar *)mat_X, (const CeedScalar *)mat_G, mat_C, n, n, n));  in CeedSimultaneousDiagonalization()
1486   CeedCall(CeedSymmetricSchurDecomposition(ceed, mat_C, lambda, n));  in CeedSimultaneousDiagonalization()
1493         for (CeedInt k = 0; k < n; k++) CeedScalarSwap(mat_C[k * n + j], mat_C[k * n + j + 1]);  in CeedSimultaneousDiagonalization()
1500 …CeedCall(CeedMatrixMatrixMultiply(ceed, (const CeedScalar *)mat_G, (const CeedScalar *)mat_C, mat_…  in CeedSimultaneousDiagonalization()
1503   CeedCall(CeedFree(&mat_C));  in CeedSimultaneousDiagonalization()