]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 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 |