From mboxrd@z Thu Jan  1 00:00:00 1970
From: Dave Love <d.love@dl.ac.uk>
To: egcs@cygnus.com
Subject: Re: That annoying dnrm2.f failure on i386
Date: Wed, 18 Feb 1998 14:54:00 -0000
Message-id: <rzqbtw4vfap.fsf@djlvig.dl.ac.uk>
References: <Pine.SOL.3.90.980218113529.12315K-100000@bond.informatik.rwth-aachen.de>
X-SW-Source: 1998-02/msg00894.html

>>>>> "Bernd" == Bernd Schmidt <crux@pool.informatik.rwth-aachen.de> writes:

 Bernd> Unfortunately, dnrm2.f still fails. This appears to be due to
 Bernd> the fact that it reads uninitialized variables. 

Indeed.  It's certainly wrong.  I've checked in a fixed version as I
think I sent it in the first place.  Apologies if the lossage was
mine.

I haven't yet been able to test the fixed version with your patch.

 Bernd> Can someone who can actually parse Fortran tell me how the
 Bernd> test case should be fixed?

Actually, I can't parse that stuff, and this infamous bit of BLAS was
replaced some time ago, thank heavens, but just won't die...

Many thanks for the work!

--- dnrm2.f     1997/11/19 16:20:12     1.2
+++ dnrm2.f     1998/02/18 12:00:51     1.3
@@ -15,7 +15,7 @@ C     (This is actually an obsolete vers
 c     current Netlib BLAS.)

       integer i
-      double precision a(1:100)
+      double precision a(1:100), dnrm2
       do i=1,100
          a(i)=0.D0
       enddo
@@ -26,6 +26,8 @@ c     current Netlib BLAS.)
       integer i, incx, ix, j, n, next
       double precision   dx(1), cutlo, cuthi, hitest, sum, xmax,zero,one
       data   zero, one /0.0d0, 1.0d0/
+      data cutlo, cuthi / 8.232d-11,  1.304d19 /
+      j = 0
       if(n .gt. 0 .and. incx.gt.0) go to 10
          dnrm2  = zero