44{
47 size_t k=m;
48 const size_t m2=2*m;
49 const size_t m3=3*m;
50
51
53
54 do {
56
57 C_MUL(scratch[0],Fout[m] , *tw1 );
58 C_MUL(scratch[1],Fout[m2] , *tw2 );
59 C_MUL(scratch[2],Fout[m3] , *tw3 );
60
61 C_SUB( scratch[5] , *Fout, scratch[1] );
63 C_ADD( scratch[3] , scratch[0] , scratch[2] );
64 C_SUB( scratch[4] , scratch[0] , scratch[2] );
65 C_SUB( Fout[m2], *Fout, scratch[3] );
66 tw1 += fstride;
67 tw2 += fstride*2;
68 tw3 += fstride*3;
70
72 Fout[m].
r = scratch[5].
r - scratch[4].
i;
73 Fout[m].
i = scratch[5].
i + scratch[4].
r;
74 Fout[m3].
r = scratch[5].
r + scratch[4].
i;
75 Fout[m3].
i = scratch[5].
i - scratch[4].
r;
76 }else{
77 Fout[m].
r = scratch[5].
r + scratch[4].
i;
78 Fout[m].
i = scratch[5].
i - scratch[4].
r;
79 Fout[m3].
r = scratch[5].
r - scratch[4].
i;
80 Fout[m3].
i = scratch[5].
i + scratch[4].
r;
81 }
82 ++Fout;
83 }while(--k);
84}