This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / d / s64.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1996/04/01 15:02:23  mclareni
6 * Mathlib gen
7 *
8 *
9 #include "gen/pilot.h"
10       SUBROUTINE S64(A,B)
11       REAL A(32),B(32),C(32)
12       COMMON /TAB/D(126),E(127)
13       CALL D701BD
14       DO 10 I=3,15,2
15       A1 = A(I-1)-A(I+1)
16       A2 = A(33-I)-A(35-I)
17       C1R = A1+A2
18       C2R = A1-A2
19       C1I = A(I)-A(34-I)
20       C2I = A(I)+A(34-I)
21       C1 = D(8*I-9)
22       C2 = D(8*I-8)
23       C3R = C1*C2R-C2*C2I
24       C3I = C1*C2I+C2*C2R
25       C(I) = C1R+C3R
26       C(I+1) = C1I+C3I
27       C(34-I) = C1R-C3R
28    10 C(35-I) = C3I-C1I
29       C(1) = 2.*(A(32)-A(2))
30       C(2) = -2.*(A(32)+A(2))
31       C(17) = 2.*(A(16)-A(18))
32       C(18) = 2.*A(17)
33       CALL I16(C,C(2),2,C,C(2),2)
34       B(17) = 0.5*C(17)
35       DO 20 I=2,16
36       E1 = 0.5*(C(I)-C(34-I))
37       E2 = E(8*I-8)*(C(I)+C(34-I))
38       B(I) = E2-E1
39    20 B(34-I) = E2+E1
40       RETURN
41       END