🚀 go-pugleaf

In article <1993Jan7.235543.25287@murdoch.acc.Virginia.EDU> ss4i@kelvin.seas.Virginia.EDU (Sriram  Srinivasan) writes:
#Hi,
#
#   I had the following question about the BALANC routine from EISPACK.
#
#   Given a matrix A, BALANC gives the balanced matrix
#
#          Ab = T*A*inv(T)
#
#   and also some information in the variables LOW, IGH, and SCALE which
#   determine the similarity transformation T.
#
#   Does anybody know how to construct a T explicity, for all A?
#
#   Thanks.                       
#
#   Sriram (spoiled by MATLAB)
#
#   (ss4i@kelvin.seas.virginia.edu)

The matrix T is constructed iteratively by balancing successively
each row and its corresponding column until the process converges.
Details can be found in

   @article{pare:69,
      author = "B. N. Parlett and C. Reinsch",
      year = "1969",
      title = "Balancing a Matrix for Calculation of Eigenvalues
               and Eigenvectors",
      journal = "Numerische Mathematik",
      volume = "13",
      pages = "293--304",
      kwds = "nla, eig, scaling, software"
   }

Pete Stewart