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