That annoying dnrm2.f failure on i386
Dave Love
d.love@dl.ac.uk
Wed Feb 18 14:54:00 GMT 1998
>>>>> "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
More information about the Gcc
mailing list