Mercurial > octave-nkf
diff liboctave/CMatrix.cc @ 4153:6b96ce9f5743
[project @ 2002-11-06 20:38:49 by jwe]
author | jwe |
---|---|
date | Wed, 06 Nov 2002 20:38:50 +0000 |
parents | 7d9bda865012 |
children | 5719210fff4c |
line wrap: on
line diff
--- a/liboctave/CMatrix.cc Wed Nov 06 18:41:50 2002 +0000 +++ b/liboctave/CMatrix.cc Wed Nov 06 20:38:50 2002 +0000 @@ -986,6 +986,8 @@ for (size_t i = 0; i < nsamples; i++) { + OCTAVE_QUIT; + octave_fftw::fft (&in[npts * i], &out[npts * i], npts); } @@ -1018,6 +1020,8 @@ for (size_t i = 0; i < nsamples; i++) { + OCTAVE_QUIT; + octave_fftw::ifft (&in[npts * i], &out[npts * i], npts); } @@ -1086,7 +1090,11 @@ F77_FUNC (cffti, CFFTI) (npts, pwsave); for (int j = 0; j < nsamples; j++) - F77_FUNC (cfftf, CFFTF) (npts, &tmp_data[npts*j], pwsave); + { + OCTAVE_QUIT; + + F77_FUNC (cfftf, CFFTF) (npts, &tmp_data[npts*j], pwsave); + } return retval; } @@ -1123,7 +1131,11 @@ F77_FUNC (cffti, CFFTI) (npts, pwsave); for (int j = 0; j < nsamples; j++) - F77_FUNC (cfftb, CFFTB) (npts, &tmp_data[npts*j], pwsave); + { + OCTAVE_QUIT; + + F77_FUNC (cfftb, CFFTB) (npts, &tmp_data[npts*j], pwsave); + } for (int j = 0; j < npts*nsamples; j++) tmp_data[j] = tmp_data[j] / static_cast<double> (npts); @@ -1163,7 +1175,11 @@ F77_FUNC (cffti, CFFTI) (npts, pwsave); for (int j = 0; j < nsamples; j++) - F77_FUNC (cfftf, CFFTF) (npts, &tmp_data[npts*j], pwsave); + { + OCTAVE_QUIT; + + F77_FUNC (cfftf, CFFTF) (npts, &tmp_data[npts*j], pwsave); + } npts = nc; nsamples = nr; @@ -1179,6 +1195,8 @@ for (int j = 0; j < nsamples; j++) { + OCTAVE_QUIT; + for (int i = 0; i < npts; i++) prow[i] = tmp_data[i*nr + j]; @@ -1223,7 +1241,11 @@ F77_FUNC (cffti, CFFTI) (npts, pwsave); for (int j = 0; j < nsamples; j++) - F77_FUNC (cfftb, CFFTB) (npts, &tmp_data[npts*j], pwsave); + { + OCTAVE_QUIT; + + F77_FUNC (cfftb, CFFTB) (npts, &tmp_data[npts*j], pwsave); + } for (int j = 0; j < npts*nsamples; j++) tmp_data[j] = tmp_data[j] / static_cast<double> (npts); @@ -1242,6 +1264,8 @@ for (int j = 0; j < nsamples; j++) { + OCTAVE_QUIT; + for (int i = 0; i < npts; i++) prow[i] = tmp_data[i*nr + j]; @@ -1937,6 +1961,8 @@ for (int j = 0; j < nc; j++) retval(i,j) *= dscale(i) / dscale(j); + OCTAVE_QUIT; + // construct balancing permutation vector Array<int> ipermute (nc); for (int i = 0; i < nc; i++) @@ -1965,6 +1991,8 @@ for (int i = 0; i < nc; i++) invpvec(ipermute(i)) = i; // Thanks to R. A. Lippert for this method + OCTAVE_QUIT; + ComplexMatrix tmpMat = retval; for (int i = 0; i < nc; i++) for (int j = 0; j < nc; j++)