This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / d / c256.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 C256(A,B)
11       REAL A(129),B(129),C(128)
12       COMMON /TAB/D(126),E(127)
13       CALL D701BD
14       S2 = A(64)+A(66)
15       DO 10 I=3,63,2
16       S2 = S2+A(I-1)+A(131-I)
17       A1 = A(I+1)-A(I-1)
18       A2 = A(131-I)-A(129-I)
19       C1R = A(I)+A(130-I)
20       C2R = A(I)-A(130-I)
21       C1I = A1-A2
22       C2I = A1+A2
23       C1 = -D(2*I-3)
24       C2 = D(2*I-2)
25       C3R = C1*C2R-C2*C2I
26       C3I = C1*C2I+C2*C2R
27       C(I) = C1R+C3R
28       C(I+1) = C1I+C3I
29       C(130-I) = C1R-C3R
30    10 C(131-I) = C3I-C1I
31       C(1) = A(1)+A(129)
32       C(2) = A(1)-A(129)
33       C(65) = 2.*A(65)
34       C(66) = 2.*(A(64)-A(66))
35       CALL D64(C,C(2),2,C,C(2),2)
36       S2 = 2.*S2
37       B(1) = C(1)+S2
38       B(129) = C(1)-S2
39       B(65) = C(65)
40       DO 20 I=2,64
41       E1 = 0.5*(C(I)+C(130-I))
42       E2 = E(2*I-2)*(C(I)-C(130-I))
43       B(I) = E1+E2
44    20 B(130-I) = E1-E2
45       RETURN
46       END