Low level routine
eng -
fr
convolr_c - finite impulse response filter by fft and overlap-add method computational routine
- n : size of the original vector
- nb_coef : length of the impulse response vector
- m1 : size of the fft vector (2x)
- [z1_r;z1_i] : address of the input complex vector 1 (signal)
- [z2_r;z2_i] : address of the input complex vector 2 (impulse reponse)
- [y_r;y_i] : address of the output complex vector
- z_r : address of the overlaped word
/* convolr_c subroutine
* FIR computation
* with FFT convolution
* and overlap-ad method
* IRCOM GROUP - Author : A.Layec
*/
/* REVISION HISTORY :
* $Log$
*/
#include "mod_num_lib.h"
/* convolr_c routine qui realise un filtre FIR
* par multiplication dans le domaine frequentiel
* avec mise mémoire .
*
* n : taille du mot original
* nb_coef : longueur de la réponse impulsionnelle
* m1 : taille de la fft
* [z1_r;z1_i] : adresses de départ du vecteur complexe 1
* [z2_r;z2_i] : adresses de départ du vecteur complexe 2
* [y_r;y_i] : adresses de départ du vecteur complexe résultat
* z_r : adresse de départ du mot mémoire
*
* utilise : dset (BLAS)
* convolfft_c
* overlapadr_c
*/
void convolr_c(int *n,int *nb_coef,int *m1,double *z1_r,double *z1_i,double *z2_r,double *z2_i,double *y_r,double *y_i,double *z_r)
{
/*déclaration*/
int i,l,k;
/*ajoute les zéros au vecteur z1_r*/
F2C(dset)((k=(*m1)-(*n),&k),(l=0,&l),&z1_r[(*n)],(i=1,&i));
/*place valeur img z1_i a zero*/
F2C(dset)((k=(*m1),&k),(l=0,&l),&z1_i[0],(i=1,&i));
/*ajoute les zéros au vecteur z2_r*/
F2C(dset)((k=(*m1)-(*nb_coef),&k),(l=0,&l),&z2_r[(*nb_coef)],(i=1,&i));
/*place valeur img z2_i a zero*/
F2C(dset)((k=(*m1),&k),(l=0,&l),&z2_i[0],(i=1,&i));
/*Appel convolfft_c*/
convolfft_c(m1,&z1_r[0],&z1_i[0],&z2_r[0],&z2_i[0],&y_r[0],&y_i[0]);
/*Appel overlapad_c*/
overlapadr_c(m1,n,nb_coef,&y_r[0],&z_r[0]);
return;
}
IRCOM Group
Alan Layec