Lines Matching refs:rmap
16 PetscInt j, k, n = A->rmap->n; in MatSeqDenseSymmetrize_Private()
20 …PetscCheck(A->rmap->n == A->cmap->n, PetscObjectComm((PetscObject)A), PETSC_ERR_SUP, "Cannot symme… in MatSeqDenseSymmetrize_Private()
41 if (!A->rmap->n || !A->cmap->n) PetscFunctionReturn(PETSC_SUCCESS); in MatSeqDenseInvertFactors_Private()
94 PetscInt m = l->lda, n = A->cmap->n, r = A->rmap->n, i, j; in MatZeroRowsColumns_SeqDense()
102 …rmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Row %" PetscInt_FMT " requested to be zeroed … in MatZeroRowsColumns_SeqDense()
112 …PetscCheck(A->rmap->n == A->cmap->n, PETSC_COMM_SELF, PETSC_ERR_SUP, "Only coded for square matric… in MatZeroRowsColumns_SeqDense()
138 …PetscCheck(A->rmap->n == A->cmap->n, PETSC_COMM_SELF, PETSC_ERR_SUP, "Only coded for square matric… in MatZeroRowsColumns_SeqDense()
153 PetscInt *ai = a->i, *aj = a->j, m = A->rmap->N, n = A->cmap->N, i; in MatConvert_SeqAIJ_SeqDense()
205 PetscCall(PetscCalloc3(A->rmap->n, &rows, A->rmap->n, &nnz, A->rmap->n, &vals)); in MatConvert_SeqDense_SeqAIJ()
208 PetscCall(MatSetSizes(B, A->rmap->n, A->cmap->n, A->rmap->N, A->cmap->N)); in MatConvert_SeqDense_SeqAIJ()
211 for (i = 0; i < A->rmap->n; i++) in MatConvert_SeqDense_SeqAIJ()
212 if (aa[i] != 0.0 || (i == j && A->cmap->n == A->rmap->n)) ++nnz[i]; in MatConvert_SeqDense_SeqAIJ()
220 for (i = 0; i < A->rmap->n; i++) in MatConvert_SeqDense_SeqAIJ()
221 if (aa[i] != 0.0 || (i == j && A->cmap->n == A->rmap->n)) { in MatConvert_SeqDense_SeqAIJ()
247 PetscCall(PetscBLASIntCast(X->rmap->n * X->cmap->n, &N)); in MatAXPY_SeqDense()
248 PetscCall(PetscBLASIntCast(X->rmap->n, &m)); in MatAXPY_SeqDense()
264 PetscLogDouble N = A->rmap->n * A->cmap->n; in MatGetInfo_SeqDense()
289 if (lda > A->rmap->n) { in MatScale_SeqDense()
290 PetscCall(PetscBLASIntCast(A->rmap->n, &nz)); in MatScale_SeqDense()
293 PetscCall(PetscBLASIntCast(A->rmap->n * A->cmap->n, &nz)); in MatScale_SeqDense()
296 PetscCall(PetscLogFlops(A->rmap->n * A->cmap->n)); in MatScale_SeqDense()
309 k = PetscMin(A->rmap->n, A->cmap->n); in MatShift_SeqDense()
319 PetscInt i, j, m = A->rmap->n, N = a->lda; in MatIsHermitian_SeqDense()
324 if (A->rmap->n != A->cmap->n) PetscFunctionReturn(PETSC_SUCCESS); in MatIsHermitian_SeqDense()
340 PetscInt i, j, m = A->rmap->n, N = a->lda; in MatIsSymmetric_SeqDense()
345 if (A->rmap->n != A->cmap->n) PetscFunctionReturn(PETSC_SUCCESS); in MatIsSymmetric_SeqDense()
365 PetscCall(PetscLayoutReference(A->rmap, &newi->rmap)); in MatDuplicateNoCreate_SeqDense()
379 m = A->rmap->n; in MatDuplicateNoCreate_SeqDense()
383 PetscCall(PetscArraycpy(v, av, A->rmap->n * A->cmap->n)); in MatDuplicateNoCreate_SeqDense()
396 PetscCall(MatSetSizes(*newmat, A->rmap->n, A->cmap->n, A->rmap->n, A->cmap->n)); in MatDuplicate_SeqDense()
493 if (A->rmap->n == A->cmap->n && mat->rank == A->rmap->n) { in MatSolveTranspose_SeqDense_Internal_QR()
528 PetscCall(PetscBLASIntCast(A->rmap->n, &m)); in MatSolve_SeqDense_SetUp()
651 PetscCall(PetscBLASIntCast(A->rmap->n, &m)); in MatMatSolve_SeqDense_SetUp()
796 PetscCall(PetscBLASIntCast(A->rmap->n, &m)); in MatLUFactor_SeqDense()
797 if (!mat->pivots) PetscCall(PetscMalloc1(A->rmap->n, &mat->pivots)); in MatLUFactor_SeqDense()
798 if (!A->rmap->n || !A->cmap->n) PetscFunctionReturn(PETSC_SUCCESS); in MatLUFactor_SeqDense()
844 if (!A->rmap->n || !A->cmap->n) PetscFunctionReturn(PETSC_SUCCESS); in MatCholeskyFactor_SeqDense()
851 if (!mat->pivots) PetscCall(PetscMalloc1(A->rmap->n, &mat->pivots)); in MatCholeskyFactor_SeqDense()
867 if (!mat->pivots) PetscCall(PetscMalloc1(A->rmap->n, &mat->pivots)); in MatCholeskyFactor_SeqDense()
921 PetscCall(PetscBLASIntCast(A->rmap->n, &m)); in MatQRFactor_SeqDense()
927 if (!A->rmap->n || !A->cmap->n) PetscFunctionReturn(PETSC_SUCCESS); in MatQRFactor_SeqDense()
982 PetscCall(MatSetSizes(*fact, A->rmap->n, A->cmap->n, A->rmap->n, A->cmap->n)); in MatGetFactor_seqdense_petsc()
1009 PetscInt m = A->rmap->n, i; in MatSOR_SeqDense()
1053 PetscCall(PetscBLASIntCast(A->rmap->n, &m)); in MatMultColumnRangeKernel_SeqDense()
1113 PetscCall(PetscBLASIntCast(A->rmap->n, &m)); in MatMultAddColumnRangeKernel_SeqDense()
1233 …exm[i] < A->rmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Row too large: row %" PetscInt_FM… in MatSetValues_SeqDense()
1249 …exm[i] < A->rmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Row too large: row %" PetscInt_FM… in MatSetValues_SeqDense()
1261 …exm[i] < A->rmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Row too large: row %" PetscInt_FM… in MatSetValues_SeqDense()
1277 …exm[i] < A->rmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Row too large: row %" PetscInt_FM… in MatSetValues_SeqDense()
1315 … A->rmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Row %" PetscInt_FMT " requested larger th… in MatGetValues_SeqDense()
1405 if (mat->rmap->N < 0) mat->rmap->N = M; in MatLoad_Dense_Binary()
1490 for (i = 0; i < A->rmap->n; i++) { in MatView_SeqDense_ASCII()
1514 for (i = 0; i < A->rmap->n; i++) { in MatView_SeqDense_ASCII()
1524 …erASCIIPrintf(viewer, "%% Size = %" PetscInt_FMT " %" PetscInt_FMT " \n", A->rmap->n, A->cmap->n)); in MatView_SeqDense_ASCII()
1525 …rintf(viewer, "%s = zeros(%" PetscInt_FMT ",%" PetscInt_FMT ");\n", name, A->rmap->n, A->cmap->n)); in MatView_SeqDense_ASCII()
1529 for (i = 0; i < A->rmap->n; i++) { in MatView_SeqDense_ASCII()
1557 PetscInt m = A->rmap->n, n = A->cmap->n, i, j; in MatView_SeqDense_Draw_Zoom()
1629 yr = A->rmap->n; in MatView_SeqDense_Draw()
1713 …State((PetscObject)mat, "Rows %" PetscInt_FMT " Cols %" PetscInt_FMT, mat->rmap->n, mat->cmap->n)); in MatDestroy_SeqDense()
1786 PetscInt k, j, m = A->rmap->n, M = mat->lda, n = A->cmap->n; in MatTranspose_SeqDense()
1821 tmplayout = A->rmap; in MatTranspose_SeqDense()
1822 A->rmap = A->cmap; in MatTranspose_SeqDense()
1833 PetscCall(MatSetSizes(tmat, A->cmap->n, A->rmap->n, A->cmap->n, A->rmap->n)); in MatTranspose_SeqDense()
1862 if (A1->rmap->n != A2->rmap->n) { in MatEqual_SeqDense()
1873 PetscCall(PetscArraycmp(v1, v2, A1->rmap->n, flg)); in MatEqual_SeqDense()
1894 len = PetscMin(A->rmap->n, A->cmap->n); in MatGetDiagonal_SeqDense()
1896 PetscCheck(n == A->rmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Nonconforming mat and vec"); in MatGetDiagonal_SeqDense()
1908 PetscInt i, j, m = A->rmap->n, n = A->cmap->n; in MatDiagonalScale_SeqDense()
1915 PetscCheck(m == A->rmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Left scaling vec wrong size"); in MatDiagonalScale_SeqDense()
1948 PetscInt lda, m = A->rmap->n, i, j; in MatNorm_SeqDense()
1965 PetscBLASInt one = 1, cnt = A->cmap->n * A->rmap->n; in MatNorm_SeqDense()
1969 for (i = 0; i < A->cmap->n * A->rmap->n; i++) { in MatNorm_SeqDense()
1976 PetscCall(PetscLogFlops(2.0 * A->cmap->n * A->rmap->n)); in MatNorm_SeqDense()
1982 for (i = 0; i < A->rmap->n; i++) {
1988 PetscCall(PetscLogFlops(1.0 * A->cmap->n * A->rmap->n));
1991 for (j = 0; j < A->rmap->n; j++) {
2000 PetscCall(PetscLogFlops(1.0 * A->cmap->n * A->rmap->n));
2024 PetscInt lda = l->lda, m = A->rmap->n, n = A->cmap->n, j; in MatZeroEntries_SeqDense()
2049 …rmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Row %" PetscInt_FMT " requested to be zeroed … in MatZeroRows_SeqDense()
2072 …PetscCheck(A->rmap->n == A->cmap->n, PETSC_COMM_SELF, PETSC_ERR_SUP, "Only coded for square matric… in MatZeroRows_SeqDense()
2620 PetscInt lda1 = a->lda, lda2 = b->lda, m = A->rmap->n, n = A->cmap->n, j; in MatCopy_SeqDense()
2628 …PetscCheck(m == B->rmap->n && n == B->cmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "size(B) != siz… in MatCopy_SeqDense()
2634 PetscCall(PetscArraycpy(vb, va, A->rmap->n * A->cmap->n)); in MatCopy_SeqDense()
2646 PetscCall(PetscLayoutSetUp(A->rmap)); in MatSetUp_SeqDense()
2656 PetscInt min = PetscMin(A->rmap->n, A->cmap->n); in MatConjugate_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()
2701 PetscInt m = A->rmap->n, n = B->cmap->n; in MatMatMultSymbolic_SeqDense_SeqDense()
2731 PetscCall(PetscBLASIntCast(C->rmap->n, &m)); in MatMatMultNumeric_SeqDense_SeqDense()
2748 PetscInt m = A->rmap->n, n = B->rmap->n; in MatMatTransposeMultSymbolic_SeqDense_SeqDense()
2780 PetscCall(PetscBLASIntCast(C->rmap->n, &m)); in MatMatTransposeMultNumeric_SeqDense_SeqDense()
2829 PetscCall(PetscBLASIntCast(C->rmap->n, &m)); in MatTransposeMatMultNumeric_SeqDense_SeqDense()
2831 PetscCall(PetscBLASIntCast(A->rmap->n, &k)); in MatTransposeMatMultNumeric_SeqDense_SeqDense()
2892 PetscInt i, j, m = A->rmap->n, n = A->cmap->n, p; in MatGetRowMax_SeqDense()
2901 …PetscCheck(p == A->rmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Nonconforming matrix and vector"); in MatGetRowMax_SeqDense()
2920 PetscInt i, j, m = A->rmap->n, n = A->cmap->n, p; in MatGetRowMaxAbs_SeqDense()
2930 …PetscCheck(p == A->rmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Nonconforming matrix and vector"); in MatGetRowMaxAbs_SeqDense()
2949 PetscInt i, j, m = A->rmap->n, n = A->cmap->n, p; in MatGetRowMin_SeqDense()
2958 …PetscCheck(p == A->rmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Nonconforming matrix and vector"); in MatGetRowMin_SeqDense()
2984 PetscCall(PetscArraycpy(x, aa + col * a->lda, A->rmap->n)); in MatGetColumnVector_SeqDense()
3288 PetscCall(PetscLayoutSetUp(B->rmap)); in MatSeqDenseSetPreallocation_SeqDense()
3291 if (b->lda <= 0) PetscCall(PetscBLASIntCast(B->rmap->n, &b->lda)); in MatSeqDenseSetPreallocation_SeqDense()
3313 PetscInt M = A->rmap->N, N = A->cmap->N, i, j, k, *rows, *cols; in MatConvert_SeqDense_Elemental()
3353 data = (B->rmap->n > 0 && B->cmap->n > 0) ? (b->v ? PETSC_TRUE : PETSC_FALSE) : PETSC_FALSE; in MatDenseSetLDA_SeqDense()
3355 …(lda >= B->rmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "LDA %" PetscInt_FMT " must be at least ma… in MatDenseSetLDA_SeqDense()
3384 …all(VecCreateMPIWithArray(PetscObjectComm((PetscObject)A), A->rmap->bs, A->rmap->n, A->rmap->N, NU… in MatDenseCreateColumnVec_Private()
3385 …else PetscCall(VecCreateSeqWithArray(PetscObjectComm((PetscObject)A), A->rmap->bs, A->rmap->n, NUL… in MatDenseCreateColumnVec_Private()
3389 …cCreateMPIKokkosWithArray(PetscObjectComm((PetscObject)A), A->rmap->bs, A->rmap->n, A->rmap->N, NU… in MatDenseCreateColumnVec_Private()
3390 …all(VecCreateSeqKokkosWithArray(PetscObjectComm((PetscObject)A), A->rmap->bs, A->rmap->n, NULL, v)… in MatDenseCreateColumnVec_Private()
3396 …ecCreateMPICUDAWithArrays(PetscObjectComm((PetscObject)A), A->rmap->bs, A->rmap->n, A->rmap->N, NU… in MatDenseCreateColumnVec_Private()
3397 …else PetscCall(VecCreateSeqCUDAWithArrays(PetscObjectComm((PetscObject)A), A->rmap->bs, A->rmap->n… in MatDenseCreateColumnVec_Private()
3403 …VecCreateMPIHIPWithArrays(PetscObjectComm((PetscObject)A), A->rmap->bs, A->rmap->n, A->rmap->N, NU… in MatDenseCreateColumnVec_Private()
3404 …else PetscCall(VecCreateSeqHIPWithArrays(PetscObjectComm((PetscObject)A), A->rmap->bs, A->rmap->n,… in MatDenseCreateColumnVec_Private()
3510 …if (a->cmat && (cend - cbegin != a->cmat->cmap->N || rend - rbegin != a->cmat->rmap->N)) PetscCall… in MatDenseGetSubMatrix_SeqDense()
3893 if (rend == PETSC_DECIDE) rend = A->rmap->N; in MatDenseGetSubMatrix()
3897 …>rmap->N, PetscObjectComm((PetscObject)A), PETSC_ERR_ARG_WRONG, "Invalid rbegin %" PetscInt_FMT ",… in MatDenseGetSubMatrix()
3898 …rmap->N, PetscObjectComm((PetscObject)A), PETSC_ERR_ARG_WRONG, "Invalid rend %" PetscInt_FMT ", sh… in MatDenseGetSubMatrix()