| SETERRQ(PETSC_COMM_SELF,PETSC_ERR_SUP,"STEVR - Lapack routine is unavailable."); |
| #else |
| PetscErrorCode ierr; |
| PetscInt i,j; |
| PetscInt i; |
| PetscBLASInt n1,n2,n3,lwork,liwork,info,l,n,m,ld,off,il,iu,*isuppz; |
| PetscScalar *A,*Q,*W,one=1.0,zero=0.0; |
| PetscReal *d,*e,abstol=0.0,vl,vu; |
| #if defined(PETSC_USE_COMPLEX) |
| PetscInt j; |
| PetscReal *ritz; |
| #endif |
| #endif |
| norm = BLASnrm2_(&n,Y,&inc); |
| #if !defined(PETSC_USE_COMPLEX) |
| tmp = BLASnrm2_(&n,Y+ld,&inc); |
| norm = SlepcAbsEigenvalue(norm,tmp); |
| if (iscomplex) { |
| tmp = BLASnrm2_(&n,Y+ld,&inc); |
| norm = SlepcAbsEigenvalue(norm,tmp); |
| } |
| #endif |
| tmp = 1.0 / norm; |
| BLASscal_(&n,&tmp,Y,&inc); |
| #if !defined(PETSC_USE_COMPLEX) |
| BLASscal_(&n,&tmp,Y+ld,&inc); |
| if (iscomplex) BLASscal_(&n,&tmp,Y+ld,&inc); |
| #endif |
| } |