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

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1996/02/15 17:49:54  mclareni
* Kernlib
*
*
#include "kerngen/pilot.h"
      SUBROUTINE TRLTA (U,A,B,M,N)
C
C CERN PROGLIB# F112    TRLTA           .VERSION KERNFOR  4.15  861204
C ORIG. 18/12/74 WH
C
#if defined(CERNLIB_INTDOUBL)
      DOUBLE PRECISION SUM
#endif
      DIMENSION A(*),U(*),B(*)
C
      IPIV = 0
      MX = M*N
      MXPN = MX + N
      IB = 0
C
      I = 0
   10 I = I + 1
      IPIV = IPIV + I
C
   20 IU = IPIV
      NSTEP = I
      IB = IB + 1
      IA = IB
C
      SUM = 0.
   30 SUM = SUM + A(IA)*U(IU)
      IA = IA + N
      IU = IU + NSTEP
      NSTEP = NSTEP + 1
      IF (IA.LE.MX)           GO TO 30
C
      B(IB) = SUM
      IF (IA.LT.MXPN)         GO TO 20
C
      IF (I.LT.M)             GO TO 10
C
      RETURN
      END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.1 1996/02/15 17:49:54 mclareni
	6	* Kernlib
	7	*
	8	*
	9	#include "kerngen/pilot.h"
	10	SUBROUTINE TRLTA (U,A,B,M,N)
	11	C
	12	C CERN PROGLIB# F112 TRLTA .VERSION KERNFOR 4.15 861204
	13	C ORIG. 18/12/74 WH
	14	C
	15	#if defined(CERNLIB_INTDOUBL)
	16	DOUBLE PRECISION SUM
	17	#endif
	18	DIMENSION A(),U(),B(*)
	19	C
	20	IPIV = 0
	21	MX = M*N
	22	MXPN = MX + N
	23	IB = 0
	24	C
	25	I = 0
	26	10 I = I + 1
	27	IPIV = IPIV + I
	28	C
	29	20 IU = IPIV
	30	NSTEP = I
	31	IB = IB + 1
	32	IA = IB
	33	C
	34	SUM = 0.
	35	30 SUM = SUM + A(IA)*U(IU)
	36	IA = IA + N
	37	IU = IU + NSTEP
	38	NSTEP = NSTEP + 1
	39	IF (IA.LE.MX) GO TO 30
	40	C
	41	B(IB) = SUM
	42	IF (IA.LT.MXPN) GO TO 20
	43	C
	44	IF (I.LT.M) GO TO 10
	45	C
	46	RETURN
	47	END