Actual source code: slepceps.h

slepc-3.15.1 2021-05-28
Report Typos and Errors
  1: /*
  2:    - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
  3:    SLEPc - Scalable Library for Eigenvalue Problem Computations
  4:    Copyright (c) 2002-2021, Universitat Politecnica de Valencia, Spain

  6:    This file is part of SLEPc.
  7:    SLEPc is distributed under a 2-clause BSD license (see LICENSE).
  8:    - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
  9: */
 10: /*
 11:    User interface for the SLEPc linear eigenvalue solvers
 12: */

 14: #if !defined(SLEPCEPS_H)
 15: #define SLEPCEPS_H
 16: #include <slepcst.h>
 17: #include <slepcbv.h>
 18: #include <slepcds.h>
 19: #include <slepcrg.h>
 20: #include <slepclme.h>
 21: #include <petscsnes.h>

 23: SLEPC_EXTERN PetscErrorCode EPSInitializePackage(void);

 25: /*S
 26:     EPS - Abstract SLEPc object that manages all the eigenvalue
 27:     problem solvers.

 29:     Level: beginner

 31: .seealso:  EPSCreate(), ST
 32: S*/
 33: typedef struct _p_EPS* EPS;

 35: /*J
 36:     EPSType - String with the name of a SLEPc eigensolver

 38:     Level: beginner

 40: .seealso: EPSSetType(), EPS
 41: J*/
 42: typedef const char* EPSType;
 43: #define EPSPOWER       "power"
 44: #define EPSSUBSPACE    "subspace"
 45: #define EPSARNOLDI     "arnoldi"
 46: #define EPSLANCZOS     "lanczos"
 47: #define EPSKRYLOVSCHUR "krylovschur"
 48: #define EPSGD          "gd"
 49: #define EPSJD          "jd"
 50: #define EPSRQCG        "rqcg"
 51: #define EPSLOBPCG      "lobpcg"
 52: #define EPSCISS        "ciss"
 53: #define EPSLYAPII      "lyapii"
 54: #define EPSLAPACK      "lapack"
 55: #define EPSARPACK      "arpack"
 56: #define EPSBLZPACK     "blzpack"
 57: #define EPSTRLAN       "trlan"
 58: #define EPSBLOPEX      "blopex"
 59: #define EPSPRIMME      "primme"
 60: #define EPSFEAST       "feast"
 61: #define EPSSCALAPACK   "scalapack"
 62: #define EPSELPA        "elpa"
 63: #define EPSELEMENTAL   "elemental"
 64: #define EPSEVSL        "evsl"

 66: /* Logging support */
 67: SLEPC_EXTERN PetscClassId EPS_CLASSID;

 69: /*E
 70:     EPSProblemType - Determines the type of eigenvalue problem

 72:     Level: beginner

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

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

 87:     Level: advanced

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

 99: /*E
100:     EPSWhich - Determines which part of the spectrum is requested

102:     Level: intermediate

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

118: /*E
119:     EPSBalance - The type of balancing used for non-Hermitian problems

121:     Level: intermediate

123: .seealso: EPSSetBalance()
124: E*/
125: typedef enum { EPS_BALANCE_NONE,
126:                EPS_BALANCE_ONESIDE,
127:                EPS_BALANCE_TWOSIDE,
128:                EPS_BALANCE_USER } EPSBalance;
129: SLEPC_EXTERN const char *EPSBalanceTypes[];

131: /*E
132:     EPSErrorType - The error type used to assess accuracy of computed solutions

134:     Level: intermediate

136: .seealso: EPSComputeError()
137: E*/
138: typedef enum { EPS_ERROR_ABSOLUTE,
139:                EPS_ERROR_RELATIVE,
140:                EPS_ERROR_BACKWARD } EPSErrorType;
141: SLEPC_EXTERN const char *EPSErrorTypes[];

143: /*E
144:     EPSConv - Determines the convergence test

146:     Level: intermediate

148: .seealso: EPSSetConvergenceTest(), EPSSetConvergenceTestFunction()
149: E*/
150: typedef enum { EPS_CONV_ABS,
151:                EPS_CONV_REL,
152:                EPS_CONV_NORM,
153:                EPS_CONV_USER } EPSConv;

155: /*E
156:     EPSStop - Determines the stopping test

158:     Level: advanced

160: .seealso: EPSSetStoppingTest(), EPSSetStoppingTestFunction()
161: E*/
162: typedef enum { EPS_STOP_BASIC,
163:                EPS_STOP_USER } EPSStop;

165: /*E
166:     EPSConvergedReason - Reason an eigensolver was said to
167:          have converged or diverged

169:     Level: intermediate

171: .seealso: EPSSolve(), EPSGetConvergedReason(), EPSSetTolerances()
172: E*/
173: typedef enum {/* converged */
174:               EPS_CONVERGED_TOL                =  1,
175:               EPS_CONVERGED_USER               =  2,
176:               /* diverged */
177:               EPS_DIVERGED_ITS                 = -1,
178:               EPS_DIVERGED_BREAKDOWN           = -2,
179:               EPS_DIVERGED_SYMMETRY_LOST       = -3,
180:               EPS_CONVERGED_ITERATING          =  0} EPSConvergedReason;
181: SLEPC_EXTERN const char *const*EPSConvergedReasons;

183: SLEPC_EXTERN PetscErrorCode EPSCreate(MPI_Comm,EPS*);
184: SLEPC_EXTERN PetscErrorCode EPSDestroy(EPS*);
185: SLEPC_EXTERN PetscErrorCode EPSReset(EPS);
186: SLEPC_EXTERN PetscErrorCode EPSSetType(EPS,EPSType);
187: SLEPC_EXTERN PetscErrorCode EPSGetType(EPS,EPSType*);
188: SLEPC_EXTERN PetscErrorCode EPSSetProblemType(EPS,EPSProblemType);
189: SLEPC_EXTERN PetscErrorCode EPSGetProblemType(EPS,EPSProblemType*);
190: SLEPC_EXTERN PetscErrorCode EPSSetExtraction(EPS,EPSExtraction);
191: SLEPC_EXTERN PetscErrorCode EPSGetExtraction(EPS,EPSExtraction*);
192: SLEPC_EXTERN PetscErrorCode EPSSetBalance(EPS,EPSBalance,PetscInt,PetscReal);
193: SLEPC_EXTERN PetscErrorCode EPSGetBalance(EPS,EPSBalance*,PetscInt*,PetscReal*);
194: SLEPC_EXTERN PetscErrorCode EPSSetOperators(EPS,Mat,Mat);
195: SLEPC_EXTERN PetscErrorCode EPSGetOperators(EPS,Mat*,Mat*);
196: SLEPC_EXTERN PetscErrorCode EPSSetFromOptions(EPS);
197: SLEPC_EXTERN PetscErrorCode EPSSetUp(EPS);
198: SLEPC_EXTERN PetscErrorCode EPSSolve(EPS);
199: SLEPC_EXTERN PetscErrorCode EPSView(EPS,PetscViewer);
200: SLEPC_EXTERN PetscErrorCode EPSViewFromOptions(EPS,PetscObject,const char[]);
201: SLEPC_EXTERN PetscErrorCode EPSErrorView(EPS,EPSErrorType,PetscViewer);
202: PETSC_DEPRECATED_FUNCTION("Use EPSErrorView()") PETSC_STATIC_INLINE PetscErrorCode EPSPrintSolution(EPS eps,PetscViewer v) {return EPSErrorView(eps,EPS_ERROR_RELATIVE,v);}
203: SLEPC_EXTERN PetscErrorCode EPSErrorViewFromOptions(EPS);
204: SLEPC_EXTERN PetscErrorCode EPSConvergedReasonView(EPS,PetscViewer);
205: SLEPC_EXTERN PetscErrorCode EPSConvergedReasonViewFromOptions(EPS);
206: PETSC_DEPRECATED_FUNCTION("Use EPSConvergedReasonView() (since version 3.14)") PETSC_STATIC_INLINE PetscErrorCode EPSReasonView(EPS eps,PetscViewer v) {return EPSConvergedReasonView(eps,v);}
207: PETSC_DEPRECATED_FUNCTION("Use EPSConvergedReasonViewFromOptions() (since version 3.14)") PETSC_STATIC_INLINE PetscErrorCode EPSReasonViewFromOptions(EPS eps) {return EPSConvergedReasonViewFromOptions(eps);}
208: SLEPC_EXTERN PetscErrorCode EPSValuesView(EPS,PetscViewer);
209: SLEPC_EXTERN PetscErrorCode EPSValuesViewFromOptions(EPS);
210: SLEPC_EXTERN PetscErrorCode EPSVectorsView(EPS,PetscViewer);
211: SLEPC_EXTERN PetscErrorCode EPSVectorsViewFromOptions(EPS);

213: SLEPC_EXTERN PetscErrorCode EPSSetTarget(EPS,PetscScalar);
214: SLEPC_EXTERN PetscErrorCode EPSGetTarget(EPS,PetscScalar*);
215: SLEPC_EXTERN PetscErrorCode EPSSetInterval(EPS,PetscReal,PetscReal);
216: SLEPC_EXTERN PetscErrorCode EPSGetInterval(EPS,PetscReal*,PetscReal*);
217: SLEPC_EXTERN PetscErrorCode EPSSetST(EPS,ST);
218: SLEPC_EXTERN PetscErrorCode EPSGetST(EPS,ST*);
219: SLEPC_EXTERN PetscErrorCode EPSSetBV(EPS,BV);
220: SLEPC_EXTERN PetscErrorCode EPSGetBV(EPS,BV*);
221: SLEPC_EXTERN PetscErrorCode EPSSetRG(EPS,RG);
222: SLEPC_EXTERN PetscErrorCode EPSGetRG(EPS,RG*);
223: SLEPC_EXTERN PetscErrorCode EPSSetDS(EPS,DS);
224: SLEPC_EXTERN PetscErrorCode EPSGetDS(EPS,DS*);
225: SLEPC_EXTERN PetscErrorCode EPSSetTolerances(EPS,PetscReal,PetscInt);
226: SLEPC_EXTERN PetscErrorCode EPSGetTolerances(EPS,PetscReal*,PetscInt*);
227: SLEPC_EXTERN PetscErrorCode EPSSetConvergenceTestFunction(EPS,PetscErrorCode (*)(EPS,PetscScalar,PetscScalar,PetscReal,PetscReal*,void*),void*,PetscErrorCode (*)(void*));
228: SLEPC_EXTERN PetscErrorCode EPSSetConvergenceTest(EPS,EPSConv);
229: SLEPC_EXTERN PetscErrorCode EPSGetConvergenceTest(EPS,EPSConv*);
230: SLEPC_EXTERN PetscErrorCode EPSConvergedAbsolute(EPS,PetscScalar,PetscScalar,PetscReal,PetscReal*,void*);
231: SLEPC_EXTERN PetscErrorCode EPSConvergedRelative(EPS,PetscScalar,PetscScalar,PetscReal,PetscReal*,void*);
232: SLEPC_EXTERN PetscErrorCode EPSConvergedNorm(EPS,PetscScalar,PetscScalar,PetscReal,PetscReal*,void*);
233: SLEPC_EXTERN PetscErrorCode EPSSetStoppingTestFunction(EPS,PetscErrorCode (*)(EPS,PetscInt,PetscInt,PetscInt,PetscInt,EPSConvergedReason*,void*),void*,PetscErrorCode (*)(void*));
234: SLEPC_EXTERN PetscErrorCode EPSSetStoppingTest(EPS,EPSStop);
235: SLEPC_EXTERN PetscErrorCode EPSGetStoppingTest(EPS,EPSStop*);
236: SLEPC_EXTERN PetscErrorCode EPSStoppingBasic(EPS,PetscInt,PetscInt,PetscInt,PetscInt,EPSConvergedReason*,void*);
237: SLEPC_EXTERN PetscErrorCode EPSGetConvergedReason(EPS,EPSConvergedReason*);

239: SLEPC_EXTERN PetscErrorCode EPSSetDimensions(EPS,PetscInt,PetscInt,PetscInt);
240: SLEPC_EXTERN PetscErrorCode EPSGetDimensions(EPS,PetscInt*,PetscInt*,PetscInt*);

242: SLEPC_EXTERN PetscErrorCode EPSGetConverged(EPS,PetscInt*);
243: SLEPC_EXTERN PetscErrorCode EPSGetEigenpair(EPS,PetscInt,PetscScalar*,PetscScalar*,Vec,Vec);
244: SLEPC_EXTERN PetscErrorCode EPSGetEigenvalue(EPS,PetscInt,PetscScalar*,PetscScalar*);
245: SLEPC_EXTERN PetscErrorCode EPSGetEigenvector(EPS,PetscInt,Vec,Vec);
246: SLEPC_EXTERN PetscErrorCode EPSGetLeftEigenvector(EPS,PetscInt,Vec,Vec);

248: SLEPC_EXTERN PetscErrorCode EPSComputeError(EPS,PetscInt,EPSErrorType,PetscReal*);
249: PETSC_DEPRECATED_FUNCTION("Use EPSComputeError()") PETSC_STATIC_INLINE PetscErrorCode EPSComputeRelativeError(EPS eps,PetscInt i,PetscReal *r) {return EPSComputeError(eps,i,EPS_ERROR_RELATIVE,r);}
250: PETSC_DEPRECATED_FUNCTION("Use EPSComputeError() with EPS_ERROR_ABSOLUTE") PETSC_STATIC_INLINE PetscErrorCode EPSComputeResidualNorm(EPS eps,PetscInt i,PetscReal *r) {return EPSComputeError(eps,i,EPS_ERROR_ABSOLUTE,r);}
251: SLEPC_EXTERN PetscErrorCode EPSGetInvariantSubspace(EPS,Vec[]);
252: SLEPC_EXTERN PetscErrorCode EPSGetErrorEstimate(EPS,PetscInt,PetscReal*);
253: SLEPC_EXTERN PetscErrorCode EPSGetIterationNumber(EPS,PetscInt*);

255: SLEPC_EXTERN PetscErrorCode EPSSetWhichEigenpairs(EPS,EPSWhich);
256: SLEPC_EXTERN PetscErrorCode EPSGetWhichEigenpairs(EPS,EPSWhich*);
257: SLEPC_EXTERN PetscErrorCode EPSSetTwoSided(EPS,PetscBool);
258: SLEPC_EXTERN PetscErrorCode EPSGetTwoSided(EPS,PetscBool*);
259: SLEPC_EXTERN PetscErrorCode EPSSetTrueResidual(EPS,PetscBool);
260: SLEPC_EXTERN PetscErrorCode EPSGetTrueResidual(EPS,PetscBool*);
261: SLEPC_EXTERN PetscErrorCode EPSSetPurify(EPS,PetscBool);
262: SLEPC_EXTERN PetscErrorCode EPSGetPurify(EPS,PetscBool*);
263: SLEPC_EXTERN PetscErrorCode EPSSetEigenvalueComparison(EPS,PetscErrorCode (*func)(PetscScalar,PetscScalar,PetscScalar,PetscScalar,PetscInt*,void*),void*);
264: SLEPC_EXTERN PetscErrorCode EPSSetArbitrarySelection(EPS,PetscErrorCode (*func)(PetscScalar,PetscScalar,Vec,Vec,PetscScalar*,PetscScalar*,void*),void*);
265: SLEPC_EXTERN PetscErrorCode EPSIsGeneralized(EPS,PetscBool*);
266: SLEPC_EXTERN PetscErrorCode EPSIsHermitian(EPS,PetscBool*);
267: SLEPC_EXTERN PetscErrorCode EPSIsPositive(EPS,PetscBool*);

269: SLEPC_EXTERN PetscErrorCode EPSSetTrackAll(EPS,PetscBool);
270: SLEPC_EXTERN PetscErrorCode EPSGetTrackAll(EPS,PetscBool*);

272: SLEPC_EXTERN PetscErrorCode EPSSetDeflationSpace(EPS,PetscInt,Vec[]);
273: SLEPC_EXTERN PetscErrorCode EPSSetInitialSpace(EPS,PetscInt,Vec[]);
274: SLEPC_EXTERN PetscErrorCode EPSSetLeftInitialSpace(EPS,PetscInt,Vec[]);

276: SLEPC_EXTERN PetscErrorCode EPSMonitor(EPS,PetscInt,PetscInt,PetscScalar*,PetscScalar*,PetscReal*,PetscInt);
277: SLEPC_EXTERN PetscErrorCode EPSMonitorSet(EPS,PetscErrorCode (*)(EPS,PetscInt,PetscInt,PetscScalar*,PetscScalar*,PetscReal*,PetscInt,void*),void*,PetscErrorCode (*)(void**));
278: SLEPC_EXTERN PetscErrorCode EPSMonitorCancel(EPS);
279: SLEPC_EXTERN PetscErrorCode EPSGetMonitorContext(EPS,void**);

281: SLEPC_EXTERN PetscErrorCode EPSMonitorSetFromOptions(EPS,const char[],const char[],void*,PetscBool);
282: SLEPC_EXTERN PetscErrorCode EPSMonitorLGCreate(MPI_Comm,const char[],const char[],const char[],PetscInt,const char*[],int,int,int,int,PetscDrawLG*);
283: SLEPC_EXTERN PetscErrorCode EPSMonitorFirst(EPS,PetscInt,PetscInt,PetscScalar*,PetscScalar*,PetscReal*,PetscInt,PetscViewerAndFormat*);
284: SLEPC_EXTERN PetscErrorCode EPSMonitorFirstDrawLG(EPS,PetscInt,PetscInt,PetscScalar*,PetscScalar*,PetscReal*,PetscInt,PetscViewerAndFormat*);
285: SLEPC_EXTERN PetscErrorCode EPSMonitorFirstDrawLGCreate(PetscViewer,PetscViewerFormat,void *,PetscViewerAndFormat**);
286: SLEPC_EXTERN PetscErrorCode EPSMonitorAll(EPS,PetscInt,PetscInt,PetscScalar*,PetscScalar*,PetscReal*,PetscInt,PetscViewerAndFormat*);
287: SLEPC_EXTERN PetscErrorCode EPSMonitorAllDrawLG(EPS,PetscInt,PetscInt,PetscScalar*,PetscScalar*,PetscReal*,PetscInt,PetscViewerAndFormat*);
288: SLEPC_EXTERN PetscErrorCode EPSMonitorAllDrawLGCreate(PetscViewer,PetscViewerFormat,void *,PetscViewerAndFormat**);
289: SLEPC_EXTERN PetscErrorCode EPSMonitorConverged(EPS,PetscInt,PetscInt,PetscScalar*,PetscScalar*,PetscReal*,PetscInt,PetscViewerAndFormat*);
290: SLEPC_EXTERN PetscErrorCode EPSMonitorConvergedCreate(PetscViewer,PetscViewerFormat,void *,PetscViewerAndFormat**);
291: SLEPC_EXTERN PetscErrorCode EPSMonitorConvergedDrawLG(EPS,PetscInt,PetscInt,PetscScalar*,PetscScalar*,PetscReal*,PetscInt,PetscViewerAndFormat*);
292: SLEPC_EXTERN PetscErrorCode EPSMonitorConvergedDrawLGCreate(PetscViewer,PetscViewerFormat,void *,PetscViewerAndFormat**);
293: SLEPC_EXTERN PetscErrorCode EPSMonitorConvergedDestroy(PetscViewerAndFormat**);

295: SLEPC_EXTERN PetscErrorCode EPSSetOptionsPrefix(EPS,const char*);
296: SLEPC_EXTERN PetscErrorCode EPSAppendOptionsPrefix(EPS,const char*);
297: SLEPC_EXTERN PetscErrorCode EPSGetOptionsPrefix(EPS,const char*[]);

299: SLEPC_EXTERN PetscFunctionList EPSList;
300: SLEPC_EXTERN PetscFunctionList EPSMonitorList;
301: SLEPC_EXTERN PetscFunctionList EPSMonitorCreateList;
302: SLEPC_EXTERN PetscFunctionList EPSMonitorDestroyList;
303: SLEPC_EXTERN PetscErrorCode EPSRegister(const char[],PetscErrorCode(*)(EPS));
304: SLEPC_EXTERN PetscErrorCode EPSMonitorRegister(const char[],PetscViewerType,PetscViewerFormat,PetscErrorCode(*)(EPS,PetscInt,PetscInt,PetscScalar*,PetscScalar*,PetscReal*,PetscInt,PetscViewerAndFormat*),PetscErrorCode(*)(PetscViewer,PetscViewerFormat,void*,PetscViewerAndFormat**),PetscErrorCode(*)(PetscViewerAndFormat**));

306: SLEPC_EXTERN PetscErrorCode EPSSetWorkVecs(EPS,PetscInt);
307: SLEPC_EXTERN PetscErrorCode EPSAllocateSolution(EPS,PetscInt);

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

311: /*E
312:     EPSPowerShiftType - determines the type of shift used in the Power iteration

314:     Level: advanced

316: .seealso: EPSPowerSetShiftType(), EPSPowerGetShiftType()
317: E*/
318: typedef enum { EPS_POWER_SHIFT_CONSTANT,
319:                EPS_POWER_SHIFT_RAYLEIGH,
320:                EPS_POWER_SHIFT_WILKINSON } EPSPowerShiftType;
321: SLEPC_EXTERN const char *EPSPowerShiftTypes[];

323: SLEPC_EXTERN PetscErrorCode EPSPowerSetShiftType(EPS,EPSPowerShiftType);
324: SLEPC_EXTERN PetscErrorCode EPSPowerGetShiftType(EPS,EPSPowerShiftType*);
325: SLEPC_EXTERN PetscErrorCode EPSPowerSetNonlinear(EPS,PetscBool);
326: SLEPC_EXTERN PetscErrorCode EPSPowerGetNonlinear(EPS,PetscBool*);
327: SLEPC_EXTERN PetscErrorCode EPSPowerSetUpdate(EPS,PetscBool);
328: SLEPC_EXTERN PetscErrorCode EPSPowerGetUpdate(EPS,PetscBool*);
329: SLEPC_EXTERN PetscErrorCode EPSPowerSetSNES(EPS,SNES);
330: SLEPC_EXTERN PetscErrorCode EPSPowerGetSNES(EPS,SNES*);

332: SLEPC_EXTERN PetscErrorCode EPSArnoldiSetDelayed(EPS,PetscBool);
333: SLEPC_EXTERN PetscErrorCode EPSArnoldiGetDelayed(EPS,PetscBool*);

335: SLEPC_EXTERN PetscErrorCode EPSKrylovSchurSetRestart(EPS,PetscReal);
336: SLEPC_EXTERN PetscErrorCode EPSKrylovSchurGetRestart(EPS,PetscReal*);
337: SLEPC_EXTERN PetscErrorCode EPSKrylovSchurSetLocking(EPS,PetscBool);
338: SLEPC_EXTERN PetscErrorCode EPSKrylovSchurGetLocking(EPS,PetscBool*);
339: SLEPC_EXTERN PetscErrorCode EPSKrylovSchurSetPartitions(EPS,PetscInt);
340: SLEPC_EXTERN PetscErrorCode EPSKrylovSchurGetPartitions(EPS,PetscInt*);
341: SLEPC_EXTERN PetscErrorCode EPSKrylovSchurSetDetectZeros(EPS,PetscBool);
342: SLEPC_EXTERN PetscErrorCode EPSKrylovSchurGetDetectZeros(EPS,PetscBool*);
343: SLEPC_EXTERN PetscErrorCode EPSKrylovSchurSetDimensions(EPS,PetscInt,PetscInt,PetscInt);
344: SLEPC_EXTERN PetscErrorCode EPSKrylovSchurGetDimensions(EPS,PetscInt*,PetscInt*,PetscInt*);
345: SLEPC_EXTERN PetscErrorCode EPSKrylovSchurSetSubintervals(EPS,PetscReal*);
346: SLEPC_EXTERN PetscErrorCode EPSKrylovSchurGetSubintervals(EPS,PetscReal**);
347: SLEPC_EXTERN PetscErrorCode EPSKrylovSchurGetInertias(EPS,PetscInt*,PetscReal**,PetscInt**);
348: SLEPC_EXTERN PetscErrorCode EPSKrylovSchurGetSubcommInfo(EPS,PetscInt*,PetscInt*,Vec*);
349: SLEPC_EXTERN PetscErrorCode EPSKrylovSchurGetSubcommPairs(EPS,PetscInt,PetscScalar*,Vec);
350: SLEPC_EXTERN PetscErrorCode EPSKrylovSchurGetSubcommMats(EPS,Mat*,Mat*);
351: SLEPC_EXTERN PetscErrorCode EPSKrylovSchurUpdateSubcommMats(EPS,PetscScalar,PetscScalar,Mat,PetscScalar,PetscScalar, Mat,MatStructure,PetscBool);
352: SLEPC_EXTERN PetscErrorCode EPSKrylovSchurGetKSP(EPS,KSP*);

354: /*E
355:     EPSLanczosReorthogType - determines the type of reorthogonalization
356:     used in the Lanczos method

358:     Level: advanced

360: .seealso: EPSLanczosSetReorthog(), EPSLanczosGetReorthog()
361: E*/
362: typedef enum { EPS_LANCZOS_REORTHOG_LOCAL,
363:                EPS_LANCZOS_REORTHOG_FULL,
364:                EPS_LANCZOS_REORTHOG_SELECTIVE,
365:                EPS_LANCZOS_REORTHOG_PERIODIC,
366:                EPS_LANCZOS_REORTHOG_PARTIAL,
367:                EPS_LANCZOS_REORTHOG_DELAYED } EPSLanczosReorthogType;
368: SLEPC_EXTERN const char *EPSLanczosReorthogTypes[];

370: SLEPC_EXTERN PetscErrorCode EPSLanczosSetReorthog(EPS,EPSLanczosReorthogType);
371: SLEPC_EXTERN PetscErrorCode EPSLanczosGetReorthog(EPS,EPSLanczosReorthogType*);

373: SLEPC_EXTERN PetscErrorCode EPSBlzpackSetBlockSize(EPS,PetscInt);
374: SLEPC_EXTERN PetscErrorCode EPSBlzpackGetBlockSize(EPS,PetscInt*);
375: SLEPC_EXTERN PetscErrorCode EPSBlzpackSetNSteps(EPS,PetscInt);
376: SLEPC_EXTERN PetscErrorCode EPSBlzpackGetNSteps(EPS,PetscInt*);

378: /*E
379:     EPSPRIMMEMethod - determines the method selected in the PRIMME library

381:     Level: advanced

383: .seealso: EPSPRIMMESetMethod(), EPSPRIMMEGetMethod()
384: E*/
385: typedef enum { EPS_PRIMME_DYNAMIC=1,
386:                EPS_PRIMME_DEFAULT_MIN_TIME,
387:                EPS_PRIMME_DEFAULT_MIN_MATVECS,
388:                EPS_PRIMME_ARNOLDI,
389:                EPS_PRIMME_GD,
390:                EPS_PRIMME_GD_PLUSK,
391:                EPS_PRIMME_GD_OLSEN_PLUSK,
392:                EPS_PRIMME_JD_OLSEN_PLUSK,
393:                EPS_PRIMME_RQI,
394:                EPS_PRIMME_JDQR,
395:                EPS_PRIMME_JDQMR,
396:                EPS_PRIMME_JDQMR_ETOL,
397:                EPS_PRIMME_SUBSPACE_ITERATION,
398:                EPS_PRIMME_LOBPCG_ORTHOBASIS,
399:                EPS_PRIMME_LOBPCG_ORTHOBASISW } EPSPRIMMEMethod;
400: SLEPC_EXTERN const char *EPSPRIMMEMethods[];

402: SLEPC_EXTERN PetscErrorCode EPSPRIMMESetBlockSize(EPS,PetscInt);
403: SLEPC_EXTERN PetscErrorCode EPSPRIMMEGetBlockSize(EPS,PetscInt*);
404: SLEPC_EXTERN PetscErrorCode EPSPRIMMESetMethod(EPS,EPSPRIMMEMethod);
405: SLEPC_EXTERN PetscErrorCode EPSPRIMMEGetMethod(EPS,EPSPRIMMEMethod*);

407: SLEPC_EXTERN PetscErrorCode EPSGDSetKrylovStart(EPS,PetscBool);
408: SLEPC_EXTERN PetscErrorCode EPSGDGetKrylovStart(EPS,PetscBool*);
409: SLEPC_EXTERN PetscErrorCode EPSGDSetBlockSize(EPS,PetscInt);
410: SLEPC_EXTERN PetscErrorCode EPSGDGetBlockSize(EPS,PetscInt*);
411: SLEPC_EXTERN PetscErrorCode EPSGDSetRestart(EPS,PetscInt,PetscInt);
412: SLEPC_EXTERN PetscErrorCode EPSGDGetRestart(EPS,PetscInt*,PetscInt*);
413: SLEPC_EXTERN PetscErrorCode EPSGDSetInitialSize(EPS,PetscInt);
414: SLEPC_EXTERN PetscErrorCode EPSGDGetInitialSize(EPS,PetscInt*);
415: SLEPC_EXTERN PetscErrorCode EPSGDSetBOrth(EPS,PetscBool);
416: SLEPC_EXTERN PetscErrorCode EPSGDGetBOrth(EPS,PetscBool*);
417: SLEPC_EXTERN PetscErrorCode EPSGDSetDoubleExpansion(EPS,PetscBool);
418: SLEPC_EXTERN PetscErrorCode EPSGDGetDoubleExpansion(EPS,PetscBool*);

420: SLEPC_EXTERN PetscErrorCode EPSJDSetKrylovStart(EPS,PetscBool);
421: SLEPC_EXTERN PetscErrorCode EPSJDGetKrylovStart(EPS,PetscBool*);
422: SLEPC_EXTERN PetscErrorCode EPSJDSetBlockSize(EPS,PetscInt);
423: SLEPC_EXTERN PetscErrorCode EPSJDGetBlockSize(EPS,PetscInt*);
424: SLEPC_EXTERN PetscErrorCode EPSJDSetRestart(EPS,PetscInt,PetscInt);
425: SLEPC_EXTERN PetscErrorCode EPSJDGetRestart(EPS,PetscInt*,PetscInt*);
426: SLEPC_EXTERN PetscErrorCode EPSJDSetInitialSize(EPS,PetscInt);
427: SLEPC_EXTERN PetscErrorCode EPSJDGetInitialSize(EPS,PetscInt*);
428: SLEPC_EXTERN PetscErrorCode EPSJDSetFix(EPS,PetscReal);
429: SLEPC_EXTERN PetscErrorCode EPSJDGetFix(EPS,PetscReal*);
430: SLEPC_EXTERN PetscErrorCode EPSJDSetConstCorrectionTol(EPS,PetscBool);
431: SLEPC_EXTERN PetscErrorCode EPSJDGetConstCorrectionTol(EPS,PetscBool*);
432: SLEPC_EXTERN PetscErrorCode EPSJDSetBOrth(EPS,PetscBool);
433: SLEPC_EXTERN PetscErrorCode EPSJDGetBOrth(EPS,PetscBool*);

435: SLEPC_EXTERN PetscErrorCode EPSRQCGSetReset(EPS,PetscInt);
436: SLEPC_EXTERN PetscErrorCode EPSRQCGGetReset(EPS,PetscInt*);

438: SLEPC_EXTERN PetscErrorCode EPSLOBPCGSetBlockSize(EPS,PetscInt);
439: SLEPC_EXTERN PetscErrorCode EPSLOBPCGGetBlockSize(EPS,PetscInt*);
440: SLEPC_EXTERN PetscErrorCode EPSLOBPCGSetRestart(EPS,PetscReal);
441: SLEPC_EXTERN PetscErrorCode EPSLOBPCGGetRestart(EPS,PetscReal*);
442: SLEPC_EXTERN PetscErrorCode EPSLOBPCGSetLocking(EPS,PetscBool);
443: SLEPC_EXTERN PetscErrorCode EPSLOBPCGGetLocking(EPS,PetscBool*);

445: /*E
446:     EPSCISSQuadRule - determines the quadrature rule in the CISS solver

448:     Level: advanced

450: .seealso: EPSCISSSetQuadRule(), EPSCISSGetQuadRule()
451: E*/
452: typedef enum { EPS_CISS_QUADRULE_TRAPEZOIDAL=1,
453:                EPS_CISS_QUADRULE_CHEBYSHEV } EPSCISSQuadRule;
454: SLEPC_EXTERN const char *EPSCISSQuadRules[];

456: /*E
457:     EPSCISSExtraction - determines the extraction technique in the CISS solver

459:     Level: advanced

461: .seealso: EPSCISSSetExtraction(), EPSCISSGetExtraction()
462: E*/
463: typedef enum { EPS_CISS_EXTRACTION_RITZ,
464:                EPS_CISS_EXTRACTION_HANKEL } EPSCISSExtraction;
465: SLEPC_EXTERN const char *EPSCISSExtractions[];

467: SLEPC_EXTERN PetscErrorCode EPSCISSSetExtraction(EPS,EPSCISSExtraction);
468: SLEPC_EXTERN PetscErrorCode EPSCISSGetExtraction(EPS,EPSCISSExtraction*);
469: SLEPC_EXTERN PetscErrorCode EPSCISSSetQuadRule(EPS,EPSCISSQuadRule);
470: SLEPC_EXTERN PetscErrorCode EPSCISSGetQuadRule(EPS,EPSCISSQuadRule*);
471: SLEPC_EXTERN PetscErrorCode EPSCISSSetSizes(EPS,PetscInt,PetscInt,PetscInt,PetscInt,PetscInt,PetscBool);
472: SLEPC_EXTERN PetscErrorCode EPSCISSGetSizes(EPS,PetscInt*,PetscInt*,PetscInt*,PetscInt*,PetscInt*,PetscBool*);
473: SLEPC_EXTERN PetscErrorCode EPSCISSSetThreshold(EPS,PetscReal,PetscReal);
474: SLEPC_EXTERN PetscErrorCode EPSCISSGetThreshold(EPS,PetscReal*,PetscReal*);
475: SLEPC_EXTERN PetscErrorCode EPSCISSSetRefinement(EPS,PetscInt,PetscInt);
476: SLEPC_EXTERN PetscErrorCode EPSCISSGetRefinement(EPS,PetscInt*,PetscInt*);
477: SLEPC_EXTERN PetscErrorCode EPSCISSSetUseST(EPS,PetscBool);
478: SLEPC_EXTERN PetscErrorCode EPSCISSGetUseST(EPS,PetscBool*);
479: SLEPC_EXTERN PetscErrorCode EPSCISSGetKSPs(EPS,PetscInt*,KSP**);

481: SLEPC_EXTERN PetscErrorCode EPSLyapIISetLME(EPS,LME);
482: SLEPC_EXTERN PetscErrorCode EPSLyapIIGetLME(EPS,LME*);
483: SLEPC_EXTERN PetscErrorCode EPSLyapIISetRanks(EPS,PetscInt,PetscInt);
484: SLEPC_EXTERN PetscErrorCode EPSLyapIIGetRanks(EPS,PetscInt*,PetscInt*);

486: SLEPC_EXTERN PetscErrorCode EPSBLOPEXSetBlockSize(EPS,PetscInt);
487: SLEPC_EXTERN PetscErrorCode EPSBLOPEXGetBlockSize(EPS,PetscInt*);

489: /*E
490:     EPSEVSLDOSMethod - the method to approximate the density of states (DOS) in the EVSL solver

492:     Level: advanced

494: .seealso: EPSEVSLSetDOSParameters(), EPSEVSLGetDOSParameters()
495: E*/
496: typedef enum { EPS_EVSL_DOS_KPM,
497:                EPS_EVSL_DOS_LANCZOS } EPSEVSLDOSMethod;
498: SLEPC_EXTERN const char *EPSEVSLDOSMethods[];

500: /*E
501:     EPSEVSLDamping - the damping type used in the EVSL solver

503:     Level: advanced

505: .seealso: EPSEVSLSetDOSParameters(), EPSEVSLGetDOSParameters()
506: E*/
507: typedef enum { EPS_EVSL_DAMPING_NONE,
508:                EPS_EVSL_DAMPING_JACKSON,
509:                EPS_EVSL_DAMPING_SIGMA } EPSEVSLDamping;
510: SLEPC_EXTERN const char *EPSEVSLDampings[];

512: SLEPC_EXTERN PetscErrorCode EPSEVSLSetRange(EPS,PetscReal,PetscReal);
513: SLEPC_EXTERN PetscErrorCode EPSEVSLGetRange(EPS,PetscReal*,PetscReal*);
514: SLEPC_EXTERN PetscErrorCode EPSEVSLSetSlices(EPS,PetscInt);
515: SLEPC_EXTERN PetscErrorCode EPSEVSLGetSlices(EPS,PetscInt*);
516: SLEPC_EXTERN PetscErrorCode EPSEVSLSetDOSParameters(EPS,EPSEVSLDOSMethod,PetscInt,PetscInt,PetscInt,PetscInt);
517: SLEPC_EXTERN PetscErrorCode EPSEVSLGetDOSParameters(EPS,EPSEVSLDOSMethod*,PetscInt*,PetscInt*,PetscInt*,PetscInt*);
518: SLEPC_EXTERN PetscErrorCode EPSEVSLSetPolParameters(EPS,PetscInt,PetscReal);
519: SLEPC_EXTERN PetscErrorCode EPSEVSLGetPolParameters(EPS,PetscInt*,PetscReal*);
520: SLEPC_EXTERN PetscErrorCode EPSEVSLSetDamping(EPS,EPSEVSLDamping);
521: SLEPC_EXTERN PetscErrorCode EPSEVSLGetDamping(EPS,EPSEVSLDamping*);

523: SLEPC_EXTERN PetscErrorCode EPSFEASTSetNumPoints(EPS,PetscInt);
524: SLEPC_EXTERN PetscErrorCode EPSFEASTGetNumPoints(EPS,PetscInt*);

526: #endif