📄 fft.cpp
字号:
ch[*ido + ((k << 2) + 2) * ch_dim1] = cc[(k + cc_dim2) * cc_dim1 + 1]
- cc[(k + cc_dim2 * 3) * cc_dim1 + 1];
ch[((k << 2) + 3) * ch_dim1 + 1] = cc[(k + (cc_dim2 << 2)) * cc_dim1
+ 1] - cc[(k + (cc_dim2 << 1)) * cc_dim1 + 1];
/* L101: */
}
if ((i__1 = *ido - 2) < 0) {
goto L107;
} else if (i__1 == 0) {
goto L105;
} else {
goto L102;
}
L102:
idp2 = *ido + 2;
i__1 = *l1;
for (k = 1; k <= i__1; ++k) {
i__2 = *ido;
for (i__ = 3; i__ <= i__2; i__ += 2) {
ic = idp2 - i__;
cr2 = wa1[i__ - 2] * cc[i__ - 1 + (k + (cc_dim2 << 1)) * cc_dim1]
+ wa1[i__ - 1] * cc[i__ + (k + (cc_dim2 << 1)) * cc_dim1];
ci2 = wa1[i__ - 2] * cc[i__ + (k + (cc_dim2 << 1)) * cc_dim1] -
wa1[i__ - 1] * cc[i__ - 1 + (k + (cc_dim2 << 1)) *
cc_dim1];
cr3 = wa2[i__ - 2] * cc[i__ - 1 + (k + cc_dim2 * 3) * cc_dim1] +
wa2[i__ - 1] * cc[i__ + (k + cc_dim2 * 3) * cc_dim1];
ci3 = wa2[i__ - 2] * cc[i__ + (k + cc_dim2 * 3) * cc_dim1] - wa2[
i__ - 1] * cc[i__ - 1 + (k + cc_dim2 * 3) * cc_dim1];
cr4 = wa3[i__ - 2] * cc[i__ - 1 + (k + (cc_dim2 << 2)) * cc_dim1]
+ wa3[i__ - 1] * cc[i__ + (k + (cc_dim2 << 2)) * cc_dim1];
ci4 = wa3[i__ - 2] * cc[i__ + (k + (cc_dim2 << 2)) * cc_dim1] -
wa3[i__ - 1] * cc[i__ - 1 + (k + (cc_dim2 << 2)) *
cc_dim1];
tr1 = cr2 + cr4;
tr4 = cr4 - cr2;
ti1 = ci2 + ci4;
ti4 = ci2 - ci4;
ti2 = cc[i__ + (k + cc_dim2) * cc_dim1] + ci3;
ti3 = cc[i__ + (k + cc_dim2) * cc_dim1] - ci3;
tr2 = cc[i__ - 1 + (k + cc_dim2) * cc_dim1] + cr3;
tr3 = cc[i__ - 1 + (k + cc_dim2) * cc_dim1] - cr3;
ch[i__ - 1 + ((k << 2) + 1) * ch_dim1] = tr1 + tr2;
ch[ic - 1 + ((k << 2) + 4) * ch_dim1] = tr2 - tr1;
ch[i__ + ((k << 2) + 1) * ch_dim1] = ti1 + ti2;
ch[ic + ((k << 2) + 4) * ch_dim1] = ti1 - ti2;
ch[i__ - 1 + ((k << 2) + 3) * ch_dim1] = ti4 + tr3;
ch[ic - 1 + ((k << 2) + 2) * ch_dim1] = tr3 - ti4;
ch[i__ + ((k << 2) + 3) * ch_dim1] = tr4 + ti3;
ch[ic + ((k << 2) + 2) * ch_dim1] = tr4 - ti3;
/* L103: */
}
/* L104: */
}
if (*ido % 2 == 1) {
return 0;
}
L105:
i__1 = *l1;
for (k = 1; k <= i__1; ++k) {
ti1 = -hsqt2 * (cc[*ido + (k + (cc_dim2 << 1)) * cc_dim1] + cc[*ido +
(k + (cc_dim2 << 2)) * cc_dim1]);
tr1 = hsqt2 * (cc[*ido + (k + (cc_dim2 << 1)) * cc_dim1] - cc[*ido + (
k + (cc_dim2 << 2)) * cc_dim1]);
ch[*ido + ((k << 2) + 1) * ch_dim1] = tr1 + cc[*ido + (k + cc_dim2) *
cc_dim1];
ch[*ido + ((k << 2) + 3) * ch_dim1] = cc[*ido + (k + cc_dim2) *
cc_dim1] - tr1;
ch[((k << 2) + 2) * ch_dim1 + 1] = ti1 - cc[*ido + (k + cc_dim2 * 3) *
cc_dim1];
ch[((k << 2) + 4) * ch_dim1 + 1] = ti1 + cc[*ido + (k + cc_dim2 * 3) *
cc_dim1];
/* L106: */
}
L107:
return 0;
}
int FFT::radf5_(int *ido, int *l1, double *cc, double *ch, double *wa1, double *wa2, double *wa3, double *wa4)
{
/* Initialized data */
static double tr11 = .309016994374947f;
static double ti11 = .951056516295154f;
static double tr12 = -.809016994374947f;
static double ti12 = .587785252292473f;
/* System generated locals */
int cc_dim1, cc_dim2, cc_offset, ch_dim1, ch_offset, i__1, i__2;
/* Local variables */
static int i__, k, ic;
static double ci2, di2, ci4, ci5, di3, di4, di5, ci3, cr2, cr3, dr2, dr3,
dr4, dr5, cr5, cr4, ti2, ti3, ti5, ti4, tr2, tr3, tr4, tr5;
static int idp2;
/* Parameter adjustments */
ch_dim1 = *ido;
ch_offset = 1 + ch_dim1 * 6;
ch -= ch_offset;
cc_dim1 = *ido;
cc_dim2 = *l1;
cc_offset = 1 + cc_dim1 * (1 + cc_dim2);
cc -= cc_offset;
--wa1;
--wa2;
--wa3;
--wa4;
/* Function Body */
i__1 = *l1;
for (k = 1; k <= i__1; ++k) {
cr2 = cc[(k + cc_dim2 * 5) * cc_dim1 + 1] + cc[(k + (cc_dim2 << 1)) *
cc_dim1 + 1];
ci5 = cc[(k + cc_dim2 * 5) * cc_dim1 + 1] - cc[(k + (cc_dim2 << 1)) *
cc_dim1 + 1];
cr3 = cc[(k + (cc_dim2 << 2)) * cc_dim1 + 1] + cc[(k + cc_dim2 * 3) *
cc_dim1 + 1];
ci4 = cc[(k + (cc_dim2 << 2)) * cc_dim1 + 1] - cc[(k + cc_dim2 * 3) *
cc_dim1 + 1];
ch[(k * 5 + 1) * ch_dim1 + 1] = cc[(k + cc_dim2) * cc_dim1 + 1] + cr2
+ cr3;
ch[*ido + (k * 5 + 2) * ch_dim1] = cc[(k + cc_dim2) * cc_dim1 + 1] +
tr11 * cr2 + tr12 * cr3;
ch[(k * 5 + 3) * ch_dim1 + 1] = ti11 * ci5 + ti12 * ci4;
ch[*ido + (k * 5 + 4) * ch_dim1] = cc[(k + cc_dim2) * cc_dim1 + 1] +
tr12 * cr2 + tr11 * cr3;
ch[(k * 5 + 5) * ch_dim1 + 1] = ti12 * ci5 - ti11 * ci4;
/* L101: */
}
if (*ido == 1) {
return 0;
}
idp2 = *ido + 2;
i__1 = *l1;
for (k = 1; k <= i__1; ++k) {
i__2 = *ido;
for (i__ = 3; i__ <= i__2; i__ += 2) {
ic = idp2 - i__;
dr2 = wa1[i__ - 2] * cc[i__ - 1 + (k + (cc_dim2 << 1)) * cc_dim1]
+ wa1[i__ - 1] * cc[i__ + (k + (cc_dim2 << 1)) * cc_dim1];
di2 = wa1[i__ - 2] * cc[i__ + (k + (cc_dim2 << 1)) * cc_dim1] -
wa1[i__ - 1] * cc[i__ - 1 + (k + (cc_dim2 << 1)) *
cc_dim1];
dr3 = wa2[i__ - 2] * cc[i__ - 1 + (k + cc_dim2 * 3) * cc_dim1] +
wa2[i__ - 1] * cc[i__ + (k + cc_dim2 * 3) * cc_dim1];
di3 = wa2[i__ - 2] * cc[i__ + (k + cc_dim2 * 3) * cc_dim1] - wa2[
i__ - 1] * cc[i__ - 1 + (k + cc_dim2 * 3) * cc_dim1];
dr4 = wa3[i__ - 2] * cc[i__ - 1 + (k + (cc_dim2 << 2)) * cc_dim1]
+ wa3[i__ - 1] * cc[i__ + (k + (cc_dim2 << 2)) * cc_dim1];
di4 = wa3[i__ - 2] * cc[i__ + (k + (cc_dim2 << 2)) * cc_dim1] -
wa3[i__ - 1] * cc[i__ - 1 + (k + (cc_dim2 << 2)) *
cc_dim1];
dr5 = wa4[i__ - 2] * cc[i__ - 1 + (k + cc_dim2 * 5) * cc_dim1] +
wa4[i__ - 1] * cc[i__ + (k + cc_dim2 * 5) * cc_dim1];
di5 = wa4[i__ - 2] * cc[i__ + (k + cc_dim2 * 5) * cc_dim1] - wa4[
i__ - 1] * cc[i__ - 1 + (k + cc_dim2 * 5) * cc_dim1];
cr2 = dr2 + dr5;
ci5 = dr5 - dr2;
cr5 = di2 - di5;
ci2 = di2 + di5;
cr3 = dr3 + dr4;
ci4 = dr4 - dr3;
cr4 = di3 - di4;
ci3 = di3 + di4;
ch[i__ - 1 + (k * 5 + 1) * ch_dim1] = cc[i__ - 1 + (k + cc_dim2) *
cc_dim1] + cr2 + cr3;
ch[i__ + (k * 5 + 1) * ch_dim1] = cc[i__ + (k + cc_dim2) *
cc_dim1] + ci2 + ci3;
tr2 = cc[i__ - 1 + (k + cc_dim2) * cc_dim1] + tr11 * cr2 + tr12 *
cr3;
ti2 = cc[i__ + (k + cc_dim2) * cc_dim1] + tr11 * ci2 + tr12 * ci3;
tr3 = cc[i__ - 1 + (k + cc_dim2) * cc_dim1] + tr12 * cr2 + tr11 *
cr3;
ti3 = cc[i__ + (k + cc_dim2) * cc_dim1] + tr12 * ci2 + tr11 * ci3;
tr5 = ti11 * cr5 + ti12 * cr4;
ti5 = ti11 * ci5 + ti12 * ci4;
tr4 = ti12 * cr5 - ti11 * cr4;
ti4 = ti12 * ci5 - ti11 * ci4;
ch[i__ - 1 + (k * 5 + 3) * ch_dim1] = tr2 + tr5;
ch[ic - 1 + (k * 5 + 2) * ch_dim1] = tr2 - tr5;
ch[i__ + (k * 5 + 3) * ch_dim1] = ti2 + ti5;
ch[ic + (k * 5 + 2) * ch_dim1] = ti5 - ti2;
ch[i__ - 1 + (k * 5 + 5) * ch_dim1] = tr3 + tr4;
ch[ic - 1 + (k * 5 + 4) * ch_dim1] = tr3 - tr4;
ch[i__ + (k * 5 + 5) * ch_dim1] = ti3 + ti4;
ch[ic + (k * 5 + 4) * ch_dim1] = ti4 - ti3;
/* L102: */
}
/* L103: */
}
return 0;
}
int FFT::radfg_(int *ido, int *ip, int *l1, int *idl1, double *cc, double *c1, double *c2, double *ch, double *ch2, double *wa)
{
/* System generated locals */
int ch_dim1, ch_dim2, ch_offset, cc_dim1, cc_dim2, cc_offset, c1_dim1,
c1_dim2, c1_offset, c2_dim1, c2_offset, ch2_dim1, ch2_offset,
i__1, i__2, i__3;
/* Builtin functions */
// double cos(doublereal), sin(doublereal);
/* Local variables */
static int i__, j, k, l, j2, ic, jc, lc, ik, is;
static double dc2, ai1, ai2, ar1, ar2, ds2;
static int nbd;
static double dcp, arg, dum, dsp, tpi, ar1h, ar2h;
static int idp2, ipp2, idij, ipph;
// extern double pimach_(real *);
/* Parameter adjustments */
ch_dim1 = *ido;
ch_dim2 = *l1;
ch_offset = 1 + ch_dim1 * (1 + ch_dim2);
ch -= ch_offset;
c1_dim1 = *ido;
c1_dim2 = *l1;
c1_offset = 1 + c1_dim1 * (1 + c1_dim2);
c1 -= c1_offset;
cc_dim1 = *ido;
cc_dim2 = *ip;
cc_offset = 1 + cc_dim1 * (1 + cc_dim2);
cc -= cc_offset;
ch2_dim1 = *idl1;
ch2_offset = 1 + ch2_dim1;
ch2 -= ch2_offset;
c2_dim1 = *idl1;
c2_offset = 1 + c2_dim1;
c2 -= c2_offset;
--wa;
/* Function Body */
tpi = pimach_(&dum) * 2.f;
arg = tpi / (double) (*ip);
dcp = cos(arg);
dsp = sin(arg);
ipph = (*ip + 1) / 2;
ipp2 = *ip + 2;
idp2 = *ido + 2;
nbd = (*ido - 1) / 2;
if (*ido == 1) {
goto L119;
}
i__1 = *idl1;
for (ik = 1; ik <= i__1; ++ik) {
ch2[ik + ch2_dim1] = c2[ik + c2_dim1];
/* L101: */
}
i__1 = *ip;
for (j = 2; j <= i__1; ++j) {
i__2 = *l1;
for (k = 1; k <= i__2; ++k) {
ch[(k + j * ch_dim2) * ch_dim1 + 1] = c1[(k + j * c1_dim2) *
c1_dim1 + 1];
/* L102: */
}
/* L103: */
}
if (nbd > *l1) {
goto L107;
}
is = -(*ido);
i__1 = *ip;
for (j = 2; j <= i__1; ++j) {
is += *ido;
idij = is;
i__2 = *ido;
for (i__ = 3; i__ <= i__2; i__ += 2) {
idij += 2;
i__3 = *l1;
for (k = 1; k <= i__3; ++k) {
ch[i__ - 1 + (k + j * ch_dim2) * ch_dim1] = wa[idij - 1] * c1[
i__ - 1 + (k + j * c1_dim2) * c1_dim1] + wa[idij] *
c1[i__ + (k + j * c1_dim2) * c1_dim1];
ch[i__ + (k + j * ch_dim2) * ch_dim1] = wa[idij - 1] * c1[i__
+ (k + j * c1_dim2) * c1_dim1] - wa[idij] * c1[i__ -
1 + (k + j * c1_dim2) * c1_dim1];
/* L104: */
}
/* L105: */
}
/* L106: */
}
goto L111;
L107:
is = -(*ido);
i__1 = *ip;
for (j = 2; j <= i__1; ++j) {
is += *ido;
i__2 = *l1;
for (k = 1; k <= i__2; ++k) {
idij = is;
i__3 = *ido;
for (i__ = 3; i__ <= i__3; i__ += 2) {
idij += 2;
ch[i__ - 1 + (k + j * ch_dim2) * ch_dim1] = wa[idij - 1] * c1[
i__ - 1 + (k + j * c1_dim2) * c1_dim1] + wa[idij] *
c1[i__ + (k + j * c1_dim2) * c1_dim1];
ch[i__ + (k + j * ch_dim2) * ch_dim1] = wa[idij - 1] * c1[i__
+ (k + j * c1_dim2) * c1_dim1] - wa[idij] * c1[i__ -
1 + (k + j * c1_dim2) * c1_dim1];
/* L108: */
}
/* L109: */
}
/* L110: */
}
L111:
if (nbd < *l1) {
goto L115;
}
i__1 = ipph;
for (j = 2; j <= i__1; ++j) {
jc = ipp2 - j;
i__2 = *l1;
for (k = 1; k <= i__2; ++k) {
i__3 = *ido;
for (i__ = 3; i__ <= i__3; i__ += 2) {
c1[i__ - 1 + (k + j * c1_dim2) * c1_dim1] = ch[i__ - 1 + (k +
j * ch_dim2) * ch_dim1] + ch[i__ - 1 + (k + jc *
ch_dim2) * ch_dim1];
c1[i__ - 1 + (k + jc * c1_dim2) * c1_dim1] = ch[i__ + (k + j *
ch_dim2) * ch_dim1] - ch[i__ + (k + jc * ch_dim2) *
ch_dim1];
c1[i__ + (k + j * c1_dim2) * c1_dim1] = ch[i__ + (k + j *
ch_dim2) * ch_dim1] + ch[i__ + (k + jc * ch_dim2) *
ch_dim1];
c1[i__ + (k + jc * c1_dim2) * c1_dim1] = ch[i__ - 1 + (k + jc
* ch_dim2) * ch_dim1] - ch[i__ - 1 + (k + j * ch_dim2)
* ch_dim1];
/* L112: */
}
/* L113: */
}
/* L114: */
}
goto L121;
L115:
i__1 = ipph;
for (j = 2; j <= i__1; ++j) {
jc = ipp2 - j;
i__2 = *ido;
for (i__ = 3; i__ <= i__2; i__ += 2) {
i__3 = *l1;
for (k = 1; k <= i__3; ++k) {
c1[i__ - 1 + (k + j * c1_dim2) * c1_dim1] = ch[i__ - 1 + (k +
j * ch_dim2) * ch_dim1] + ch[i__ - 1 + (k + jc *
ch_dim2) * ch_dim1];
c1[i__ - 1 + (k + jc * c1_dim2) * c1_dim1] = ch[i__ + (k + j *
ch_dim2) * ch_dim1] - ch[i__ + (k + jc * ch_dim2) *
ch_dim1];
c1[i__ + (k + j * c1_dim2) * c1_dim1] = ch[i__ + (k + j *
ch_dim2) * ch_dim1] + ch[i__ + (k + jc * ch_dim2) *
ch_dim1];
c1[i__ + (k + jc * c1_dim2) * c1_dim1] = ch[i__ - 1 + (k + jc
* ch_dim2) * ch_dim1] - ch[i__ - 1 + (k + j * ch_dim2)
* ch_dim1];
/* L116: */
}
/* L117: */
}
/* L118: */
}
goto L121;
L119:
i__1 = *idl1;
for (ik = 1; ik <= i__1; ++ik) {
c2[ik + c2_dim1] = ch2[ik + ch2_dim1];
/* L120: */
}
L121:
i__1 = ipph;
for (j = 2; j <= i__1; ++j) {
jc = ipp2 - j;
i__2 = *l1;
for (k = 1; k <= i__2; ++k) {
c1[(k + j * c1_dim2) * c1_dim1 + 1] = ch[(k + j * ch_dim2) *
ch_dim1 + 1] + ch[(k + jc * ch_dim2) * ch_dim1 + 1];
c1[(k + jc * c1_dim2) * c1_dim1 + 1] = ch[(k + jc * ch_dim2) *
ch_dim1 + 1] - ch[(k + j * ch_dim2) * ch_dim1 + 1];
/* L122: */
}
/* L123: */
}
ar1 = 1.f;
ai1 = 0.f;
i__1 = ipph;
for (l = 2; l <= i__1; ++l) {
lc = ipp2 - l;
ar1h = dcp * ar1 - dsp * ai1;
ai1 = dcp * ai1 + dsp * ar1;
ar1 = ar1h;
i__2 = *idl1;
for (ik = 1; ik <= i__2; ++ik) {
ch2[ik + l * ch2_dim1] = c2[ik + c2_dim1] + ar1 * c2[ik + (
c2_dim1 << 1)];
ch2[ik + lc * ch2_dim1] = ai1 * c2[ik + *ip * c2_dim1];
/* L124: */
}
dc2 = ar1;
ds2 = ai1;
ar2 = ar1;
ai2 = ai1;
i__2 = ipph;
for (j = 3; j <= i__2; ++j) {
jc = ipp2 - j;
ar2h = dc2 * ar2 - ds2 * ai2;
ai2 = dc2 * ai2 + ds2 * ar2;
ar2 = ar2h;
i__3 = *idl1;
for (ik = 1; ik <= i__3; ++ik) {
ch2[ik + l * ch2_dim1] += ar2 * c2[ik + j * c2_dim1];
ch2[ik + lc * ch2_dim1] += ai2 * c2[ik + jc * c2_dim1];
/* L125: */
}
/* L126: */
}
/* L127: */
}
i__1 = ipph;
for (j = 2; j <= i__1; ++j) {
i__2 = *idl1;
for (ik = 1; ik <= i__2; ++ik) {
ch2[ik + ch2_dim1] += c2[ik + j * c2_dim1];
/* L128: */
}
/* L129: */
}
if (*ido < *l1) {
goto L132;
}
i__1 = *l1;
for (k = 1; k <= i__1; ++k) {
i__2 = *ido;
for (i__ = 1; i__ <= i__2; ++i__) {
cc[i__ + (k * cc_dim2 + 1) * cc_dim1] = ch[i__ + (k + ch_dim2) *
ch_dim1];
/* L130: */
}
/* L131: */
}
goto L135;
L132:
i__1 = *ido;
for (i__ = 1; i__ <= i__1; ++i__) {
i__2 = *l1;
for (k = 1; k <= i__2; ++k) {
cc[i__ + (k * cc_dim2 + 1) * cc_dim1] = ch[i__ + (k + ch_dim2) *
ch_dim1];
/* L133: */
}
/* L134: */
}
L135:
i__1 = ipph;
for (j = 2; j <= i__1; ++j) {
jc = ipp2 - j;
j2 = j + j;
i__2 = *l1;
for (k = 1; k <= i__2; ++k) {
cc[*ido + (j2 - 2 + k * cc_dim2) * cc_dim1] = ch[(k + j * ch_dim2)
* ch_dim1 + 1];
cc[(j2 - 1 + k * cc_dim2) * cc_dim1 + 1] = ch[(k + jc * ch_dim2) *
ch_dim1 + 1];
/* L136: */
}
/* L137: */
}
if (*ido == 1) {
return 0;
}
if (nbd < *l1) {
goto L141;
}
i__1 = ipph;
for (j = 2; j <= i__1; ++j) {
jc = ipp2 - j;
j2 = j + j;
i__2 = *l1;
for (k = 1; k <= i__2; ++k) {
i__3 = *ido;
for (i__ = 3; i__ <= i__3; i__ += 2) {
ic = idp2 - i__;
cc[i__ - 1 + (j2 - 1 + k * cc_dim2) * cc_dim1] = ch[i__ - 1 +
(k + j * ch_dim2) * ch_dim1] + ch[i__ - 1 + (k + jc *
ch_dim2) * ch_dim1];
cc[ic - 1 + (j2 - 2 + k * cc_dim2) * cc_dim1] = ch[i__ - 1 + (
k + j * ch_dim2) * ch_dim1] - ch[i__ - 1 + (k + jc *
ch_dim2) * ch_dim1];
cc[i__ + (j2 - 1 + k * cc_dim2) * cc_dim1] = ch[i__ + (k + j *
ch_dim2) * ch_dim1] + ch[i__ + (k + jc * ch_dim2) *
ch_dim1];
cc[ic + (j2 - 2 + k * cc_dim2) * cc_dim1] = ch[i__ + (k + jc *
ch_dim2) * ch_dim1] - ch[i__ + (k + j * ch_dim2) *
ch_dim1];
/* L138: */
}
/* L139: */
}
/* L140: */
}
return 0;
L141:
i__1 = ipph;
for (j = 2; j <= i__1; ++j) {
jc = ipp2 - j;
j2 = j + j;
i__2 = *ido;
for (i__ = 3; i__ <= i__2; i__ += 2) {
ic = idp2 - i__;
i__3 = *l1;
for (k = 1; k <= i__3; ++k) {
cc[i__ - 1 + (j2 - 1 + k * cc_dim2) * cc_dim1] = ch[i__ - 1 +
(k + j * ch_dim2) * ch_dim1] + ch[i__ - 1 + (k + jc *
ch_dim2) * ch_dim1];
cc[ic - 1 + (j2 - 2 + k * cc_dim2) * cc_dim1] = ch[i__ - 1 + (
k + j * ch_dim2) * ch_dim1] - ch[i__ - 1 + (k + jc *
ch_dim2) * ch_dim1];
cc[i__ + (j2 - 1 + k * cc_dim2) * cc_dim1] = ch[i__ + (k + j *
ch_dim2) * ch_dim1] + ch[i__ + (k + jc * ch_dim2) *
ch_dim1];
cc[ic + (j2 - 2 + k * cc_dim2) * cc_dim1] = ch[i__ + (k + jc *
ch_dim2) * ch_dim1] - ch[i__ + (k + j * ch_dim2) *
ch_dim1];
/* L142: */
}
/* L143: */
}
/* L144: */
}
return 0;
}
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -