diff options
Diffstat (limited to 'recipes-demo/de10-nano-fftsw-apps/files/c16_256x32x128.c')
-rw-r--r-- | recipes-demo/de10-nano-fftsw-apps/files/c16_256x32x128.c | 113 |
1 files changed, 113 insertions, 0 deletions
diff --git a/recipes-demo/de10-nano-fftsw-apps/files/c16_256x32x128.c b/recipes-demo/de10-nano-fftsw-apps/files/c16_256x32x128.c new file mode 100644 index 0000000..c015f80 --- /dev/null +++ b/recipes-demo/de10-nano-fftsw-apps/files/c16_256x32x128.c @@ -0,0 +1,113 @@ +#include <stdio.h> +#include <stdlib.h> +#include <sys/types.h> +#include <sys/stat.h> +#include <fcntl.h> +#include <unistd.h> +#include <sys/mman.h> +#include <error.h> +#include <errno.h> + +#include "NE10.h" + +/* customize these macros and variables for a particular implementation */ +#define FFT_IN_ROUNDS (128) +#define FFT_CALC_POINTS (256) +#define FFT_CALC_ROUNDS (32) +#define FFT_POINTS (FFT_CALC_POINTS * FFT_CALC_ROUNDS) +#define IN_SAMPLE_TYPE ne10_fft_cpx_int16_t +#define OUT_SAMPLE_TYPE ne10_fft_cpx_int16_t +#define CFG_TYPE ne10_fft_cfg_int16_t +#define CFG_ALLOC_FUNC ne10_fft_alloc_c2c_int16 +#define FFT_FUNC ne10_fft_c2c_1d_int16_c + +const char *g_test = "c16_256x32x128"; +const char *g_description = "256x32x128 point FFT, C implementation, 16-bit precision"; + +/* these functions and variables are implemented in overhead.c as a library */ +extern void initialize_everything(int argc, char **argv); +extern uint64_t get_gt_value(void); +extern void release_everything(void); +extern void print_results(void); + +extern char *g_input_filename; +extern char *g_output_filename; +extern uint64_t g_start_time; +extern uint64_t g_end_time; + +int main(int argc, char **argv) +{ + int input_fd; + int output_fd; + int result; + IN_SAMPLE_TYPE *in_buf; + OUT_SAMPLE_TYPE *out_buf; + CFG_TYPE cfg; + int i; + int j; + + initialize_everything(argc, argv); + + /* allocate storage for input, output and config buffers */ + in_buf = (IN_SAMPLE_TYPE*) NE10_MALLOC (FFT_POINTS * sizeof(IN_SAMPLE_TYPE)); + if(in_buf == NULL) + error(1, errno, "in_buf allocation"); + + out_buf = (OUT_SAMPLE_TYPE*) NE10_MALLOC (FFT_POINTS * sizeof(OUT_SAMPLE_TYPE)); + if(out_buf == NULL) + error(1, errno, "out_buf allocation"); + + cfg = CFG_ALLOC_FUNC(FFT_CALC_POINTS); + if(cfg == NULL) + error(1, errno, "cfg allocation"); + + /* open the input and output files */ + input_fd = open(g_input_filename, O_RDONLY); + if(input_fd < 0) + error(1, errno, "opening input file '%s'", g_input_filename); + + output_fd = open(g_output_filename, O_WRONLY | O_CREAT, ~((mode_t)umask)); + if(output_fd < 0) + error(1, errno, "opening output file '%s'", g_output_filename); + + /* capture the start value of the GT */ + g_start_time = get_gt_value(); + + for(j = 0 ; j < FFT_IN_ROUNDS ; j++) { + /* read the input data */ + result = read(input_fd, in_buf, FFT_POINTS * sizeof(IN_SAMPLE_TYPE)); + if(result < 0) + error(1, errno, "read input file"); + if(result != (FFT_POINTS * sizeof(IN_SAMPLE_TYPE))) + error(1, 0, "input data size, expected %d but got %d", FFT_POINTS * sizeof(IN_SAMPLE_TYPE), result); + + /* compute FFT */ + for (i = 0; i < FFT_CALC_ROUNDS ; i++) { + FFT_FUNC(out_buf + (i * FFT_CALC_POINTS), in_buf + (i * FFT_CALC_POINTS), cfg, 0, 1); + } + + /* write the output data */ + result = write(output_fd, out_buf, FFT_POINTS * sizeof(OUT_SAMPLE_TYPE)); + if(result < 0) + error(1, errno, "write output file"); + if(result != (FFT_POINTS * sizeof(OUT_SAMPLE_TYPE))) + error(1, 0, "output data size, expected %d but got %d", FFT_POINTS * sizeof(OUT_SAMPLE_TYPE), result); + } + + /* capture the end value of the GT */ + g_end_time = get_gt_value(); + + /* close the input and output files */ + close(output_fd); + close(input_fd); + + /* free storage for input, output and config buffers */ + NE10_FREE (cfg); + NE10_FREE (out_buf); + NE10_FREE (in_buf); + + print_results(); + release_everything(); + return 0; +} + |