Actual source code: slepceps.h

  1: /*
  2:    User interface for the SLEPC eigenproblem solvers. 

  4:    - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
  5:    SLEPc - Scalable Library for Eigenvalue Problem Computations
  6:    Copyright (c) 2002-2012, Universitat Politecnica de Valencia, Spain

  8:    This file is part of SLEPc.
  9:       
 10:    SLEPc is free software: you can redistribute it and/or modify it under  the
 11:    terms of version 3 of the GNU Lesser General Public License as published by
 12:    the Free Software Foundation.

 14:    SLEPc  is  distributed in the hope that it will be useful, but WITHOUT  ANY 
 15:    WARRANTY;  without even the implied warranty of MERCHANTABILITY or  FITNESS 
 16:    FOR  A  PARTICULAR PURPOSE. See the GNU Lesser General Public  License  for 
 17:    more details.

 19:    You  should have received a copy of the GNU Lesser General  Public  License
 20:    along with SLEPc. If not, see <http://www.gnu.org/licenses/>.
 21:    - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
 22: */

 26:  #include slepcsys.h
 27:  #include slepcst.h
 28:  #include slepcip.h
 29:  #include slepcds.h

 31: PETSC_EXTERN PetscErrorCode EPSInitializePackage(const char[]);

 33: /*S
 34:     EPS - Abstract SLEPc object that manages all the eigenvalue 
 35:     problem solvers.

 37:     Level: beginner

 39: .seealso:  EPSCreate(), ST
 40: S*/
 41: typedef struct _p_EPS* EPS;

 43: /*E
 44:     EPSType - String with the name of a SLEPc eigensolver

 46:     Level: beginner

 48: .seealso: EPSSetType(), EPS
 49: E*/
 50: #define EPSType        char*
 51: #define EPSPOWER       "power"
 52: #define EPSSUBSPACE    "subspace"
 53: #define EPSARNOLDI     "arnoldi"
 54: #define EPSLANCZOS     "lanczos"
 55: #define EPSKRYLOVSCHUR "krylovschur"
 56: #define EPSGD          "gd"
 57: #define EPSJD          "jd"
 58: #define EPSRQCG        "rqcg"
 59: #define EPSLAPACK      "lapack"
 60: /* the next ones are interfaces to external libraries */
 61: #define EPSARPACK      "arpack"
 62: #define EPSBLZPACK     "blzpack"
 63: #define EPSTRLAN       "trlan"
 64: #define EPSBLOPEX      "blopex"
 65: #define EPSPRIMME      "primme"

 67: /* Logging support */
 68: PETSC_EXTERN PetscClassId EPS_CLASSID;

 70: /*E
 71:     EPSProblemType - determines the type of eigenvalue problem

 73:     Level: beginner

 75: .seealso: EPSSetProblemType(), EPSGetProblemType()
 76: E*/
 77: typedef enum { EPS_HEP=1,
 78:                EPS_GHEP,
 79:                EPS_NHEP,
 80:                EPS_GNHEP,
 81:                EPS_PGNHEP,
 82:                EPS_GHIEP } EPSProblemType;

 84: /*E
 85:     EPSExtraction - determines the type of extraction technique employed
 86:     by the eigensolver

 88:     Level: beginner

 90: .seealso: EPSSetExtraction(), EPSGetExtraction()
 91: E*/
 92: typedef enum { EPS_RITZ=1,
 93:                EPS_HARMONIC,
 94:                EPS_HARMONIC_RELATIVE,
 95:                EPS_HARMONIC_RIGHT,
 96:                EPS_HARMONIC_LARGEST,
 97:                EPS_REFINED,
 98:                EPS_REFINED_HARMONIC } EPSExtraction;

100: /*E
101:     EPSWhich - determines which part of the spectrum is requested

103:     Level: intermediate

105: .seealso: EPSSetWhichEigenpairs(), EPSGetWhichEigenpairs()
106: E*/
107: typedef enum { EPS_LARGEST_MAGNITUDE=1,
108:                EPS_SMALLEST_MAGNITUDE,
109:                EPS_LARGEST_REAL,
110:                EPS_SMALLEST_REAL,
111:                EPS_LARGEST_IMAGINARY,
112:                EPS_SMALLEST_IMAGINARY,
113:                EPS_TARGET_MAGNITUDE,
114:                EPS_TARGET_REAL,
115:                EPS_TARGET_IMAGINARY,
116:                EPS_ALL,
117:                EPS_WHICH_USER } EPSWhich;

119: /*E
120:     EPSBalance - the type of balancing used for non-Hermitian problems

122:     Level: intermediate

124: .seealso: EPSSetBalance()
125: E*/
126: typedef enum { EPS_BALANCE_NONE=1,
127:                EPS_BALANCE_ONESIDE,
128:                EPS_BALANCE_TWOSIDE,
129:                EPS_BALANCE_USER } EPSBalance;

131: /*E
132:     EPSConv - determines the convergence test

134:     Level: intermediate

136: .seealso: EPSSetConvergenceTest(), EPSSetConvergenceTestFunction()
137: E*/
138: typedef enum { EPS_CONV_ABS=1,
139:                EPS_CONV_EIG,
140:                EPS_CONV_NORM,
141:                EPS_CONV_USER } EPSConv;

143: PETSC_EXTERN PetscErrorCode EPSCreate(MPI_Comm,EPS *);
144: PETSC_EXTERN PetscErrorCode EPSDestroy(EPS*);
145: PETSC_EXTERN PetscErrorCode EPSReset(EPS);
146: PETSC_EXTERN PetscErrorCode EPSSetType(EPS,const EPSType);
147: PETSC_EXTERN PetscErrorCode EPSGetType(EPS,const EPSType*);
148: PETSC_EXTERN PetscErrorCode EPSSetProblemType(EPS,EPSProblemType);
149: PETSC_EXTERN PetscErrorCode EPSGetProblemType(EPS,EPSProblemType*);
150: PETSC_EXTERN PetscErrorCode EPSSetExtraction(EPS,EPSExtraction);
151: PETSC_EXTERN PetscErrorCode EPSGetExtraction(EPS,EPSExtraction*);
152: PETSC_EXTERN PetscErrorCode EPSSetBalance(EPS,EPSBalance,PetscInt,PetscReal);
153: PETSC_EXTERN PetscErrorCode EPSGetBalance(EPS,EPSBalance*,PetscInt*,PetscReal*);
154: PETSC_EXTERN PetscErrorCode EPSSetOperators(EPS,Mat,Mat);
155: PETSC_EXTERN PetscErrorCode EPSGetOperators(EPS,Mat*,Mat*);
156: PETSC_EXTERN PetscErrorCode EPSSetFromOptions(EPS);
157: PETSC_EXTERN PetscErrorCode EPSSetUp(EPS);
158: PETSC_EXTERN PetscErrorCode EPSSolve(EPS);
159: PETSC_EXTERN PetscErrorCode EPSView(EPS,PetscViewer);
160: PETSC_EXTERN PetscErrorCode EPSPrintSolution(EPS,PetscViewer);

162: PETSC_EXTERN PetscErrorCode EPSSetTarget(EPS,PetscScalar);
163: PETSC_EXTERN PetscErrorCode EPSGetTarget(EPS,PetscScalar*);
164: PETSC_EXTERN PetscErrorCode EPSSetInterval(EPS,PetscReal,PetscReal);
165: PETSC_EXTERN PetscErrorCode EPSGetInterval(EPS,PetscReal*,PetscReal*);
166: PETSC_EXTERN PetscErrorCode EPSSetST(EPS,ST);
167: PETSC_EXTERN PetscErrorCode EPSGetST(EPS,ST*);
168: PETSC_EXTERN PetscErrorCode EPSSetIP(EPS,IP);
169: PETSC_EXTERN PetscErrorCode EPSGetIP(EPS,IP*);
170: PETSC_EXTERN PetscErrorCode EPSSetDS(EPS,DS);
171: PETSC_EXTERN PetscErrorCode EPSGetDS(EPS,DS*);
172: PETSC_EXTERN PetscErrorCode EPSSetTolerances(EPS,PetscReal,PetscInt);
173: PETSC_EXTERN PetscErrorCode EPSGetTolerances(EPS,PetscReal*,PetscInt*);
174: PETSC_EXTERN PetscErrorCode EPSSetConvergenceTestFunction(EPS,PetscErrorCode (*)(EPS,PetscScalar,PetscScalar,PetscReal,PetscReal*,void*),void*);
175: PETSC_EXTERN PetscErrorCode EPSSetConvergenceTest(EPS eps,EPSConv conv);
176: PETSC_EXTERN PetscErrorCode EPSGetConvergenceTest(EPS eps,EPSConv *conv);
177: PETSC_EXTERN PetscErrorCode EPSConvergedEigRelative(EPS,PetscScalar,PetscScalar,PetscReal,PetscReal*,void*);
178: PETSC_EXTERN PetscErrorCode EPSConvergedAbsolute(EPS,PetscScalar,PetscScalar,PetscReal,PetscReal*,void*);
179: PETSC_EXTERN PetscErrorCode EPSConvergedNormRelative(EPS,PetscScalar,PetscScalar,PetscReal,PetscReal*,void*);
180: PETSC_EXTERN PetscErrorCode EPSSetDimensions(EPS,PetscInt,PetscInt,PetscInt);
181: PETSC_EXTERN PetscErrorCode EPSGetDimensions(EPS,PetscInt*,PetscInt*,PetscInt*);

183: PETSC_EXTERN PetscErrorCode EPSGetConverged(EPS,PetscInt*);
184: PETSC_EXTERN PetscErrorCode EPSGetEigenpair(EPS,PetscInt,PetscScalar*,PetscScalar*,Vec,Vec);
185: PETSC_EXTERN PetscErrorCode EPSGetEigenvalue(EPS,PetscInt,PetscScalar*,PetscScalar*);
186: PETSC_EXTERN PetscErrorCode EPSGetEigenvector(EPS,PetscInt,Vec,Vec);
187: PETSC_EXTERN PetscErrorCode EPSGetEigenvectorLeft(EPS,PetscInt,Vec,Vec);
188: PETSC_EXTERN PetscErrorCode EPSComputeRelativeError(EPS,PetscInt,PetscReal*);
189: PETSC_EXTERN PetscErrorCode EPSComputeRelativeErrorLeft(EPS,PetscInt,PetscReal*);
190: PETSC_EXTERN PetscErrorCode EPSComputeResidualNorm(EPS,PetscInt,PetscReal*);
191: PETSC_EXTERN PetscErrorCode EPSComputeResidualNormLeft(EPS,PetscInt,PetscReal*);
192: PETSC_EXTERN PetscErrorCode EPSGetInvariantSubspace(EPS,Vec*);
193: PETSC_EXTERN PetscErrorCode EPSGetInvariantSubspaceLeft(EPS,Vec*);
194: PETSC_EXTERN PetscErrorCode EPSGetErrorEstimate(EPS,PetscInt,PetscReal*);
195: PETSC_EXTERN PetscErrorCode EPSGetErrorEstimateLeft(EPS,PetscInt,PetscReal*);

197: PETSC_EXTERN PetscErrorCode EPSMonitorSet(EPS,PetscErrorCode (*)(EPS,PetscInt,PetscInt,PetscScalar*,PetscScalar*,PetscReal*,PetscInt,void*),void*,PetscErrorCode (*)(void**));
198: PETSC_EXTERN PetscErrorCode EPSMonitorCancel(EPS);
199: PETSC_EXTERN PetscErrorCode EPSGetMonitorContext(EPS,void **);
200: PETSC_EXTERN PetscErrorCode EPSGetIterationNumber(EPS,PetscInt*);
201: PETSC_EXTERN PetscErrorCode EPSGetOperationCounters(EPS,PetscInt*,PetscInt*,PetscInt*);

203: PETSC_EXTERN PetscErrorCode EPSSetWhichEigenpairs(EPS,EPSWhich);
204: PETSC_EXTERN PetscErrorCode EPSGetWhichEigenpairs(EPS,EPSWhich*);
205: PETSC_EXTERN PetscErrorCode EPSSetLeftVectorsWanted(EPS,PetscBool);
206: PETSC_EXTERN PetscErrorCode EPSGetLeftVectorsWanted(EPS,PetscBool*);
207: PETSC_EXTERN PetscErrorCode EPSSetMatrixNorms(EPS,PetscReal,PetscReal,PetscBool);
208: PETSC_EXTERN PetscErrorCode EPSGetMatrixNorms(EPS,PetscReal*,PetscReal*,PetscBool*);
209: PETSC_EXTERN PetscErrorCode EPSSetTrueResidual(EPS,PetscBool);
210: PETSC_EXTERN PetscErrorCode EPSGetTrueResidual(EPS,PetscBool*);
211: PETSC_EXTERN PetscErrorCode EPSSetEigenvalueComparison(EPS,PetscErrorCode (*func)(PetscScalar,PetscScalar,PetscScalar,PetscScalar,PetscInt*,void*),void*);
212: PETSC_EXTERN PetscErrorCode EPSSetArbitrarySelection(EPS,PetscErrorCode (*func)(PetscScalar,PetscScalar,Vec,Vec,PetscScalar*,PetscScalar*,void*),void*);
213: PETSC_EXTERN PetscErrorCode EPSIsGeneralized(EPS,PetscBool*);
214: PETSC_EXTERN PetscErrorCode EPSIsHermitian(EPS,PetscBool*);
215: PETSC_EXTERN PetscErrorCode EPSIsPositive(EPS,PetscBool*);

217: PETSC_EXTERN PetscErrorCode EPSMonitorFirst(EPS,PetscInt,PetscInt,PetscScalar*,PetscScalar*,PetscReal*,PetscInt,void*);
218: PETSC_EXTERN PetscErrorCode EPSMonitorAll(EPS,PetscInt,PetscInt,PetscScalar*,PetscScalar*,PetscReal*,PetscInt,void*);
219: PETSC_EXTERN PetscErrorCode EPSMonitorConverged(EPS,PetscInt,PetscInt,PetscScalar*,PetscScalar*,PetscReal*,PetscInt,void*);
220: PETSC_EXTERN PetscErrorCode EPSMonitorLG(EPS,PetscInt,PetscInt,PetscScalar*,PetscScalar*,PetscReal*,PetscInt,void*);
221: PETSC_EXTERN PetscErrorCode EPSMonitorLGAll(EPS,PetscInt,PetscInt,PetscScalar*,PetscScalar*,PetscReal*,PetscInt,void*);

223: PETSC_EXTERN PetscErrorCode EPSSetTrackAll(EPS,PetscBool);
224: PETSC_EXTERN PetscErrorCode EPSGetTrackAll(EPS,PetscBool*);

226: PETSC_EXTERN PetscErrorCode EPSSetDeflationSpace(EPS,PetscInt,Vec*);
227: PETSC_EXTERN PetscErrorCode EPSRemoveDeflationSpace(EPS);
228: PETSC_EXTERN PetscErrorCode EPSSetInitialSpace(EPS,PetscInt,Vec*);
229: PETSC_EXTERN PetscErrorCode EPSSetInitialSpaceLeft(EPS,PetscInt,Vec*);

231: PETSC_EXTERN PetscErrorCode EPSSetOptionsPrefix(EPS,const char*);
232: PETSC_EXTERN PetscErrorCode EPSAppendOptionsPrefix(EPS,const char*);
233: PETSC_EXTERN PetscErrorCode EPSGetOptionsPrefix(EPS,const char*[]);

235: /*E
236:     EPSConvergedReason - reason an eigensolver was said to 
237:          have converged or diverged

239:     Level: beginner

241: .seealso: EPSSolve(), EPSGetConvergedReason(), EPSSetTolerances()
242: E*/
243: typedef enum {/* converged */
244:               EPS_CONVERGED_TOL                =  2,
245:               /* diverged */
246:               EPS_DIVERGED_ITS                 = -3,
247:               EPS_DIVERGED_BREAKDOWN           = -4,
248:               EPS_CONVERGED_ITERATING          =  0} EPSConvergedReason;

250: PETSC_EXTERN PetscErrorCode EPSGetConvergedReason(EPS,EPSConvergedReason *);

252: PETSC_EXTERN PetscErrorCode EPSSortEigenvalues(EPS,PetscInt,PetscScalar*,PetscScalar*,PetscInt*);
253: PETSC_EXTERN PetscErrorCode EPSCompareEigenvalues(EPS,PetscScalar,PetscScalar,PetscScalar,PetscScalar,PetscInt*);

255: PETSC_EXTERN PetscErrorCode EPSGetStartVector(EPS,PetscInt,Vec,PetscBool*);
256: PETSC_EXTERN PetscErrorCode EPSGetStartVectorLeft(EPS,PetscInt,Vec,PetscBool*);

258: PETSC_EXTERN PetscFList EPSList;
259: PETSC_EXTERN PetscBool  EPSRegisterAllCalled;
260: PETSC_EXTERN PetscErrorCode EPSRegisterAll(const char[]);
261: PETSC_EXTERN PetscErrorCode EPSRegisterDestroy(void);
262: PETSC_EXTERN PetscErrorCode EPSRegister(const char[],const char[],const char[],PetscErrorCode(*)(EPS));

264: /*MC
265:    EPSRegisterDynamic - Adds a method to the eigenproblem solver package.

267:    Synopsis:
268:    PetscErrorCode EPSRegisterDynamic(const char *name_solver,const char *path,const char *name_create,PetscErrorCode (*routine_create)(EPS))

270:    Not Collective

272:    Input Parameters:
273: +  name_solver - name of a new user-defined solver
274: .  path - path (either absolute or relative) the library containing this solver
275: .  name_create - name of routine to create the solver context
276: -  routine_create - routine to create the solver context

278:    Notes:
279:    EPSRegisterDynamic() may be called multiple times to add several user-defined solvers.

281:    If dynamic libraries are used, then the fourth input argument (routine_create)
282:    is ignored.

284:    Sample usage:
285: .vb
286:    EPSRegisterDynamic("my_solver",/home/username/my_lib/lib/libO/solaris/mylib.a,
287:                "MySolverCreate",MySolverCreate);
288: .ve

290:    Then, your solver can be chosen with the procedural interface via
291: $     EPSSetType(eps,"my_solver")
292:    or at runtime via the option
293: $     -eps_type my_solver

295:    Level: advanced

297: .seealso: EPSRegisterDestroy(), EPSRegisterAll()
298: M*/
299: #if defined(PETSC_USE_DYNAMIC_LIBRARIES)
300: #define EPSRegisterDynamic(a,b,c,d) EPSRegister(a,b,c,0)
301: #else
302: #define EPSRegisterDynamic(a,b,c,d) EPSRegister(a,b,c,d)
303: #endif

305: /* --------- options specific to particular eigensolvers -------- */

307: /*E
308:     EPSPowerShiftType - determines the type of shift used in the Power iteration

310:     Level: advanced

312: .seealso: EPSPowerSetShiftType(), EPSPowerGetShiftType()
313: E*/
314: typedef enum { EPS_POWER_SHIFT_CONSTANT,
315:                EPS_POWER_SHIFT_RAYLEIGH,
316:                EPS_POWER_SHIFT_WILKINSON } EPSPowerShiftType;
317: PETSC_EXTERN const char *EPSPowerShiftTypes[];

319: PETSC_EXTERN PetscErrorCode EPSPowerSetShiftType(EPS,EPSPowerShiftType);
320: PETSC_EXTERN PetscErrorCode EPSPowerGetShiftType(EPS,EPSPowerShiftType*);

322: PETSC_EXTERN PetscErrorCode EPSArnoldiSetDelayed(EPS,PetscBool);
323: PETSC_EXTERN PetscErrorCode EPSArnoldiGetDelayed(EPS,PetscBool*);

325: PETSC_EXTERN PetscErrorCode EPSKrylovSchurSetRestart(EPS,PetscReal);
326: PETSC_EXTERN PetscErrorCode EPSKrylovSchurGetRestart(EPS,PetscReal*);

328: /*E
329:     EPSLanczosReorthogType - determines the type of reorthogonalization
330:     used in the Lanczos method

332:     Level: advanced

334: .seealso: EPSLanczosSetReorthog(), EPSLanczosGetReorthog()
335: E*/
336: typedef enum { EPS_LANCZOS_REORTHOG_LOCAL,
337:                EPS_LANCZOS_REORTHOG_FULL,
338:                EPS_LANCZOS_REORTHOG_SELECTIVE,
339:                EPS_LANCZOS_REORTHOG_PERIODIC,
340:                EPS_LANCZOS_REORTHOG_PARTIAL,
341:                EPS_LANCZOS_REORTHOG_DELAYED } EPSLanczosReorthogType;
342: PETSC_EXTERN const char *EPSLanczosReorthogTypes[];

344: PETSC_EXTERN PetscErrorCode EPSLanczosSetReorthog(EPS,EPSLanczosReorthogType);
345: PETSC_EXTERN PetscErrorCode EPSLanczosGetReorthog(EPS,EPSLanczosReorthogType*);

347: PETSC_EXTERN PetscErrorCode EPSBlzpackSetBlockSize(EPS,PetscInt);
348: PETSC_EXTERN PetscErrorCode EPSBlzpackSetNSteps(EPS,PetscInt);

350: /*E
351:     EPSPRIMMEMethod - determines the method selected in the PRIMME library

353:     Level: advanced

355: .seealso: EPSPRIMMESetMethod(), EPSPRIMMEGetMethod()
356: E*/
357: typedef enum { EPS_PRIMME_DYNAMIC,
358:                EPS_PRIMME_DEFAULT_MIN_TIME,
359:                EPS_PRIMME_DEFAULT_MIN_MATVECS,
360:                EPS_PRIMME_ARNOLDI,
361:                EPS_PRIMME_GD,
362:                EPS_PRIMME_GD_PLUSK,
363:                EPS_PRIMME_GD_OLSEN_PLUSK,
364:                EPS_PRIMME_JD_OLSEN_PLUSK,
365:                EPS_PRIMME_RQI,
366:                EPS_PRIMME_JDQR,
367:                EPS_PRIMME_JDQMR,
368:                EPS_PRIMME_JDQMR_ETOL,
369:                EPS_PRIMME_SUBSPACE_ITERATION,
370:                EPS_PRIMME_LOBPCG_ORTHOBASIS,
371:                EPS_PRIMME_LOBPCG_ORTHOBASISW } EPSPRIMMEMethod;
372: PETSC_EXTERN const char *EPSPRIMMEMethods[];

374: PETSC_EXTERN PetscErrorCode EPSPRIMMESetBlockSize(EPS eps,PetscInt bs);
375: PETSC_EXTERN PetscErrorCode EPSPRIMMESetMethod(EPS eps, EPSPRIMMEMethod method);
376: PETSC_EXTERN PetscErrorCode EPSPRIMMEGetBlockSize(EPS eps,PetscInt *bs);
377: PETSC_EXTERN PetscErrorCode EPSPRIMMEGetMethod(EPS eps, EPSPRIMMEMethod *method);

379: /*E
380:     EPSOrthType - determines the orthogonalization used in the search subspace

382:     Level: advanced

384: .seealso: EPSGDSetBOrth(), EPSJDSetBOrth()
385: E*/
386: typedef enum { EPS_ORTH_I=1,
387:                EPS_ORTH_B,
388:                EPS_ORTH_BOPT } EPSOrthType;

390: PETSC_EXTERN PetscErrorCode EPSGDSetKrylovStart(EPS eps,PetscBool krylovstart);
391: PETSC_EXTERN PetscErrorCode EPSGDGetKrylovStart(EPS eps,PetscBool *krylovstart);
392: PETSC_EXTERN PetscErrorCode EPSGDSetBlockSize(EPS eps,PetscInt blocksize);
393: PETSC_EXTERN PetscErrorCode EPSGDGetBlockSize(EPS eps,PetscInt *blocksize);
394: PETSC_EXTERN PetscErrorCode EPSGDSetRestart(EPS eps,PetscInt minv,PetscInt plusk);
395: PETSC_EXTERN PetscErrorCode EPSGDGetRestart(EPS eps,PetscInt *minv,PetscInt *plusk);
396: PETSC_EXTERN PetscErrorCode EPSGDSetInitialSize(EPS eps,PetscInt initialsize);
397: PETSC_EXTERN PetscErrorCode EPSGDGetInitialSize(EPS eps,PetscInt *initialsize);
398: PETSC_EXTERN PetscErrorCode EPSGDSetBOrth(EPS eps,EPSOrthType borth);
399: PETSC_EXTERN PetscErrorCode EPSGDGetBOrth(EPS eps,EPSOrthType *borth);
400: PETSC_EXTERN PetscErrorCode EPSGDGetWindowSizes(EPS eps,PetscInt *pwindow,PetscInt *qwindow);
401: PETSC_EXTERN PetscErrorCode EPSGDSetWindowSizes(EPS eps,PetscInt pwindow,PetscInt qwindow);
402: PETSC_EXTERN PetscErrorCode EPSGDSetDoubleExpansion(EPS eps,PetscBool use_gd2);
403: PETSC_EXTERN PetscErrorCode EPSGDGetDoubleExpansion(EPS eps,PetscBool *flg);

405: PETSC_EXTERN PetscErrorCode EPSJDSetKrylovStart(EPS eps,PetscBool krylovstart);
406: PETSC_EXTERN PetscErrorCode EPSJDGetKrylovStart(EPS eps,PetscBool *krylovstart);
407: PETSC_EXTERN PetscErrorCode EPSJDSetBlockSize(EPS eps,PetscInt blocksize);
408: PETSC_EXTERN PetscErrorCode EPSJDGetBlockSize(EPS eps,PetscInt *blocksize);
409: PETSC_EXTERN PetscErrorCode EPSJDSetRestart(EPS eps,PetscInt minv,PetscInt plusk);
410: PETSC_EXTERN PetscErrorCode EPSJDGetRestart(EPS eps,PetscInt *minv,PetscInt *plusk);
411: PETSC_EXTERN PetscErrorCode EPSJDSetInitialSize(EPS eps,PetscInt initialsize);
412: PETSC_EXTERN PetscErrorCode EPSJDGetInitialSize(EPS eps,PetscInt *initialsize);
413: PETSC_EXTERN PetscErrorCode EPSJDSetFix(EPS eps,PetscReal fix);
414: PETSC_EXTERN PetscErrorCode EPSJDGetFix(EPS eps,PetscReal *fix);
415: PETSC_EXTERN PetscErrorCode EPSJDSetConstantCorrectionTolerance(EPS eps,PetscBool dynamic);
416: PETSC_EXTERN PetscErrorCode EPSJDGetConstantCorrectionTolerance(EPS eps,PetscBool *dynamic);
417: PETSC_EXTERN PetscErrorCode EPSJDSetBOrth(EPS eps,EPSOrthType borth);
418: PETSC_EXTERN PetscErrorCode EPSJDGetBOrth(EPS eps,EPSOrthType *borth);
419: PETSC_EXTERN PetscErrorCode EPSJDGetWindowSizes(EPS eps,PetscInt *pwindow,PetscInt *qwindow);
420: PETSC_EXTERN PetscErrorCode EPSJDSetWindowSizes(EPS eps,PetscInt pwindow,PetscInt qwindow);

422: PETSC_EXTERN PetscErrorCode EPSRQCGSetReset(EPS,PetscInt);
423: PETSC_EXTERN PetscErrorCode EPSRQCGGetReset(EPS,PetscInt*);

425: #endif