| Line 20... |
Line 20... |
along with SLEPc. If not, see <http://www.gnu.org/licenses/>.
|
along with SLEPc. If not, see <http://www.gnu.org/licenses/>.
|
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
*/
|
*/
|
|
|
#include <private/epsimpl.h> /*I "slepceps.h" I*/
|
#include <private/epsimpl.h> /*I "slepceps.h" I*/
|
#include "../src/sys/viewer/impls/ascii/asciiimpl.h"
|
|
|
|
PetscFList EPSList = 0;
|
PetscFList EPSList = 0;
|
PetscBool EPSRegisterAllCalled = PETSC_FALSE;
|
PetscBool EPSRegisterAllCalled = PETSC_FALSE;
|
PetscClassId EPS_CLASSID = 0;
|
PetscClassId EPS_CLASSID = 0;
|
PetscLogEvent EPS_SetUp = 0,EPS_Solve = 0,EPS_Dense = 0;
|
PetscLogEvent EPS_SetUp = 0,EPS_Solve = 0,EPS_Dense = 0;
|
| Line 325... |
Line 324... |
PetscBool terse,errok,isascii;
|
PetscBool terse,errok,isascii;
|
PetscReal error,re,im;
|
PetscReal error,re,im;
|
PetscScalar kr,ki;
|
PetscScalar kr,ki;
|
PetscInt i,j;
|
PetscInt i,j;
|
PetscErrorCode ierr;
|
PetscErrorCode ierr;
|
PetscViewer_ASCII *vascii;
|
|
int err;
|
|
|
|
PetscFunctionBegin;
|
PetscFunctionBegin;
|
PetscValidHeaderSpecific(eps,EPS_CLASSID,1);
|
PetscValidHeaderSpecific(eps,EPS_CLASSID,1);
|
if (!viewer) viewer = PETSC_VIEWER_STDOUT_(((PetscObject)eps)->comm);
|
if (!viewer) viewer = PETSC_VIEWER_STDOUT_(((PetscObject)eps)->comm);
|
PetscValidHeaderSpecific(viewer,PETSC_VIEWER_CLASSID,2);
|
PetscValidHeaderSpecific(viewer,PETSC_VIEWER_CLASSID,2);
|
| Line 338... |
Line 335... |
if (!eps->eigr || !eps->eigi || !eps->V) {
|
if (!eps->eigr || !eps->eigi || !eps->V) {
|
SETERRQ(((PetscObject)eps)->comm,PETSC_ERR_ARG_WRONGSTATE,"EPSSolve must be called first");
|
SETERRQ(((PetscObject)eps)->comm,PETSC_ERR_ARG_WRONGSTATE,"EPSSolve must be called first");
|
}
|
}
|
ierr = PetscTypeCompare((PetscObject)viewer,PETSCVIEWERASCII,&isascii);CHKERRQ(ierr);
|
ierr = PetscTypeCompare((PetscObject)viewer,PETSCVIEWERASCII,&isascii);CHKERRQ(ierr);
|
if (!isascii) PetscFunctionReturn(0);
|
if (!isascii) PetscFunctionReturn(0);
|
vascii = (PetscViewer_ASCII *)viewer->data;
|
|
err = fflush(vascii->fd);
|
|
//ierr = PetscViewerFlush(viewer);CHKERRQ(ierr);
|
|
|
|
ierr = PetscOptionsHasName(PETSC_NULL,"-eps_terse",&terse);CHKERRQ(ierr);
|
ierr = PetscOptionsHasName(PETSC_NULL,"-eps_terse",&terse);CHKERRQ(ierr);
|
if (terse) {
|
if (terse) {
|
if (eps->nconv<eps->nev) {
|
if (eps->nconv<eps->nev) {
|
ierr = PetscViewerASCIIPrintf(viewer," Problem: less than %D eigenvalues converged\n\n",eps->nev);CHKERRQ(ierr);
|
ierr = PetscViewerASCIIPrintf(viewer," Problem: less than %D eigenvalues converged\n\n",eps->nev);CHKERRQ(ierr);
|