]>
Commit | Line | Data |
---|---|---|
3820ca8e | 1 | |
2 | CDECK ID>, HWULOR. | |
3 | ||
4 | *CMZ :- -26/04/91 11.11.56 by Bryan Webber | |
5 | ||
6 | *-- Author : Giovanni Abbiendi & Luca Stanco | |
7 | ||
8 | C----------------------------------------------------------------------- | |
9 | ||
10 | SUBROUTINE HWULOR (TRANSF,PI,PF) | |
11 | ||
12 | C----------------------------------------------------------------------- | |
13 | ||
14 | C Makes the HWULOR transformation specified by TRANSF on the | |
15 | ||
16 | C quadrivector PI(5), giving PF(5). | |
17 | ||
18 | C----------------------------------------------------------------------- | |
19 | ||
20 | DOUBLE PRECISION TRANSF(4,4),PI(5),PF(5) | |
21 | ||
22 | INTEGER I,J | |
23 | ||
24 | DO 1 I=1,5 | |
25 | ||
26 | PF(I)=0.D0 | |
27 | ||
28 | 1 CONTINUE | |
29 | ||
30 | DO 3 I=1,4 | |
31 | ||
32 | DO 2 J=1,4 | |
33 | ||
34 | PF(I) = PF(I) + TRANSF(I,J) * PI(J) | |
35 | ||
36 | 2 CONTINUE | |
37 | ||
38 | 3 CONTINUE | |
39 | ||
40 | PF(5) = PI(5) | |
41 | ||
42 | RETURN | |
43 | ||
44 | END |