Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
| Rev 1073 | Rev 1087 |
|---|---|
| Line 28... | Line 28... |
! eps eigenproblem solver context |
! eps eigenproblem solver context |
Mat A |
Mat A |
EPS eps |
EPS eps |
EPSType type |
EPSType type |
PetscReal tol, error, re, im |
PetscReal tol, error |
PetscScalar kr, ki |
PetscScalar kr, ki |
integer rank, n, nev, ierr, maxit, i, its, nconv |
integer rank, n, nev, ierr, maxit, i, its, nconv |
integer col(3), Istart, Iend |
integer col(3), Istart, Iend |
PetscTruth flg |
PetscTruth flg |
PetscScalar value(3) |
PetscScalar value(3) |
| Line 154... | Line 154... |
call EPSGetEigenpair(eps,i,kr,ki,PETSC_NULL,PETSC_NULL,ierr) |
call EPSGetEigenpair(eps,i,kr,ki,PETSC_NULL,PETSC_NULL,ierr) |
! ** Compute the relative error associated to each eigenpair |
! ** Compute the relative error associated to each eigenpair |
call EPSComputeRelativeError(eps,i,error,ierr) |
call EPSComputeRelativeError(eps,i,error,ierr) |
#if defined(PETSC_USE_COMPLEX) |
write(*,160) PetscRealPart(kr), error |
re = PetscRealPart(kr); |
160 format (1P,' ',E12.4,' ',E12.4) |
im = PetscImaginaryPart(kr); |
|
#else |
|
re = kr; |
|
im = ki; |
|
#endif |
|
if (im.ne.0.D0) then |
|
write(*,160) re, im, error |
|
160 format (1P,E11.4,E11.4,' j ',E12.4) |
|
else |
|
write(*,170) re, error |
|
170 format (1P,' ',E12.4,' ',E12.4) |
|
endif |
|
enddo |
enddo |
write(*,*) |
write(*,*) |
endif |
endif |