Searched refs:mat_copy (Results 1 – 1 of 1) sorted by relevance
| /libCEED/interface/ |
| H A D | ceed-basis.c | 1257 CeedScalar *tau, *I, *mat_copy; in CeedMatrixPseudoinverse() local
1261 CeedCall(CeedCalloc(m * n, &mat_copy)); in CeedMatrixPseudoinverse()
1262 memcpy(mat_copy, mat, m * n * sizeof mat[0]); 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()
1272 mat_pinv[j + m * (n - 1)] = I[j + m * (n - 1)] / mat_copy[n * n - 1]; in CeedMatrixPseudoinverse()
1275 …for (CeedInt k = i + 1; k < n; k++) mat_pinv[j + m * i] -= mat_copy[k + n * i] * mat_pinv[j + m * … in CeedMatrixPseudoinverse()
1276 mat_pinv[j + m * i] /= mat_copy[i + n * i]; in CeedMatrixPseudoinverse()
1283 CeedCall(CeedFree(&mat_copy)); in CeedMatrixPseudoinverse()
|