Changeset 515 in Sophya for trunk/SophyaLib/NTools/fftserver.cc


Ignore:
Timestamp:
Oct 26, 1999, 2:44:44 PM (26 years ago)
Author:
ansari
Message:

Portage SGI-CC Reza 26/10/99

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/SophyaLib/NTools/fftserver.cc

    r514 r515  
    9292  checkint_cfft(l);
    9393  float* foo = new float[2*l];
    94   for (int i=0;i<l;i++){
     94  int i;
     95  for (i=0;i<l;i++){
    9596    foo[2*i]=inout[i].real();
    9697    foo[2*i+1]=inout[i].imag();
     
    9899  cfftf_(&l, foo, ws_cfft);
    99100  inout[0]=complex<float> (foo[0],foo[1]);
    100   for (int i=1;i<l;i++) inout[l-i]= complex<float> (foo[2*i], foo[2*i+1]);
     101  for (i=1;i<l;i++) inout[l-i]= complex<float> (foo[2*i], foo[2*i+1]);
    101102  delete[] foo;
    102103}
     
    106107  checkint_cdfft(l);
    107108  double* foo=new double[2*l];
    108   for (int i=0;i<l;i++){
     109  int i;
     110  for (i=0;i<l;i++){
    109111    foo[2*i]=inout[i].real();
    110112    foo[2*i+1]=inout[i].imag();
     
    112114  cdfftf_(&l, foo, ws_cdfft);
    113115  inout[0]=complex<double> (foo[0],foo[1]);
    114   for (int i=1;i<l;i++) {
     116  for (i=1;i<l;i++) {
    115117    inout[l-i]= complex<double> (foo[2*i],foo[2*i+1]);
    116118  }
     
    134136  checkint_cfft(l);
    135137  float* foo = new float[2*l];
    136   for (int i=0;i<l;i++){
     138  int i;
     139  for (i=0;i<l;i++){
    137140    foo[2*i]=inout[i].real();
    138141    foo[2*i+1]=inout[i].imag();
    139142  }
    140143  cfftf_(&l, foo, ws_cfft);
    141   for (int i=0;i<l;i++) inout[i]=complex<float> (foo[2*i],foo[2*i+1]);
     144  for (i=0;i<l;i++) inout[i]=complex<float> (foo[2*i],foo[2*i+1]);
    142145  delete[] foo;
    143146}
     
    147150  checkint_cdfft(l);
    148151  double* foo = new double[2*l];
    149   for (int i=0;i<l;i++){
     152  int i;
     153  for (i=0;i<l;i++){
    150154    foo[2*i]=inout[i].real();
    151155    foo[2*i+1]=inout[i].imag();
    152156  }
    153157  cdfftf_(&l, foo, ws_cdfft);
    154   for (int i=0;i<l;i++) inout[i]=complex<double> (foo[2*i],foo[2*i+1]);
     158  for (i=0;i<l;i++) inout[i]=complex<double> (foo[2*i],foo[2*i+1]);
    155159  delete[] foo;
    156160}
Note: See TracChangeset for help on using the changeset viewer.