home *** CD-ROM | disk | FTP | other *** search
/ Geek Gadgets 1 / ADE-1.bin / ade-dist / eispack-1.0-src.tgz / tar.out / contrib / eispack / rt.f < prev    next >
Text File  |  1996-09-28  |  2KB  |  65 lines
  1.       subroutine rt(nm,n,a,w,matz,z,fv1,ierr)
  2. c
  3.       integer n,nm,ierr,matz
  4.       double precision a(nm,3),w(n),z(nm,n),fv1(n)
  5. c
  6. c     this subroutine calls the recommended sequence of
  7. c     subroutines from the eigensystem subroutine package (eispack)
  8. c     to find the eigenvalues and eigenvectors (if desired)
  9. c     of a special real tridiagonal matrix.
  10. c
  11. c     on input
  12. c
  13. c        nm  must be set to the row dimension of the two-dimensional
  14. c        array parameters as declared in the calling program
  15. c        dimension statement.
  16. c
  17. c        n  is the order of the matrix  a.
  18. c
  19. c        a  contains the special real tridiagonal matrix in its
  20. c        first three columns.  the subdiagonal elements are stored
  21. c        in the last  n-1  positions of the first column, the
  22. c        diagonal elements in the second column, and the superdiagonal
  23. c        elements in the first  n-1  positions of the third column.
  24. c        elements  a(1,1)  and  a(n,3)  are arbitrary.
  25. c
  26. c        matz  is an integer variable set equal to zero if
  27. c        only eigenvalues are desired.  otherwise it is set to
  28. c        any non-zero integer for both eigenvalues and eigenvectors.
  29. c
  30. c     on output
  31. c
  32. c        w  contains the eigenvalues in ascending order.
  33. c
  34. c        z  contains the eigenvectors if matz is not zero.
  35. c
  36. c        ierr  is an integer output variable set equal to an error
  37. c           completion code described in the documentation for imtql1
  38. c           and imtql2.  the normal completion code is zero.
  39. c
  40. c        fv1  is a temporary storage array.
  41. c
  42. c     questions and comments should be directed to burton s. garbow,
  43. c     mathematics and computer science div, argonne national laboratory
  44. c
  45. c     this version dated august 1983.
  46. c
  47. c     ------------------------------------------------------------------
  48. c
  49.       if (n .le. nm) go to 10
  50.       ierr = 10 * n
  51.       go to 50
  52. c
  53.    10 if (matz .ne. 0) go to 20
  54. c     .......... find eigenvalues only ..........
  55.       call  figi(nm,n,a,w,fv1,fv1,ierr)
  56.       if (ierr .gt. 0) go to 50
  57.       call  imtql1(n,w,fv1,ierr)
  58.       go to 50
  59. c     .......... find both eigenvalues and eigenvectors ..........
  60.    20 call  figi2(nm,n,a,w,fv1,z,ierr)
  61.       if (ierr .ne. 0) go to 50
  62.       call  imtql2(nm,n,w,fv1,z,ierr)
  63.    50 return
  64.       end
  65.