499     double xr, xi, yr, yi;  in dfct()  local
515     yi = a[m];  in dfct()
518     t[0] = xi - yi;  in dfct()
519     t[m] = xi + yi;  in dfct()
527             yi = a[k] + a[n - k];  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
615             yi = a[k] - a[n - k];  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
752                 yi = a[k1 + 1];  in bitrv2()
754                 a[j1 + 1] = yi;  in bitrv2()
762                 yi = a[k1 + 1];  in bitrv2()
764                 a[j1 + 1] = yi;  in bitrv2()
772                 yi = a[k1 + 1];  in bitrv2()
774                 a[j1 + 1] = yi;  in bitrv2()
782                 yi = a[k1 + 1];  in bitrv2()
784                 a[j1 + 1] = yi;  in bitrv2()
793             yi = a[k1 + 1];  in bitrv2()
795             a[j1 + 1] = yi;  in bitrv2()
807                 yi = a[k1 + 1];  in bitrv2()
809                 a[j1 + 1] = yi;  in bitrv2()
817                 yi = a[k1 + 1];  in bitrv2()
819                 a[j1 + 1] = yi;  in bitrv2()
831     double xr, xi, yr, yi;  in bitrv2conj()  local
853                 yi = -a[k1 + 1];  in bitrv2conj()
855                 a[j1 + 1] = yi;  in bitrv2conj()
863                 yi = -a[k1 + 1];  in bitrv2conj()
865                 a[j1 + 1] = yi;  in bitrv2conj()
873                 yi = -a[k1 + 1];  in bitrv2conj()
875                 a[j1 + 1] = yi;  in bitrv2conj()
883                 yi = -a[k1 + 1];  in bitrv2conj()
885                 a[j1 + 1] = yi;  in bitrv2conj()
896             yi = -a[k1 + 1];  in bitrv2conj()
898             a[j1 + 1] = yi;  in bitrv2conj()
914                 yi = -a[k1 + 1];  in bitrv2conj()
916                 a[j1 + 1] = yi;  in bitrv2conj()
924                 yi = -a[k1 + 1];  in bitrv2conj()
926                 a[j1 + 1] = yi;  in bitrv2conj()
1273     double wkr, wki, xr, xi, yr, yi;  in rftfsub()  local
1286         yi = wkr * xi + wki * xr;  in rftfsub()
1288         a[j + 1] -= yi;  in rftfsub()
1290         a[k + 1] -= yi;  in rftfsub()
1298     double wkr, wki, xr, xi, yr, yi;  in rftbsub()  local
1312         yi = wkr * xi - wki * xr;  in rftbsub()
1314         a[j + 1] = yi - a[j + 1];  in rftbsub()
1316         a[k + 1] = yi - a[k + 1];  in rftbsub()