MINICERN/packlib/kernlib/kerngen/tcgen/trsat.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 TRSAT (S,A,B,M,N)
  11 C
  12 C CERN PROGLIB# F112    TRSAT           .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       MN = M*N
  21       INDS = 0
  22       IB = 0
  23       I = 0
  24 C
  25     5 INDS = INDS + I
  26       IA = 0
  27 C
  28       DO 40 J=1,N
  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 IA = IA + 1
  38       SUM = SUM + S(IS)*A(IA)
  39       K = K + 1
  40       IF (K.LT.M)             GO TO 10
  41       IB = IB + 1
  42    40 B(IB) = SUM
  43       I = I + 1
  44       IF (I.LT.M)             GO TO 5
  45 C
  46       RETURN
  47       END