MINICERN/packlib/kernlib/kerngen/tcgen/trata.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 TRATA  (A,R,M,N)
  11 C
  12 C CERN PROGLIB# F112    TRATA           .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(*),R(*)
  19 C
  20       MN = M*N
  21       IR = 0
  22 C
  23       DO 30 I=1,M
  24 C
  25       DO 30 J=1,I
  26       IA = I
  27       IAT = J
  28 C
  29       SUM = 0.
  30    20 SUM = SUM + A(IA)*A(IAT)
  31       IA = IA + M
  32       IAT = IAT + M
  33       IF (IA.LE.MN)           GO TO 20
  34 C
  35       IR = IR + 1
  36    30 R(IR) = SUM
  37 C
  38       RETURN
  39       END