#ifndef __KSP_PACKAGE
#define __KSP_PACKAGE
#include "petsc.h"
#include "vec.h"

typedef struct _KSP*     KSP;

/*  Possible Krylov Space Methods */
typedef enum { KSPRICHARDSON, KSPCHEBYCHEV, KSPCG, KSPGMRES, 
               KSPTCQMR, KSPBCGS, KSPCGS, KSPTFQMR, KSPCR, KSPLSQR,
               KSPPREONLY } KSPMETHOD;

extern int KSPPrintMethods(char *);

extern int KSPCreate(KSP *);
extern int KSPSetMethod(KSP,KSPMETHOD);
extern int KSPSetUp(KSP);
extern int KSPSolve(KSP,int *);
extern int KSPDestroy(KSP);

extern int KSPRegisterAll();
extern int KSPRegister(KSPMETHOD,char *,int (*)(KSP));

extern int KSPGetMethodName(KSPMETHOD,char **);
extern int KSPGetWorkCounts(KSP,int*,int*,int*,int*,int*);
extern int KSPClearWorkCounts(KSP);

extern int KSPSetIterations(KSP,int);
extern int KSPSetRightPreconditioner(KSP);
extern int KSPGetPreconditionerSide(KSP,int *);
extern int KSPGetMethodFromContext(KSP,KSPMETHOD *);
extern int KSPGetMethodFromOptions(int,char*,KSPMETHOD *);
extern int KSPSetRelativeTolerance(KSP,double);
extern int KSPSetAbsoluteTolerance(KSP,double);
extern int KSPSetDivergenceTolerance(KSP,double);
extern int KSPSetCalculateResidual(KSP);
extern int KSPSetDoNotCalculateResidual(KSP);
extern int KSPSetUsePreconditionedResidual(KSP);
extern int KSPSetInitialGuessZero(KSP);
extern int KSPSetCalculateEigenvalues(KSP);
extern int KSPSetRhs(KSP,Vec);
extern int KSPGetRhs(KSP,Vec *);
extern int KSPSetSolution(KSP,Vec);
extern int KSPGetSolution(KSP,Vec *);
extern int KSPSetAmult(KSP,int (*)(void *,Vec,Vec),void *);

extern int KSPGetAmultContext(KSP,void **);
extern int KSPSetAmultTranspose(KSP,int (*)(void *,Vec,Vec));
extern int KSPSetBinv(KSP,int(*)(void *,Vec,Vec),void *);
extern int KSPGetBinvContext(KSP,void **);
extern int KSPSetBinvTranspose(KSP,int (*)(void *,Vec,Vec));
extern int KSPSetMatop(KSP,int (*)(void *,void *,Vec,Vec));
extern int KSPSetMatopTranspose(KSP,int (*)(void *,void *,Vec,Vec));
extern int KSPSetMonitor(KSP,int (*)(KSP,int,double, void*), void *);
extern int KSPGetMonitorContext(KSP,void **);
extern int KSPSetResidualHistory(KSP, double *,int);
extern int KSPSetConvergenceTest(KSP,int (*)(KSP,int,double, void*), void *);
extern int KSPGetConvergenceContext(KSP,void **);

extern int KSPBuildSolution(KSP, Vec,Vec *);
extern int KSPBuildResidual(KSP, Vec, Vec,Vec *);

extern int KSPRichardsonSetScale(KSP , double);
extern int KSPChebychevSetEigenvalues(KSP , double, double);
extern int KSPCGGetEigenvalues(KSP, int n,double*,double*);
extern int KSPGMRESSetRestart(KSP, int);
extern int KSPSetFromOptions(KSP);

extern int KSPDefaultCGMonitor(KSP,int,double, void * );
extern int KSPDefaultCGConverged(KSP,int,double, void *);
extern int KSPDefaultMonitor(KSP,int,double, void *);
extern int KSPLineGraphMonitor(KSP,int,double, void *);
extern int KSPDefaultConverged(KSP,int,double, void *);

extern int KSPPrintHelp(KSP);

extern int KSPRichardsonSetFast(KSP,int (*)(void*,Vec,Vec,Vec,int),void*);

extern int KSPView(KSP,Viewer);
#endif