seq/aijperm/aijperm.c

120   Mat_SeqAIJPERM *aijperm_dest;
127     aijperm_dest = (Mat_SeqAIJPERM *)(*M)->spptr;
128     PetscCall(PetscFree(aijperm_dest->xgroup));
129     PetscCall(PetscFree(aijperm_dest->nzgroup));
130     PetscCall(PetscFree(aijperm_dest->iperm));
132     PetscCall(PetscNew(&aijperm_dest));
133     (*M)->spptr = (void *)aijperm_dest;
137   PetscCall(PetscArraycpy(aijperm_dest, aijperm, 1));
144   PetscCall(PetscMalloc1(A->rmap->n, &aijperm_dest->iperm));
145   PetscCall(PetscMalloc1(aijperm->ngroup + 1, &aijperm_dest->xgroup));
146   PetscCall(PetscMalloc1(aijperm->ngroup, &aijperm_dest->nzgroup));
147   PetscCall(PetscArraycpy(aijperm_dest->iperm, aijperm->iperm, A->rmap->n));
148   PetscCall(PetscArraycpy(aijperm_dest->xgroup, aijperm->xgroup, aijperm->ngroup + 1));
149   PetscCall(PetscArraycpy(aijperm_dest->nzgroup, aijperm->nzgroup, aijperm->ngroup));