main.c

ADSP-TS101S and ADSP-TS201S Real and Complex radix-2 C-callable FFT

// Example C usage of fft_flp32

//************************* Includes ************************************

#include <stdio.h>
#include <sysreg.h>
#include <builtins.h>
#include "FFTDef.h"

//************************* Externs *************************************
extern FFT32( float (*)[], float (*)[], float (*)[], float (*)[], int, int );
extern init();

#pragma align 4
section ("data1ab")
float output[outbuff_size];
#if N==32
  #ifndef FFT_Real
      float input[2*N] = {
                  #include "inputs/input64.dat"
                  };
  #endif
#elif N==64
  #ifdef FFT_Real
      float input[N] = {
                  #include "inputs/input64.dat"
                  };
  #else
      float input[2*N] = {
                  #include "inputs/input128.dat"
                  };
  #endif
#elif N==128
  #ifdef FFT_Real
      float input[N] = {
                  #include "inputs/input128.dat"
                  };
  #else
      float input[2*N] = {
                  #include "inputs/input256.dat"
                  };
  #endif
#elif N==256
  #ifdef FFT_Real
      float input[N] = {
                  #include "inputs/input256.dat"
                  };
  #else
      float input[2*N] = {
                  #include "inputs/input512.dat"
                  };
  #endif
#elif N==512
  #ifdef FFT_Real
      float input[N] = {
                  #include "inputs/input512.dat"
                  };
  #else
      float input[2*N] = {
                  #include "inputs/input1024.dat"
                  };
  #endif
#elif N==1024
  #ifdef FFT_Real
      float input[N] = {
                  #include "inputs/input1024.dat"
                  };
  #else
      float input[2*N] = {
                  #include "inputs/input2048.dat"
                  };
  #endif
#elif N==2048
  #ifdef FFT_Real
      float input[N] = {
                  #include "inputs/input2048.dat"
                  };
  #else
      float input[2*N] = {
                  #include "inputs/input4096.dat"
                  };
  #endif
#elif N==4096
  #ifdef FFT_Real
      float input[N] = {
                  #include "inputs/input4096.dat"
                  };
  #else
      float input[2*N] = {
                  #include "inputs/input8192.dat"
                  };
  #endif
#elif N==8192
  #ifdef FFT_Real
      float input[N] = {
                  #include "inputs/input8192.dat"
                  };
  #else
      float input[2*N] = {
                  #include "inputs/input16384.dat"
                  };
  #endif
#elif N==16384
  #ifdef FFT_Real
      float input[N] = {
                  #include "inputs/input16384.dat"
                  };
  #else
      #ifdef __ADSPTS201__
          float input[2*N] = {
                  #include "inputs/input32768.dat"
                  };
      #endif
  #endif
#elif N==32768
  #ifdef __ADSPTS201__
      #ifdef FFT_Real
          float input[N] = {
                  #include "inputs/input32768.dat"
                  };
      #endif
  #endif
#endif

#ifdef initialization_error
      float input[16384];
#endif


#pragma align 4
section ("data2ab")
	float ping_pong_buffer1[outbuff_size];

#pragma align 4
section ("data3ab")
  float ping_pong_buffer2[outbuff_size];

volatile int
	i,
	tmp_i0,
	tmp_i1;


//*********************************** main ********************************************
void main( void )
{

#ifdef initialization_error

      printf("ERROR: This FFT function works only in the following conditions:\n");
      printf("      -N must be a power of 2\n");
  #ifdef __ADSPTS201__
      printf("      -for real inputs, 64<=N<=32768\n");
      printf("      -for complex inputs, 32<=N<=16384\n");
  #else
      printf("      -for real inputs, 64<=N<=8192\n");
      printf("      -for complex inputs, 32<=N<=4096\n");
  #endif
      printf("Please change the settings in the file FFTDef.h\n");
#else

/*in the case of TS201, at the beginning of the program the
cache must be enabled. The procedure is contained in the
cache_enable macro that uses the refresh rate as input parameter
      -if CCLK=500MHz, refresh_rate=750
      -if CCLK=400MHz, refresh_rate=600
      -if CCLK=300MHz, refresh_rate=450
      -if CCLK=250MHz, refresh_rate=375
*/

#ifdef __ADSPTS201__
  asm("#include <defts201.h>");
  asm("#include <cache_macros.h>" );
  asm("cache_enable(750);");
  asm("#include <ini_cache.h>");
  asm("#include <fftdef.h>");
  asm("preload_cache;");
#endif


	tmp_i0 = __builtin_sysreg_read( __CCNT0 );         // read initial cycle count

    #ifdef FFT_Real
        FFT32(&(input), &(ping_pong_buffer1), &(ping_pong_buffer2), &(output), N, REAL);
    #else
        FFT32(&(input), &(ping_pong_buffer1), &(ping_pong_buffer2), &(output), N, COMPLEX);
    #endif

	tmp_i1 = __builtin_sysreg_read( __CCNT0 );         // read final cycle count
	printf("cycle count = %d", tmp_i1 - tmp_i0);       // print the cycle count

#endif

}