*
* $Id: tlstep.F,v 1.1.1.1 1996/02/15 17:49:53 mclareni Exp $
*
* $Log: tlstep.F,v $
* Revision 1.1.1.1  1996/02/15 17:49:53  mclareni
* Kernlib
*
*
#include "kerngen/pilot.h"
      SUBROUTINE TLSTEP (A,B,IASEP,IBSEP,NR,NC,BETA)
C
C CERN PROGLIB# E230    TLSTEP          .VERSION KERNFOR  2.06  740511
C ORIG. 11/05/74 WH+WM
C
C.  SUBROUTINE TLSTEP (A,B,IASEP,IBSEP,NR,NC,BETA)
C.
C.       PERFORMS AN ORTHOGONAL TRANSFORMATION.
C.
C.-------------------------------------------------------------------
C
      COMMON /SLATE/ DUM(34),H,I,IB,J,JA,JB
      DIMENSION      A(*), B(*)
C
      IB = 0
      DO           30        J=1,NC
      IB = IB + 1
      H  = 0.
      JA = 1
      JB = IB
C
      DO           10        I=1,NR
      H  = H + A(JA) * B(JB)
      JA = JA +IASEP
   10 JB = JB + IBSEP
C
      H  = H * BETA
      JA = 1
      JB = IB
      DO           20        I=1,NR
      B(JB) = B(JB) - A(JA) * H
      JA = JA +IASEP
   20 JB = JB + IBSEP
   30 CONTINUE
C
      RETURN
      END