| Line 98... |
Line 98... |
/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
Solve for interior eigenvalues (target=2.1)
|
Solve for interior eigenvalues (target=2.1)
|
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
|
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
|
ierr = EPSSetWhichEigenpairs(eps,EPS_TARGET_MAGNITUDE);CHKERRQ(ierr);
|
ierr = EPSSetWhichEigenpairs(eps,EPS_TARGET_MAGNITUDE);CHKERRQ(ierr);
|
ierr = EPSSetTarget(eps,2.1);CHKERRQ(ierr);
|
ierr = EPSSetTarget(eps,2.1);CHKERRQ(ierr);
|
ierr = PetscTypeCompare((PetscObject)eps,EPSLANCZOS,&flg);CHKERRQ(ierr);
|
ierr = PetscObjectTypeCompare((PetscObject)eps,EPSLANCZOS,&flg);CHKERRQ(ierr);
|
if (flg) {
|
if (flg) {
|
ierr = EPSGetST(eps,&st);CHKERRQ(ierr);
|
ierr = EPSGetST(eps,&st);CHKERRQ(ierr);
|
ierr = STSetType(st,STSINVERT);CHKERRQ(ierr);
|
ierr = STSetType(st,STSINVERT);CHKERRQ(ierr);
|
} else {
|
} else {
|
ierr = PetscTypeCompare((PetscObject)eps,EPSKRYLOVSCHUR,&flg);CHKERRQ(ierr);
|
ierr = PetscObjectTypeCompare((PetscObject)eps,EPSKRYLOVSCHUR,&flg);CHKERRQ(ierr);
|
if (!flg) {
|
if (!flg) {
|
ierr = PetscTypeCompare((PetscObject)eps,EPSARNOLDI,&flg);CHKERRQ(ierr);
|
ierr = PetscObjectTypeCompare((PetscObject)eps,EPSARNOLDI,&flg);CHKERRQ(ierr);
|
}
|
}
|
if (flg) {
|
if (flg) {
|
ierr = EPSSetExtraction(eps,EPS_HARMONIC);CHKERRQ(ierr);
|
ierr = EPSSetExtraction(eps,EPS_HARMONIC);CHKERRQ(ierr);
|
}
|
}
|
}
|
}
|