1*b4fd4287SBarry Smith /* $Id: pc.h,v 1.40 1995/12/21 16:39:01 bsmith Exp bsmith $ */ 2d03aef70SBarry Smith 3d03aef70SBarry Smith /* 4f26ada1bSBarry Smith Preconditioner module. Defines the preconditioner routines. 5d03aef70SBarry Smith */ 6d03aef70SBarry Smith #if !defined(__PC_PACKAGE) 7d03aef70SBarry Smith #define __PC_PACKAGE 8d03aef70SBarry Smith #include "petsc.h" 9d03aef70SBarry Smith #include "mat.h" 10d03aef70SBarry Smith 11070d5079SBarry Smith typedef enum { PCNONE, PCJACOBI, PCSOR, PCLU, PCSHELL, PCBJACOBI, PCMG, 1221e95762SBarry Smith PCEISENSTAT, PCILU, PCICC, PCASM, PCBGS } PCMethod; 13d03aef70SBarry Smith 14d03aef70SBarry Smith typedef struct _PC* PC; 159e25ed09SBarry Smith #define PC_COOKIE PETSC_COOKIE+9 16d03aef70SBarry Smith 1721e95762SBarry Smith typedef struct _PCNullSpace* PCNullSpace; 1821e95762SBarry Smith #define PCNULLSPACE_COOKIE PETSC_COOKIE+17 1921e95762SBarry Smith 2021e95762SBarry Smith extern int PCNullSpaceCreate(MPI_Comm,int,int,Vec *,PCNullSpace*); 2121e95762SBarry Smith extern int PCNullSpaceDestroy(PCNullSpace); 2221e95762SBarry Smith extern int PCNullSpaceRemove(PCNullSpace,Vec); 2321e95762SBarry Smith 246b5873e3SBarry Smith extern int PCCreate(MPI_Comm,PC*); 25ca9b4cbeSLois Curfman McInnes extern int PCSetMethod(PC,PCMethod); 268ed539a5SBarry Smith extern int PCSetUp(PC); 27aabeff55SBarry Smith extern int PCApply(PC,Vec,Vec); 28aabeff55SBarry Smith extern int PCApplyBAorAB(PC,int,Vec,Vec,Vec); 29aabeff55SBarry Smith extern int PCApplyTrans(PC,Vec,Vec); 30aabeff55SBarry Smith extern int PCApplyBAorABTrans(PC,int,Vec,Vec,Vec); 31aabeff55SBarry Smith extern int PCApplyRichardson(PC,Vec,Vec,Vec,int); 328ed539a5SBarry Smith extern int PCApplyRichardsonExists(PC); 338ed539a5SBarry Smith extern int PCRegisterAll(); 34d6dfbf8fSBarry Smith extern int PCRegisterDestroy(); 35ca9b4cbeSLois Curfman McInnes extern int PCRegister(PCMethod,char *,int (*)(PC)); 368ed539a5SBarry Smith extern int PCDestroy(PC); 378ed539a5SBarry Smith extern int PCSetFromOptions(PC); 38ca9b4cbeSLois Curfman McInnes extern int PCGetMethodFromContext(PC,PCMethod*); 39ca9b4cbeSLois Curfman McInnes extern int PCGetMethodName(PCMethod,char **); 40291315efSLois Curfman McInnes extern int PCGetFactoredMatrix(PC,Mat*); 415b116368SBarry Smith 425b116368SBarry Smith /* Flags for PCSetOperators */ 4396e03844SLois Curfman McInnes typedef enum {ALLMAT_DIFFERENT_NONZERO_PATTERN=0,MAT_SAME_NONZERO_PATTERN=1, 44f26ada1bSBarry Smith PMAT_SAME_NONZERO_PATTERN=2,ALLMAT_SAME_NONZERO_PATTERN=3} MatStructure; 455b116368SBarry Smith 4696e03844SLois Curfman McInnes extern int PCSetOperators(PC,Mat,Mat,MatStructure); 478ed539a5SBarry Smith extern int PCSetVector(PC,Vec); 488ed539a5SBarry Smith extern int PCPrintHelp(PC); 49023758ceSLois Curfman McInnes extern int PCView(PC,Viewer); 501eb62cbbSBarry Smith extern int PCSetOptionsPrefix(PC,char*); 518ed539a5SBarry Smith 52ddda3453SLois Curfman McInnes extern int PCSORSetSymmetric(PC, MatSORType); 538ed539a5SBarry Smith extern int PCSORSetOmega(PC, double); 545b116368SBarry Smith extern int PCEisenstatSetOmega(PC, double); 5520563c6bSBarry Smith extern int PCSORSetIterations(PC, int); 56d03aef70SBarry Smith 57*b4fd4287SBarry Smith typedef enum {BGS_FORWARD_SWEEP=1,BGS_BACKWARD_SWEEP=2,BGS_SYMMETRIC_SWEEP=3} PCBGSType; 58*b4fd4287SBarry Smith extern int PCBGSSetSymmetric(PC, PCBGSType); 59*b4fd4287SBarry Smith 6049c46530SBarry Smith #define USE_PRECONDITIONER_MATRIX 0 6149c46530SBarry Smith #define USE_TRUE_MATRIX 1 6249c46530SBarry Smith 6349c46530SBarry Smith extern int PCBJacobiSetUseTrueLocal(PC); 6449c46530SBarry Smith extern int PCBJacobiSetTotalBlocks(PC, int, int*,int*); 6549c46530SBarry Smith extern int PCBJacobiSetLocalBlocks(PC, int, int*,int*); 661eb62cbbSBarry Smith 67e3e2a4cfSLois Curfman McInnes extern int PCBSIterSetBlockSolve(PC); 687d5e31e2SLois Curfman McInnes extern int PCBSIterSetFromOptions(PC); 69e3e2a4cfSLois Curfman McInnes extern int PCBSIterSolve(PC,Vec,Vec,int*); 70c06847eeSLois Curfman McInnes 71aabeff55SBarry Smith extern int PCShellSetApply(PC, int (*)(void*,Vec,Vec), void*); 72f0479e8cSBarry Smith extern int PCShellSetApplyRichardson(PC,int (*)(void*,Vec,Vec,Vec,int),void*); 73aabeff55SBarry Smith 74722fc279SLois Curfman McInnes extern int PCGetOperators(PC,Mat*,Mat*,MatStructure*); 75da3a660dSBarry Smith 760213d9f2SLois Curfman McInnes extern int PCLUSetUseInPlace(PC); 770213d9f2SLois Curfman McInnes extern int PCILUSetUseInPlace(PC); 789e25ed09SBarry Smith extern int PCILUSetLevels(PC,int); 7924a74354SBarry Smith extern int PCEisenstatUseDiagonalScaling(PC); 809e25ed09SBarry Smith 81d7e8b826SBarry Smith extern int PCASMCreateSubdomains2D(int,int,int,int,int,int,int *,IS **); 82d7e8b826SBarry Smith extern int PCASMSetSubdomains(PC, int, IS *); 8360907fceSBarry Smith 84d03aef70SBarry Smith #endif 85da3a660dSBarry Smith 86