MINICERN/packlib/kernlib/kerngen/tcgen/trsa.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 TRSA (S,A,B,M,N)
  11 C
  12 C CERN PROGLIB# F112    TRSA            .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 S(*),A(*),B(*)
  19 C
  20       INDS = 0
  21       IB = 0
  22       I = 0
  23 C
  24     5 INDS = INDS + I
  25 C
  26       DO 40 J=1,N
  27       IA = J
  28       IS = INDS
  29       SUM = 0.
  30       K = 0
  31 C
  32    10 IF (K.GT.I)             GO TO 20
  33       IS = IS + 1
  34         GO  TO                30
  35    20 IS = IS + K
  36    30 SUM = SUM + S(IS)*A(IA)
  37       IA = IA + N
  38       K = K + 1
  39       IF (K.LT.M)             GO TO 10
  40       IB = IB + 1
  41    40 B(IB) = SUM
  42       I = I + 1
  43       IF (I.LT.M)             GO TO 5
  44 C
  45       RETURN
  46       END