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

   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 TRAL (A,U,B,M,N)
  11 C
  12 C CERN PROGLIB# F112    TRAL            .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       IB = 1
  21 C
  22       DO 40 I=1,M
  23       INDU = 0
  24 C
  25       DO 40 J=1,N
  26       INDU = INDU + J
  27       IA = IB
  28       IU = INDU
  29       SUM = 0.
  30 C
  31       DO 30 K=J,N
  32       SUM = SUM + A(IA)*U(IU)
  33       IA = IA + 1
  34    30 IU = IU + K
  35       B(IB) = SUM
  36    40 IB = IB + 1
  37 C
  38       RETURN
  39       END