391     double xr;  in ddct()  local
407         xr = a[n - 1];  in ddct()
412         a[1] = a[0] - xr;  in ddct()
413         a[0] += xr;  in ddct()
431         xr = a[0] - a[1];  in ddct()
437         a[n - 1] = xr;  in ddct()
445     double xr;  in ddst()  local
461         xr = a[n - 1];  in ddst()
466         a[1] = a[0] + xr;  in ddst()
467         a[0] -= xr;  in ddst()
485         xr = a[0] - a[1];  in ddst()
491         a[n - 1] = -xr;  in ddst()
499     double xr, xi, yr, yi;  in dfct()  local
524             xr = a[j] - a[n - j];  in dfct()
528             a[j] = xr;  in dfct()
592     double xr, xi, yr, yi;  in dfst()  local
612             xr = a[j] + a[n - j];  in dfst()
616             a[j] = xr;  in dfst()
730     double xr, xi, yr, yi;  in bitrv2()  local
749                 xr = a[j1];  in bitrv2()
755                 a[k1] = xr;  in bitrv2()
759                 xr = a[j1];  in bitrv2()
765                 a[k1] = xr;  in bitrv2()
769                 xr = a[j1];  in bitrv2()
775                 a[k1] = xr;  in bitrv2()
779                 xr = a[j1];  in bitrv2()
785                 a[k1] = xr;  in bitrv2()
790             xr = a[j1];  in bitrv2()
796             a[k1] = xr;  in bitrv2()
804                 xr = a[j1];  in bitrv2()
810                 a[k1] = xr;  in bitrv2()
814                 xr = a[j1];  in bitrv2()
820                 a[k1] = xr;  in bitrv2()
831     double xr, xi, yr, yi;  in bitrv2conj()  local
850                 xr = a[j1];  in bitrv2conj()
856                 a[k1] = xr;  in bitrv2conj()
860                 xr = a[j1];  in bitrv2conj()
866                 a[k1] = xr;  in bitrv2conj()
870                 xr = a[j1];  in bitrv2conj()
876                 a[k1] = xr;  in bitrv2conj()
880                 xr = a[j1];  in bitrv2conj()
886                 a[k1] = xr;  in bitrv2conj()
893             xr = a[j1];  in bitrv2conj()
899             a[k1] = xr;  in bitrv2conj()
911                 xr = a[j1];  in bitrv2conj()
917                 a[k1] = xr;  in bitrv2conj()
921                 xr = a[j1];  in bitrv2conj()
927                 a[k1] = xr;  in bitrv2conj()
1273     double wkr, wki, xr, xi, yr, yi;  in rftfsub()  local
1283         xr = a[j] - a[k];  in rftfsub()
1285         yr = wkr * xr - wki * xi;  in rftfsub()
1286         yi = wkr * xi + wki * xr;  in rftfsub()
1298     double wkr, wki, xr, xi, yr, yi;  in rftbsub()  local
1309         xr = a[j] - a[k];  in rftbsub()
1311         yr = wkr * xr + wki * xi;  in rftbsub()
1312         yi = wkr * xi - wki * xr;  in rftbsub()
1325     double wkr, wki, xr;  in dctsub()  local
1335         xr = wki * a[j] - wkr * a[k];  in dctsub()
1337         a[k] = xr;  in dctsub()
1346     double wkr, wki, xr;  in dstsub()  local
1356         xr = wki * a[k] - wkr * a[j];  in dstsub()
1358         a[j] = xr;  in dstsub()