315 static void dctsub(int n, double *a, int nc, double const *c);
316 static void dstsub(int n, double *a, int nc, double const *c);
317 static void makect(int nc, int *ip, double *c);
319 static void rftbsub(int n, double *a, int nc, double const *c);
320 static void rftfsub(int n, double *a, int nc, double const *c);
347     int nw, nc;  in rdft()  local
358     nc = ip[1];  in rdft()
359     if (n > (nc << 2)) {  in rdft()
360         nc = n >> 2;  in rdft()
361         makect(nc, ip, w + nw);  in rdft()
367             rftfsub(n, a, nc, w + nw);  in rdft()
378             rftbsub(n, a, nc, w + nw);  in rdft()
390     int j, nw, nc;  in ddct()  local
401     nc = ip[1];  in ddct()
402     if (n > nc) {  in ddct()
403         nc = n;  in ddct()
404         makect(nc, ip, w + nw);  in ddct()
415             rftbsub(n, a, nc, w + nw);  in ddct()
422     dctsub(n, a, nc, w + nw);  in ddct()
427             rftfsub(n, a, nc, w + nw);  in ddct()
444     int j, nw, nc;  in ddst()  local
455     nc = ip[1];  in ddst()
456     if (n > nc) {  in ddst()
457         nc = n;  in ddst()
458         makect(nc, ip, w + nw);  in ddst()
469             rftbsub(n, a, nc, w + nw);  in ddst()
476     dstsub(n, a, nc, w + nw);  in ddst()
481             rftfsub(n, a, nc, w + nw);  in ddst()
498     int j, k, l, m, mh, nw, nc;  in dfct()  local
509     nc = ip[1];  in dfct()
510     if (n > (nc << 1)) {  in dfct()
511         nc = n >> 1;  in dfct()
512         makect(nc, ip, w + nw);  in dfct()
535         dctsub(m, a, nc, w + nw);  in dfct()
539             rftfsub(m, a, nc, w + nw);  in dfct()
552             dctsub(m, t, nc, w + nw);  in dfct()
556                 rftfsub(m, t, nc, w + nw);  in dfct()
591     int j, k, l, m, mh, nw, nc;  in dfst()  local
602     nc = ip[1];  in dfst()
603     if (n > (nc << 1)) {  in dfst()
604         nc = n >> 1;  in dfst()
605         makect(nc, ip, w + nw);  in dfst()
624         dstsub(m, a, nc, w + nw);  in dfst()
628             rftfsub(m, a, nc, w + nw);  in dfst()
641             dstsub(m, t, nc, w + nw);  in dfst()
645                 rftfsub(m, t, nc, w + nw);  in dfst()
705 static void makect(int nc, int *ip, double *c)  in makect()  argument
710     ip[1] = nc;  in makect()
711     if (nc > 1) {  in makect()
712         nch = nc >> 1;  in makect()
718             c[nc - j] = 0.5 * sin(delta * j);  in makect()
1270 static void rftfsub(int n, double *a, int nc, double const *c)  in rftfsub()  argument
1276     ks = 2 * nc / m;  in rftfsub()
1281         wkr = 0.5 - c[nc - kk];  in rftfsub()
1295 static void rftbsub(int n, double *a, int nc, double const *c)  in rftbsub()  argument
1302     ks = 2 * nc / m;  in rftbsub()
1307         wkr = 0.5 - c[nc - kk];  in rftbsub()
1322 static void dctsub(int n, double *a, int nc, double const *c)  in dctsub()  argument
1328     ks = nc / n;  in dctsub()
1333         wkr = c[kk] - c[nc - kk];  in dctsub()
1334         wki = c[kk] + c[nc - kk];  in dctsub()
1343 static void dstsub(int n, double *a, int nc, double const *c)  in dstsub()  argument
1349     ks = nc / n;  in dstsub()
1354         wkr = c[kk] - c[nc - kk];  in dstsub()
1355         wki = c[kk] + c[nc - kk];  in dstsub()