bug report

Posted on : May 11, 2003 (Sun) 16:00:13

by h. akai

A bug is found in "getnum.f".
The corrected code (complete subroutine) is attached below.
The correction should be made on the data statement:

"scale/214748364d0/" -> scale/2147483647d0/

This bug causes no serious problems as far as it is
used in the KKR-CPA package.

getnum.f

subroutine getnum(ise,r,n)
c-----------------------------------------------------------------------
implicit real*8 (a-h,o-z)
dimension r(n)
data a/16807d0/, scale/2147483647d0/
se=ise
do 10 i=1,n
se=se*a
in=se/scale
se=se-scale*dble(in)
10 r(i)=se/scale
return
end