/////////////////////////////////////////////////////////////////
const int OFFSET_ENTREE=0;
const int OFFSET_SORTIE=0;
const int VALEUR_MAX=10000;
const int VALEUR_MIN=-10000;

int NB_COEFF_B;
int NB_COEFF_A;
int MEMORYSIZE;
double *b;           //Coefficients du numérateur de la fonction de transfert du filtre
double *a;           //Coefficients du dénominateur de la fonction de transfert du filtre
                     //le premier élement est a1
                     //le tableau denCoeffs est de taille 1 pour un RIF
double *memoireVk;  //buffer rotatif pour stocker les valeurs de vk
int indice_ecr;              //indice d'écriture dans le buffer rotatif memoireVk
/////////////////////////////////////////////////////////////////
void initFiltreUnEchantillon(sFiltre * pFiltre){
    int i;
    //met à jour les variables globales
    NB_COEFF_B=  pFiltre->nbCoeffB;
    b= pFiltre->pCoeffB;
    NB_COEFF_A= pFiltre->nbCoeffA;
    a= pFiltre->pCoeffA;
    MEMORYSIZE=pFiltre->nbMemoireVk;
    memoireVk=pFiltre->memoireVk;
    for (i=0;i<MEMORYSIZE;i++)
        memoireVk[i]=0;
    indice_ecr=0;
}
/////////////////////////////////////////////////////////////////
double filtreUnEchantillon(double ek){
    double vktemp; // calcul intermediaire pour vk
    int i;  //indice de lecture pour les valeurs des coefficients du filtre
    int j; //indice de lecture pour les valeurs des coefficients du filtre
    int indice_lec; //indice de lecture pour les échantillons d'entrée
    double skout; //valeur pour la sortie calculée
    //suppression offset numérique
    vktemp=ek-OFFSET_ENTREE;
    //pour un RIF, NB_COEFF_A=1 donc vktemp sera inchangé
    indice_lec = indice_ecr;
    //  le premier element utile du tableau denCoeffs est a1
    for (j=1;j<NB_COEFF_A;j++) {  //M itérations
        //ici le premier vk pris en compte est vk-1, donc on décremente indice_lec avant utilisation
        indice_lec=indice_lec-1;
        if (indice_lec < 0)
            indice_lec = MEMORYSIZE - 1;
        vktemp=vktemp-a[j] * memoireVk[indice_lec];
    }
    //rangement de la valeur calculée dans tableau vk à l'indice indice_ecr
    memoireVk[indice_ecr]=vktemp;
    //calcul de sk
    skout = 0; //valeur par défaut pour le résultat
    indice_lec = indice_ecr;
    for (i = 0; i <NB_COEFF_B; i++){ //N+1 itérations
        skout = skout + b[i] * memoireVk[indice_lec];
        //ici le premier vk pris en compte est vk, donc on décremente indice_lec avant après utilisation
        indice_lec=indice_lec-1;
        if (indice_lec < 0)
            indice_lec = MEMORYSIZE - 1;
    }
    //incrémentation du pointeur d'écriture pour la prochaine itération
    indice_ecr = (indice_ecr + 1);
    if (indice_ecr>=MEMORYSIZE)
        indice_ecr=0;
    if (skout >  VALEUR_MAX )	 // gestion des saturations du résultat
        skout = VALEUR_MAX ;
    else if (skout < VALEUR_MIN )
        skout = VALEUR_MIN ;
    //ajout offset numérique
    return  skout + OFFSET_SORTIE;
}
/////////////////////////////////////////////////////////////////

void calculReponseFiltreEnLigne(sFiltre * pFiltre, double * ek, double * sk, int nbPts)
{
   initFiltreUnEchantillon(pFiltre);
    for (int k=0; k < nbPts; k++)    {
        sk[k] = filtreUnEchantillon(ek[k]);
    }
}

/////////////////////////////////////////////////////////////////