home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1993 #1 / NN_1993_1.iso / spool / comp / lang / fortran / 4939 < prev    next >
Encoding:
Text File  |  1993-01-08  |  1.5 KB  |  34 lines
  1. Newsgroups: comp.lang.fortran
  2. Path: sparky!uunet!cis.ohio-state.edu!zaphod.mps.ohio-state.edu!sol.ctr.columbia.edu!spool.mu.edu!umn.edu!csus.edu!netcom.com!walt
  3. From: walt@netcom.com (Walt Brainerd)
  4. Subject: Re: inverse matrix
  5. Message-ID: <1993Jan8.171121.13855@netcom.com>
  6. Summary: No
  7. Organization: Netcom - Online Communication Services  (408 241-9760 guest) 
  8. References:  <C0I49C.Jrr@athena.cs.uga.edu> <1993Jan8.100315.5851@edf.fr> <1993Jan8.134612.18468@news.unige.ch>
  9. Date: Fri, 8 Jan 1993 17:11:21 GMT
  10. Lines: 22
  11.  
  12. In article <1993Jan8.134612.18468@news.unige.ch>, pfennige@scsun.unige.ch (PFENNIGER Daniel) writes:
  13. > 
  14. > Is it true that a Fortran-90 programme inverting matrix A would be:
  15. >       B = 1./A
  16. > ?
  17.  
  18. No; it is important to keep in mind the model that Fortran 90
  19. does _array processing_, not matrix arithmetic.  1./A produces
  20. an array with the same shape of A in which each element is
  21. the reciprocal of the corresponding element of A.  Similarly,
  22. A * B is element-by-element multiplication, not matrix
  23. multiplication.  In Fortran 90, MATMUL(A, B) is an intrinsic
  24. function, but there is no inversion intrinsic.
  25.  
  26. However, in Fortran 90 it is possible to define a new data type,
  27. called "MATRIX", say, and define A * B to be matrix multiplication
  28. and perhaps a unary operation .INVERSE.
  29. -- 
  30. Walt Brainerd                walt@netcom.com
  31. Unicomp, Inc.                +1-415-949-4052
  32. 235 Mt. Hamilton Ave.        +1-415-949-4058 (fax)
  33. Los Altos, CA 94022 USA      "F90--Fortran for the '90s"
  34.