home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1993 #1 / NN_1993_1.iso / spool / comp / lang / fortran / 4964 < prev    next >
Encoding:
Text File  |  1993-01-11  |  1.8 KB  |  38 lines

  1. Newsgroups: comp.lang.fortran
  2. Path: sparky!uunet!zaphod.mps.ohio-state.edu!darwin.sura.net!news.udel.edu!perelandra.cms.udel.edu!mccalpin
  3. From: mccalpin@perelandra.cms.udel.edu (John D. McCalpin)
  4. Subject: Re: inverse matrix
  5. Message-ID: <C0p5JD.M0L@news.udel.edu>
  6. Sender: usenet@news.udel.edu
  7. Nntp-Posting-Host: perelandra.cms.udel.edu
  8. Organization: College of Marine Studies, U. Del.
  9. References: <C0I49C.Jrr@athena.cs.uga.edu> <93008.125409HDK@psuvm.psu.edu> <1993Jan8.201645.14915@news.eng.convex.com>
  10. Date: Mon, 11 Jan 1993 15:54:48 GMT
  11. Lines: 25
  12.  
  13. In article <1993Jan8.201645.14915@news.eng.convex.com> dodson@convex.COM (Dave Dodson) writes:
  14. >I'd like to point out that it is almost never required or desirable to
  15. >compute the inverse of a matrix.  Almost without exception, you can do
  16. >any computation in which you would use the inverse in a better way that
  17. >does not use the inverse.  By 'better' I mean faster, uses less memory,
  18. >more accurate, etc.
  19.  
  20. The direct use of the inverse matrix is generally the fastest way to
  21. solve a dense system of equations with multiple, consecutive right-hand-sides
  22. (as in a time-dependent fluid dynamics problem).
  23.  
  24. On the Cray Y series, for example, multiplying by the inverse matrix is
  25. always (?) faster than performing the forward-backward substitution
  26. step for each new RHS.  This is true even though both operations have
  27. the same computational complexity in terms of operation counts (=2 N^2).
  28. The difference is that matrix multiplication is a significantly
  29. simpler algorithm with less overhead than forward-backward substitution.   
  30.  
  31. I suspect that the same sorts of speedups would be observed on Dodson's
  32. Convex machines....
  33. -- 
  34. --
  35. John D. McCalpin                        mccalpin@perelandra.cms.udel.edu
  36. Assistant Professor                     mccalpin@brahms.udel.edu
  37. College of Marine Studies, U. Del.      John.McCalpin@mvs.udel.edu
  38.