117   int t0,t1,t2,t3,t4,t5,t6;  in dradf2()  local
121   t3=ido<<1;  in dradf2()
124     ch[(t1<<1)+t3-1]=cc[t1]-cc[t2];  in dradf2()
135     t3=t2;  in dradf2()
140       t3+=2;  in dradf2()
144       tr2=wa1[i-2]*cc[t3-1]+wa1[i-1]*cc[t3];  in dradf2()
145       ti2=wa1[i-2]*cc[t3]-wa1[i-1]*cc[t3-1];  in dradf2()
158   t3=(t2=(t1=ido)-1);  in dradf2()
162     ch[t1-1]=cc[t3];  in dradf2()
165     t3+=ido;  in dradf2()
172   int i,k,t0,t1,t2,t3,t4,t5,t6;  in dradf4()  local
179   t3=0;  in dradf4()
183     tr2=cc[t3]+cc[t4];  in dradf4()
185     ch[t5=t3<<2]=tr1+tr2;  in dradf4()
187     ch[(t5+=(ido<<1))-1]=cc[t3]-cc[t4];  in dradf4()
192     t3+=ido;  in dradf4()
206       t3=(t2+=2);  in dradf4()
210       t3+=t0;  in dradf4()
211       cr2=wa1[i-2]*cc[t3-1]+wa1[i-1]*cc[t3];  in dradf4()
212       ci2=wa1[i-2]*cc[t3]-wa1[i-1]*cc[t3-1];  in dradf4()
213       t3+=t0;  in dradf4()
214       cr3=wa2[i-2]*cc[t3-1]+wa2[i-1]*cc[t3];  in dradf4()
215       ci3=wa2[i-2]*cc[t3]-wa2[i-1]*cc[t3-1];  in dradf4()
216       t3+=t0;  in dradf4()
217       cr4=wa3[i-2]*cc[t3-1]+wa3[i-1]*cc[t3];  in dradf4()
218       ci4=wa3[i-2]*cc[t3]-wa3[i-1]*cc[t3-1];  in dradf4()
249   t3=ido<<2;  in dradf4()
266     t4+=t3;  in dradf4()
276   int t0,t1,t2,t3,t4,t5,t6,t7,t8,t9,t10;  in dradfg()  local
315         t3=t2;  in dradfg()
318           t3+=2;  in dradfg()
319           ch[t3-1]=wa[idij-1]*c1[t3-1]+wa[idij]*c1[t3];  in dradfg()
320           ch[t3]=wa[idij-1]*c1[t3]-wa[idij]*c1[t3-1];  in dradfg()
334         t3=t2;  in dradfg()
336           ch[t3-1]=wa[idij-1]*c1[t3-1]+wa[idij]*c1[t3];  in dradfg()
337           ch[t3]=wa[idij-1]*c1[t3]-wa[idij]*c1[t3-1];  in dradfg()
338           t3+=ido;  in dradfg()
350       t3=t1;  in dradfg()
353         t3+=2;  in dradfg()
355         t5=t3-ido;  in dradfg()
371       t3=t1;  in dradfg()
374         t5=t3;  in dradfg()
384         t3+=ido;  in dradfg()
398     t3=t1-ido;  in dradfg()
401       t3+=ido;  in dradfg()
403       c1[t3]=ch[t3]+ch[t4];  in dradfg()
404       c1[t4]=ch[t4]-ch[t3];  in dradfg()
412   t3=(ip-1)*idl1;  in dradfg()
421     t6=t3;  in dradfg()
467     t3=t1;  in dradfg()
469     for(i=0;i<ido;i++)cc[t4++]=ch[t3++];  in dradfg()
490   t3=0;  in dradfg()
495     t3+=t0;  in dradfg()
499     t6=t3;  in dradfg()
515   t3=0;  in dradfg()
520     t3+=t2;  in dradfg()
524     t7=t3;  in dradfg()
546   t3=0;  in dradfg()
551     t3+=t2;  in dradfg()
556       t7=i+t3;  in dradfg()
635   int i,k,t0,t1,t2,t3,t4,t5,t6;  in dradb2()  local
642   t3=(ido<<1)-1;  in dradb2()
644     ch[t1]=cc[t2]+cc[t3+t2];  in dradb2()
645     ch[t1+t0]=cc[t2]-cc[t3+t2];  in dradb2()
655     t3=t1;  in dradb2()
659       t3+=2;  in dradb2()
663       ch[t3-1]=cc[t4-1]+cc[t5-1];  in dradb2()
665       ch[t3]=cc[t4]-cc[t5];  in dradb2()
690   int i,k,t0,t1,t2,t3,t4,t5,t6,t7,t8,t9,t10;  in dradb3()  local
696   t3=ido<<1;  in dradb3()
700     tr2=cc[t3-1]+cc[t3-1];  in dradb3()
703     ci3=taui*(cc[t3]+cc[t3]);  in dradb3()
707     t3+=t4;  in dradb3()
714   t3=ido<<1;  in dradb3()
717     t6=(t5=t7+t3);  in dradb3()
752   int i,k,t0,t1,t2,t3,t4,t5,t6,t7,t8;  in dradb4()  local
758   t3=0;  in dradb4()
761     t4=t3+t6;  in dradb4()
765     tr1=cc[t3]-cc[(t4+=t6)-1];  in dradb4()
766     tr2=cc[t3]+cc[t4-1];  in dradb4()
772     t3+=t2;  in dradb4()
780     t5=(t4=(t3=(t2=t1<<2)+t6))+t6;  in dradb4()
784       t3+=2;  in dradb4()
790       ti3=cc[t3]-cc[t4];  in dradb4()
791       tr4=cc[t3]+cc[t4];  in dradb4()
794       ti4=cc[t3-1]-cc[t4-1];  in dradb4()
795       tr3=cc[t3-1]+cc[t4-1];  in dradb4()
821   t3=ido-1;  in dradb4()
824     t5=t3;  in dradb4()
834     t3+=ido;  in dradb4()
843   int idij,ipph,i,j,k,l,ik,is,t0,t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,  in dradbg()  local
863     t3=t1;  in dradbg()
866       ch[t3]=cc[t4];  in dradbg()
867       t3++;  in dradbg()
879     t3=t1;  in dradbg()
881       ch[t2]=cc[t3];  in dradbg()
883       t3+=t10;  in dradbg()
895     t3=t1;  in dradbg()
899       ch[t3]=cc[t6-1]+cc[t6-1];  in dradbg()
901       t3+=ido;  in dradbg()
917     t3=t1;  in dradbg()
923       t5=t3;  in dradbg()
937       t3+=ido;  in dradbg()
951     t3=t1;  in dradbg()
957       t3+=2;  in dradbg()
961       t5=t3;  in dradbg()
983   t3=(ip-1)*idl1;  in dradbg()
995     t8=t3;  in dradbg()
1036     t3=t1;  in dradbg()
1039       ch[t3]=c1[t3]-c1[t4];  in dradbg()
1040       ch[t4]=c1[t3]+c1[t4];  in dradbg()
1041       t3+=ido;  in dradbg()
1054     t3=t1;  in dradbg()
1057       t5=t3;  in dradbg()
1067       t3+=ido;  in dradbg()
1079     t3=t1;  in dradbg()
1082       t3+=2;  in dradbg()
1084       t5=t3;  in dradbg()
1123       t3=t2;  in dradbg()
1125         c1[t3-1]=wa[idij-1]*ch[t3-1]-wa[idij]*ch[t3];  in dradbg()
1126         c1[t3]=wa[idij-1]*ch[t3]+wa[idij]*ch[t3-1];  in dradbg()
1127         t3+=ido;  in dradbg()
1142       t3=t2;  in dradbg()
1145         t3+=2;  in dradbg()
1146         c1[t3-1]=wa[idij-1]*ch[t3-1]-wa[idij]*ch[t3];  in dradbg()
1147         c1[t3]=wa[idij-1]*ch[t3]+wa[idij]*ch[t3-1];  in dradbg()