seq/aijsell/aijsell.c

22   Mat_SeqAIJSELL *aijsell = (Mat_SeqAIJSELL *)A->spptr;
39   if (reuse == MAT_INITIAL_MATRIX) aijsell = (Mat_SeqAIJSELL *)B->spptr;
43   PetscCall(MatDestroy(&aijsell->S));
55   Mat_SeqAIJSELL *aijsell = (Mat_SeqAIJSELL *)A->spptr;
60   if (aijsell) {
62     PetscCall(MatDestroy(&aijsell->S));
81   Mat_SeqAIJSELL  *aijsell = (Mat_SeqAIJSELL *)A->spptr;
86   if (aijsell->S && aijsell->state == state) {
92   if (aijsell->S) {
93     PetscCall(MatConvert_SeqAIJ_SeqSELL(A, MATSEQSELL, MAT_REUSE_MATRIX, &aijsell->S));
95     PetscCall(MatConvert_SeqAIJ_SeqSELL(A, MATSEQSELL, MAT_INITIAL_MATRIX, &aijsell->S));
100   PetscCall(PetscObjectStateGet((PetscObject)A, &aijsell->state));
106   Mat_SeqAIJSELL *aijsell;
111   aijsell      = (Mat_SeqAIJSELL *)A->spptr;
113   PetscCall(PetscArraycpy(aijsell_dest, aijsell, 1));
116   if (aijsell->eager_shadow) PetscCall(MatSeqAIJSELL_build_shadow(A));
123   Mat_SeqAIJSELL *aijsell = (Mat_SeqAIJSELL *)A->spptr;
143   if (aijsell->eager_shadow) PetscCall(MatSeqAIJSELL_build_shadow(A));
149   Mat_SeqAIJSELL *aijsell = (Mat_SeqAIJSELL *)A->spptr;
153   PetscCall(MatMult_SeqSELL(aijsell->S, xx, yy));
159   Mat_SeqAIJSELL *aijsell = (Mat_SeqAIJSELL *)A->spptr;
163   PetscCall(MatMultTranspose_SeqSELL(aijsell->S, xx, yy));
169   Mat_SeqAIJSELL *aijsell = (Mat_SeqAIJSELL *)A->spptr;
173   PetscCall(MatMultAdd_SeqSELL(aijsell->S, xx, yy, zz));
179   Mat_SeqAIJSELL *aijsell = (Mat_SeqAIJSELL *)A->spptr;
183   PetscCall(MatMultTransposeAdd_SeqSELL(aijsell->S, xx, yy, zz));
189   Mat_SeqAIJSELL *aijsell = (Mat_SeqAIJSELL *)A->spptr;
193   PetscCall(MatSOR_SeqSELL(aijsell->S, bb, omega, flag, fshift, its, lits, xx));
205   Mat_SeqAIJSELL *aijsell;
215   PetscCall(PetscNew(&aijsell));
217   B->spptr = (void *)aijsell;
230   aijsell->S            = NULL;
231   aijsell->eager_shadow = PETSC_FALSE;
235   PetscCall(PetscOptionsBool("-mat_aijsell_eager_shadow", "Eager Shadowing", "None", aijsell->eager_shadow, &aijsell->eager_shadow, &set));
239   if (A->assembled && aijsell->eager_shadow) PetscCall(MatSeqAIJSELL_build_shadow(A));