The R Project SVN R

Rev

Rev 12778 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 12778 Rev 42390
Line 68... Line 68...
68 
c     begin block with ...exits to 40
 68 
c     begin block with ...exits to 40
69 
c
 69 
c
70 
c
 70 
c
71 
         do 30 j = 1, n
 71 
         do 30 j = 1, n
72 
            info = j
 72 
            info = j
73 
            s = 0.0e0
 73 
            s = 0.0d0
74 
            ik = m + 1
 74 
            ik = m + 1
75 
            jk = max0(j-m,1)
 75 
            jk = max0(j-m,1)
76 
            mu = max0(m+2-j,1)
 76 
            mu = max0(m+2-j,1)
77 
            if (m .lt. mu) go to 20
 77 
            if (m .lt. mu) go to 20
78 
            do 10 k = mu, m
 78 
            do 10 k = mu, m
Line 86... Line 86...
86 
   10       continue
 86 
   10       continue
87 
   20       continue
 87 
   20       continue
88 
 
 88 
 
89 
            s = abd(m+1,j) - s
 89 
            s = abd(m+1,j) - s
90 
c     ......exit
 90 
c     ......exit
91 
            if (s .le. 0.0e0) go to 40
 91 
            if (s .le. 0.0d0) go to 40
92 
 
 92 
 
93 
            abd(m+1,j) = sqrt(s)
 93 
            abd(m+1,j) = sqrt(s)
94 
 
 94 
 
95 
   30    continue
 95 
   30    continue
96 
         info = 0
 96 
         info = 0