[u/mrichter/AliRoot.git] / MINICERN / packlib / kernlib / kerngen / tcgen / tlstep.F

*
* $Id$
*
* $Log$
* 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
Commit	Line	Data
fe4da5cc	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