Lines Matching refs:ch
114 static void dradf2(int ido,int l1,float *cc,float *ch,float *wa1){
123 ch[t1<<1]=cc[t1]+cc[t2];
124 ch[(t1<<1)+t3-1]=cc[t1]-cc[t2];
146 ch[t6]=cc[t5]+ti2;
147 ch[t4]=ti2-cc[t5];
148 ch[t6-1]=cc[t5-1]+tr2;
149 ch[t4-1]=cc[t5-1]-tr2;
161 ch[t1]=-cc[t2];
162 ch[t1-1]=cc[t3];
169 static void dradf4(int ido,int l1,float *cc,float *ch,float *wa1,
185 ch[t5=t3<<2]=tr1+tr2;
186 ch[(ido<<2)+t5-1]=tr2-tr1;
187 ch[(t5+=(ido<<1))-1]=cc[t3]-cc[t4];
188 ch[t5]=cc[t2]-cc[t1];
230 ch[t4-1]=tr1+tr2;
231 ch[t4]=ti1+ti2;
233 ch[t5-1]=tr3-ti4;
234 ch[t5]=tr4-ti3;
236 ch[t4+t6-1]=ti4+tr3;
237 ch[t4+t6]=tr4+ti3;
239 ch[t5+t6-1]=tr2-tr1;
240 ch[t5+t6]=ti1-ti2;
258 ch[t4-1]=tr1+cc[t6-1];
259 ch[t4+t5-1]=cc[t6-1]-tr1;
261 ch[t4]=ti1-cc[t1+t0];
262 ch[t4+t5]=ti1+cc[t1+t0];
272 float *c2,float *ch,float *ch2,float *wa){
300 ch[t2]=c1[t2];
319 ch[t3-1]=wa[idij-1]*c1[t3-1]+wa[idij]*c1[t3];
320 ch[t3]=wa[idij-1]*c1[t3]-wa[idij]*c1[t3-1];
336 ch[t3-1]=wa[idij-1]*c1[t3-1]+wa[idij]*c1[t3];
337 ch[t3]=wa[idij-1]*c1[t3]-wa[idij]*c1[t3-1];
360 c1[t5-1]=ch[t5-1]+ch[t6-1];
361 c1[t6-1]=ch[t5]-ch[t6];
362 c1[t5]=ch[t5]+ch[t6];
363 c1[t6]=ch[t6-1]-ch[t5-1];
379 c1[t5-1]=ch[t5-1]+ch[t6-1];
380 c1[t6-1]=ch[t5]-ch[t6];
381 c1[t5]=ch[t5]+ch[t6];
382 c1[t6]=ch[t6-1]-ch[t5-1];
403 c1[t3]=ch[t3]+ch[t4];
404 c1[t4]=ch[t4]-ch[t3];
469 for(i=0;i<ido;i++)cc[t4++]=ch[t3++];
481 cc[t2]=ch[t1];
503 cc[t5-1]=ch[t6];
504 cc[t5]=ch[t7];
530 cc[i+t7-1]=ch[i+t8-1]+ch[i+t9-1];
531 cc[ic+t6-1]=ch[i+t8-1]-ch[i+t9-1];
532 cc[i+t7]=ch[i+t8]+ch[i+t9];
533 cc[ic+t6]=ch[i+t9]-ch[i+t8];
560 cc[t7-1]=ch[t8-1]+ch[t9-1];
561 cc[t6-1]=ch[t8-1]-ch[t9-1];
562 cc[t7]=ch[t8]+ch[t9];
563 cc[t6]=ch[t9]-ch[t8];
573 static void drftf1(int n,float *c,float *ch,float *wa,int *ifac){
597 dradf4(ido,l1,ch,c,wa+iw-1,wa+ix2-1,wa+ix3-1);
599 dradf4(ido,l1,c,ch,wa+iw-1,wa+ix2-1,wa+ix3-1);
606 dradf2(ido,l1,c,ch,wa+iw-1);
610 dradf2(ido,l1,ch,c,wa+iw-1);
617 dradfg(ido,ip,l1,idl1,c,c,c,ch,ch,wa+iw-1);
622 dradfg(ido,ip,l1,idl1,ch,ch,ch,c,c,wa+iw-1);
631 for(i=0;i<n;i++)c[i]=ch[i];
634 static void dradb2(int ido,int l1,float *cc,float *ch,float *wa1){
644 ch[t1]=cc[t2]+cc[t3+t2];
645 ch[t1+t0]=cc[t2]-cc[t3+t2];
663 ch[t3-1]=cc[t4-1]+cc[t5-1];
665 ch[t3]=cc[t4]-cc[t5];
667 ch[t6-1]=wa1[i-2]*tr2-wa1[i-1]*ti2;
668 ch[t6]=wa1[i-2]*ti2+wa1[i-1]*tr2;
679 ch[t1]=cc[t2]+cc[t2];
680 ch[t1+t0]=-(cc[t2+1]+cc[t2+1]);
686 static void dradb3(int ido,int l1,float *cc,float *ch,float *wa1,
702 ch[t1]=cc[t5]+tr2;
704 ch[t1+t0]=cr2-ci3;
705 ch[t1+t2]=cr2+ci3;
730 ch[t8-1]=cc[t7-1]+tr2;
733 ch[t8]=cc[t7]+ti2;
740 ch[t9-1]=wa1[i-2]*dr2-wa1[i-1]*di2;
741 ch[t9]=wa1[i-2]*di2+wa1[i-1]*dr2;
742 ch[t10-1]=wa2[i-2]*dr3-wa2[i-1]*di3;
743 ch[t10]=wa2[i-2]*di3+wa2[i-1]*dr3;
749 static void dradb4(int ido,int l1,float *cc,float *ch,float *wa1,
767 ch[t5]=tr2+tr3;
768 ch[t5+=t0]=tr1-tr4;
769 ch[t5+=t0]=tr2-tr3;
770 ch[t5+=t0]=tr1+tr4;
796 ch[t7-1]=tr2+tr3;
798 ch[t7]=ti2+ti3;
805 ch[(t8=t7+t0)-1]=wa1[i-2]*cr2-wa1[i-1]*ci2;
806 ch[t8]=wa1[i-2]*ci2+wa1[i-1]*cr2;
807 ch[(t8+=t0)-1]=wa2[i-2]*cr3-wa2[i-1]*ci3;
808 ch[t8]=wa2[i-2]*ci3+wa2[i-1]*cr3;
809 ch[(t8+=t0)-1]=wa3[i-2]*cr4-wa3[i-1]*ci4;
810 ch[t8]=wa3[i-2]*ci4+wa3[i-1]*cr4;
829 ch[t5]=tr2+tr2;
830 ch[t5+=t0]=sqrt2*(tr1-ti1);
831 ch[t5+=t0]=ti2+ti2;
832 ch[t5+=t0]=-sqrt2*(tr1+ti1);
841 float *c2,float *ch,float *ch2,float *wa){
866 ch[t3]=cc[t4];
881 ch[t2]=cc[t3];
899 ch[t3]=cc[t6-1]+cc[t6-1];
900 ch[t4]=cc[t6]+cc[t6];
932 ch[t5-1]=cc[t9-1]+cc[t11-1];
933 ch[t6-1]=cc[t9-1]-cc[t11-1];
934 ch[t5]=cc[t9]-cc[t11];
935 ch[t6]=cc[t9]+cc[t11];
966 ch[t5-1]=cc[t11-1]+cc[t12-1];
967 ch[t6-1]=cc[t11-1]-cc[t12-1];
968 ch[t5]=cc[t11]-cc[t12];
969 ch[t6]=cc[t11]+cc[t12];
1039 ch[t3]=c1[t3]-c1[t4];
1040 ch[t4]=c1[t3]+c1[t4];
1062 ch[t5-1]=c1[t5-1]-c1[t6];
1063 ch[t6-1]=c1[t5-1]+c1[t6];
1064 ch[t5]=c1[t5]+c1[t6-1];
1065 ch[t6]=c1[t5]-c1[t6-1];
1087 ch[t5-1]=c1[t5-1]-c1[t6];
1088 ch[t6-1]=c1[t5-1]+c1[t6];
1089 ch[t5]=c1[t5]+c1[t6-1];
1090 ch[t6]=c1[t5]-c1[t6-1];
1106 c1[t2]=ch[t2];
1125 c1[t3-1]=wa[idij-1]*ch[t3-1]-wa[idij]*ch[t3];
1126 c1[t3]=wa[idij-1]*ch[t3]+wa[idij]*ch[t3-1];
1146 c1[t3-1]=wa[idij-1]*ch[t3-1]-wa[idij]*ch[t3];
1147 c1[t3]=wa[idij-1]*ch[t3]+wa[idij]*ch[t3-1];
1154 static void drftb1(int n, float *c, float *ch, float *wa, int *ifac){
1174 dradb4(ido,l1,ch,c,wa+iw-1,wa+ix2-1,wa+ix3-1);
1176 dradb4(ido,l1,c,ch,wa+iw-1,wa+ix2-1,wa+ix3-1);
1184 dradb2(ido,l1,ch,c,wa+iw-1);
1186 dradb2(ido,l1,c,ch,wa+iw-1);
1195 dradb3(ido,l1,ch,c,wa+iw-1,wa+ix2-1);
1197 dradb3(ido,l1,c,ch,wa+iw-1,wa+ix2-1);
1209 dradb5(ido,l1,ch,c,wa+iw-1,wa+ix2-1,wa+ix3-1,wa+ix4-1);
1211 dradb5(ido,l1,c,ch,wa+iw-1,wa+ix2-1,wa+ix3-1,wa+ix4-1);
1217 dradbg(ido,ip,l1,idl1,ch,ch,ch,c,c,wa+iw-1);
1219 dradbg(ido,ip,l1,idl1,c,c,c,ch,ch,wa+iw-1);
1229 for(i=0;i<n;i++)c[i]=ch[i];
