MINICERN/packlib/kernlib/kerngen/tcgen/trlta.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 TRLTA (U,A,B,M,N)
  11 C
  12 C CERN PROGLIB# F112    TRLTA           .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       IPIV = 0
  21       MX = M*N
  22       MXPN = MX + N
  23       IB = 0
  24 C
  25       I = 0
  26    10 I = I + 1
  27       IPIV = IPIV + I
  28 C
  29    20 IU = IPIV
  30       NSTEP = I
  31       IB = IB + 1
  32       IA = IB
  33 C
  34       SUM = 0.
  35    30 SUM = SUM + A(IA)*U(IU)
  36       IA = IA + N
  37       IU = IU + NSTEP
  38       NSTEP = NSTEP + 1
  39       IF (IA.LE.MX)           GO TO 30
  40 C
  41       B(IB) = SUM
  42       IF (IA.LT.MXPN)         GO TO 20
  43 C
  44       IF (I.LT.M)             GO TO 10
  45 C
  46       RETURN
  47       END