SLESSetOperators

Sets the matrix associated with the linear system and a (possibly) different one associated with the preconditioner.

Synopsis

#include "petscsles.h"    
int SLESSetOperators(SLES sles,Mat Amat,Mat Pmat,MatStructure flag)
Collective on SLES and Mat

Input Parameters

sles - the SLES context
Amat - the matrix associated with the linear system
Pmat - the matrix to be used in constructing the preconditioner, usually the same as Amat.
flag - flag indicating information about the preconditioner matrix structure during successive linear solves. This flag is ignored the first time a linear system is solved, and thus is irrelevant when solving just one linear system.

Notes

The flag can be used to eliminate unnecessary work in the preconditioner during the repeated solution of linear systems of the same size. The available options are
   SAME_PRECONDITIONER -
     Pmat is identical during successive linear solves.
     This option is intended for folks who are using
     different Amat and Pmat matrices and want to reuse the
     same preconditioner matrix.  For example, this option
     saves work by not recomputing incomplete factorization
     for ILU/ICC preconditioners.
   SAME_NONZERO_PATTERN -
     Pmat has the same nonzero structure during
     successive linear solves. 
   DIFFERENT_NONZERO_PATTERN -
     Pmat does not have the same nonzero structure.

Caution

If you specify SAME_NONZERO_PATTERN, PETSc believes your assertion and does not check the structure of the matrix. If you erroneously claim that the structure is the same when it actually is not, the new preconditioner will not function correctly. Thus, use this optimization feature carefully!

If in doubt about whether your preconditioner matrix has changed structure or not, use the flag DIFFERENT_NONZERO_PATTERN.

Keywords

SLES, set, operators, matrix, preconditioner, linear system

See Also

SLESSolve(), SLESGetPC(), PCGetOperators()

Level:beginner
Location:
src/sles/interface/sles.c
Index of all SLES routines
Table of Contents for all manual pages
Index of all manual pages

Examples

src/sles/examples/tutorials/ex1.c.html
src/sles/examples/tutorials/ex2.c.html
src/sles/examples/tutorials/ex3.c.html
src/sles/examples/tutorials/ex4.c.html
src/sles/examples/tutorials/ex5.c.html
src/sles/examples/tutorials/ex7.c.html
src/sles/examples/tutorials/ex8.c.html
src/sles/examples/tutorials/ex9.c.html
src/sles/examples/tutorials/ex10.c.html
src/sles/examples/tutorials/ex11.c.html
src/sles/examples/tutorials/ex12.c.html