45   int k1, l1, l2, ib;  in drfti1()  local
82   l1=1;  in drfti1()
89     l2=l1*ip;  in drfti1()
94       ld+=l1;  in drfti1()
106     l1=l2;  in drfti1()
116 static void dradf2(int ido,int l1,float *cc,float *ch,float *wa1){  in dradf2()  argument
122   t0=(t2=l1*ido);  in dradf2()
124   for(k=0;k<l1;k++){  in dradf2()
136   for(k=0;k<l1;k++){  in dradf2()
162   for(k=0;k<l1;k++){  in dradf2()
171 static void dradf4(int ido,int l1,float *cc,float *ch,float *wa1,  in dradf4()  argument
176   t0=l1*ido;  in dradf4()
183   for(k=0;k<l1;k++){  in dradf4()
203   for(k=0;k<l1;k++){  in dradf4()
256   for(k=0;k<l1;k++){  in dradf4()
273 static void dradfg(int ido,int ip,int l1,int idl1,float *cc,float *c1,  in dradfg()  argument
291   t0=l1*ido;  in dradfg()
301     for(k=0;k<l1;k++){  in dradfg()
309   if(nbd>l1){  in dradfg()
314       for(k=0;k<l1;k++){  in dradfg()
337         for(k=0;k<l1;k++){  in dradfg()
348   if(nbd<l1){  in dradfg()
359         for(k=0;k<l1;k++){  in dradfg()
375       for(k=0;k<l1;k++){  in dradfg()
402     for(k=0;k<l1;k++){  in dradfg()
464   if(ido<l1)goto L132;  in dradfg()
468   for(k=0;k<l1;k++){  in dradfg()
482     for(k=0;k<l1;k++){  in dradfg()
504     for(k=0;k<l1;k++){  in dradfg()
514   if(nbd<l1)goto L141;  in dradfg()
529     for(k=0;k<l1;k++){  in dradfg()
561       for(k=0;k<l1;k++){  in dradfg()
576   int i,k1,l1,l2;  in drftf1()  local
588     l1=l2/ip;  in drftf1()
590     idl1=ido*l1;  in drftf1()
599       dradf4(ido,l1,ch,c,wa+iw-1,wa+ix2-1,wa+ix3-1);  in drftf1()
601       dradf4(ido,l1,c,ch,wa+iw-1,wa+ix2-1,wa+ix3-1);  in drftf1()
608     dradf2(ido,l1,c,ch,wa+iw-1);  in drftf1()
612     dradf2(ido,l1,ch,c,wa+iw-1);  in drftf1()
619     dradfg(ido,ip,l1,idl1,c,c,c,ch,ch,wa+iw-1);  in drftf1()
624     dradfg(ido,ip,l1,idl1,ch,ch,ch,c,c,wa+iw-1);  in drftf1()
628     l2=l1;  in drftf1()
636 static void dradb2(int ido,int l1,float *cc,float *ch,float *wa1){  in dradb2()  argument
640   t0=l1*ido;  in dradb2()
645   for(k=0;k<l1;k++){  in dradb2()
656   for(k=0;k<l1;k++){  in dradb2()
680   for(k=0;k<l1;k++){  in dradb2()
688 static void dradb3(int ido,int l1,float *cc,float *ch,float *wa1,  in dradb3()  argument
694   t0=l1*ido;  in dradb3()
701   for(k=0;k<l1;k++){  in dradb3()
717   for(k=0;k<l1;k++){  in dradb3()
751 static void dradb4(int ido,int l1,float *cc,float *ch,float *wa1,  in dradb4()  argument
756   t0=l1*ido;  in dradb4()
762   for(k=0;k<l1;k++){  in dradb4()
781   for(k=0;k<l1;k++){  in dradb4()
825   for(k=0;k<l1;k++){  in dradb4()
842 static void dradbg(int ido,int ip,int l1,int idl1,float *cc,float *c1,  in dradbg()  argument
853   t0=l1*ido;  in dradbg()
860   if(ido<l1)goto L103;  in dradbg()
864   for(k=0;k<l1;k++){  in dradbg()
882     for(k=0;k<l1;k++){  in dradbg()
900     for(k=0;k<l1;k++){  in dradbg()
911   if(nbd<l1)goto L112;  in dradbg()
924     for(k=0;k<l1;k++){  in dradbg()
967       for(k=0;k<l1;k++){  in dradbg()
1040     for(k=0;k<l1;k++){  in dradbg()
1049   if(nbd<l1)goto L128;  in dradbg()
1058     for(k=0;k<l1;k++){  in dradbg()
1088       for(k=0;k<l1;k++){  in dradbg()
1107     for(k=0;k<l1;k++){  in dradbg()
1113   if(nbd>l1)goto L139;  in dradbg()
1126       for(k=0;k<l1;k++){  in dradbg()
1142     for(k=0;k<l1;k++){  in dradbg()
1157   int i,k1,l1,l2;  in drftb1()  local
1163   l1=1;  in drftb1()
1168     l2=ip*l1;  in drftb1()
1170     idl1=ido*l1;  in drftb1()
1176       dradb4(ido,l1,ch,c,wa+iw-1,wa+ix2-1,wa+ix3-1);  in drftb1()
1178       dradb4(ido,l1,c,ch,wa+iw-1,wa+ix2-1,wa+ix3-1);  in drftb1()
1186       dradb2(ido,l1,ch,c,wa+iw-1);  in drftb1()
1188       dradb2(ido,l1,c,ch,wa+iw-1);  in drftb1()
1197       dradb3(ido,l1,ch,c,wa+iw-1,wa+ix2-1);  in drftb1()
1199       dradb3(ido,l1,c,ch,wa+iw-1,wa+ix2-1);  in drftb1()
1219       dradbg(ido,ip,l1,idl1,ch,ch,ch,c,c,wa+iw-1);  in drftb1()
1221       dradbg(ido,ip,l1,idl1,c,c,c,ch,ch,wa+iw-1);  in drftb1()
1225     l1=l2;  in drftb1()