Changeset 515 in Sophya for trunk/SophyaLib/NTools/fftserver.cc
- Timestamp:
- Oct 26, 1999, 2:44:44 PM (26 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaLib/NTools/fftserver.cc
r514 r515 92 92 checkint_cfft(l); 93 93 float* foo = new float[2*l]; 94 for (int i=0;i<l;i++){ 94 int i; 95 for (i=0;i<l;i++){ 95 96 foo[2*i]=inout[i].real(); 96 97 foo[2*i+1]=inout[i].imag(); … … 98 99 cfftf_(&l, foo, ws_cfft); 99 100 inout[0]=complex<float> (foo[0],foo[1]); 100 for (i nt 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]); 101 102 delete[] foo; 102 103 } … … 106 107 checkint_cdfft(l); 107 108 double* foo=new double[2*l]; 108 for (int i=0;i<l;i++){ 109 int i; 110 for (i=0;i<l;i++){ 109 111 foo[2*i]=inout[i].real(); 110 112 foo[2*i+1]=inout[i].imag(); … … 112 114 cdfftf_(&l, foo, ws_cdfft); 113 115 inout[0]=complex<double> (foo[0],foo[1]); 114 for (i nt i=1;i<l;i++) {116 for (i=1;i<l;i++) { 115 117 inout[l-i]= complex<double> (foo[2*i],foo[2*i+1]); 116 118 } … … 134 136 checkint_cfft(l); 135 137 float* foo = new float[2*l]; 136 for (int i=0;i<l;i++){ 138 int i; 139 for (i=0;i<l;i++){ 137 140 foo[2*i]=inout[i].real(); 138 141 foo[2*i+1]=inout[i].imag(); 139 142 } 140 143 cfftf_(&l, foo, ws_cfft); 141 for (i nt 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]); 142 145 delete[] foo; 143 146 } … … 147 150 checkint_cdfft(l); 148 151 double* foo = new double[2*l]; 149 for (int i=0;i<l;i++){ 152 int i; 153 for (i=0;i<l;i++){ 150 154 foo[2*i]=inout[i].real(); 151 155 foo[2*i+1]=inout[i].imag(); 152 156 } 153 157 cdfftf_(&l, foo, ws_cdfft); 154 for (i nt 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]); 155 159 delete[] foo; 156 160 }
Note:
See TracChangeset
for help on using the changeset viewer.