| Line 149... |
Line 149... |
switch (qep->problem_type) {
|
switch (qep->problem_type) {
|
case QEP_GENERAL: type = "general quadratic eigenvalue problem"; break;
|
case QEP_GENERAL: type = "general quadratic eigenvalue problem"; break;
|
case QEP_HERMITIAN: type = HERM " quadratic eigenvalue problem"; break;
|
case QEP_HERMITIAN: type = HERM " quadratic eigenvalue problem"; break;
|
case QEP_GYROSCOPIC: type = "gyroscopic quadratic eigenvalue problem"; break;
|
case QEP_GYROSCOPIC: type = "gyroscopic quadratic eigenvalue problem"; break;
|
case 0: type = "not yet set"; break;
|
case 0: type = "not yet set"; break;
|
default: SETERRQ(1,"Wrong value of qep->problem_type");
|
default: SETERRQ(((PetscObject)qep)->comm,1,"Wrong value of qep->problem_type");
|
}
|
}
|
ierr = PetscViewerASCIIPrintf(viewer," problem type: %s\n",type);CHKERRQ(ierr);
|
ierr = PetscViewerASCIIPrintf(viewer," problem type: %s\n",type);CHKERRQ(ierr);
|
ierr = QEPGetType(qep,&type);CHKERRQ(ierr);
|
ierr = QEPGetType(qep,&type);CHKERRQ(ierr);
|
if (type) {
|
if (type) {
|
ierr = PetscViewerASCIIPrintf(viewer," method: %s\n",type);CHKERRQ(ierr);
|
ierr = PetscViewerASCIIPrintf(viewer," method: %s\n",type);CHKERRQ(ierr);
|
| Line 180... |
Line 180... |
ierr = PetscViewerASCIIPrintf(viewer,"largest imaginary parts\n");CHKERRQ(ierr);
|
ierr = PetscViewerASCIIPrintf(viewer,"largest imaginary parts\n");CHKERRQ(ierr);
|
break;
|
break;
|
case QEP_SMALLEST_IMAGINARY:
|
case QEP_SMALLEST_IMAGINARY:
|
ierr = PetscViewerASCIIPrintf(viewer,"smallest imaginary parts\n");CHKERRQ(ierr);
|
ierr = PetscViewerASCIIPrintf(viewer,"smallest imaginary parts\n");CHKERRQ(ierr);
|
break;
|
break;
|
default: SETERRQ(1,"Wrong value of qep->which");
|
default: SETERRQ(((PetscObject)qep)->comm,1,"Wrong value of qep->which");
|
}
|
}
|
if (qep->leftvecs) {
|
if (qep->leftvecs) {
|
ierr = PetscViewerASCIIPrintf(viewer," computing left eigenvectors also\n");CHKERRQ(ierr);
|
ierr = PetscViewerASCIIPrintf(viewer," computing left eigenvectors also\n");CHKERRQ(ierr);
|
}
|
}
|
ierr = PetscViewerASCIIPrintf(viewer," number of eigenvalues (nev): %d\n",qep->nev);CHKERRQ(ierr);
|
ierr = PetscViewerASCIIPrintf(viewer," number of eigenvalues (nev): %d\n",qep->nev);CHKERRQ(ierr);
|
| Line 341... |
Line 341... |
qep->data = 0;
|
qep->data = 0;
|
}
|
}
|
|
|
ierr = PetscFListFind(QEPList,((PetscObject)qep)->comm,type,(void (**)(void))&r);CHKERRQ(ierr);
|
ierr = PetscFListFind(QEPList,((PetscObject)qep)->comm,type,(void (**)(void))&r);CHKERRQ(ierr);
|
|
|
if (!r) SETERRQ1(1,"Unknown QEP type given: %s",type);
|
if (!r) SETERRQ1(((PetscObject)qep)->comm,1,"Unknown QEP type given: %s",type);
|
|
|
qep->setupcalled = 0;
|
qep->setupcalled = 0;
|
ierr = PetscMemzero(qep->ops,sizeof(struct _QEPOps));CHKERRQ(ierr);
|
ierr = PetscMemzero(qep->ops,sizeof(struct _QEPOps));CHKERRQ(ierr);
|
ierr = (*r)(qep);CHKERRQ(ierr);
|
ierr = (*r)(qep);CHKERRQ(ierr);
|
|
|