1171 int CeedQRFactorization(Ceed ceed, CeedScalar *mat, CeedScalar *tau, CeedInt m, CeedInt n) {  in CeedQRFactorization()  argument
1181       tau[i] = 0.;  in CeedQRFactorization()
1197     tau[i] = 2 * v[i] * v[i] / (v[i] * v[i] + sigma);  in CeedQRFactorization()
1201     CeedHouseholderReflect(&mat[i * n + i + 1], &v[i], tau[i], m - i, n - i - 1, n, 1);  in CeedQRFactorization()
1228 int CeedHouseholderApplyQ(CeedScalar *mat_A, const CeedScalar *mat_Q, const CeedScalar *tau, CeedTr…  in CeedHouseholderApplyQ()  argument
1237     CeedCall(CeedHouseholderReflect(&mat_A[i * row], &v[i], tau[i], m - i, n, row, col));  in CeedHouseholderApplyQ()
1257   CeedScalar *tau, *I, *mat_copy;  in CeedMatrixPseudoinverse()  local
1259   CeedCall(CeedCalloc(m, &tau));  in CeedMatrixPseudoinverse()
1265   CeedCall(CeedQRFactorization(ceed, mat_copy, tau, m, n));  in CeedMatrixPseudoinverse()
1269   CeedCall(CeedHouseholderApplyQ(I, mat_copy, tau, CEED_TRANSPOSE, m, m, n, m, 1));  in CeedMatrixPseudoinverse()
1282   CeedCall(CeedFree(&tau));  in CeedMatrixPseudoinverse()
1304   CeedScalar v[n - 1], tau[n - 1], mat_T[n * n];  in CeedSymmetricSchurDecomposition()  local
1329     tau[i] = i == n - 2 ? 2 : 2 * v[i] * v[i] / (v[i] * v[i] + sigma);  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()
1350     if (tau[i] > 0.0) {  in CeedSymmetricSchurDecomposition()
1356 …CeedHouseholderReflect(&mat[(i + 1) + n * (i + 1)], &v[i], tau[i], n - (i + 1), n - (i + 1), n, 1);  in CeedSymmetricSchurDecomposition()
1391           const CeedScalar tau = -x / z;  in CeedSymmetricSchurDecomposition()  local
1393           s = 1 / sqrt(1 + tau * tau);  in CeedSymmetricSchurDecomposition()
1394           c = s * tau;  in CeedSymmetricSchurDecomposition()
1396           const CeedScalar tau = -z / x;  in CeedSymmetricSchurDecomposition()  local
1398           c = 1 / sqrt(1 + tau * tau);  in CeedSymmetricSchurDecomposition()
1399           s = c * tau;  in CeedSymmetricSchurDecomposition()