baij/mpi/mpiaijbaij.c

10   Mat         M;
11   Mat_MPIAIJ *mpimat = (Mat_MPIAIJ *)A->data;
13   PetscInt    m, n, lm, ln, bs = A->rmap->bs;
18     PetscBool3 sym = A->symmetric, hermitian = A->hermitian, structurally_symmetric = A->structurally_symmetric, spd = A->spd;
20     PetscCall(MatGetSize(A, &m, &n));
22     PetscCall(MatConvert_SeqAIJ_SeqBAIJ_Preallocate(mpimat->A, &d_nnz));
23     PetscCall(MatConvert_SeqAIJ_SeqBAIJ_Preallocate(mpimat->B, &o_nnz));
24     PetscCall(MatCreate(PetscObjectComm((PetscObject)A), &M));
25     PetscCall(MatSetSizes(M, lm, ln, m, n));
26     PetscCall(MatSetType(M, MATMPIBAIJ));
27     PetscCall(MatSeqBAIJSetPreallocation(M, bs, 0, d_nnz));
28     PetscCall(MatMPIBAIJSetPreallocation(M, bs, 0, d_nnz, 0, o_nnz));
35     A->symmetric              = sym;
36     A->hermitian              = hermitian;
37     A->structurally_symmetric = structurally_symmetric;
38     A->spd                    = spd;
40     M = *newmat;
41     PetscCall(MatGetOption(M, MAT_NO_OFF_PROC_ENTRIES, &flg));
47   PetscCall(MatSetOption(M, MAT_NO_OFF_PROC_ENTRIES, PETSC_TRUE));
48   PetscCall(MatConvert_Basic(A, newtype, MAT_REUSE_MATRIX, &M));
49   PetscCall(MatSetOption(M, MAT_NO_OFF_PROC_ENTRIES, flg));
51   if (reuse == MAT_INPLACE_MATRIX) PetscCall(MatHeaderReplace(A, &M));
52   else *newmat = M;