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

   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