[Bug tree-optimization/31280] segfault in remove_referenced_var
pinskia at gcc dot gnu dot org
gcc-bugzilla@gcc.gnu.org
Thu Mar 29 00:41:00 GMT 2007
------- Comment #5 from pinskia at gcc dot gnu dot org 2007-03-29 01:41 -------
Reduced self contained testcase:
Module NsLin
Real*8, Dimension(54,54) :: VXX, VXY,
& VYX, VYY, VXX2, VXY2, VYX2, VYY2
End Module NsLin
Subroutine SysNSL (N, NX, NY)
Use NsLin
Implicit Real*8 (A-H,O-Z)
Real*8, Dimension(1000,1000) :: EnsLap, EnsArh
Real*8, Dimension(1000,1000) :: EnsTim, HtTim
Integer, Allocatable, Dimension(:) :: IPIV
Allocate ( IPIV(N), stat=jerr )
Do L=1,NX
Do M=1,NY
Do i=1,NX
Do j=1,NY
EnsTim(NY+ M, NY+ j ) =
& ( VXX(i,L)* VXY(j,M) + VYX(i,L)* VYY(j,M) )
& * W
EnsLap(NY*(L-1) + M,NY*(i-1) + j ) =
& (VXX(i,L) * VXY2(j,M) + VYX(i,L) * VYY2(j,M) )
& * W
& + (VXX2(i,L) * VXY(j,M) + VYX2(i,L) * VYY(j,M) )
& / W
End Do
End Do
End Do
End Do
Call Dgetrf (EnsTim(1:N,1:N), info)
If (info .NE. 0) then
Write (*,*) 'sysnsL(Dgetrf): info=',info
Stop
End If
Call Dgetri ( EnsTim(1:N,1:N), IPIV, info)
If (info .NE. 0) then
Stop
End If
EnsLap(1:N,1:N) = MATMUL( EnsTim(1:N,1:N), EnsLap(1:N,1:N) )
Deallocate ( IPIV, stat=jerr)
Return
End
--
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=31280
More information about the Gcc-bugs
mailing list