home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1993 #1 / NN_1993_1.iso / spool / comp / lang / fortran / 4941 < prev    next >
Encoding:
Internet Message Format  |  1993-01-08  |  1.2 KB
  1. Path: sparky!uunet!paladin.american.edu!gatech!rpi!uwm.edu!psuvax1!psuvm!hdk
  2. Organization: Penn State University
  3. Date: Fri, 8 Jan 1993 12:54:09 EST
  4. From: H. D. Knoble <HDK@psuvm.psu.edu>
  5. Message-ID: <93008.125409HDK@psuvm.psu.edu>
  6. Newsgroups: comp.lang.fortran
  7. Subject: Re: inverse matrix
  8. References: <C0I49C.Jrr@athena.cs.uga.edu>
  9. Lines: 15
  10.  
  11. Regardless of the kind of matrix that is being dealt with, the code
  12. (e.g., IMSL, NAG) should check for computational singularity (or in
  13. the case of generalized inverse check that the generalized inverse is
  14. also a least squares inverse and set an appropriate flag). The code
  15. should also effect as accurate an inverse as feasible on each platform
  16. it is run on. Else one can get an "inverse" that is just plain wrong
  17. too easily.
  18.  
  19. If this matrix inversion is for an important project, it might be wise
  20. to stay away from Numerical Recipes algorithms (and textbook code
  21. in general with a few exceptions) as they are not usually robust in
  22. the above sense. What this amounts to is that good algorithms are most often
  23. more than taking the algebra and applying FORmula TRANslation to it.
  24. So know WHO you're getting code from, and what algorithm(s) are being
  25. used.
  26.