This is the mail archive of the
fortran@gcc.gnu.org
mailing list for the GNU Fortran project.
Re: matmul intrinsic procedure
- From: Richard E Maine <Richard dot Maine at nasa dot gov>
- To: Jan Stanisław Owoc <jsowoc at ucalgary dot ca>
- Cc: fortran at gcc dot gnu dot org
- Date: Thu, 30 Jun 2005 13:31:29 -0700
- Subject: Re: matmul intrinsic procedure
- References: <Pine.LNX.4.44.0506301400030.27330-100000@c1.ucalgary.ca>
On Jun 30, 2005, at 1:11 PM, Jan Stanisław Owoc wrote:
If there was some way I could get an exact specification for the
intrinsic
procedure "matmul",
It is in the standard, or pretty much any textbook. As specs go, it is
pretty short. My informal summary is
1. You presumably know what matrix multiplication means. If not, it
spells it out.
2. It must handle rank2*rank2, rank1*rank2, and rank2*rank1.
3. It must work with all types and kinds of numeric arguments (integer,
real, and complex, with all KIND values), with the "usual" rules for
handling type mixing. (You really ought to know the "usual" rules to do
any kind of implementation work... or even much programming).
4. It also must work with all kinds of logicals, with .or. playing the
role of addition and .and. playing the role of multiplication.
P.S. Is it really not yet implemented in gfortran? That sort of
surprises me, but I haven't actually checked.
--
Richard Maine | Good judgment comes from experience;
Richard.Maine@nasa.gov | experience comes from bad judgment.
| -- Mark Twain