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

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1996/02/15 17:49:55  mclareni
* Kernlib
*
*
#include "kerngen/pilot.h"
      SUBROUTINE TRATS (A,S,B,M,N)
C
C CERN PROGLIB# F112    TRATS           .VERSION KERNFOR  4.15  861204
C ORIG. 18/12/74 WH
C
#if defined(CERNLIB_INTDOUBL)
      DOUBLE PRECISION SUM
#endif
      DIMENSION A(*),S(*),B(*)
C
      IB = 0
      INDS = 0
      I = 0
C
    5 INDS = INDS + I
      IB = I + 1
C
      DO 40 J=1,M
      IA = J
      IS = INDS
      SUM = 0.
      K = 0
C
   10 IF (K.GT.I)             GO TO 20
      IS = IS + 1
        GO  TO                30
   20 IS = IS + K
   30 SUM = SUM + A(IA)*S(IS)
      IA = IA + M
      K = K + 1
      IF (K.LT.N)             GO TO 10
C
      B(IB) = SUM
   40 IB = IB + N
      I = I + 1
      IF (I.LT.N)             GO TO 5
C
      RETURN
      END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.1 1996/02/15 17:49:55 mclareni
	6	* Kernlib
	7	*
	8	*
	9	#include "kerngen/pilot.h"
	10	SUBROUTINE TRATS (A,S,B,M,N)
	11	C
	12	C CERN PROGLIB# F112 TRATS .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(),S(),B(*)
	19	C
	20	IB = 0
	21	INDS = 0
	22	I = 0
	23	C
	24	5 INDS = INDS + I
	25	IB = I + 1
	26	C
	27	DO 40 J=1,M
	28	IA = J
	29	IS = INDS
	30	SUM = 0.
	31	K = 0
	32	C
	33	10 IF (K.GT.I) GO TO 20
	34	IS = IS + 1
	35	GO TO 30
	36	20 IS = IS + K
	37	30 SUM = SUM + A(IA)*S(IS)
	38	IA = IA + M
	39	K = K + 1
	40	IF (K.LT.N) GO TO 10
	41	C
	42	B(IB) = SUM
	43	40 IB = IB + N
	44	I = I + 1
	45	IF (I.LT.N) GO TO 5
	46	C
	47	RETURN
	48	END