Searched refs:atfill (Results 1 – 2 of 2) sorted by relevance
| /petsc/src/mat/impls/aij/seq/ |
| H A D | symtranspose.c | 16 PetscInt *ati, *atj, *atfill, *ai = a->i, *aj = a->j; in MatTransposeSymbolic_SeqAIJ() local 30 PetscCall(PetscMalloc1(an, &atfill)); in MatTransposeSymbolic_SeqAIJ() 31 PetscCall(PetscArraycpy(atfill, ati, an)); in MatTransposeSymbolic_SeqAIJ() 37 atj[atfill[*aj]] = i; in MatTransposeSymbolic_SeqAIJ() 38 atfill[*aj++] += 1; in MatTransposeSymbolic_SeqAIJ() 41 PetscCall(PetscFree(atfill)); in MatTransposeSymbolic_SeqAIJ() 61 PetscInt *ati, *atj, *atfill, *ai = a->i, *aj = a->j; in MatTranspose_SeqAIJ() local 97 PetscCall(PetscMalloc1(an, &atfill)); in MatTranspose_SeqAIJ() 98 PetscCall(PetscArraycpy(atfill, ati, an)); in MatTranspose_SeqAIJ() 105 atj[atfill[*aj]] = i; in MatTranspose_SeqAIJ() [all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baij.c | 1647 …etscInt i, j, k, *aj = a->j, *ai = a->i, bs = A->rmap->bs, mbs = a->mbs, nbs = a->nbs, *atfill; in MatTranspose_SeqBAIJ() local 1653 PetscCall(PetscCalloc1(1 + nbs, &atfill)); in MatTranspose_SeqBAIJ() 1655 for (i = 0; i < ai[mbs]; i++) atfill[aj[i]] += 1; /* count num of non-zeros in row aj[i] */ in MatTranspose_SeqBAIJ() 1660 PetscCall(MatSeqBAIJSetPreallocation(C, bs, 0, atfill)); in MatTranspose_SeqBAIJ() 1675 PetscCall(PetscArraycpy(atfill, ati, nbs)); in MatTranspose_SeqBAIJ() 1681 atj[atfill[*aj]] = i; in MatTranspose_SeqBAIJ() 1683 for (k = 0; k < bs; k++) ata[bs2 * atfill[*aj] + k * bs + kr] = *aa++; in MatTranspose_SeqBAIJ() 1685 atfill[*aj++] += 1; in MatTranspose_SeqBAIJ() 1692 PetscCall(PetscFree(atfill)); in MatTranspose_SeqBAIJ()
|