URL
https://opencores.org/ocsvn/openrisc/openrisc/trunk
Subversion Repositories openrisc
[/] [openrisc/] [tags/] [gnu-src/] [gcc-4.5.1/] [gcc-4.5.1-or32-1.0rc3/] [gcc/] [testsuite/] [gfortran.dg/] [g77/] [20000511-2.f] - Rev 516
Compare with Previous | Blame | View Log
c { dg-do compile } subroutine sgbcon(norm,n,kl,ku,ab,ldab,ipiv,anorm,rcond,work,iwork &,info) C C -- LAPACK routine (version 3.0) -- C Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd., C Courant Institute, Argonne National Lab, and Rice University C September 30, 1994 C C .. Scalar Arguments .. character norm integer info,kl,ku,ldab,n real anorm,rcond C .. C .. Array Arguments .. integer ipiv(n),iwork(n) real ab(ldab,n),work(n) C .. C C Purpose C ======= C demonstrate g77 bug at -O -funroll-loops C ===================================================================== C C .. Parameters .. real one,zero parameter(one= 1.0e+0,zero= 0.0e+0) C .. C .. Local Scalars .. logical lnoti,onenrm character normin integer ix,j,jp,kase,kase1,kd,lm real ainvnm,scale,smlnum,t C .. C .. External Functions .. logical lsame integer isamax real sdot,slamch externallsame,isamax,sdot,slamch C .. C .. External Subroutines .. externalsaxpy,slacon,slatbs,srscl,xerbla C .. C .. Executable Statements .. C C Multiply by inv(L). C do j= 1,n-1 C the following min() intrinsic provokes this bug lm= min(kl,n-j) jp= ipiv(j) t= work(jp) if(jp.ne.j)then C but only when combined with this if block work(jp)= work(j) work(j)= t endif C and this subroutine call call saxpy(lm,-t,ab(kd+1,j),1,work(j+1),1) enddo return end