Lines Matching refs:lda
24 for (j = k; j < n; j++) v[j * mat->lda + k] = v[k * mat->lda + j]; in MatSeqDenseSymmetrize_Private()
28 for (j = k; j < n; j++) v[j * mat->lda + k] = PetscConj(v[k * mat->lda + j]); in MatSeqDenseSymmetrize_Private()
50 …PetscCallBLAS("LAPACKgetri", LAPACKgetri_(&n, mat->v, &mat->lda, mat->pivots, mat->fwork, &mat->lf… in MatSeqDenseInvertFactors_Private()
56 PetscCallBLAS("LAPACKpotri", LAPACKpotri_("L", &n, mat->v, &mat->lda, &info)); in MatSeqDenseInvertFactors_Private()
64 …PetscCallBLAS("LAPACKhetri", LAPACKhetri_("L", &n, mat->v, &mat->lda, mat->pivots, mat->fwork, &in… in MatSeqDenseInvertFactors_Private()
72 …PetscCallBLAS("LAPACKsytri", LAPACKsytri_("L", &n, mat->v, &mat->lda, mat->pivots, mat->fwork, &in… in MatSeqDenseInvertFactors_Private()
94 PetscInt m = l->lda, n = A->cmap->n, r = A->rmap->n, i, j; in MatZeroRowsColumns_SeqDense()
170 for (i = 0; i < n; i++) PetscCall(PetscArrayzero(b->v + i * b->lda, m)); in MatConvert_SeqAIJ_SeqDense()
176 b->v[*aj * b->lda + i] = *av; in MatConvert_SeqAIJ_SeqDense()
213 aa += a->lda; in MatConvert_SeqDense_SeqAIJ()
226 aa += a->lda; in MatConvert_SeqDense_SeqAIJ()
249 PetscCall(PetscBLASIntCast(x->lda, &ldax)); in MatAXPY_SeqDense()
250 PetscCall(PetscBLASIntCast(y->lda, &lday)); in MatAXPY_SeqDense()
284 PetscBLASInt one = 1, j, nz, lda = 0; in MatScale_SeqDense() local
288 PetscCall(PetscBLASIntCast(a->lda, &lda)); in MatScale_SeqDense()
289 if (lda > A->rmap->n) { in MatScale_SeqDense()
291 …for (j = 0; j < A->cmap->n; j++) PetscCallBLAS("BLASscal", BLASscal_(&nz, &alpha, v + j * lda, &on… in MatScale_SeqDense()
310 for (j = 0; j < k; j++) v[j + j * a->lda] += alpha; in MatShift_SeqDense()
319 PetscInt i, j, m = A->rmap->n, N = a->lda; in MatIsHermitian_SeqDense()
340 PetscInt i, j, m = A->rmap->n, N = a->lda; in MatIsSymmetric_SeqDense()
361 PetscInt lda = mat->lda, j, m, nlda = lda; in MatDuplicateNoCreate_SeqDense() local
368 PetscCall(MatDenseSetLDA(newi, lda)); in MatDuplicateNoCreate_SeqDense()
380 if (lda > m || nlda > m) { in MatDuplicateNoCreate_SeqDense()
381 …etscArraycpy(PetscSafePointerPlusOffset(v, j * nlda), PetscSafePointerPlusOffset(av, j * lda), m)); in MatDuplicateNoCreate_SeqDense()
409 …PetscCallBLAS("LAPACKgetrs", LAPACKgetrs_(T ? "T" : "N", &m, &nrhs, mat->v, &mat->lda, mat->pivots… in MatSolve_SeqDense_Internal_LU()
425 PetscCallBLAS("LAPACKpotrs", LAPACKpotrs_("L", &m, &nrhs, mat->v, &mat->lda, x, &m, &info)); in MatSolve_SeqDense_Internal_Cholesky()
433 …PetscCallBLAS("LAPACKhetrs", LAPACKhetrs_("L", &m, &nrhs, mat->v, &mat->lda, mat->pivots, x, &m, &… in MatSolve_SeqDense_Internal_Cholesky()
440 …PetscCallBLAS("LAPACKsytrs", LAPACKsytrs_("L", &m, &nrhs, mat->v, &mat->lda, mat->pivots, x, &m, &… in MatSolve_SeqDense_Internal_Cholesky()
465 …PetscCallBLAS("LAPACKormqr", LAPACKormqr_("L", &trans, &m, &nrhs, &mat->rank, mat->v, &mat->lda, m… in MatSolve_SeqDense_Internal_QR()
473 …PetscCallBLAS("LAPACKormqr", LAPACKormqr_("L", &trans, &m, &nrhs, &mat->rank, mat->v, &mat->lda, m… in MatSolve_SeqDense_Internal_QR()
477 …PetscCallBLAS("LAPACKtrtrs", LAPACKtrtrs_("U", "N", "N", &mat->rank, &nrhs, mat->v, &mat->lda, x, … in MatSolve_SeqDense_Internal_QR()
495 …PetscCallBLAS("LAPACKtrtrs", LAPACKtrtrs_("U", "T", "N", &m, &nrhs, mat->v, &mat->lda, x, &ldx, &i… in MatSolveTranspose_SeqDense_Internal_QR()
503 …PetscCallBLAS("LAPACKormqr", LAPACKormqr_("L", "N", &m, &nrhs, &mat->rank, mat->v, &mat->lda, mat-… in MatSolveTranspose_SeqDense_Internal_QR()
512 …PetscCallBLAS("LAPACKormqr", LAPACKormqr_("L", "N", &m, &nrhs, &mat->rank, mat->v, &mat->lda, mat-… in MatSolveTranspose_SeqDense_Internal_QR()
800 PetscCallBLAS("LAPACKgetrf", LAPACKgetrf_(&m, &n, mat->v, &mat->lda, mat->pivots, &info)); in MatLUFactor_SeqDense()
847 PetscCallBLAS("LAPACKpotrf", LAPACKpotrf_("L", &n, mat->v, &mat->lda, &info)); in MatCholeskyFactor_SeqDense()
857 …PetscCallBLAS("LAPACKhetrf", LAPACKhetrf_("L", &n, mat->v, &mat->lda, mat->pivots, &dummy, &mat->l… in MatCholeskyFactor_SeqDense()
863 …PetscCallBLAS("LAPACKhetrf", LAPACKhetrf_("L", &n, mat->v, &mat->lda, mat->pivots, mat->fwork, &ma… in MatCholeskyFactor_SeqDense()
873 …PetscCallBLAS("LAPACKsytrf", LAPACKsytrf_("L", &n, mat->v, &mat->lda, mat->pivots, &dummy, &mat->l… in MatCholeskyFactor_SeqDense()
879 …PetscCallBLAS("LAPACKsytrf", LAPACKsytrf_("L", &n, mat->v, &mat->lda, mat->pivots, mat->fwork, &ma… in MatCholeskyFactor_SeqDense()
933 …PetscCallBLAS("LAPACKgeqrf", LAPACKgeqrf_(&m, &n, mat->v, &mat->lda, mat->tau, &dummy, &mat->lfwor… in MatQRFactor_SeqDense()
939 …PetscCallBLAS("LAPACKgeqrf", LAPACKgeqrf_(&m, &n, mat->v, &mat->lda, mat->tau, mat->fwork, &mat->l… in MatQRFactor_SeqDense()
1065 …llBLAS("BLASgemv", BLASgemv_("C", &m, &n, &_DOne, v + c_start * mat->lda, &mat->lda, x, &_One, &_D… in MatMultColumnRangeKernel_SeqDense()
1066 …llBLAS("BLASgemv", BLASgemv_("T", &m, &n, &_DOne, v + c_start * mat->lda, &mat->lda, x, &_One, &_D… in MatMultColumnRangeKernel_SeqDense()
1068 …PetscCallBLAS("BLASgemv", BLASgemv_("N", &m, &n, &_DOne, v + c_start * mat->lda, &mat->lda, x + c_… in MatMultColumnRangeKernel_SeqDense()
1120 …llBLAS("BLASgemv", BLASgemv_("C", &m, &n, &_DOne, v + c_start * mat->lda, &mat->lda, x, &_One, &_D… in MatMultAddColumnRangeKernel_SeqDense()
1121 …else PetscCallBLAS("BLASgemv", BLASgemv_("T", &m, &n, &_DOne, v + c_start * mat->lda, &mat->lda, x… in MatMultAddColumnRangeKernel_SeqDense()
1123 …PetscCallBLAS("BLASgemv", BLASgemv_("N", &m, &n, &_DOne, v + c_start * mat->lda, &mat->lda, x + c_… in MatMultAddColumnRangeKernel_SeqDense()
1194 v += mat->lda; in MatGetRow_SeqDense()
1234 av[indexn[j] * mat->lda + indexm[i]] = v ? v[idx++] : (idx++, 0.0); in MatSetValues_SeqDense()
1250 av[indexn[j] * mat->lda + indexm[i]] += v ? v[idx++] : (idx++, 0.0); in MatSetValues_SeqDense()
1268 av[indexn[j] * mat->lda + indexm[i]] = v ? v[idx++] : (idx++, 0.0); in MatSetValues_SeqDense()
1284 av[indexn[j] * mat->lda + indexm[i]] += v ? v[idx++] : (idx++, 0.0); in MatSetValues_SeqDense()
1322 *v++ = vv[indexn[j] * mat->lda + indexm[i]]; in MatGetValues_SeqDense()
1333 PetscInt header[4], M, N, m, lda, i, j; in MatView_Dense_Binary() local
1369 PetscCall(MatDenseGetLDA(mat, &lda)); in MatView_Dense_Binary()
1371 for (j = 0; j < N; j++, k++) vwork[k] = v[i + (size_t)lda * j]; in MatView_Dense_Binary()
1381 PetscInt header[4], M, N, m, nz, lda, i, j, k; in MatLoad_Dense_Binary() local
1415 PetscCall(MatDenseGetLDA(mat, &lda)); in MatLoad_Dense_Binary()
1423 for (i = 0; i < m; i++) v[i + (size_t)lda * j] = vwork[(size_t)i * N + j]; in MatLoad_Dense_Binary()
1437 for (j = 0; j < rlens[i]; j++, k++) v[i + lda * icols[k]] = vwork[k]; in MatLoad_Dense_Binary()
1503 v += a->lda; in MatView_SeqDense_ASCII()
1513 v = av + j * a->lda; in MatView_SeqDense_ASCII()
1541 v += a->lda; in MatView_SeqDense_ASCII()
1786 PetscInt k, j, m = A->rmap->n, M = mat->lda, n = A->cmap->n; in MatTranspose_SeqDense()
1820 PetscCall(PetscBLASIntCast(n, &mat->lda)); in MatTranspose_SeqDense()
1841 M2 = tmatd->lda; in MatTranspose_SeqDense()
1875 v1 += mat1->lda; in MatEqual_SeqDense()
1876 v2 += mat2->lda; in MatEqual_SeqDense()
1897 for (i = 0; i < len; i++) x[i] = vv[i * mat->lda + i]; in MatGetDiagonal_SeqDense()
1921 v += mat->lda; in MatDiagonalScale_SeqDense()
1933 v = vv + i * mat->lda; in MatDiagonalScale_SeqDense()
1948 PetscInt lda, m = A->rmap->n, i, j; in MatNorm_SeqDense() local
1952 PetscCall(MatDenseGetLDA(A, &lda)); in MatNorm_SeqDense()
1955 if (lda > m) { in MatNorm_SeqDense()
1957 v = vv + j * lda; in MatNorm_SeqDense()
1980 v = vv + j * mat->lda;
1996 v += mat->lda;
2024 PetscInt lda = l->lda, m = A->rmap->n, n = A->cmap->n, j; in MatZeroEntries_SeqDense() local
2029 if (lda > m) { in MatZeroEntries_SeqDense()
2030 for (j = 0; j < n; j++) PetscCall(PetscArrayzero(v + j * lda, m)); in MatZeroEntries_SeqDense()
2041 PetscInt m = l->lda, n = A->cmap->n, i, j; in MatZeroRows_SeqDense()
2082 static PetscErrorCode MatDenseGetLDA_SeqDense(Mat A, PetscInt *lda) in MatDenseGetLDA_SeqDense() argument
2087 *lda = mat->lda; in MatDenseGetLDA_SeqDense()
2123 PetscErrorCode MatDenseGetLDA(Mat A, PetscInt *lda) in MatDenseGetLDA() argument
2127 PetscAssertPointer(lda, 2); in MatDenseGetLDA()
2129 PetscUseMethod(A, "MatDenseGetLDA_C", (Mat, PetscInt *), (A, lda)); in MatDenseGetLDA()
2146 PetscErrorCode MatDenseSetLDA(Mat A, PetscInt lda) in MatDenseSetLDA() argument
2150 PetscTryMethod(A, "MatDenseSetLDA_C", (Mat, PetscInt), (A, lda)); in MatDenseSetLDA()
2587 av = v + mat->lda * icol[i]; in MatCreateSubMatrix_SeqDense()
2620 PetscInt lda1 = a->lda, lda2 = b->lda, m = A->rmap->n, n = A->cmap->n, j; in MatCopy_SeqDense()
2662 for (i = 0; i < A->rmap->n; i++) aa[i + j * mat->lda] = PetscConj(aa[i + j * mat->lda]); in MatConjugate_SeqDense()
2678 for (i = 0; i < A->rmap->n; i++) aa[i + j * mat->lda] = PetscRealPart(aa[i + j * mat->lda]); in MatRealPart_SeqDense()
2693 … for (i = 0; i < A->rmap->n; i++) aa[i + j * mat->lda] = PetscImaginaryPart(aa[i + j * mat->lda]); in MatImaginaryPart_SeqDense()
2738 …LASgemm", BLASgemm_("N", "N", &m, &n, &k, &_DOne, av, &a->lda, bv, &b->lda, &_DZero, cv, &c->lda)); in MatMatMultNumeric_SeqDense_SeqDense()
2787 …LASgemm", BLASgemm_("N", "T", &m, &n, &k, &_DOne, av, &a->lda, bv, &b->lda, &_DZero, cv, &c->lda)); in MatMatTransposeMultNumeric_SeqDense_SeqDense()
2836 …LASgemm", BLASgemm_("T", "N", &m, &n, &k, &_DOne, av, &a->lda, bv, &b->lda, &_DZero, cv, &c->lda)); in MatTransposeMatMultNumeric_SeqDense_SeqDense()
2906 if (PetscRealPart(x[i]) < PetscRealPart(aa[i + a->lda * j])) { in MatGetRowMax_SeqDense()
2907 x[i] = aa[i + a->lda * j]; in MatGetRowMax_SeqDense()
2934 atmp = PetscAbsScalar(aa[i + a->lda * j]); in MatGetRowMaxAbs_SeqDense()
2963 if (PetscRealPart(x[i]) > PetscRealPart(aa[i + a->lda * j])) { in MatGetRowMin_SeqDense()
2964 x[i] = aa[i + a->lda * j]; in MatGetRowMin_SeqDense()
2984 PetscCall(PetscArraycpy(x, aa + col * a->lda, A->rmap->n)); in MatGetColumnVector_SeqDense()
3037 PetscInt lda, m, n, i, j; in MatSetRandom_SeqDense() local
3041 PetscCall(MatDenseGetLDA(x, &lda)); in MatSetRandom_SeqDense()
3044 for (i = 0; i < m; i++) PetscCall(PetscRandomGetValue(rctx, a + j * lda + i)); in MatSetRandom_SeqDense()
3059 *vals = v + col * a->lda; in MatDenseGetColumn_SeqDense()
3291 if (b->lda <= 0) PetscCall(PetscBLASIntCast(B->rmap->n, &b->lda)); in MatSeqDenseSetPreallocation_SeqDense()
3295 PetscCall(PetscCalloc1((size_t)b->lda * B->cmap->n, &b->v)); in MatSeqDenseSetPreallocation_SeqDense()
3347 PetscErrorCode MatDenseSetLDA_SeqDense(Mat B, PetscInt lda) in MatDenseSetLDA_SeqDense() argument
3354 …PetscCheck(b->user_alloc || !data || b->lda == lda, PETSC_COMM_SELF, PETSC_ERR_ORDER, "LDA cannot … in MatDenseSetLDA_SeqDense()
3355 …etscCheck(lda >= B->rmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "LDA %" PetscInt_FMT " must be at… in MatDenseSetLDA_SeqDense()
3356 PetscCall(PetscBLASIntCast(lda, &b->lda)); in MatDenseSetLDA_SeqDense()
3421 PetscCall(VecPlaceArray(a->cvec, a->ptrinuse + (size_t)col * (size_t)a->lda)); in MatDenseGetColumnVec_SeqDense()
3451 …all(VecPlaceArray(a->cvec, PetscSafePointerPlusOffset(a->ptrinuse, (size_t)col * (size_t)a->lda))); in MatDenseGetColumnVecRead_SeqDense()
3483 …all(VecPlaceArray(a->cvec, PetscSafePointerPlusOffset(a->ptrinuse, (size_t)col * (size_t)a->lda))); in MatDenseGetColumnVecWrite_SeqDense()
3512 …gin, cend - cbegin, PetscSafePointerPlusOffset(a->v, rbegin + (size_t)cbegin * a->lda), &a->cmat)); in MatDenseGetSubMatrix_SeqDense()
3514 …l(MatDensePlaceArray(a->cmat, PetscSafePointerPlusOffset(a->v, rbegin + (size_t)cbegin * a->lda))); in MatDenseGetSubMatrix_SeqDense()
3516 PetscCall(MatDenseSetLDA(a->cmat, a->lda)); in MatDenseGetSubMatrix_SeqDense()