File: HLSR.FT of Disk: Disks/MyPDP/m8-2-rka1-rkb1

(Source file text) 

---

      SUBROUTINE HLSR(A,B,M,N,L,IER,AUX)
C 
C        PURPOSE
C           TO CALCULATE THE FITTED RESIDUALS AFTER A CALL TO HLS 
C 
C        USAGE
C           CALL HLSR (A,B,M,N,L,IER,AUX) 
C 
C        DESCRIPTION OF PARAMETERS
C           AS BEFORE WITH MATRIX B OUTPUTING THE FITTED RESIDUALS. 
C 
C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
C           HHSTEP
C 
C        METHOD 
C           APPLY INVERSE HH TRANSFORMATION TO RESIDUAL VECTOR
C 
C     ..................................................................
C 
      DIMENSION A(2), B(2), AUX(2)
C 
      IF(IER.EQ.0) GO TO 100
      IER = IABS(IER) 
C 
C     SET FIRST IER ELEMENTS OF EACH B VECTOR TO ZERO 
      IST = - M 
      DO 10 J=1,L 
      IST = IST + M 
      I = IST 
      DO 10 K=1,IER 
      I = I + 1 
   10 B(I) = 0. 
C 
C      TRANSFORMATION LOOP
      K1 = MAX0(N,L)
      IST = (IER-1)*(M+1) + 1 
      DO 30 KK=1,IER
      LV = M - IER + KK 
      KN = K1 + IER - KK + 1
      BETA = -1./(AUX(KN)*A(IST)) 
      IB = IER - KK + 1 
      CALL HHSTEP(A(IST),B(IB),1,M,LV,L,BETA) 
   30 IST = IST - M - 1 
  100 RETURN
      END 


---