/*@C
MatCreateSeqAIJPERM - Creates a sparse matrix of type SEQAIJPERM.
This type inherits from AIJ, but calculates some additional permutation
information that is used to allow better vectorization of some
operations. At the cost of increased storage, the AIJ formatted
matrix can be copied to a format in which pieces of the matrix are
stored in ELLPACK format, allowing the vectorized matrix multiply
routine to use stride-1 memory accesses. As with the AIJ type, it is
important to preallocate matrix storage in order to get good assembly
performance.
Collective on MPI_Comm
Input Parameters:
+ comm - MPI communicator, set to PETSC_COMM_SELF
. m - number of rows
. n - number of columns
. nz - number of nonzeros per row (same for all rows)
- nnz - array containing the number of nonzeros in the various rows
(possibly different for each row) or NULL
Output Parameter:
. A - the matrix
Notes:
If nnz is given then nz is ignored
Level: intermediate
.keywords: matrix, cray, sparse, parallel
.seealso: MatCreate(), MatCreateMPIAIJPERM(), MatSetValues()
@*/
PetscErrorCode MatCreateSeqAIJPERM(MPI_Comm comm,PetscInt m,PetscInt n,PetscInt nz,const PetscInt nnz[],Mat *A)
{
PetscErrorCode ierr;
PetscFunctionBegin;
ierr = MatCreate(comm,A);CHKERRQ(ierr);
ierr = MatSetSizes(*A,m,n,m,n);CHKERRQ(ierr);
ierr = MatSetType(*A,MATSEQAIJPERM);CHKERRQ(ierr);
ierr = MatSeqAIJSetPreallocation_SeqAIJ(*A,nz,nnz);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
EXTERN_C_END
#undef __FUNCT__
#define __FUNCT__ "MatCreateSeqCRL"
/*@C
MatCreateSeqCRL - Creates a sparse matrix of type SEQCRL.
This type inherits from AIJ, but stores some additional
information that is used to allow better vectorization of
the matrix-vector product. At the cost of increased storage, the AIJ formatted
matrix can be copied to a format in which pieces of the matrix are
stored in ELLPACK format, allowing the vectorized matrix multiply
routine to use stride-1 memory accesses. As with the AIJ type, it is
important to preallocate matrix storage in order to get good assembly
performance.
Collective on MPI_Comm
Input Parameters:
+ comm - MPI communicator, set to PETSC_COMM_SELF
. m - number of rows
. n - number of columns
. nz - number of nonzeros per row (same for all rows)
- nnz - array containing the number of nonzeros in the various rows
(possibly different for each row) or PETSC_NULL
Output Parameter:
. A - the matrix
Notes:
If nnz is given then nz is ignored
Level: intermediate
.keywords: matrix, cray, sparse, parallel
.seealso: MatCreate(), MatCreateMPICSRPERM(), MatSetValues()
@*/
PetscErrorCode PETSCMAT_DLLEXPORT MatCreateSeqCRL(MPI_Comm comm,PetscInt m,PetscInt n,PetscInt nz,const PetscInt nnz[],Mat *A)
{
PetscErrorCode ierr;
PetscFunctionBegin;
ierr = MatCreate(comm,A);CHKERRQ(ierr);
ierr = MatSetSizes(*A,m,n,m,n);CHKERRQ(ierr);
ierr = MatSetType(*A,MATSEQCRL);CHKERRQ(ierr);
ierr = MatSeqAIJSetPreallocation_SeqAIJ(*A,nz,(PetscInt*)nnz);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
