Lines Matching refs:cmap
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()
42 PetscCall(PetscBLASIntCast(A->cmap->n, &n)); in MatSeqDenseInvertFactors_Private()
52 PetscCall(PetscLogFlops((1.0 * A->cmap->n * A->cmap->n * A->cmap->n) / 3.0)); in MatSeqDenseInvertFactors_Private()
77 PetscCall(PetscLogFlops((1.0 * A->cmap->n * A->cmap->n * A->cmap->n) / 3.0)); in MatSeqDenseInvertFactors_Private()
94 PetscInt m = l->lda, n = A->cmap->n, r = A->rmap->n, i, j; in MatZeroRowsColumns_SeqDense()
103 …cmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Col %" 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()
208 PetscCall(MatSetSizes(B, A->rmap->n, A->cmap->n, A->rmap->N, A->cmap->N)); in MatConvert_SeqDense_SeqAIJ()
210 for (j = 0; j < A->cmap->n; j++) { in MatConvert_SeqDense_SeqAIJ()
212 if (aa[i] != 0.0 || (i == j && A->cmap->n == A->rmap->n)) ++nnz[i]; in MatConvert_SeqDense_SeqAIJ()
218 for (j = 0; j < A->cmap->n; j++) { 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()
252 …for (PetscInt j = 0; j < X->cmap->n; j++) PetscCallBLAS("BLASaxpy", BLASaxpy_(&m, &alpha, PetscSaf… in MatAXPY_SeqDense()
264 PetscLogDouble N = A->rmap->n * A->cmap->n; in MatGetInfo_SeqDense()
291 …for (j = 0; j < A->cmap->n; j++) PetscCallBLAS("BLASscal", BLASscal_(&nz, &alpha, v + j * lda, &on… 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()
324 if (A->rmap->n != A->cmap->n) PetscFunctionReturn(PETSC_SUCCESS); in MatIsHermitian_SeqDense()
345 if (A->rmap->n != A->cmap->n) PetscFunctionReturn(PETSC_SUCCESS); in MatIsSymmetric_SeqDense()
366 PetscCall(PetscLayoutReference(A->cmap, &newi->cmap)); in MatDuplicateNoCreate_SeqDense()
381 …for (j = 0; j < A->cmap->n; j++) PetscCall(PetscArraycpy(PetscSafePointerPlusOffset(v, j * nlda), … 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()
529 PetscCall(PetscBLASIntCast(A->cmap->n, &k)); in MatSolve_SeqDense_SetUp()
652 PetscCall(PetscBLASIntCast(A->cmap->n, &k)); in MatMatSolve_SeqDense_SetUp()
795 PetscCall(PetscBLASIntCast(A->cmap->n, &n)); in MatLUFactor_SeqDense()
798 if (!A->rmap->n || !A->cmap->n) PetscFunctionReturn(PETSC_SUCCESS); in MatLUFactor_SeqDense()
815 PetscCall(PetscLogFlops((2.0 * A->cmap->n * A->cmap->n * A->cmap->n) / 3)); in MatLUFactor_SeqDense()
843 PetscCall(PetscBLASIntCast(A->cmap->n, &n)); in MatCholeskyFactor_SeqDense()
844 if (!A->rmap->n || !A->cmap->n) PetscFunctionReturn(PETSC_SUCCESS); in MatCholeskyFactor_SeqDense()
893 PetscCall(PetscLogFlops((1.0 * A->cmap->n * A->cmap->n * A->cmap->n) / 3.0)); in MatCholeskyFactor_SeqDense()
920 PetscCall(PetscBLASIntCast(A->cmap->n, &n)); 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()
1087 PetscCall(MatMultColumnRangeKernel_SeqDense(A, xx, yy, 0, A->cmap->n, PETSC_FALSE, PETSC_FALSE)); in MatMult_SeqDense()
1094 PetscCall(MatMultColumnRangeKernel_SeqDense(A, xx, yy, 0, A->cmap->n, PETSC_TRUE, PETSC_FALSE)); in MatMultTranspose_SeqDense()
1101 PetscCall(MatMultColumnRangeKernel_SeqDense(A, xx, yy, 0, A->cmap->n, PETSC_TRUE, PETSC_TRUE)); in MatMultHermitianTranspose_SeqDense()
1157 …PetscCall(MatMultAddColumnRangeKernel_SeqDense(A, xx, zz, yy, 0, A->cmap->n, PETSC_FALSE, PETSC_FA… in MatMultAdd_SeqDense()
1164 …PetscCall(MatMultAddColumnRangeKernel_SeqDense(A, xx, zz, yy, 0, A->cmap->n, PETSC_TRUE, PETSC_FAL… in MatMultTransposeAdd_SeqDense()
1171 …PetscCall(MatMultAddColumnRangeKernel_SeqDense(A, xx, zz, yy, 0, A->cmap->n, PETSC_TRUE, PETSC_TRU… in MatMultHermitianTransposeAdd_SeqDense()
1181 if (ncols) *ncols = A->cmap->n; in MatGetRow_SeqDense()
1183 PetscCall(PetscMalloc1(A->cmap->n, cols)); in MatGetRow_SeqDense()
1184 for (i = 0; i < A->cmap->n; i++) (*cols)[i] = i; in MatGetRow_SeqDense()
1190 PetscCall(PetscMalloc1(A->cmap->n, vals)); in MatGetRow_SeqDense()
1192 for (i = 0; i < A->cmap->n; i++) { in MatGetRow_SeqDense()
1227 …n[j] < A->cmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Column too large: col %" PetscInt_F… in MatSetValues_SeqDense()
1243 …n[j] < A->cmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Column too large: col %" PetscInt_F… in MatSetValues_SeqDense()
1267 …n[j] < A->cmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Column too large: col %" PetscInt_F… in MatSetValues_SeqDense()
1283 …n[j] < A->cmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Column too large: col %" PetscInt_F… in MatSetValues_SeqDense()
1321 …>cmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Column %" PetscInt_FMT " requested larger th… in MatGetValues_SeqDense()
1406 if (mat->cmap->N < 0) mat->cmap->N = N; in MatLoad_Dense_Binary()
1493 for (j = 0; j < A->cmap->n; j++) { in MatView_SeqDense_ASCII()
1512 for (j = 0; j < A->cmap->n; j++) { 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()
1531 for (j = 0; j < A->cmap->n; j++) { in MatView_SeqDense_ASCII()
1557 PetscInt m = A->rmap->n, n = A->cmap->n, i, j; in MatView_SeqDense_Draw_Zoom()
1628 xr = A->cmap->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()
1822 A->rmap = A->cmap; in MatTranspose_SeqDense()
1823 A->cmap = tmplayout; in MatTranspose_SeqDense()
1833 PetscCall(MatSetSizes(tmat, A->cmap->n, A->rmap->n, A->cmap->n, A->rmap->n)); in MatTranspose_SeqDense()
1866 if (A1->cmap->n != A2->cmap->n) { in MatEqual_SeqDense()
1872 for (i = 0; i < A1->cmap->n; i++) { in MatEqual_SeqDense()
1894 len = PetscMin(A->rmap->n, A->cmap->n); in MatGetDiagonal_SeqDense()
1908 PetscInt i, j, m = A->rmap->n, n = A->cmap->n; in MatDiagonalScale_SeqDense()
1930 PetscCheck(n == A->cmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Right scaling vec wrong size"); in MatDiagonalScale_SeqDense()
1956 for (j = 0; j < A->cmap->n; 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()
1979 for (j = 0; j < A->cmap->n; j++) {
1988 PetscCall(PetscLogFlops(1.0 * A->cmap->n * A->rmap->n));
1994 for (i = 0; i < A->cmap->n; i++) {
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()
2041 PetscInt m = l->lda, n = A->cmap->n, i, j; 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()
2647 PetscCall(PetscLayoutSetUp(A->cmap)); in MatSetUp_SeqDense()
2656 PetscInt min = PetscMin(A->rmap->n, A->cmap->n); in MatConjugate_SeqDense()
2661 for (j = 0; j < A->cmap->n; j++) in MatConjugate_SeqDense()
2677 for (j = 0; j < A->cmap->n; j++) { in MatRealPart_SeqDense()
2692 for (j = 0; j < A->cmap->n; j++) { in MatImaginaryPart_SeqDense()
2701 PetscInt m = A->rmap->n, n = B->cmap->n; in MatMatMultSymbolic_SeqDense_SeqDense()
2732 PetscCall(PetscBLASIntCast(C->cmap->n, &n)); in MatMatMultNumeric_SeqDense_SeqDense()
2733 PetscCall(PetscBLASIntCast(A->cmap->n, &k)); in MatMatMultNumeric_SeqDense_SeqDense()
2781 PetscCall(PetscBLASIntCast(C->cmap->n, &n)); in MatMatTransposeMultNumeric_SeqDense_SeqDense()
2782 PetscCall(PetscBLASIntCast(A->cmap->n, &k)); in MatMatTransposeMultNumeric_SeqDense_SeqDense()
2797 PetscInt m = A->cmap->n, n = B->cmap->n; in MatTransposeMatMultSymbolic_SeqDense_SeqDense()
2830 PetscCall(PetscBLASIntCast(C->cmap->n, &n)); in MatTransposeMatMultNumeric_SeqDense_SeqDense()
2892 PetscInt i, j, m = A->rmap->n, n = A->cmap->n, p; in MatGetRowMax_SeqDense()
2920 PetscInt i, j, m = A->rmap->n, n = A->cmap->n, p; in MatGetRowMaxAbs_SeqDense()
2949 PetscInt i, j, m = A->rmap->n, n = A->cmap->n, p; in MatGetRowMin_SeqDense()
3289 PetscCall(PetscLayoutSetUp(B->cmap)); in MatSeqDenseSetPreallocation_SeqDense()
3295 PetscCall(PetscCalloc1((size_t)b->lda * B->cmap->n, &b->v)); 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()
3510 …if (a->cmat && (cend - cbegin != a->cmat->cmap->N || rend - rbegin != a->cmat->rmap->N)) PetscCall… in MatDenseGetSubMatrix_SeqDense()
3702 … A->cmap->N, PetscObjectComm((PetscObject)A), PETSC_ERR_ARG_WRONG, "Invalid col %" PetscInt_FMT ",… in MatDenseGetColumnVec()
3729 … A->cmap->N, PetscObjectComm((PetscObject)A), PETSC_ERR_ARG_WRONG, "Invalid col %" PetscInt_FMT ",… in MatDenseRestoreColumnVec()
3765 … A->cmap->N, PetscObjectComm((PetscObject)A), PETSC_ERR_ARG_WRONG, "Invalid col %" PetscInt_FMT ",… in MatDenseGetColumnVecRead()
3792 … A->cmap->N, PetscObjectComm((PetscObject)A), PETSC_ERR_ARG_WRONG, "Invalid col %" PetscInt_FMT ",… in MatDenseRestoreColumnVecRead()
3826 … A->cmap->N, PetscObjectComm((PetscObject)A), PETSC_ERR_ARG_WRONG, "Invalid col %" PetscInt_FMT ",… in MatDenseGetColumnVecWrite()
3853 … A->cmap->N, PetscObjectComm((PetscObject)A), PETSC_ERR_ARG_WRONG, "Invalid col %" PetscInt_FMT ",… in MatDenseRestoreColumnVecWrite()
3895 if (cend == PETSC_DECIDE) cend = A->cmap->N; in MatDenseGetSubMatrix()
3899 …>cmap->N, PetscObjectComm((PetscObject)A), PETSC_ERR_ARG_WRONG, "Invalid cbegin %" PetscInt_FMT ",… in MatDenseGetSubMatrix()
3900 …cmap->N, PetscObjectComm((PetscObject)A), PETSC_ERR_ARG_WRONG, "Invalid cend %" PetscInt_FMT ", sh… in MatDenseGetSubMatrix()