]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.1 1995/10/24 10:20:51 cernlib | |
6 | * Geant | |
7 | * | |
8 | * | |
9 | #include "geant321/pilot.h" | |
10 | *CMZ : 3.21/02 29/03/94 15.41.29 by S.Giani | |
11 | *-- Author : | |
12 | SUBROUTINE GINROT(X,RMAT,XNEW) | |
13 | C. | |
14 | C. ****************************************************************** | |
15 | C * * | |
16 | C * ROUTINE TO ROTATE THE VECTOR X BY THE INVERSE OF THE * | |
17 | C * ROTATION MATRIX RMAT AND STORE THE RESULT AS THE * | |
18 | C * VECTOR XNEW * | |
19 | C. * * | |
20 | C. * ==>Called by : GDFR3D,GDSSUR,GD3D3D,GFLCAR,GFLPHI,GFLRAD * | |
21 | C. * GFLTHE * | |
22 | C. * Author A.McPherson ********* * | |
23 | C. * * | |
24 | C. ****************************************************************** | |
25 | C. | |
26 | DIMENSION X(3),XNEW(3),RMAT(10) | |
27 | C. | |
28 | C. ------------------------------------------------------------------ | |
29 | C. | |
30 | IF(RMAT(10).NE.0.0) THEN | |
31 | XNEW(1) = X(1)*RMAT(1) + X(2)*RMAT(4) + X(3)*RMAT(7) | |
32 | XNEW(2) = X(1)*RMAT(2) + X(2)*RMAT(5) + X(3)*RMAT(8) | |
33 | XNEW(3) = X(1)*RMAT(3) + X(2)*RMAT(6) + X(3)*RMAT(9) | |
34 | * | |
35 | ELSE | |
36 | XNEW(1) = X(1) | |
37 | XNEW(2) = X(2) | |
38 | XNEW(3) = X(3) | |
39 | * | |
40 | ENDIF | |
41 | END |