seq/crl/crl.c

15   Mat_AIJCRL *aijcrl = (Mat_AIJCRL *)A->spptr;
19   if (aijcrl) PetscCall(PetscFree2(aijcrl->acols, aijcrl->icols));
34   Mat_AIJCRL  *aijcrl = (Mat_AIJCRL *)A->spptr;
42   aijcrl->nz   = a->nz;
43   aijcrl->m    = A->rmap->n;
44   aijcrl->rmax = rmax;
46   PetscCall(PetscFree2(aijcrl->acols, aijcrl->icols));
47   PetscCall(PetscMalloc2(rmax * m, &aijcrl->acols, rmax * m, &aijcrl->icols));
48   acols = aijcrl->acols;
49   icols = aijcrl->icols;
88   Mat_AIJCRL        *aijcrl = (Mat_AIJCRL *)A->spptr;
89   PetscInt           m      = aijcrl->m; /* Number of rows in the matrix. */
90   PetscInt           rmax = aijcrl->rmax, *icols = aijcrl->icols;
91   PetscScalar       *acols = aijcrl->acols;
103   if (aijcrl->xscat) {
104     PetscCall(VecCopy(xx, aijcrl->xwork));
106     PetscCall(VecScatterBegin(aijcrl->xscat, xx, aijcrl->fwork, INSERT_VALUES, SCATTER_FORWARD));
107     PetscCall(VecScatterEnd(aijcrl->xscat, xx, aijcrl->fwork, INSERT_VALUES, SCATTER_FORWARD));
108     xx = aijcrl->xwork;
133   PetscCall(PetscLogFlops(2.0 * aijcrl->nz - m));
146   Mat_AIJCRL *aijcrl;
154   PetscCall(PetscNew(&aijcrl));
155   B->spptr = (void *)aijcrl;