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

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1996/02/15 17:49:53  mclareni
* Kernlib
*
*
#include "kerngen/pilot.h"
      SUBROUTINE TLPIV(A,B,IASEP,IBSEP,M,H,G)
C
C CERN PROGLIB# E230    TLPIV           .VERSION KERNFOR  2.06  740511
C ORIG. 11/05/74 WH+WM
C
C.  SUBROUTINE TLPIV (A,B,IASEP,IBSEP,M,H,G)
C.
C.        CALCULATES PIVOTING FACTORS FOR TLS
C.
C.-------------------------------------------------------------------
C
      COMMON /SLATE/ DUM(37),I,JA,JB
      DIMENSION A(*),B(*)
C
      H=0.0
      G=0.0
      JA=1
      JB=1
      DO 2 I=1,M
      G=G+A(JA)*B(JB)
      H=H+A(JA)*A(JA)
      JA=JA+IASEP
    2 JB=JB+IBSEP
      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 TLPIV(A,B,IASEP,IBSEP,M,H,G)
	11	C
	12	C CERN PROGLIB# E230 TLPIV .VERSION KERNFOR 2.06 740511
	13	C ORIG. 11/05/74 WH+WM
	14	C
	15	C. SUBROUTINE TLPIV (A,B,IASEP,IBSEP,M,H,G)
	16	C.
	17	C. CALCULATES PIVOTING FACTORS FOR TLS
	18	C.
	19	C.-------------------------------------------------------------------
	20	C
	21	COMMON /SLATE/ DUM(37),I,JA,JB
	22	DIMENSION A(),B()
	23	C
	24	H=0.0
	25	G=0.0
	26	JA=1
	27	JB=1
	28	DO 2 I=1,M
	29	G=G+A(JA)*B(JB)
	30	H=H+A(JA)*A(JA)
	31	JA=JA+IASEP
	32	2 JB=JB+IBSEP
	33	RETURN
	34	END