| Line 173... |
Line 173... |
PetscErrorCode ierr;
|
PetscErrorCode ierr;
|
PetscInt i,k,l,lwork,nv;
|
PetscInt i,k,l,lwork,nv;
|
Vec u=eps->work[0];
|
Vec u=eps->work[0];
|
PetscScalar *S=eps->T,*Q,*g,*work;
|
PetscScalar *S=eps->T,*Q,*g,*work;
|
PetscReal beta,gnorm;
|
PetscReal beta,gnorm;
|
PetscTruth breakdown;
|
PetscBool breakdown;
|
|
|
PetscFunctionBegin;
|
PetscFunctionBegin;
|
ierr = PetscMemzero(S,eps->ncv*eps->ncv*sizeof(PetscScalar));CHKERRQ(ierr);
|
ierr = PetscMemzero(S,eps->ncv*eps->ncv*sizeof(PetscScalar));CHKERRQ(ierr);
|
ierr = PetscMalloc(eps->ncv*eps->ncv*sizeof(PetscScalar),&Q);CHKERRQ(ierr);
|
ierr = PetscMalloc(eps->ncv*eps->ncv*sizeof(PetscScalar),&Q);CHKERRQ(ierr);
|
lwork = PetscMax((eps->ncv+1)*eps->ncv,7*eps->ncv);
|
lwork = PetscMax((eps->ncv+1)*eps->ncv,7*eps->ncv);
|