25 static void SlowFT(const ColumnVector& a, const ColumnVector&b,  in SlowFT()  argument
28    int n = a.Nrows();  in SlowFT()
38 	 sumx += c * a(k) - s * b(k); sumy += s * a(k) + c * b(k);  in SlowFT()
44 static void SlowDTT_II(const ColumnVector& a, ColumnVector& c, ColumnVector& s)  in SlowDTT_II()  argument
46    int n = a.Nrows(); c.ReSize(n); s.ReSize(n);  in SlowDTT_II()
54       for (int j=1; j<=n; j++) sum += cos(k1 * (2*j-1) * f) * a(j);  in SlowDTT_II()
61       for (int j=1; j<=n; j++) sum += sin(k * (2*j-1) * f) * a(j);  in SlowDTT_II()
66 static void SlowDTT(const ColumnVector& a, ColumnVector& c, ColumnVector& s)  in SlowDTT()  argument
68    int n1 = a.Nrows(); int n = n1 - 1;  in SlowDTT()
77       for (int j=2; j<=n; j++) sum += cos((j-1) * (k-1) * f) * a(j);  in SlowDTT()
78       c(k) = sum + (a(1) + sign * a(n1)) / 2.0;  in SlowDTT()
85       for (int j=2; j<=n; j++) sum += sin((j-1) * (k-1) * f) * a(j);  in SlowDTT()
91 void SlowFT2(const Matrix& U, const Matrix& V, Matrix& X, Matrix& Y)  in SlowFT2()  argument
97    X.ReSize(U); Y.ReSize(V);  in SlowFT2()
104          Real a = -pi2 * ( (Real)k * (Real)i / (Real)m  in SlowFT2()  local
106          Real cs = cos(a); Real sn = sin(a);  in SlowFT2()
110       X(i+1,j+1) = sumr; Y(i+1,j+1) = sumi;  in SlowFT2()
120    ColumnVector A(n), B(n), X, Y;  in test()  local
123    FFT(A, B, X, Y); FFTI(X, Y, X, Y);  in test()
124    X = X - A; Y = Y - B;  in test()
125    Clean(X,0.000000001); Clean(Y,0.000000001); Print(X); Print(Y);  in test()
133    ColumnVector A(n), B(n), X, Y;  in test1()  local
136    FFT(A, B, X, Y);  in test1()
141    X2 = X1 - X.Rows(1,n2); Y2 = Y1 - Y.Rows(1,n2);  in test1()
154    ColumnVector A(n), B(n), X, Y, X1, Y1;  in test2()  local
157    FFT(A, B, X, Y);  in test2()
159    X = X - X1; Y = Y - Y1;  in test2()
160    Clean(X,0.000000001); Clean(Y,0.000000001); Print(X); Print(Y);  in test2()
168    ColumnVector A(n), X, Y, X1, Y1;  in test3()  local
170    DCT_II(A, X); DST_II(A, Y);  in test3()
172    X -= X1; Y -= Y1;  in test3()
173    Clean(X,0.000000001); Clean(Y,0.000000001); Print(X); Print(Y);  in test3()
181    ColumnVector A1(n);  in test4()  local
182    for (int i=0; i<n; i++) A1.element(i) = mwc.Next();  in test4()
185    P = 0.0; Q = 0.0; P.Rows(1,n) = A1;  in test4()
195    DCT_II(A1,B2);  in test4()
199    A1 -= A2; B1 -= B2;  in test4()
200    Clean(A1,0.000000001); Clean(B1,0.000000001); Print(A1); Print(B1);  in test4()
208    ColumnVector A1(n);  in test5()  local
209    for (int i=0; i<n; i++) A1.element(i) = mwc.Next();  in test5()
212    P = 0.0; Q = 0.0; P.Rows(1,n) = A1;  in test5()
222    DST_II(A1,B2);  in test5()
226    A1 -= A2; B1 -= B2;  in test5()
227    Clean(A1,0.000000001); Clean(B1,0.000000001); Print(A1); Print(B1);  in test5()
235    ColumnVector A1(n+1);  in test6()  local
236    A1(1) = A1(n+1) = 0;  in test6()
237    for (int i=1; i<n; i++) A1.element(i) = mwc.Next();  in test6()
240    ColumnVector P(2*n), Q(2*n); P = 0.0; Q = 0.0; P.Rows(1,n+1) = A1;  in test6()
245    DST(A1,B2);  in test6()
249    A1 -= A2; B1 -= B2;  in test6()
250    Clean(A1,0.000000001); Clean(B1,0.000000001); Print(A1); Print(B1);  in test6()
260    ColumnVector A1(n+1);  in test7()  local
261    for (int i=0; i<=n; i++) A1.element(i) = mwc.Next();  in test7()
264    ColumnVector P(2*n), Q(2*n); P = 0.0; Q = 0.0; P.Rows(1,n+1) = A1;  in test7()
270    DCT(A1,B2);  in test7()
274    A1 -= A2; B1 -= B2;  in test7()
275    Clean(A1,0.000000001); Clean(B1,0.000000001); Print(A1); Print(B1);  in test7()
283    ColumnVector A(n+1), X, Y, X1, Y1;  in test8()  local
286    DCT(A, X); DST(A, Y);  in test8()
288    X -= X1; Y -= Y1;  in test8()
289    Clean(X,0.000000001); Clean(Y,0.000000001); Print(X); Print(Y);  in test8()
297    Matrix A(m,n), B(m,n), X, Y, X1, Y1;  in test9()  local
300    FFT2(A, B, X, Y);  in test9()
302    X = X - X1; Y = Y - Y1;  in test9()
303    Clean(X,0.000000001); Clean(Y,0.000000001); Print(X); Print(Y);  in test9()