MATSEQSBAIJ#

MATSEQSBAIJ = “seqsbaij” - A matrix type to be used for sequential symmetric block sparse matrices, based on block compressed sparse row format. Only the upper triangular portion of the matrix is stored. For complex numbers by default this matrix is symmetric, NOT Hermitian symmetric. To make it Hermitian symmetric you can call MatSetOption(A, MAT_HERMITIAN, PETSC_TRUE).

Options Database Key#

-mat_type seqsbaij - sets the matrix type to “seqsbaij” during a call to MatSetFromOptions()

Notes#

Call MatSetOption(A, MAT_STRUCTURE_ONLY, PETSC_TRUE) before preallocation or MatSetUp() to store only the nonzero pattern. The assembled matrix has no numerical value array. Row and column indices supplied during insertion are retained, while numerical values are ignored. Such matrices can be used for structural operations, but not for numerical operations.

By default if you insert values into the lower triangular part of the matrix they are simply ignored (since they are not stored and they are assumed to be symmetric to the upper triangular). If you call MatSetOption(A, MAT_IGNORE_LOWER_TRIANGULAR, PETSC_FALSE) or use the options database -mat_ignore_lower_triangular false it will generate an error if you try to set a value in the lower triangular portion.

The number of rows in the matrix must be less than or equal to the number of columns

Level#

beginner

Location#

src/mat/impls/sbaij/seq/sbaij.c

Examples#

src/mat/tests/ex74.c
src/mat/tests/ex128.c
src/mat/tests/ex141.c
src/mat/tests/ex174.cxx
src/mat/tests/ex182.c
src/mat/tests/ex192.c
src/mat/tests/ex232.c
src/mat/tests/ex234.c
src/mat/tests/ex237.c
src/mat/tests/ex244.cxx

Index of all Mat routines
Table of Contents for all manual pages
Index of all manual pages