]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 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 TRSA (S,A,B,M,N) | |
11 | C | |
12 | C CERN PROGLIB# F112 TRSA .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 | INDS = 0 | |
21 | IB = 0 | |
22 | I = 0 | |
23 | C | |
24 | 5 INDS = INDS + I | |
25 | C | |
26 | DO 40 J=1,N | |
27 | IA = J | |
28 | IS = INDS | |
29 | SUM = 0. | |
30 | K = 0 | |
31 | C | |
32 | 10 IF (K.GT.I) GO TO 20 | |
33 | IS = IS + 1 | |
34 | GO TO 30 | |
35 | 20 IS = IS + K | |
36 | 30 SUM = SUM + S(IS)*A(IA) | |
37 | IA = IA + N | |
38 | K = K + 1 | |
39 | IF (K.LT.M) GO TO 10 | |
40 | IB = IB + 1 | |
41 | 40 B(IB) = SUM | |
42 | I = I + 1 | |
43 | IF (I.LT.M) GO TO 5 | |
44 | C | |
45 | RETURN | |
46 | END |