| Line 34... |
Line 34... |
EPS eps; /* eigenproblem solver context */
|
EPS eps; /* eigenproblem solver context */
|
const EPSType type;
|
const EPSType type;
|
PetscReal error,tol,re,im;
|
PetscReal error,tol,re,im;
|
PetscScalar kr,ki;
|
PetscScalar kr,ki;
|
PetscInt nev,maxit,i,its,nconv;
|
PetscInt nev,maxit,i,its,nconv;
|
char filename[256];
|
char filename[PETSC_MAX_PATH_LEN];
|
PetscViewer viewer;
|
PetscViewer viewer;
|
PetscBool flg;
|
PetscBool flg;
|
PetscErrorCode ierr;
|
PetscErrorCode ierr;
|
|
|
SlepcInitialize(&argc,&argv,(char*)0,help);
|
SlepcInitialize(&argc,&argv,(char*)0,help);
|
| Line 46... |
Line 46... |
/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
Load the operator matrix that defines the eigensystem, Ax=kx
|
Load the operator matrix that defines the eigensystem, Ax=kx
|
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
|
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
|
|
|
ierr = PetscPrintf(PETSC_COMM_WORLD,"\nEigenproblem stored in file.\n\n");CHKERRQ(ierr);
|
ierr = PetscPrintf(PETSC_COMM_WORLD,"\nEigenproblem stored in file.\n\n");CHKERRQ(ierr);
|
ierr = PetscOptionsGetString(PETSC_NULL,"-file",filename,256,&flg);CHKERRQ(ierr);
|
ierr = PetscOptionsGetString(PETSC_NULL,"-file",filename,PETSC_MAX_PATH_LEN,&flg);CHKERRQ(ierr);
|
if (!flg) {
|
if (!flg) {
|
SETERRQ(PETSC_COMM_WORLD,1,"Must indicate a file name with the -file option.");
|
SETERRQ(PETSC_COMM_WORLD,1,"Must indicate a file name with the -file option.");
|
}
|
}
|
|
|
#if defined(PETSC_USE_COMPLEX)
|
#if defined(PETSC_USE_COMPLEX)
|