Article View: comp.lang.fortran
Article #4574Re: inverse matrix
From: Walt Brainerd
Date: Fri, 08 Jan 1993 17:11
Date: Fri, 08 Jan 1993 17:11
22 lines
1533 bytes
1533 bytes
In article <1993Jan8.134612.18468@news.unige.ch>, pfennige@scsun.unige.ch (PFENNIGER Daniel) writes: > > Is it true that a Fortran-90 programme inverting matrix A would be: > B = 1./A > ? No; it is important to keep in mind the model that Fortran 90 does _array processing_, not matrix arithmetic. 1./A produces an array with the same shape of A in which each element is the reciprocal of the corresponding element of A. Similarly, A * B is element-by-element multiplication, not matrix multiplication. In Fortran 90, MATMUL(A, B) is an intrinsic function, but there is no inversion intrinsic. However, in Fortran 90 it is possible to define a new data type, called "MATRIX", say, and define A * B to be matrix multiplication and perhaps a unary operation .INVERSE. -- Walt Brainerd walt@netcom.com Unicomp, Inc. +1-415-949-4052 235 Mt. Hamilton Ave. +1-415-949-4058 (fax) Los Altos, CA 94022 USA "F90--Fortran for the '90s"
Technical Details
Message-ID:
<1993Jan8.171121.13855@netcom.com>Path:
arc1-ca.pugleaf.net!sparky!uunet!cis.ohio-state.edu!zaphod.mps.ohio-state.edu!sol.ctr.columbia.edu!spool.mu.edu!umn.edu!csus.edu!netcom.com!waltReferences:
<C0I49C.Jrr@athena.cs.uga.edu> <1993Jan8.100315.5851@edf.fr> <1993Jan8.134612.18468@news.unige.ch>Show Full Headers
Complete Article Headers
Newsgroups: comp.lang.fortran Path: arc1-ca.pugleaf.net!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 From: walt@netcom.com (Walt Brainerd) Subject: Re: inverse matrix Message-ID: <1993Jan8.171121.13855@netcom.com> Summary: No Organization: Netcom - Online Communication Services (408 241-9760 guest) References: <C0I49C.Jrr@athena.cs.uga.edu> <1993Jan8.100315.5851@edf.fr> <1993Jan8.134612.18468@news.unige.ch> Date: Fri, 8 Jan 1993 17:11:21 GMT Lines: 22