The R Project SVN R

Rev

Rev 86632 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 86632 Rev 86633
Line 180... Line 180...
180 
               do 160 j = lp1, p
 180 
               do 160 j = lp1, p
181 
                  t = -ddot(n-l+1,x(l,l),1,x(l,j),1)/x(l,l)
 181 
                  t = -ddot(n-l+1,x(l,l),1,x(l,j),1)/x(l,l)
182 
                  call daxpy(n-l+1,t,x(l,l),1,x(l,j),1)
 182 
                  call daxpy(n-l+1,t,x(l,l),1,x(l,j),1)
183 
                  if (j .lt. pl .or. j .gt. pu) go to 150
 183 
                  if (j .lt. pl .or. j .gt. pu) go to 150
184 
                  if (qraux(j) .eq. 0.0d0) go to 150
 184 
                  if (qraux(j) .eq. 0.0d0) go to 150
185 
                     tt = 1.0d0 - (aabs(x(l,j))/qraux(j))**2
 185 
                     tt = 1.0d0 - (abs(x(l,j))/qraux(j))**2
186 
                     tt = max(tt,0.0d0)
 186 
                     tt = max(tt,0.0d0)
187 
                     t = tt
 187 
                     t = tt
188 
                     tt = 1.0d0 + 0.05d0*tt*(qraux(j)/work(j))**2
 188 
                     tt = 1.0d0 + 0.05d0*tt*(qraux(j)/work(j))**2
189 
                     if (tt .eq. 1.0d0) go to 130
 189 
                     if (tt .eq. 1.0d0) go to 130
190 
                        qraux(j) = qraux(j)*sqrt(t)
 190 
                        qraux(j) = qraux(j)*sqrt(t)