petsc-3.5.4 2015-05-23
MatCreateSeqAIJCRL
Creates a sparse matrix of type SEQAIJCRL. 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.
Synopsis
PetscErrorCode MatCreateSeqAIJCRL(MPI_Comm comm,PetscInt m,PetscInt n,PetscInt nz,const PetscInt nnz[],Mat *A)
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
Keywords
matrix, cray, sparse, parallel
See Also
MatCreate(), MatCreateMPIAIJPERM(), MatSetValues()
Level:intermediate
Location:src/mat/impls/aij/seq/crl/crl.c
Index of all Mat routines
Table of Contents for all manual pages
Index of all manual pages