The R Project SVN R

Rev

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

Rev 5458 Rev 24972
Line 92... Line 92...
92 
c
 92 
c
93 
      call dqrdc2(x,n,n,p,tol,k,qraux,jpvt,work)
 93 
      call dqrdc2(x,n,n,p,tol,k,qraux,jpvt,work)
94 
c
 94 
c
95 
c     solve the truncated least squares problem for each rhs.
 95 
c     solve the truncated least squares problem for each rhs.
96 
c
 96 
c
97 
      if(k .eq. 0) go to 30
 97 
      if(k .gt. 0) then
98 
      do 20 jj=1,ny
 98 
         do 20 jj=1,ny
99 
   20    call dqrsl(x,n,n,k,qraux,y(1,jj),rsd(1,jj),qty(1,jj),
 99 
   20       call dqrsl(x,n,n,k,qraux,y(1,jj),rsd(1,jj),qty(1,jj),
100 
     1      b(1,jj),rsd(1,jj),rsd(1,jj),1110,info)
 100 
     1           b(1,jj),rsd(1,jj),rsd(1,jj),1110,info)
- 
 
 101 
      else
- 
 
 102 
         do 30 i=1,n
- 
 
 103 
            do 30 jj=1,ny
- 
 
 104 
   30           rsd(i,jj) = y(i,jj)
101 
   30 continue
 105 
      endif
102 
c
 106 
c
103 
c     set the unused components of b to zero.
 107 
c     set the unused components of b to zero.
104 
c
 108 
c
105 
      kk = k + 1
 109 
      kk = k + 1
106 
      do 50 j=kk,p
 110 
      do 50 j=kk,p