📄 gpfa5f.c
字号:
jn = jm + jstepl;
/*< if (jn.lt.istart) jn = jn + ninc >*/
if (jn < istart) {
jn += ninc;
}
/*< jo = jn + jstepl >*/
jo = jn + jstepl;
/*< if (jo.lt.istart) jo = jo + ninc >*/
if (jo < istart) {
jo += ninc;
}
/*< jp = jk + laincl >*/
jp = jk + laincl;
/*< if (jp.lt.istart) jp = jp + ninc >*/
if (jp < istart) {
jp += ninc;
}
/*< jq = jp + jstepl >*/
jq = jp + jstepl;
/*< if (jq.lt.istart) jq = jq + ninc >*/
if (jq < istart) {
jq += ninc;
}
/*< jr = jq + jstepl >*/
jr = jq + jstepl;
/*< if (jr.lt.istart) jr = jr + ninc >*/
if (jr < istart) {
jr += ninc;
}
/*< js = jr + jstepl >*/
js = jr + jstepl;
/*< if (js.lt.istart) js = js + ninc >*/
if (js < istart) {
js += ninc;
}
/*< jt = js + jstepl >*/
jt = js + jstepl;
/*< if (jt.lt.istart) jt = jt + ninc >*/
if (jt < istart) {
jt += ninc;
}
/*< ju = jp + laincl >*/
ju = jp + laincl;
/*< if (ju.lt.istart) ju = ju + ninc >*/
if (ju < istart) {
ju += ninc;
}
/*< jv = ju + jstepl >*/
jv = ju + jstepl;
/*< if (jv.lt.istart) jv = jv + ninc >*/
if (jv < istart) {
jv += ninc;
}
/*< jw = jv + jstepl >*/
jw = jv + jstepl;
/*< if (jw.lt.istart) jw = jw + ninc >*/
if (jw < istart) {
jw += ninc;
}
/*< jx = jw + jstepl >*/
jx = jw + jstepl;
/*< if (jx.lt.istart) jx = jx + ninc >*/
if (jx < istart) {
jx += ninc;
}
/*< jy = jx + jstepl >*/
jy = jx + jstepl;
/*< if (jy.lt.istart) jy = jy + ninc >*/
if (jy < istart) {
jy += ninc;
}
/*< j = 0 >*/
j = 0;
/* loop across transforms */
/* ---------------------- */
/*< if (k.eq.0) then >*/
if (k == 0) {
/* dir$ ivdep, shortloop */
/*< do 410 l = 1 , nvex >*/
i__10 = nvex;
for (l = 1; l <= i__10; ++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;
/*< ajf = a(jf+j) >*/
ajf = a[jf + j];
/*< ajb = ajf >*/
ajb = ajf;
/*< 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;
/*< ajk = a(jk+j) >*/
ajk = a[jk + j];
/*< ajc = ajk >*/
ajc = ajk;
/*< 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;
/*< bjf = b(jf+j) >*/
bjf = b[jf + j];
/*< bjb = bjf >*/
bjb = bjf;
/*< 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;
/*< bjk = b(jk+j) >*/
bjk = b[jk + j];
/*< bjc = bjk >*/
bjc = bjk;
/*< u10 = c3 * u3 - c2 * u4 >*/
u10 = c3 * u3 - c2 * u4;
/*< u11 = c2 * u3 + c3 * u4 >*/
u11 = c2 * u3 + c3 * u4;
/*< a(jf+j) = t8 - u11 >*/
a[jf + j] = t8 - u11;
/*< b(jf+j) = u8 + t11 >*/
b[jf + j] = u8 + t11;
/*< aje = t8 + u11 >*/
aje = t8 + u11;
/*< bje = u8 - t11 >*/
bje = u8 - t11;
/*< a(jk+j) = t9 - u10 >*/
a[jk + j] = t9 - u10;
/*< b(jk+j) = u9 + t10 >*/
b[jk + j] = u9 + t10;
/*< ajd = t9 + u10 >*/
ajd = t9 + u10;
/*< bjd = u9 - t10 >*/
bjd = u9 - t10;
/* ---------------------- */
/*< ajg = a(jg+j) >*/
ajg = a[jg + j];
/*< ajj = a(jj+j) >*/
ajj = a[jj + j];
/*< t1 = ajg + ajj >*/
t1 = ajg + ajj;
/*< ajh = a(jh+j) >*/
ajh = a[jh + j];
/*< aji = a(ji+j) >*/
aji = a[ji + j];
/*< t2 = ajh + aji >*/
t2 = ajh + aji;
/*< t3 = ajg - ajj >*/
t3 = ajg - ajj;
/*< t4 = ajh - aji >*/
t4 = ajh - aji;
/*< ajl = a(jl+j) >*/
ajl = a[jl + j];
/*< ajh = ajl >*/
ajh = ajl;
/*< t5 = t1 + t2 >*/
t5 = t1 + t2;
/*< t6 = c1 * ( t1 - t2 ) >*/
t6 = c1 * (t1 - t2);
/*< t7 = ajb - 0.25 * t5 >*/
t7 = ajb - t5 * (float).25;
/*< a(jb+j) = ajb + t5 >*/
a[jb + j] = ajb + t5;
/*< t8 = t7 + t6 >*/
t8 = t7 + t6;
/*< t9 = t7 - t6 >*/
t9 = t7 - t6;
/*< ajq = a(jq+j) >*/
ajq = a[jq + j];
/*< aji = ajq >*/
aji = ajq;
/*< t10 = c3 * t3 - c2 * t4 >*/
t10 = c3 * t3 - c2 * t4;
/*< t11 = c2 * t3 + c3 * t4 >*/
t11 = c2 * t3 + c3 * t4;
/*< bjg = b(jg+j) >*/
bjg = b[jg + j];
/*< bjj = b(jj+j) >*/
bjj = b[jj + j];
/*< u1 = bjg + bjj >*/
u1 = bjg + bjj;
/*< bjh = b(jh+j) >*/
bjh = b[jh + j];
/*< bji = b(ji+j) >*/
bji = b[ji + j];
/*< u2 = bjh + bji >*/
u2 = bjh + bji;
/*< u3 = bjg - bjj >*/
u3 = bjg - bjj;
/*< u4 = bjh - bji >*/
u4 = bjh - bji;
/*< bjl = b(jl+j) >*/
bjl = b[jl + j];
/*< bjh = bjl >*/
bjh = bjl;
/*< u5 = u1 + u2 >*/
u5 = u1 + u2;
/*< u6 = c1 * ( u1 - u2 ) >*/
u6 = c1 * (u1 - u2);
/*< u7 = bjb - 0.25 * u5 >*/
u7 = bjb - u5 * (float).25;
/*< b(jb+j) = bjb + u5 >*/
b[jb + j] = bjb + u5;
/*< u8 = u7 + u6 >*/
u8 = u7 + u6;
/*< u9 = u7 - u6 >*/
u9 = u7 - u6;
/*< bjq = b(jq+j) >*/
bjq = b[jq + j];
/*< bji = bjq >*/
bji = bjq;
/*< u10 = c3 * u3 - c2 * u4 >*/
u10 = c3 * u3 - c2 * u4;
/*< u11 = c2 * u3 + c3 * u4 >*/
u11 = c2 * u3 + c3 * u4;
/*< a(jg+j) = t8 - u11 >*/
a[jg + j] = t8 - u11;
/*< b(jg+j) = u8 + t11 >*/
b[jg + j] = u8 + t11;
/*< ajj = t8 + u11 >*/
ajj = t8 + u11;
/*< bjj = u8 - t11 >*/
bjj = u8 - t11;
/*< a(jl+j) = t9 - u10 >*/
a[jl + j] = t9 - u10;
/*< b(jl+j) = u9 + t10 >*/
b[jl + j] = u9 + t10;
/*< a(jq+j) = t9 + u10 >*/
a[jq + j] = t9 + u10;
/*< b(jq+j) = u9 - t10 >*/
b[jq + j] = u9 - t10;
/* ---------------------- */
/*< ajo = a(jo+j) >*/
ajo = a[jo + j];
/*< t1 = ajh + ajo >*/
t1 = ajh + ajo;
/*< ajm = a(jm+j) >*/
ajm = a[jm + j];
/*< ajn = a(jn+j) >*/
ajn = a[jn + j];
/*< t2 = ajm + ajn >*/
t2 = ajm + ajn;
/*< t3 = ajh - ajo >*/
t3 = ajh - ajo;
/*< t4 = ajm - ajn >*/
t4 = ajm - ajn;
/*< ajr = a(jr+j) >*/
ajr = a[jr + j];
/*< ajn = ajr >*/
ajn = ajr;
/*< t5 = t1 + t2 >*/
t5 = t1 + t2;
/*< t6 = c1 * ( t1 - t2 ) >*/
t6 = c1 * (t1 - t2);
/*< t7 = ajc - 0.25 * t5 >*/
t7 = ajc - t5 * (float).25;
/*< a(jc+j) = ajc + t5 >*/
a[jc + j] = ajc + t5;
/*< t8 = t7 + t6 >*/
t8 = t7 + t6;
/*< t9 = t7 - t6 >*/
t9 = t7 - t6;
/*< ajw = a(jw+j) >*/
ajw = a[jw + j];
/*< ajo = ajw >*/
ajo = ajw;
/*< t10 = c3 * t3 - c2 * t4 >*/
t10 = c3 * t3 - c2 * t4;
/*< t11 = c2 * t3 + c3 * t4 >*/
t11 = c2 * t3 + c3 * t4;
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -