Re: That annoying dnrm2.f failure on i386

[Dave Love <d.love@dl.ac.uk>]

>>>>> "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