Part Number: TMS320C6678
Tool/software: TI-RTOS
We have used the DSPF_SP_IFFT example in the documentation and combined the DSPF_SP_FFT example. We also used the linker of the former file.
We generated the data and ifft it then fft. Also, we used the individual gen. twiddle of the files.
However, we cannot recover the original data after fft. Please help us. Also, if it is possible please do not refer us to another link/forum. Thank you.
/* ======================================================================= */ /* DSPF_sp_ifftSPxSP.c -- Complex Inverse FFT with Mixed Radix */ /* Driver code; tests kernel and reports result in stdout */ /* */ /* Rev 0.0.2 */ /* */ /* Copyright (C) 2011 Texas Instruments Incorporated - http://www.ti.com/ */ /* */ /* */ /* Redistribution and use in source and binary forms, with or without */ /* modification, are permitted provided that the following conditions */ /* are met: */ /* */ /* Redistributions of source code must retain the above copyright */ /* notice, this list of conditions and the following disclaimer. */ /* */ /* Redistributions in binary form must reproduce the above copyright */ /* notice, this list of conditions and the following disclaimer in the */ /* documentation and/or other materials provided with the */ /* distribution. */ /* */ /* Neither the name of Texas Instruments Incorporated nor the names of */ /* its contributors may be used to endorse or promote products derived */ /* from this software without specific prior written permission. */ /* */ /* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS */ /* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT */ /* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR */ /* A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT */ /* OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, */ /* SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT */ /* LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, */ /* DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY */ /* THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT */ /* (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE */ /* OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. */ /* */ /* ======================================================================= */ #include <stdio.h> #include <time.h> #include <stdlib.h> #include <limits.h> #include <math.h> #include "DSPF_sp_ifftSPxSP.h" #include "DSPF_sp_ifftSPxSP_opt.h" #include "DSPF_sp_ifftSPxSP_cn.h" /* Defines */ #if defined(__TI_EABI__) #define kernel_size _kernel_size #endif extern char kernel_size; #define CYCLE_FORMULA_NX_PT1 128 #define CYCLE_FORMULA_NX_PT2 256 #define FORMULA_SIZE 2 int form_cycle [FORMULA_SIZE]; void tw_gen (float *w, int n); void tw_gen1 (float *w, int n); /* ======================================================================== */ /* Kernel-specific alignments */ /* ======================================================================== */ #pragma DATA_ALIGN(x_i, 8); #pragma DATA_ALIGN(x_sa, 8); #pragma DATA_ALIGN(x_cn, 8); #pragma DATA_ALIGN(w, 8); #pragma DATA_ALIGN(y_i, 8); #pragma DATA_ALIGN(y_sa, 8); #pragma DATA_ALIGN(y_cn, 8); #pragma DATA_ALIGN(brev, 8); /* ======================================================================== */ /* Parameters of fixed dataset. */ /* ======================================================================== */ #define MAXN (1024) #define M (2*MAXN) #define PAD (16) /* ======================================================================== */ /* Initialized arrays with fixed test data. */ /* ======================================================================== */ float x_i [M + 2 * PAD]; float x_sa[M + 2 * PAD]; float x_cn[M + 2 * PAD]; float w [M + 2 * PAD]; float y_i [M + 2 * PAD]; float y_sa[M + 2 * PAD]; float y_cn[M + 2 * PAD]; /* ======================================================================== */ /* Generate pointers to skip beyond array padding */ /* ======================================================================== */ float *const ptr_x_i = x_i + PAD; float *const ptr_x_sa = x_sa + PAD; float *const ptr_x_cn = x_cn + PAD; float *const ptr_w = w + PAD; float *const ptr_y_i = y_i + PAD; float *const ptr_y_sa = y_sa + PAD; float *const ptr_y_cn = y_cn + PAD; unsigned char brev[64] = { 0x0, 0x20, 0x10, 0x30, 0x8, 0x28, 0x18, 0x38, 0x4, 0x24, 0x14, 0x34, 0xc, 0x2c, 0x1c, 0x3c, 0x2, 0x22, 0x12, 0x32, 0xa, 0x2a, 0x1a, 0x3a, 0x6, 0x26, 0x16, 0x36, 0xe, 0x2e, 0x1e, 0x3e, 0x1, 0x21, 0x11, 0x31, 0x9, 0x29, 0x19, 0x39, 0x5, 0x25, 0x15, 0x35, 0xd, 0x2d, 0x1d, 0x3d, 0x3, 0x23, 0x13, 0x33, 0xb, 0x2b, 0x1b, 0x3b, 0x7, 0x27, 0x17, 0x37, 0xf, 0x2f, 0x1f, 0x3f }; /* ======================================================================== */ /* MAIN -- Top level driver for the test. */ /* ======================================================================== */ int main () { int i, j, rad, N, k; clock_t t_start, t_stop, t_overhead, t_opt, t_i, t_cn; float pct_diff, max_pct_diff = 0; /* for (N = 8, k = 1; N <= MAXN; N = N * 2, k++) */ N=64; k=1; memset (x_i, 0x55, sizeof (x_i)); memset (x_sa, 0x55, sizeof (x_sa)); memset (x_cn, 0x55, sizeof (x_cn)); /* -------------------------------------------------------------------- */ /* Initialize input vector temporarily. */ /* -------------------------------------------------------------------- */ for (i = 0; i < N; i++) { x_cn[PAD + 2 * i] = sin (2 * 3.1415 * 50 * i / (double) N); x_cn[PAD + 2 * i + 1] = sin (2 * 3.1415 * 100 * i / (double) N); x_sa[PAD + 2 * i] = x_cn[PAD + 2 * i]; x_sa[PAD + 2 * i + 1] = x_cn[PAD + 2 * i + 1]; x_i [PAD + 2 * i] = x_cn[PAD + 2 * i]; x_i [PAD + 2 * i + 1] = x_cn[PAD + 2 * i + 1]; printf("signal generated: %f,\n",x_cn[PAD + 2 * i] ); printf("signal generated: %f,\n", x_cn[PAD + 2 * i + 1]); } /* -------------------------------------------------------------------- */ /* Force uninitialized arrays to fixed values. */ /* -------------------------------------------------------------------- */ memset (y_i, 0xA5, sizeof (y_i)); memset (y_sa, 0xA5, sizeof (y_sa)); memset (y_cn, 0xA5, sizeof (y_cn)); /* -------------------------------------------------------------------- */ /* Compute the overhead of calling clock() twice to get timing info. */ /* -------------------------------------------------------------------- */ t_start = clock (); t_stop = clock (); t_overhead = t_stop - t_start; /* -------------------------------------------------------------------- */ /* Generate twiddle factors. */ /* -------------------------------------------------------------------- */ /* j = 0; for (i = 0; i <= 31; i++) if ((N & (1 << i)) == 0) j++; else break; if (j % 2 == 0) rad = 4; else */ rad = 2; tw_gen (ptr_w, N); printf("DSPF_sp_ifftSPxSP\tIter#: %d\t", k); t_start = clock (); DSPF_sp_ifftSPxSP_cn (N, &ptr_x_cn[0], &ptr_w[0], ptr_y_cn, brev, rad, 0, N); t_stop = clock (); t_cn = (t_stop - t_start) - t_overhead; /* t_start = clock (); DSPF_sp_ifftSPxSP_opt (N, &ptr_x_i[0], &ptr_w[0], ptr_y_cn, brev, rad, 0, N); t_stop = clock (); t_i = (t_stop - t_start) - t_overhead; t_start = clock (); DSPF_sp_ifftSPxSP (N, &ptr_x_sa[0], &ptr_w[0], ptr_y_sa, brev, rad, 0, N); t_stop = clock (); t_opt = (t_stop - t_start) - t_overhead; */ /* ---------------------------------------------------------------- */ /* compute percent difference and track max difference */ /* ---------------------------------------------------------------- */ /* pct_diff = 0; max_pct_diff = 0; for(i=0; i<N; i++) { pct_diff = (ptr_y_cn[i] - ptr_y_i[i]) / ptr_y_cn[i] * 100.0; if (pct_diff < 0) pct_diff *= -1; if (pct_diff > max_pct_diff) max_pct_diff = pct_diff; } if (max_pct_diff > 0.001) printf("Intrinsic Failure max_pct_diff = %f", max_pct_diff); else printf("Intrinsic Successful "); pct_diff = 0; max_pct_diff = 0; for(i=0; i<N; i++) { pct_diff = (ptr_y_cn[i] - ptr_y_i[i]) / ptr_y_cn[i] * 100.0; if (pct_diff < 0) pct_diff *= -1; if (pct_diff > max_pct_diff) max_pct_diff = pct_diff; } if (max_pct_diff > 0.001) printf("SA Failure max_pct_diff = %f", max_pct_diff); else printf("SA Successful "); */ printf("\tN = %d\tradix = %d\tnatC: %d\toptC: %d\tSA: %d\n", N, rad, t_cn, t_i, t_opt); if (N == CYCLE_FORMULA_NX_PT1) form_cycle[0] = t_opt; if (N == CYCLE_FORMULA_NX_PT2) form_cycle[1] = t_opt; /* Provide memory information */ printf("Memory: %d bytes\n", &kernel_size); /* Provide profiling information */ printf("Cycles: %d (N=128) %d (N=256)\n", form_cycle[0], form_cycle[1]); for(i=PAD;i<(N*2+PAD);i++){ printf("IFFT output Signal:%f\n",y_cn[i]); x_i[i]=y_cn[i]; } printf("DSPF_sp_fftSPxSP\tIter#: %d\t\n", k); printf("output Signal:%f\n",x_cn[PAD]); //for test t_start = clock (); tw_gen1 (ptr_w, N); DSPF_sp_fftSPxSP_cn (N, &ptr_x_i[0], &ptr_w[0], ptr_y_i, brev, rad, 0, N); t_stop = clock (); t_cn = (t_stop - t_start) - t_overhead; for(i=PAD;i<(N*2+PAD);i++){ printf("fft output Signal:%f\n",y_i[i]); //x_cn[PAD + i]=y_cn[i]; } } /* Function for generating Specialized sequence of twiddle factors */ void tw_gen (float *w, int n) { int i, j, k; const double PI = 3.141592654; for (j = 1, k = 0; j <= n >> 2; j = j << 2) { for (i = 0; i < n >> 2; i += j) { #ifdef _LITTLE_ENDIAN w[k] = (float) -sin (2 * PI * i / n); w[k + 1] = (float) cos (2 * PI * i / n); w[k + 2] = (float) -sin (4 * PI * i / n); w[k + 3] = (float) cos (4 * PI * i / n); w[k + 4] = (float) -sin (6 * PI * i / n); w[k + 5] = (float) cos (6 * PI * i / n); #else w[k] = (float) cos (2 * PI * i / n); w[k + 1] = (float) sin (2 * PI * i / n); w[k + 2] = (float) cos (4 * PI * i / n); w[k + 3] = (float) sin (4 * PI * i / n); w[k + 4] = (float) cos (6 * PI * i / n); w[k + 5] = (float) sin (6 * PI * i / n); #endif k += 6; } } } void tw_gen1 (float *w, int n) { int i, j, k; const double PI = 3.141592654; for (j = 1, k = 0; j <= n >> 2; j = j << 2) { for (i = 0; i < n >> 2; i += j) { #ifdef _LITTLE_ENDIAN w[k] = (float) sin (2 * PI * i / n); w[k + 1] = (float) cos (2 * PI * i / n); w[k + 2] = (float) sin (4 * PI * i / n); w[k + 3] = (float) cos (4 * PI * i / n); w[k + 4] = (float) sin (6 * PI * i / n); w[k + 5] = (float) cos (6 * PI * i / n); #else w[k] = (float) cos (2 * PI * i / n); w[k + 1] = (float) -sin (2 * PI * i / n); w[k + 2] = (float) cos (4 * PI * i / n); w[k + 3] = (float) -sin (4 * PI * i / n); w[k + 4] = (float) cos (6 * PI * i / n); w[k + 5] = (float) -sin (6 * PI * i / n); #endif k += 6; } } } /* ======================================================================= */ /* End of file: DSPF_sp_ifftSPxSP_d.c */ /* ----------------------------------------------------------------------- */ /* Copyright (c) 2011 Texas Instruments, Incorporated. */ /* All Rights Reserved. */ /* ======================================================================= */
