119   int t0,t1,t2,t3,t4,t5,t6;  in dradf2()  local
139     t5=t1;  in dradf2()
144       t5+=2;  in dradf2()
148       ch[t6]=cc[t5]+ti2;  in dradf2()
149       ch[t4]=ti2-cc[t5];  in dradf2()
150       ch[t6-1]=cc[t5-1]+tr2;  in dradf2()
151       ch[t4-1]=cc[t5-1]-tr2;  in dradf2()
174   int i,k,t0,t1,t2,t3,t4,t5,t6;  in dradf4()  local
187     ch[t5=t3<<2]=tr1+tr2;  in dradf4()
188     ch[(ido<<2)+t5-1]=tr2-tr1;  in dradf4()
189     ch[(t5+=(ido<<1))-1]=cc[t3]-cc[t4];  in dradf4()
190     ch[t5]=cc[t2]-cc[t1];  in dradf4()
206     t5=(t6=ido<<1)+t4;  in dradf4()
210       t5-=2;  in dradf4()
235       ch[t5-1]=tr3-ti4;  in dradf4()
236       ch[t5]=tr4-ti3;  in dradf4()
241       ch[t5+t6-1]=tr2-tr1;  in dradf4()
242       ch[t5+t6]=ti1-ti2;  in dradf4()
253   t5=ido<<1;  in dradf4()
261     ch[t4+t5-1]=cc[t6-1]-tr1;  in dradf4()
264     ch[t4+t5]=ti1+cc[t1+t0];  in dradf4()
278   int t0,t1,t2,t3,t4,t5,t6,t7,t8,t9,t10;  in dradfg()  local
357         t5=t3-ido;  in dradfg()
360           t5+=ido;  in dradfg()
362           c1[t5-1]=ch[t5-1]+ch[t6-1];  in dradfg()
363           c1[t6-1]=ch[t5]-ch[t6];  in dradfg()
364           c1[t5]=ch[t5]+ch[t6];  in dradfg()
365           c1[t6]=ch[t6-1]-ch[t5-1];  in dradfg()
376         t5=t3;  in dradfg()
379           t5+=2;  in dradfg()
381           c1[t5-1]=ch[t5-1]+ch[t6-1];  in dradfg()
382           c1[t6-1]=ch[t5]-ch[t6];  in dradfg()
383           c1[t5]=ch[t5]+ch[t6];  in dradfg()
384           c1[t6]=ch[t6-1]-ch[t5-1];  in dradfg()
422     t5=t2;  in dradfg()
428       ch2[t5++]=ai1*c2[t6++];  in dradfg()
437     t5=(ipp2-1)*idl1;  in dradfg()
440       t5-=idl1;  in dradfg()
449       t9=t5;  in dradfg()
500     t5=t1;  in dradfg()
505       cc[t5-1]=ch[t6];  in dradfg()
506       cc[t5]=ch[t7];  in dradfg()
507       t5+=t10;  in dradfg()
519   t5=ipp2*t0;  in dradfg()
524     t5-=t0;  in dradfg()
528     t9=t5;  in dradfg()
550   t5=ipp2*t0;  in dradfg()
555     t5-=t0;  in dradfg()
560       t9=i+t5;  in dradfg()
637   int i,k,t0,t1,t2,t3,t4,t5,t6;  in dradb2()  local
658     t5=(t4=t2)+(ido<<1);  in dradb2()
663       t5-=2;  in dradb2()
665       ch[t3-1]=cc[t4-1]+cc[t5-1];  in dradb2()
666       tr2=cc[t4-1]-cc[t5-1];  in dradb2()
667       ch[t3]=cc[t4]-cc[t5];  in dradb2()
668       ti2=cc[t4]+cc[t5];  in dradb2()
692   int i,k,t0,t1,t2,t3,t4,t5,t6,t7,t8,t9,t10;  in dradb3()  local
700   t5=0;  in dradb3()
703     cr2=cc[t5]+(taur*tr2);  in dradb3()
704     ch[t1]=cc[t5]+tr2;  in dradb3()
710     t5+=t4;  in dradb3()
719     t6=(t5=t7+t3);  in dradb3()
724       t5+=2;  in dradb3()
730       tr2=cc[t5-1]+cc[t6-1];  in dradb3()
733       ti2=cc[t5]-cc[t6];  in dradb3()
736       cr3=taui*(cc[t5-1]-cc[t6-1]);  in dradb3()
737       ci3=taui*(cc[t5]+cc[t6]);  in dradb3()
754   int i,k,t0,t1,t2,t3,t4,t5,t6,t7,t8;  in dradb4()  local
764     t5=t1;  in dradb4()
769     ch[t5]=tr2+tr3;  in dradb4()
770     ch[t5+=t0]=tr1-tr4;  in dradb4()
771     ch[t5+=t0]=tr2-tr3;  in dradb4()
772     ch[t5+=t0]=tr1+tr4;  in dradb4()
782     t5=(t4=(t3=(t2=t1<<2)+t6))+t6;  in dradb4()
788       t5-=2;  in dradb4()
790       ti1=cc[t2]+cc[t5];  in dradb4()
791       ti2=cc[t2]-cc[t5];  in dradb4()
794       tr1=cc[t2-1]-cc[t5-1];  in dradb4()
795       tr2=cc[t2-1]+cc[t5-1];  in dradb4()
826     t5=t3;  in dradb4()
831     ch[t5]=tr2+tr2;  in dradb4()
832     ch[t5+=t0]=sqrt2*(tr1-ti1);  in dradb4()
833     ch[t5+=t0]=ti2+ti2;  in dradb4()
834     ch[t5+=t0]=-sqrt2*(tr1+ti1);  in dradb4()
845   int idij,ipph,i,j,k,l,ik,is,t0,t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,  in dradbg()  local
893   t7=(t5=ido<<1);  in dradbg()
899     t6=t5;  in dradbg()
907     t5+=t7;  in dradbg()
925       t5=t3;  in dradbg()
930         t5+=2;  in dradbg()
934         ch[t5-1]=cc[t9-1]+cc[t11-1];  in dradbg()
936         ch[t5]=cc[t9]-cc[t11];  in dradbg()
963       t5=t3;  in dradbg()
968         ch[t5-1]=cc[t11-1]+cc[t12-1];  in dradbg()
970         ch[t5]=cc[t11]-cc[t12];  in dradbg()
972         t5+=ido;  in dradbg()
994     t5=t2;  in dradbg()
1000       c2[t5++]=ai1*ch2[t8++];  in dradbg()
1016       t5=t2;  in dradbg()
1021         c2[t5++]+=ai2*ch2[t12++];  in dradbg()
1059       t5=t3;  in dradbg()
1062         t5+=2;  in dradbg()
1064         ch[t5-1]=c1[t5-1]-c1[t6];  in dradbg()
1065         ch[t6-1]=c1[t5-1]+c1[t6];  in dradbg()
1066         ch[t5]=c1[t5]+c1[t6-1];  in dradbg()
1067         ch[t6]=c1[t5]-c1[t6-1];  in dradbg()
1086       t5=t3;  in dradbg()
1089         ch[t5-1]=c1[t5-1]-c1[t6];  in dradbg()
1090         ch[t6-1]=c1[t5-1]+c1[t6];  in dradbg()
1091         ch[t5]=c1[t5]+c1[t6-1];  in dradbg()
1092         ch[t6]=c1[t5]-c1[t6-1];  in dradbg()
1093         t5+=ido;  in dradbg()