1304   CeedScalar v[n - 1], tau[n - 1], mat_T[n * n];  in CeedSymmetricSchurDecomposition()  local
1307   memcpy(mat_T, mat, n * n * sizeof(mat[0]));  in CeedSymmetricSchurDecomposition()
1317     v[i] = mat_T[i + n * (i + 1)];  in CeedSymmetricSchurDecomposition()
1319       v[j] = mat_T[i + n * (j + 1)];  in CeedSymmetricSchurDecomposition()
1334       mat_T[i + n * j] = 0;  in CeedSymmetricSchurDecomposition()
1335       mat_T[j + n * i] = 0;  in CeedSymmetricSchurDecomposition()
1338 …CeedHouseholderReflect(&mat_T[(i + 1) + n * (i + 1)], &v[i], tau[i], n - (i + 1), n - (i + 1), n, …  in CeedSymmetricSchurDecomposition()
1339 …CeedHouseholderReflect(&mat_T[(i + 1) + n * (i + 1)], &v[i], tau[i], n - (i + 1), n - (i + 1), 1, …  in CeedSymmetricSchurDecomposition()
1342     mat_T[i + n * (i + 1)] = R_ii;  in CeedSymmetricSchurDecomposition()
1343     mat_T[(i + 1) + n * i] = R_ii;  in CeedSymmetricSchurDecomposition()
1345       mat_T[i + n * (j + 1)] = v[j];  in CeedSymmetricSchurDecomposition()
1353         v[j]                   = mat_T[i + n * (j + 1)];  in CeedSymmetricSchurDecomposition()
1354         mat_T[i + n * (j + 1)] = 0;  in CeedSymmetricSchurDecomposition()
1369       if (fabs(mat_T[i + n * (i + 1)]) < tol) q += 1;  in CeedSymmetricSchurDecomposition()
1373       if (fabs(mat_T[i + n * (i + 1)]) < tol) p += 1;  in CeedSymmetricSchurDecomposition()
1379 …CeedScalar t_nn = mat_T[(n - 1 - q) + n * (n - 1 - q)], t_nnm1 = mat_T[(n - 2 - q) + n * (n - 1 - …  in CeedSymmetricSchurDecomposition()
1380     CeedScalar d  = (mat_T[(n - 2 - q) + n * (n - 2 - q)] - t_nn) / 2;  in CeedSymmetricSchurDecomposition()
1382     CeedScalar x  = mat_T[p + n * p] - mu;  in CeedSymmetricSchurDecomposition()
1383     CeedScalar z  = mat_T[p + n * (p + 1)];  in CeedSymmetricSchurDecomposition()
1404       CeedGivensRotation(mat_T, c, s, CEED_NOTRANSPOSE, k, k + 1, n, n);  in CeedSymmetricSchurDecomposition()
1405       CeedGivensRotation(mat_T, c, s, CEED_TRANSPOSE, k, k + 1, n, n);  in CeedSymmetricSchurDecomposition()
1412         x = mat_T[k + n * (k + 1)];  in CeedSymmetricSchurDecomposition()
1413         z = mat_T[k + n * (k + 2)];  in CeedSymmetricSchurDecomposition()
1420   for (CeedInt i = 0; i < n; i++) lambda[i] = mat_T[i + n * i];  in CeedSymmetricSchurDecomposition()