------- Comment #6 from hjl at lucon dot org 2006-11-22 16:56 ------- A simpler testcase:
[EMAIL PROTECTED] 29921]$ cat bar.f90 SUBROUTINE foo(DPRS,FPK,LBOT,LTOP,KTS,KTE) IMPLICIT NONE INTEGER,INTENT(IN) :: KTS,KTE INTEGER,INTENT(OUT) :: LBOT,LTOP REAL,DIMENSION(KTS:KTE),INTENT(IN) :: DPRS,FPK REAL,DIMENSION(KTS:KTE) :: TREFK REAL :: PSUM, QOTSUM, SUMDP, TRFKL,DPKL,DST REAL :: FPTK, RDPSUM, RTBAR, DEN, DQREF INTEGER :: L SUMDP=0. DO L=LTOP,LBOT SUMDP=SUMDP+DPRS(L) ENDDO RDPSUM=1./SUMDP DO L=LTOP,LBOT TREFK(L)=TRFKL ENDDO PSUM =0. QOTSUM=0. DST =0. FPTK =FPK(LTOP) DO L=LTOP,LBOT DPKL =FPK(L)-FPTK PSUM =DPKL *DPRS(L)+PSUM ENDDO PSUM =PSUM*RDPSUM IF(-DEN/PSUM.LT.5.E-5)THEN LBOT=0 ENDIF END SUBROUTINE foo [EMAIL PROTECTED] 29921]$ -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=29921