petsc-3.10.5 2019-03-28
Report Typos and Errors

# MatGetOrdering

Gets a reordering for a matrix to reduce fill or to improve numerical stability of LU factorization.

### Synopsis

    #include "petscmat.h"
    PetscErrorCode  MatGetOrdering(Mat mat,MatOrderingType type,IS *rperm,IS *cperm)

Collective on Mat

### Input Parameters

| mat | - the matrix
| type | - type of reordering, one of the following:

         MATORDERINGNATURAL - Natural

         MATORDERINGND - Nested Dissection

         MATORDERING1WD - One-way Dissection

         MATORDERINGRCM - Reverse Cuthill-McKee

         MATORDERINGQMD - Quotient Minimum Degree

### Output Parameters

| rperm | - row permutation indices
| cperm | - column permutation indices

### Options Database Key

-mat_view_ordering draw -plots matrix nonzero structure in new ordering

### Notes

This DOES NOT actually reorder the matrix; it merely returns two index sets that define a reordering. This is usually not used directly, rather use the options PCFactorSetMatOrderingType()

The user can define additional orderings; see MatOrderingRegister().

These are generally only implemented for sequential sparse matrices.

The external packages that PETSc can use for direct factorization such as SuperLU do not accept orderings provided by this call.

### Keywords

matrix, set, ordering, factorization, direct, ILU, LU,
fill, reordering, natural, Nested Dissection, One-way Dissection, Cholesky, Reverse Cuthill-McKee, Quotient Minimum Degree

### See Also

MatOrderingRegister(), PCFactorSetMatOrderingType()

### Level

intermediate

### Location

src/mat/order/sorder.c

### Examples

src/mat/examples/tutorials/ex1.c.html
src/ksp/ksp/examples/tutorials/ex18.c.html

### Implementations

MatGetOrdering_Flow_SeqAIJ in src/mat/impls/aij/seq/aijfact.c
Index of all MatOrderings routines
Table of Contents for all manual pages
Index of all manual pages