MINICERN/packlib/kernlib/kerngen/tcgen/tlstep.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/02/15 17:49:53  mclareni
   6 * Kernlib
   7 *
   8 *
   9 #include "kerngen/pilot.h"
  10       SUBROUTINE TLSTEP (A,B,IASEP,IBSEP,NR,NC,BETA)
  11 C
  12 C CERN PROGLIB# E230    TLSTEP          .VERSION KERNFOR  2.06  740511
  13 C ORIG. 11/05/74 WH+WM
  14 C
  15 C.  SUBROUTINE TLSTEP (A,B,IASEP,IBSEP,NR,NC,BETA)
  16 C.
  17 C.       PERFORMS AN ORTHOGONAL TRANSFORMATION.
  18 C.
  19 C.-------------------------------------------------------------------
  20 C
  21       COMMON /SLATE/ DUM(34),H,I,IB,J,JA,JB
  22       DIMENSION      A(*), B(*)
  23 C
  24       IB = 0
  25       DO           30        J=1,NC
  26       IB = IB + 1
  27       H  = 0.
  28       JA = 1
  29       JB = IB
  30 C
  31       DO           10        I=1,NR
  32       H  = H + A(JA) * B(JB)
  33       JA = JA +IASEP
  34    10 JB = JB + IBSEP
  35 C
  36       H  = H * BETA
  37       JA = 1
  38       JB = IB
  39       DO           20        I=1,NR
  40       B(JB) = B(JB) - A(JA) * H
  41       JA = JA +IASEP
  42    20 JB = JB + IBSEP
  43    30 CONTINUE
  44 C
  45       RETURN
  46       END