| Line 10... |
Line 10... |
|
|
#undef __FUNCT__
|
#undef __FUNCT__
|
#define __FUNCT__ "EPSSolve_TS_ARNOLDI"
|
#define __FUNCT__ "EPSSolve_TS_ARNOLDI"
|
PetscErrorCode EPSSolve_TS_ARNOLDI(EPS eps)
|
PetscErrorCode EPSSolve_TS_ARNOLDI(EPS eps)
|
{
|
{
|
#if defined(SLEPC_MISSING_LAPACK_TREVC)
|
|
PetscFunctionBegin;
|
|
SETERRQ(PETSC_ERR_SUP,"TREVC - Lapack routine is unavailable.");
|
|
#else
|
|
PetscErrorCode ierr;
|
PetscErrorCode ierr;
|
int i,k,ncv=eps->ncv;
|
int i,k,ncv=eps->ncv;
|
Vec fr=eps->work[0];
|
Vec fr=eps->work[0];
|
Vec fl=eps->work[1];
|
Vec fl=eps->work[1];
|
Vec *Qr=eps->V, *Ql=eps->W;
|
Vec *Qr=eps->V, *Ql=eps->W;
|
| Line 102... |
Line 98... |
ierr = PetscFree(Ul);CHKERRQ(ierr);
|
ierr = PetscFree(Ul);CHKERRQ(ierr);
|
ierr = PetscFree(work);CHKERRQ(ierr);
|
ierr = PetscFree(work);CHKERRQ(ierr);
|
ierr = PetscFree(eigr);CHKERRQ(ierr);
|
ierr = PetscFree(eigr);CHKERRQ(ierr);
|
ierr = PetscFree(eigi);CHKERRQ(ierr);
|
ierr = PetscFree(eigi);CHKERRQ(ierr);
|
PetscFunctionReturn(0);
|
PetscFunctionReturn(0);
|
#endif
|
|
}
|
}
|
|
|
|
|
|
|