Lines Matching refs:aij
37 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)A->data; in MatGetColumnReductions_SeqAIJ() local
43 for (i = 0; i < aij->i[m]; i++) reductions[aij->j[i]] += PetscAbsScalar(aij->a[i] * aij->a[i]); in MatGetColumnReductions_SeqAIJ()
45 for (i = 0; i < aij->i[m]; i++) reductions[aij->j[i]] += PetscAbsScalar(aij->a[i]); in MatGetColumnReductions_SeqAIJ()
47 …for (i = 0; i < aij->i[m]; i++) reductions[aij->j[i]] = PetscMax(PetscAbsScalar(aij->a[i]), reduct… in MatGetColumnReductions_SeqAIJ()
49 for (i = 0; i < aij->i[m]; i++) reductions[aij->j[i]] += PetscRealPart(aij->a[i]); in MatGetColumnReductions_SeqAIJ()
51 for (i = 0; i < aij->i[m]; i++) reductions[aij->j[i]] += PetscImaginaryPart(aij->a[i]); in MatGetColumnReductions_SeqAIJ()
2248 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)A->data, *bij = (Mat_SeqAIJ *)B->data; in MatIsTranspose_SeqAIJ() local
2262 aii = aij->i; in MatIsTranspose_SeqAIJ()
2264 adx = aij->j; in MatIsTranspose_SeqAIJ()
2300 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)A->data, *bij = (Mat_SeqAIJ *)B->data; in MatIsHermitianTranspose_SeqAIJ() local
2312 aii = aij->i; in MatIsHermitianTranspose_SeqAIJ()
2314 adx = aij->j; in MatIsHermitianTranspose_SeqAIJ()
2316 va = aij->a; in MatIsHermitianTranspose_SeqAIJ()
2992 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)mat->data; in MatConjugate_SeqAIJ() local
2993 PetscInt i, nz = aij->nz; in MatConjugate_SeqAIJ()
3351 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)x->data; in MatSetRandom_SeqAIJ() local
3359 for (j = 0; j < aij->imax[i]; j++) { in MatSetRandom_SeqAIJ()
3367 for (i = 0; i < aij->nz; i++) PetscCall(PetscRandomGetValue(rctx, aa + i)); in MatSetRandom_SeqAIJ()
3378 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)x->data; in MatSetRandomSkipColumnRange_SeqAIJ_Private() local
3387 for (j = 0; j < aij->imax[i]; j++) { in MatSetRandomSkipColumnRange_SeqAIJ_Private()
3546 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)mat->data; in MatSeqAIJSetColumnIndices_SeqAIJ() local
3550 nz = aij->maxnz; in MatSeqAIJSetColumnIndices_SeqAIJ()
3552 for (i = 0; i < nz; i++) aij->j[i] = indices[i]; in MatSeqAIJSetColumnIndices_SeqAIJ()
3553 aij->nz = nz; in MatSeqAIJSetColumnIndices_SeqAIJ()
3554 for (i = 0; i < n; i++) aij->ilen[i] = aij->imax[i]; in MatSeqAIJSetColumnIndices_SeqAIJ()
3564 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)mat->data; in MatSeqAIJCompactOutExtraColumns_SeqAIJ() local
3567 PetscInt gid, lid, i, ec, nz = aij->nz; in MatSeqAIJCompactOutExtraColumns_SeqAIJ()
3568 PetscInt *garray, *jj = aij->j; in MatSeqAIJCompactOutExtraColumns_SeqAIJ()
3648 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)mat->data; in MatStoreValues_SeqAIJ() local
3649 size_t nz = aij->i[mat->rmap->n]; in MatStoreValues_SeqAIJ()
3652 …PetscCheck(aij->nonew, PETSC_COMM_SELF, PETSC_ERR_ORDER, "Must call MatSetOption(A,MAT_NEW_NONZERO… in MatStoreValues_SeqAIJ()
3655 if (!aij->saved_values) PetscCall(PetscMalloc1(nz + 1, &aij->saved_values)); in MatStoreValues_SeqAIJ()
3658 PetscCall(PetscArraycpy(aij->saved_values, aij->a, nz)); in MatStoreValues_SeqAIJ()
3722 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)mat->data; in MatRetrieveValues_SeqAIJ() local
3723 PetscInt nz = aij->i[mat->rmap->n]; in MatRetrieveValues_SeqAIJ()
3726 …PetscCheck(aij->nonew, PETSC_COMM_SELF, PETSC_ERR_ORDER, "Must call MatSetOption(A,MAT_NEW_NONZERO… in MatRetrieveValues_SeqAIJ()
3727 …PetscCheck(aij->saved_values, PETSC_COMM_SELF, PETSC_ERR_ORDER, "Must call MatStoreValues(A);first… in MatRetrieveValues_SeqAIJ()
3729 PetscCall(PetscArraycpy(aij->a, aij->saved_values, nz)); in MatRetrieveValues_SeqAIJ()
4325 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)A->data; in MatSeqAIJGetArray() local
4328 if (aij->ops->getarray) { in MatSeqAIJGetArray()
4329 PetscCall((*aij->ops->getarray)(A, array)); in MatSeqAIJGetArray()
4331 *array = aij->a; in MatSeqAIJGetArray()
4351 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)A->data; in MatSeqAIJRestoreArray() local
4354 if (aij->ops->restorearray) { in MatSeqAIJRestoreArray()
4355 PetscCall((*aij->ops->restorearray)(A, array)); in MatSeqAIJRestoreArray()
4380 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)A->data; in MatSeqAIJGetArrayRead() local
4383 if (aij->ops->getarrayread) { in MatSeqAIJGetArrayRead()
4384 PetscCall((*aij->ops->getarrayread)(A, array)); in MatSeqAIJGetArrayRead()
4386 *array = aij->a; in MatSeqAIJGetArrayRead()
4408 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)A->data; in MatSeqAIJRestoreArrayRead() local
4411 if (aij->ops->restorearrayread) { in MatSeqAIJRestoreArrayRead()
4412 PetscCall((*aij->ops->restorearrayread)(A, array)); in MatSeqAIJRestoreArrayRead()
4436 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)A->data; in MatSeqAIJGetArrayWrite() local
4439 if (aij->ops->getarraywrite) { in MatSeqAIJGetArrayWrite()
4440 PetscCall((*aij->ops->getarraywrite)(A, array)); in MatSeqAIJGetArrayWrite()
4442 *array = aij->a; in MatSeqAIJGetArrayWrite()
4465 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)A->data; in MatSeqAIJRestoreArrayWrite() local
4468 if (aij->ops->restorearraywrite) { in MatSeqAIJRestoreArrayWrite()
4469 PetscCall((*aij->ops->restorearraywrite)(A, array)); in MatSeqAIJRestoreArrayWrite()
4503 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)mat->data; in MatSeqAIJGetCSRAndMemType() local
4507 if (aij->ops->getcsrandmemtype) { in MatSeqAIJGetCSRAndMemType()
4508 PetscCall((*aij->ops->getcsrandmemtype)(mat, i, j, a, mtype)); in MatSeqAIJGetCSRAndMemType()
4510 if (i) *i = aij->i; in MatSeqAIJGetCSRAndMemType()
4511 if (j) *j = aij->j; in MatSeqAIJGetCSRAndMemType()
4512 if (a) *a = aij->a; in MatSeqAIJGetCSRAndMemType()
4535 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)A->data; in MatSeqAIJGetMaxRowNonzeros() local
4538 *nz = aij->rmax; in MatSeqAIJGetMaxRowNonzeros()
5128 Mat_SeqAIJ *aij; in MatCreateSeqAIJWithArrays() local
5138 aij = (Mat_SeqAIJ *)(*mat)->data; in MatCreateSeqAIJWithArrays()
5139 PetscCall(PetscMalloc1(m, &aij->imax)); in MatCreateSeqAIJWithArrays()
5140 PetscCall(PetscMalloc1(m, &aij->ilen)); in MatCreateSeqAIJWithArrays()
5142 aij->i = i; in MatCreateSeqAIJWithArrays()
5143 aij->j = j; in MatCreateSeqAIJWithArrays()
5144 aij->a = a; in MatCreateSeqAIJWithArrays()
5145 …aij->nonew = -1; /*this indicates that inserting a new value in the matrix that generates a new … in MatCreateSeqAIJWithArrays()
5146 aij->free_a = PETSC_FALSE; in MatCreateSeqAIJWithArrays()
5147 aij->free_ij = PETSC_FALSE; in MatCreateSeqAIJWithArrays()
5149 for (ii = 0, aij->nonzerorowcnt = 0, aij->rmax = 0; ii < m; ii++) { in MatCreateSeqAIJWithArrays()
5150 aij->ilen[ii] = aij->imax[ii] = i[ii + 1] - i[ii]; in MatCreateSeqAIJWithArrays()
5160 for (ii = 0; ii < aij->i[m]; ii++) { in MatCreateSeqAIJWithArrays()