📄 gpfa5f.c
字号:
/*< u5 = u1 + u2 >*/
u5 = u1 + u2;
/*< u6 = c1 * ( u1 - u2 ) >*/
u6 = c1 * (u1 - u2);
/*< bja = b(ja+j) >*/
bja = b[ja + j];
/*< u7 = bja - 0.25 * u5 >*/
u7 = bja - u5 * (float).25;
/*< b(ja+j) = bja + u5 >*/
b[ja + j] = bja + u5;
/*< u8 = u7 + u6 >*/
u8 = u7 + u6;
/*< u9 = u7 - u6 >*/
u9 = u7 - u6;
/*< u10 = c3 * u3 - c2 * u4 >*/
u10 = c3 * u3 - c2 * u4;
/*< u11 = c2 * u3 + c3 * u4 >*/
u11 = c2 * u3 + c3 * u4;
/*< a(jb+j) = t8 - u11 >*/
a[jb + j] = t8 - u11;
/*< b(jb+j) = u8 + t11 >*/
b[jb + j] = u8 + t11;
/*< a(je+j) = t8 + u11 >*/
a[je + j] = t8 + u11;
/*< b(je+j) = u8 - t11 >*/
b[je + j] = u8 - t11;
/*< a(jc+j) = t9 - u10 >*/
a[jc + j] = t9 - u10;
/*< b(jc+j) = u9 + t10 >*/
b[jc + j] = u9 + t10;
/*< a(jd+j) = t9 + u10 >*/
a[jd + j] = t9 + u10;
/*< b(jd+j) = u9 - t10 >*/
b[jd + j] = u9 - t10;
/*< j = j + jump >*/
j += *jump;
/*< 110 continue >*/
/* L110: */
}
/*< else >*/
} else {
/* dir$ ivdep,shortloop */
/*< do 130 l = 1 , nvex >*/
i__8 = nvex;
for (l = 1; l <= i__8; ++l) {
/*< ajb = a(jb+j) >*/
ajb = a[jb + j];
/*< aje = a(je+j) >*/
aje = a[je + j];
/*< t1 = ajb + aje >*/
t1 = ajb + aje;
/*< ajc = a(jc+j) >*/
ajc = a[jc + j];
/*< ajd = a(jd+j) >*/
ajd = a[jd + j];
/*< t2 = ajc + ajd >*/
t2 = ajc + ajd;
/*< t3 = ajb - aje >*/
t3 = ajb - aje;
/*< t4 = ajc - ajd >*/
t4 = ajc - ajd;
/*< t5 = t1 + t2 >*/
t5 = t1 + t2;
/*< t6 = c1 * ( t1 - t2 ) >*/
t6 = c1 * (t1 - t2);
/*< aja = a(ja+j) >*/
aja = a[ja + j];
/*< t7 = aja - 0.25 * t5 >*/
t7 = aja - t5 * (float).25;
/*< a(ja+j) = aja + t5 >*/
a[ja + j] = aja + t5;
/*< t8 = t7 + t6 >*/
t8 = t7 + t6;
/*< t9 = t7 - t6 >*/
t9 = t7 - t6;
/*< t10 = c3 * t3 - c2 * t4 >*/
t10 = c3 * t3 - c2 * t4;
/*< t11 = c2 * t3 + c3 * t4 >*/
t11 = c2 * t3 + c3 * t4;
/*< bjb = b(jb+j) >*/
bjb = b[jb + j];
/*< bje = b(je+j) >*/
bje = b[je + j];
/*< u1 = bjb + bje >*/
u1 = bjb + bje;
/*< bjc = b(jc+j) >*/
bjc = b[jc + j];
/*< bjd = b(jd+j) >*/
bjd = b[jd + j];
/*< u2 = bjc + bjd >*/
u2 = bjc + bjd;
/*< u3 = bjb - bje >*/
u3 = bjb - bje;
/*< u4 = bjc - bjd >*/
u4 = bjc - bjd;
/*< u5 = u1 + u2 >*/
u5 = u1 + u2;
/*< u6 = c1 * ( u1 - u2 ) >*/
u6 = c1 * (u1 - u2);
/*< bja = b(ja+j) >*/
bja = b[ja + j];
/*< u7 = bja - 0.25 * u5 >*/
u7 = bja - u5 * (float).25;
/*< b(ja+j) = bja + u5 >*/
b[ja + j] = bja + u5;
/*< u8 = u7 + u6 >*/
u8 = u7 + u6;
/*< u9 = u7 - u6 >*/
u9 = u7 - u6;
/*< u10 = c3 * u3 - c2 * u4 >*/
u10 = c3 * u3 - c2 * u4;
/*< u11 = c2 * u3 + c3 * u4 >*/
u11 = c2 * u3 + c3 * u4;
/*< a(jb+j) = co1*(t8-u11) - si1*(u8+t11) >*/
a[jb + j] = co1 * (t8 - u11) - si1 * (u8 +
t11);
/*< b(jb+j) = si1*(t8-u11) + co1*(u8+t11) >*/
b[jb + j] = si1 * (t8 - u11) + co1 * (u8 +
t11);
/*< a(je+j) = co4*(t8+u11) - si4*(u8-t11) >*/
a[je + j] = co4 * (t8 + u11) - si4 * (u8 -
t11);
/*< b(je+j) = si4*(t8+u11) + co4*(u8-t11) >*/
b[je + j] = si4 * (t8 + u11) + co4 * (u8 -
t11);
/*< a(jc+j) = co2*(t9-u10) - si2*(u9+t10) >*/
a[jc + j] = co2 * (t9 - u10) - si2 * (u9 +
t10);
/*< b(jc+j) = si2*(t9-u10) + co2*(u9+t10) >*/
b[jc + j] = si2 * (t9 - u10) + co2 * (u9 +
t10);
/*< a(jd+j) = co3*(t9+u10) - si3*(u9-t10) >*/
a[jd + j] = co3 * (t9 + u10) - si3 * (u9 -
t10);
/*< b(jd+j) = si3*(t9+u10) + co3*(u9-t10) >*/
b[jd + j] = si3 * (t9 + u10) + co3 * (u9 -
t10);
/*< j = j + jump >*/
j += *jump;
/*< 130 continue >*/
/* L130: */
}
/*< endif >*/
}
/* -----( end of loop across transforms ) */
/*< ja = ja + jstepx >*/
ja += jstepx;
/*< if (ja.lt.istart) ja = ja + ninc >*/
if (ja < istart) {
ja += ninc;
}
/*< 135 continue >*/
/* L135: */
}
/*< 140 continue >*/
/* L140: */
}
/* -----( end of loop along transforms ) */
/*< kk = kk + 2*la >*/
kk += la << 1;
/*< 150 continue >*/
/* L150: */
}
/* -----( end of loop on nonzero k ) */
/*< la = 5*la >*/
la *= 5;
/*< 160 continue >*/
/* L160: */
}
/* -----( end of loop on type I radix-5 passes) */
/*< if (n.eq.5) go to 490 >*/
if (*n == 5) {
goto L490;
}
/* loop on type II radix-5 passes */
/* ------------------------------ */
/*< 400 continue >*/
/* L400: */
/*< do 480 ipass = mh+1 , m >*/
i__2 = m;
for (ipass = mh + 1; ipass <= i__2; ++ipass) {
/*< jstep = (n*inc) / (5*la) >*/
jstep = *n * *inc / (la * 5);
/*< jstepl = jstep - ninc >*/
jstepl = jstep - ninc;
/*< laincl = la * ink - ninc >*/
laincl = la * ink - ninc;
/*< kk = 0 >*/
kk = 0;
/* loop on k */
/* --------- */
/*< do 470 k = 0 , jstep-ink , ink >*/
i__4 = jstep - ink;
i__3 = ink;
for (k = 0; i__3 < 0 ? k >= i__4 : k <= i__4; k += i__3) {
/*< if (k.gt.0) then >*/
if (k > 0) {
/*< co1 = trigs(kk+1) >*/
co1 = trigs[kk + 1];
/*< si1 = s*trigs(kk+2) >*/
si1 = s * trigs[kk + 2];
/*< co2 = trigs(2*kk+1) >*/
co2 = trigs[(kk << 1) + 1];
/*< si2 = s*trigs(2*kk+2) >*/
si2 = s * trigs[(kk << 1) + 2];
/*< co3 = trigs(3*kk+1) >*/
co3 = trigs[kk * 3 + 1];
/*< si3 = s*trigs(3*kk+2) >*/
si3 = s * trigs[kk * 3 + 2];
/*< co4 = trigs(4*kk+1) >*/
co4 = trigs[(kk << 2) + 1];
/*< si4 = s*trigs(4*kk+2) >*/
si4 = s * trigs[(kk << 2) + 2];
/*< endif >*/
}
/* double loop along first transform in block */
/* ------------------------------------------ */
/*< do 460 ll = k , (la-1)*ink , 5*jstep >*/
i__6 = (la - 1) * ink;
i__5 = jstep * 5;
for (ll = k; i__5 < 0 ? ll >= i__6 : ll <= i__6; ll += i__5) {
/*< do 450 jjj = ll , (n-1)*inc , 5*la*ink >*/
i__7 = (*n - 1) * *inc;
i__8 = la * 5 * ink;
for (jjj = ll; i__8 < 0 ? jjj >= i__7 : jjj <= i__7; jjj
+= i__8) {
/*< ja = istart + jjj >*/
ja = istart + jjj;
/* "transverse" loop */
/* ----------------- */
/*< do 445 nu = 1 , inq >*/
i__9 = inq;
for (nu = 1; nu <= i__9; ++nu) {
/*< jb = ja + jstepl >*/
jb = ja + jstepl;
/*< if (jb.lt.istart) jb = jb + ninc >*/
if (jb < istart) {
jb += ninc;
}
/*< jc = jb + jstepl >*/
jc = jb + jstepl;
/*< if (jc.lt.istart) jc = jc + ninc >*/
if (jc < istart) {
jc += ninc;
}
/*< jd = jc + jstepl >*/
jd = jc + jstepl;
/*< if (jd.lt.istart) jd = jd + ninc >*/
if (jd < istart) {
jd += ninc;
}
/*< je = jd + jstepl >*/
je = jd + jstepl;
/*< if (je.lt.istart) je = je + ninc >*/
if (je < istart) {
je += ninc;
}
/*< jf = ja + laincl >*/
jf = ja + laincl;
/*< if (jf.lt.istart) jf = jf + ninc >*/
if (jf < istart) {
jf += ninc;
}
/*< jg = jf + jstepl >*/
jg = jf + jstepl;
/*< if (jg.lt.istart) jg = jg + ninc >*/
if (jg < istart) {
jg += ninc;
}
/*< jh = jg + jstepl >*/
jh = jg + jstepl;
/*< if (jh.lt.istart) jh = jh + ninc >*/
if (jh < istart) {
jh += ninc;
}
/*< ji = jh + jstepl >*/
ji = jh + jstepl;
/*< if (ji.lt.istart) ji = ji + ninc >*/
if (ji < istart) {
ji += ninc;
}
/*< jj = ji + jstepl >*/
jj = ji + jstepl;
/*< if (jj.lt.istart) jj = jj + ninc >*/
if (jj < istart) {
jj += ninc;
}
/*< jk = jf + laincl >*/
jk = jf + laincl;
/*< if (jk.lt.istart) jk = jk + ninc >*/
if (jk < istart) {
jk += ninc;
}
/*< jl = jk + jstepl >*/
jl = jk + jstepl;
/*< if (jl.lt.istart) jl = jl + ninc >*/
if (jl < istart) {
jl += ninc;
}
/*< jm = jl + jstepl >*/
jm = jl + jstepl;
/*< if (jm.lt.istart) jm = jm + ninc >*/
if (jm < istart) {
jm += ninc;
}
/*< jn = jm + jstepl >*/
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -