This is the mail archive of the
fortran@gcc.gnu.org
mailing list for the GNU Fortran project.
Re: matmul intrinsic procedure
On Thu, Jun 30, 2005 at 03:13:13PM -0600, Jan Stanis?aw Owoc wrote:
> I have a test case and the CVS version of GCC 4.1 from June 27th:
>
> PROGRAM TEST
> IMPLICIT NONE
> DOUBLE PRECISION, ALLOCATABLE :: A(:,:),X(:),B(:)
> INTEGER, PARAMETER :: n = 3
> ALLOCATE ( A(n,n), X(n), B(n))
>
> A = RESHAPE( source = (/ 2.0, 5.0, 3.0, 6.0, 3.0, 7.0, 6.0, 3.0, 2.0 /),
> shape = (/3, 3/) )
> X = (/ 2.7, 5.4, 8.2 /)
>
> B = MATMUL( A, X )
>
> PRINT *, B
>
> END PROGRAM TEST
troutmask:kargl[223] ./z
86.9999995231628 54.2999999523163 62.3000004291534
> Where/how do I file this as a "bug"?
If gfortran is installed with your OS, then make sure that
the 4.1 version of gfortran that you installed is picking
up the right shared library.
>
> BTW, the list of intrinsic procedures in the docs does not include
> "matmul":
> http://gcc.gnu.org/onlinedocs/gcc-4.0.0/gfortran/Intrinsic-Procedures.html
>
Carefully read the first paragraph of that section.
"This portion of the document is incomplete and undergoing massive
expansion and editing. All contributions and corrections are strongly
encouraged."
--
Steve