filters.h

Go to the documentation of this file.

/*
  GPU multi-rate FIR filter bank example software

  Oxford e-Research Centre, Oxford University

  Centre for Digital Music, Queen Mary, University of London.

  This program is free software: you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation, either version 3 of the License, 
  or (at your option) any later version.
  See the file COPYING included with this distribution for more information.
*/


#ifndef FILTERS_H
#define FILTERS_H

#include <cuda_runtime.h> 

// Filter parameters

// use shared memory GPU kernel? 
// it appears to be faster and can be used with long filters 
#define CUDA_SHM 1 

#define SIZE_MUL 1 
#define B_SIZE 11 * SIZE_MUL 

#define A_SIZE (B_SIZE - 1) 

#define OFFSET 20 

#define MAX_RATES 10 
#define MAX_FILTERS 1024 
 
//general params
typedef struct{
  int nfilters; 
  int fsize; 
  int nrates; 
  int rnumf[MAX_RATES]; 
  int streams; 
} params;

//structs of arrays preferred to allow contiguous memory accesses.

typedef struct{
  float buf_in[MAX_FILTERS][B_SIZE + OFFSET]; 
  float buf_out[MAX_FILTERS][A_SIZE + OFFSET]; 
  int m_offb[MAX_FILTERS]; 
  float *bk; 
 } filter_arrays;

typedef struct{
  float *d_filters[MAX_RATES]; 
  float *d_in[MAX_RATES]; 
  int isize[MAX_RATES]; 
  float *d_out[MAX_RATES]; 
  int osize[MAX_RATES]; 
  float *h_in[MAX_RATES]; 
  float *h_out[MAX_RATES]; 
  cudaStream_t stream[MAX_RATES]; 
} gpu_arrays;

// command line arguments
typedef struct {
  int nrates; 
  int nf; 
  int insize; 
  int rconst; 
  int tim; 
}cmd_args;

static  int rdiv[MAX_RATES] = {2, 10, 50, 4, 16, 32, 8, 24, 36, 42}; 

// host cuda prototypes
void compute_ref( float *h_in[], float *h_reference[], gpu_arrays *gpuarrays, params *gparams, cmd_args *args, filter_arrays *farr,  int N);

void compute_omp(float *h_in[], float *h_reference[], gpu_arrays *gpuarrays, params *gparams, cmd_args *args, filter_arrays *farr, int N);

void check_results(float *h_reference[], float *h_out[], gpu_arrays *gpuarrays, params *gparams, int N);

void read_command_line(int argc, char *argv[], cmd_args *args);

void print_usage();

#ifdef __cplusplus

extern "C" {

  void cudaMultiFilterFirStreams(gpu_arrays *arrays, params *params);

  void cudaMultiFilterFirInit(float **in, float **out, float *filters, params *params, gpu_arrays *arrays);

  void cudaMultiFilterFirClose(params *params, gpu_arrays *arrays);
}

#endif

#endif // FILTERS_H