348     double xi;  in rdft()  local
371         xi = a[0] - a[1];  in rdft()
373         a[1] = xi;  in rdft()
499     double xr, xi, yr, yi;  in dfct()  local
516     xi = a[0] + a[n];  in dfct()
518     t[0] = xi - yi;  in dfct()
519     t[m] = xi + yi;  in dfct()
525             xi = a[j] + a[n - j];  in dfct()
530             t[j] = xi - yi;  in dfct()
531             t[k] = xi + yi;  in dfct()
592     double xr, xi, yr, yi;  in dfst()  local
613             xi = a[j] - a[n - j];  in dfst()
618             t[j] = xi + yi;  in dfst()
619             t[k] = xi - yi;  in dfst()
730     double xr, xi, yr, yi;  in bitrv2()  local
750                 xi = a[j1 + 1];  in bitrv2()
756                 a[k1 + 1] = xi;  in bitrv2()
760                 xi = a[j1 + 1];  in bitrv2()
766                 a[k1 + 1] = xi;  in bitrv2()
770                 xi = a[j1 + 1];  in bitrv2()
776                 a[k1 + 1] = xi;  in bitrv2()
780                 xi = a[j1 + 1];  in bitrv2()
786                 a[k1 + 1] = xi;  in bitrv2()
791             xi = a[j1 + 1];  in bitrv2()
797             a[k1 + 1] = xi;  in bitrv2()
805                 xi = a[j1 + 1];  in bitrv2()
811                 a[k1 + 1] = xi;  in bitrv2()
815                 xi = a[j1 + 1];  in bitrv2()
821                 a[k1 + 1] = xi;  in bitrv2()
831     double xr, xi, yr, yi;  in bitrv2conj()  local
851                 xi = -a[j1 + 1];  in bitrv2conj()
857                 a[k1 + 1] = xi;  in bitrv2conj()
861                 xi = -a[j1 + 1];  in bitrv2conj()
867                 a[k1 + 1] = xi;  in bitrv2conj()
871                 xi = -a[j1 + 1];  in bitrv2conj()
877                 a[k1 + 1] = xi;  in bitrv2conj()
881                 xi = -a[j1 + 1];  in bitrv2conj()
887                 a[k1 + 1] = xi;  in bitrv2conj()
894             xi = -a[j1 + 1];  in bitrv2conj()
900             a[k1 + 1] = xi;  in bitrv2conj()
912                 xi = -a[j1 + 1];  in bitrv2conj()
918                 a[k1 + 1] = xi;  in bitrv2conj()
922                 xi = -a[j1 + 1];  in bitrv2conj()
928                 a[k1 + 1] = xi;  in bitrv2conj()
1273     double wkr, wki, xr, xi, yr, yi;  in rftfsub()  local
1284         xi = a[j + 1] + a[k + 1];  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
1310         xi = a[j + 1] + a[k + 1];  in rftbsub()
1311         yr = wkr * xr + wki * xi;  in rftbsub()
1312         yi = wkr * xi - wki * xr;  in rftbsub()