home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1993 #1 / NN_1993_1.iso / spool / comp / lang / fortran / 4938 < prev    next >
Encoding:
Text File  |  1993-01-08  |  1.7 KB  |  44 lines
  1. Newsgroups: comp.lang.fortran
  2. Path: sparky!uunet!zaphod.mps.ohio-state.edu!cs.utexas.edu!sun-barr!ames!news.dfrf.nasa.gov!altair.dfrf.nasa.gov!maine
  3. From: maine@altair.dfrf.nasa.gov (Richard Maine)
  4. Subject: Re: inverse matrix
  5. In-Reply-To: pfennige@scsun.unige.ch's message of Fri, 8 Jan 1993 13:46:12 GMT
  6. Message-ID: <MAINE.93Jan8090104@altair.dfrf.nasa.gov>
  7. Sender: news@news.dfrf.nasa.gov (Usenet news)
  8. Organization: NASA Dryden, Edwards, Cal.
  9. References: <C0I49C.Jrr@athena.cs.uga.edu> <1993Jan8.100315.5851@edf.fr>
  10.     <1993Jan8.134612.18468@news.unige.ch>
  11. Date: Fri, 8 Jan 1993 17:01:10 GMT
  12. Lines: 30
  13.  
  14. On Fri, 8 Jan 1993 13:46:12 GMT, pfennige@scsun.unige.ch (PFENNIGER Daniel) said:
  15.  
  16.  
  17. PFENNIGER> Is it true that a Fortran-90 programme inverting matrix A would be:
  18. PFENNIGER>       B = 1./A
  19.  
  20. No.  That would be a legal statement, but it makes each element of B
  21. the inverse of the corresponding element of A.  The f90 arrays are
  22. better thought of as "arrays" than as "matrices".  In particular,
  23. all the operations are element-by-element.  There is a standard
  24. matmul subroutine for the usual matrix multiplication, but there is
  25. no matrix inversion.
  26.  
  27. I wasn't on the committee at the time, but I believe the rationale was
  28. that there was no single matrix inversion technique good for all
  29. situations, and that numerically-oriented users were better off
  30. choosing a particular, appropriate one than whatever the vendor
  31. happened to implement (which might be different on different
  32. platforms).
  33.  
  34. You (or a library writer) can make a "matrix" (or whatever you want to
  35. call it) user-defined type and define division for it such that the
  36. above sample line would do the matrix inversion, but such a type
  37. isn't part of the base language.
  38.  
  39. --
  40. --
  41. Richard Maine
  42. maine@altair.dfrf.nasa.gov
  43.  
  44.