public inbox for gcc-bugs@sourceware.org
help / color / mirror / Atom feed
* [Bug libfortran/21926] New: matmul does not deal with non-packed result
@ 2005-06-05 20:53 tkoenig at gcc dot gnu dot org
  2005-06-06 19:14 ` [Bug libfortran/21926] " tkoenig at gcc dot gnu dot org
                   ` (4 more replies)
  0 siblings, 5 replies; 6+ messages in thread
From: tkoenig at gcc dot gnu dot org @ 2005-06-05 20:53 UTC (permalink / raw)
  To: gcc-bugs

$ cat matmul.f90
program main
  real, dimension(2,2) :: a
  real, dimension(4,2) :: c
  a = reshape((/ 1.0, 1.0, 0.0, 1.0/), shape(a))
  c = 42.
  c(1:2,1:2) = matmul(a, transpose(a))
  print *,c(1:2,1:2)
  print *,matmul(a, transpose(a))
  print *,c(3:4,1:2)
end
$ gfortran matmul.f90
$ ./a.out
   1.000000       1.000000       43.00000       44.00000
   1.000000       1.000000       1.000000       2.000000
   0.000000       0.000000       42.00000       42.00000
$ gfortran -v
Using built-in specs.
Target: i686-pc-linux-gnu
Configured with: ../gcc-4.1/configure --prefix=/home/ig25 --enable-languages=c,f95
Thread model: posix
gcc version 4.1.0 20050603 (experimental)

The problem is in the zeroing of the result, which assumes
a packed return array.

-- 
           Summary: matmul does not deal with non-packed result
           Product: gcc
           Version: 4.1.0
            Status: UNCONFIRMED
          Keywords: wrong-code
          Severity: normal
          Priority: P2
         Component: libfortran
        AssignedTo: unassigned at gcc dot gnu dot org
        ReportedBy: tkoenig at gcc dot gnu dot org
                CC: gcc-bugs at gcc dot gnu dot org


http://gcc.gnu.org/bugzilla/show_bug.cgi?id=21926


^ permalink raw reply	[flat|nested] 6+ messages in thread
end of thread, other threads:[~2005-07-15 20:48 UTC | newest]

Thread overview: 6+ messages (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
2005-06-05 20:53 [Bug libfortran/21926] New: matmul does not deal with non-packed result tkoenig at gcc dot gnu dot org
2005-06-06 19:14 ` [Bug libfortran/21926] " tkoenig at gcc dot gnu dot org
2005-06-07 20:46 ` cvs-commit at gcc dot gnu dot org
2005-06-07 20:53 ` [Bug libfortran/21926] [4.0 only] " tkoenig at gcc dot gnu dot org
2005-07-15 20:47 ` cvs-commit at gcc dot gnu dot org
2005-07-15 20:49 ` tkoenig at gcc dot gnu dot org

This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).