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

(Source file text) 

---

      SUBROUTINE HHSTEP(A,B,IASEP,IBSEP,NR,NC,BETA) 
C 
C        PURPOSE
C           TO PERFORM THE ORTHOGONAL TRANSFORMATION
C 
      DIMENSION A(2), B(2)
C 
      IB = 1 - IBSEP
      DO 21 J = 1,NC
      IB = IB + IBSEP 
      H=0.
      JA = 1
      JB = IB 
      DO 20 I=1,NR
      H = H + A(JA) * B(JB) 
      JA = JA +IASEP
   20 JB = JB + 1 
      H=BETA*H
      JA = 1
      JB = IB 
      DO 21 I=1,NR
      B(JB) = B(JB) - A(JA) * H 
      JA = JA +IASEP
   21 JB = JB + 1 
      RETURN
      END 


---