]>
Commit | Line | Data |
---|---|---|
1 | * | |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.1 1995/10/24 10:21:07 cernlib | |
6 | * Geant | |
7 | * | |
8 | * | |
9 | #include "geant321/pilot.h" | |
10 | *CMZ : 3.21/02 29/03/94 15.41.38 by S.Giani | |
11 | *-- Author : | |
12 | C | |
13 | C *** VARIOUS VECTOR OPERATIONS *** | |
14 | C | |
15 | C COPIED FROM F14BLO.PAMLIB 23.4.82 | |
16 | C AUTHOR: V.BLOBEL (UNIVERSITY OF HAMBURG) | |
17 | C DESYLIB | |
18 | C | |
19 | C *** BLANK COMMON REPLACED BY /VECUTY/ TO MATCH GEANT/GHEISHA CODE *** | |
20 | C *** NOTE THAT P(10,100) HAS BECOME PV(10,200) DUE TO THIS *** | |
21 | C | |
22 | C UN-USED ENTRIES REMOVED : | |
23 | C "PCOP" "PEXC" "PZER" "PWRT" "DOT4" "IMPU" "IMPULI" "ADD3" | |
24 | C "SUB3" "CROSS" "DOT" "SMUL" "NORZ" "PARPER" "PUNIT" "TRAP" | |
25 | C | |
26 | C *** ALL ENTRIES RE-WRITTEN AS SUBROUTINES USING ONLY NECESSARY *** | |
27 | C *** "DOUBLE PRECISION" STMTS. AND ALL SPECIFIC FUNCTIONS HAVE *** | |
28 | C *** BEEN CHANGED TO THEIR GENERIC EQUIVALENCES *** | |
29 | C *** NVE 29-MAR-1988 CERN GENEVA *** | |
30 | C | |
31 | C ORIGIN : H.FESEFELDT (22-JUNE-1984) | |
32 | C | |
33 | SUBROUTINE LOR(K,L,M) | |
34 | C | |
35 | C --- PV-ARRAY --- | |
36 | #include "geant321/mxgkgh.inc" | |
37 | #include "geant321/s_blankp.inc" | |
38 | C | |
39 | DOUBLE PRECISION A,B,C | |
40 | C | |
41 | A=0.0 | |
42 | DO 6 I=1,3 | |
43 | A=A+PV(I,K)*PV(I,L) | |
44 | 6 CONTINUE | |
45 | A=(A/(PV(4,L)+PV(5,L))-PV(4,K))/PV(5,L) | |
46 | B=PV(5,K)*PV(5,K) | |
47 | DO 8 I=1,3 | |
48 | C=PV(I,K)+A*PV(I,L) | |
49 | B=B+C*C | |
50 | PV(I,M)=C | |
51 | 8 CONTINUE | |
52 | PV(4,M)=SQRT(B) | |
53 | PV(5,M)=PV(5,K) | |
54 | RETURN | |
55 | END |