git://git.uio.no / u / mrichter / AliRoot.git / blame
| Commit | Line | Data |
|---|---|---|
| fe4da5cc | 1 | * |
| 2 | * $Id$ | |
| 3 | * | |
| 4 | * $Log$ | |
| 5 | * Revision 1.1.1.1 1996/04/01 15:02:29 mclareni | |
| 6 | * Mathlib gen | |
| 7 | * | |
| 8 | * | |
| 9 | #include "gen/pilot.h" | |
| 10 | FUNCTION RTRGSM(X,A,N,B,M,IOP) | |
| 11 | ||
| 12 | DIMENSION A(0:*),B(*) | |
| 13 | ||
| 14 | PARAMETER (R1 = 1, HF = R1/2) | |
| 15 | ||
| 16 | C=2*COS(X) | |
| 17 | H=0 | |
| 18 | IF(IOP .NE. 3) THEN | |
| 19 | A2=0 | |
| 20 | A1=0 | |
| 21 | A0=A(N) | |
| 22 | DO 1 K = N-1,0,-1 | |
| 23 | A2=A1 | |
| 24 | A1=A0 | |
| 25 | A0=A(K)+C*A1-A2 | |
| 26 | 1 CONTINUE | |
| 27 | H=HF*(A(0)+A0-A2) | |
| 28 | ENDIF | |
| 29 | IF(IOP .NE. 2) THEN | |
| 30 | A2=0 | |
| 31 | A1=0 | |
| 32 | A0=B(M) | |
| 33 | DO 2 K = M-1,1,-1 | |
| 34 | A2=A1 | |
| 35 | A1=A0 | |
| 36 | A0=B(K)+C*A1-A2 | |
| 37 | 2 CONTINUE | |
| 38 | H=H+A0*SIN(X) | |
| 39 | ENDIF | |
| 40 | RTRGSM=H | |
| 41 | RETURN | |
| 42 | END |