1*5cd90555SBarry Smith /* $Id: ksp.h,v 1.64 1998/03/23 21:27:32 bsmith Exp bsmith $ */ 2f26ada1bSBarry Smith /* 3f26ada1bSBarry Smith Defines the interface functions for the Krylov subspace accelerators. 4f26ada1bSBarry Smith */ 52eac72dbSBarry Smith #ifndef __KSP_PACKAGE 62eac72dbSBarry Smith #define __KSP_PACKAGE 72eac72dbSBarry Smith #include "petsc.h" 82eac72dbSBarry Smith #include "vec.h" 9aabeff55SBarry Smith #include "mat.h" 10aabeff55SBarry Smith #include "pc.h" 112eac72dbSBarry Smith 129e25ed09SBarry Smith #define KSP_COOKIE PETSC_COOKIE+8 13f0479e8cSBarry Smith 14e2a1c21fSSatish Balay typedef struct _p_KSP* KSP; 152eac72dbSBarry Smith 1682bf6240SBarry Smith #define KSPRICHARDSON "richardson" 1782bf6240SBarry Smith #define KSPCHEBYCHEV "chebychev" 1882bf6240SBarry Smith #define KSPCG "cg" 1982bf6240SBarry Smith #define KSPGMRES "gmres" 2082bf6240SBarry Smith #define KSPTCQMR "tcqmr" 2182bf6240SBarry Smith #define KSPBCGS "bcgs" 2282bf6240SBarry Smith #define KSPCGS "cgs" 2382bf6240SBarry Smith #define KSPTFQMR "tfqmr" 2482bf6240SBarry Smith #define KSPCR "cr" 2582bf6240SBarry Smith #define KSPLSQR "lsqr" 2682bf6240SBarry Smith #define KSPPREONLY "preonly" 2782bf6240SBarry Smith #define KSPQCG "qcg" 2882bf6240SBarry Smith #define KSPTRLS "trls" 2982bf6240SBarry Smith typedef char * KSPType; 302eac72dbSBarry Smith 316b5873e3SBarry Smith extern int KSPCreate(MPI_Comm,KSP *); 324b0e389bSBarry Smith extern int KSPSetType(KSP,KSPType); 338ed539a5SBarry Smith extern int KSPSetUp(KSP); 348ed539a5SBarry Smith extern int KSPSolve(KSP,int *); 358ed539a5SBarry Smith extern int KSPDestroy(KSP); 362eac72dbSBarry Smith 3782bf6240SBarry Smith extern DLList KSPList; 3882bf6240SBarry Smith extern int KSPRegisterAll(char *); 39cf256101SBarry Smith extern int KSPRegisterDestroy(void); 402eac72dbSBarry Smith 4182bf6240SBarry Smith 4282bf6240SBarry Smith extern int KSPGetType(KSP, KSPType *); 43aad2872bSLois Curfman McInnes extern int KSPSetPreconditionerSide(KSP,PCSide); 44aad2872bSLois Curfman McInnes extern int KSPGetPreconditionerSide(KSP,PCSide*); 45ae622e99SLois Curfman McInnes extern int KSPGetTolerances(KSP,double*,double*,double*,int*); 46f9a5357eSLois Curfman McInnes extern int KSPSetTolerances(KSP,double,double,double,int); 47d4fbbf0eSBarry Smith extern int KSPSetComputeResidual(KSP,PetscTruth); 488ed539a5SBarry Smith extern int KSPSetUsePreconditionedResidual(KSP); 49c39b5c85SLois Curfman McInnes extern int KSPSetInitialGuessNonzero(KSP); 5090f02eecSBarry Smith extern int KSPSetComputeEigenvalues(KSP); 51d4fbbf0eSBarry Smith extern int KSPSetComputeSingularValues(KSP); 528ed539a5SBarry Smith extern int KSPSetRhs(KSP,Vec); 538ed539a5SBarry Smith extern int KSPGetRhs(KSP,Vec *); 548ed539a5SBarry Smith extern int KSPSetSolution(KSP,Vec); 558ed539a5SBarry Smith extern int KSPGetSolution(KSP,Vec *); 5625fce39dSBarry Smith extern int KSPGetResidualNorm(KSP,double*); 572eac72dbSBarry Smith 5877c4ece6SBarry Smith extern int KSPSetPC(KSP,PC); 5977c4ece6SBarry Smith extern int KSPGetPC(KSP,PC*); 60aabeff55SBarry Smith 618ed539a5SBarry Smith extern int KSPSetMonitor(KSP,int (*)(KSP,int,double, void*), void *); 62*5cd90555SBarry Smith extern int KSPClearMonitor(KSP); 638ed539a5SBarry Smith extern int KSPGetMonitorContext(KSP,void **); 648ed539a5SBarry Smith extern int KSPSetResidualHistory(KSP, double *,int); 654b0e389bSBarry Smith 668ed539a5SBarry Smith extern int KSPSetConvergenceTest(KSP,int (*)(KSP,int,double, void*), void *); 678ed539a5SBarry Smith extern int KSPGetConvergenceContext(KSP,void **); 682eac72dbSBarry Smith 698ed539a5SBarry Smith extern int KSPBuildSolution(KSP, Vec,Vec *); 708ed539a5SBarry Smith extern int KSPBuildResidual(KSP, Vec, Vec,Vec *); 712eac72dbSBarry Smith 728ed539a5SBarry Smith extern int KSPRichardsonSetScale(KSP , double); 738ed539a5SBarry Smith extern int KSPChebychevSetEigenvalues(KSP , double, double); 74d4fbbf0eSBarry Smith extern int KSPComputeExtremeSingularValues(KSP, double*,double*); 75d4fbbf0eSBarry Smith extern int KSPComputeEigenvalues(KSP,int,double*,double*); 76d4fbbf0eSBarry Smith extern int KSPComputeEigenvaluesExplicitly(KSP,int,double*,double*); 774b0e389bSBarry Smith 788ed539a5SBarry Smith extern int KSPGMRESSetRestart(KSP, int); 79b4fd4287SBarry Smith extern int KSPGMRESSetPreAllocateVectors(KSP); 8077c4ece6SBarry Smith extern int KSPGMRESSetOrthogonalization(KSP,int (*)(KSP,int)); 81bcd2baecSBarry Smith extern int KSPGMRESUnmodifiedGramSchmidtOrthogonalization(KSP,int); 82bcd2baecSBarry Smith extern int KSPGMRESModifiedGramSchmidtOrthogonalization(KSP,int); 83bcd2baecSBarry Smith extern int KSPGMRESIROrthogonalization(KSP,int); 84eed86810SBarry Smith extern int KSPGMRESDGKSOrthogonalization(KSP,int); 8508480c60SBarry Smith 868ed539a5SBarry Smith extern int KSPSetFromOptions(KSP); 87639f9d9dSBarry Smith extern int KSPAddOptionsChecker(int (*)(KSP)); 882eac72dbSBarry Smith 898f4c8dbaSBarry Smith extern int KSPSingularValueMonitor(KSP,int,double, void * ); 908ed539a5SBarry Smith extern int KSPDefaultMonitor(KSP,int,double, void *); 91af6b99e9SBarry Smith extern int KSPTrueMonitor(KSP,int,double, void *); 928c48e012SBarry Smith extern int KSPDefaultSMonitor(KSP,int,double, void *); 9384cb2905SBarry Smith 948ed539a5SBarry Smith extern int KSPDefaultConverged(KSP,int,double, void *); 952eac72dbSBarry Smith 96c01c455dSBarry Smith extern int KSPResidual(KSP,Vec,Vec,Vec,Vec,Vec,Vec); 9777c4ece6SBarry Smith extern int KSPUnwindPreconditioner(KSP,Vec,Vec); 98c01c455dSBarry Smith extern int KSPDefaultBuildSolution(KSP,Vec,Vec*); 99c01c455dSBarry Smith extern int KSPDefaultBuildResidual(KSP,Vec,Vec,Vec *); 100c01c455dSBarry Smith 1018ed539a5SBarry Smith extern int KSPPrintHelp(KSP); 1028ed539a5SBarry Smith 1031eb62cbbSBarry Smith extern int KSPSetOptionsPrefix(KSP,char*); 1040985bb82SSatish Balay extern int KSPAppendOptionsPrefix(KSP,char*); 1050985bb82SSatish Balay extern int KSPGetOptionsPrefix(KSP,char**); 1061eb62cbbSBarry Smith 1078ed539a5SBarry Smith extern int KSPView(KSP,Viewer); 1081eb62cbbSBarry Smith 109d4fbbf0eSBarry Smith extern int KSPComputeExplicitOperator(KSP,Mat *); 110d4fbbf0eSBarry Smith 1116d4a8577SBarry Smith typedef enum {KSP_CG_SYMMETRIC=1, KSP_CG_HERMITIAN=2} KSPCGType; 1126d4a8577SBarry Smith extern int KSPCGSetType(KSP,KSPCGType); 113e559a7feSLois Curfman McInnes 1143369ce9aSBarry Smith extern int PCPreSolve(PC,KSP); 1153369ce9aSBarry Smith extern int PCPostSolve(PC,KSP); 1163369ce9aSBarry Smith 117d7e8b826SBarry Smith extern int KSPLGMonitorCreate(char*,char*,int,int,int,int,DrawLG*); 1181eb62cbbSBarry Smith extern int KSPLGMonitor(KSP,int,double,void*); 119d7e8b826SBarry Smith extern int KSPLGMonitorDestroy(DrawLG); 1203b2fbd54SBarry Smith extern int KSPLGTrueMonitorCreate(MPI_Comm,char*,char*,int,int,int,int,DrawLG*); 121af6b99e9SBarry Smith extern int KSPLGTrueMonitor(KSP,int,double,void*); 122af6b99e9SBarry Smith extern int KSPLGTrueMonitorDestroy(DrawLG); 1231eb62cbbSBarry Smith 1242eac72dbSBarry Smith #endif 1252eac72dbSBarry Smith 1262eac72dbSBarry Smith 127