filter_util.h

Go to the documentation of this file.

#ifndef _FILTER_UTIL_H_
#define _FILTER_UTIL_H_

#include "num.h"

#ifdef __GNUC__
__attribute__((pure,unused))
#endif

static long extension(long N, long n)
{
    while (1)
        if (n < 0)
            n = -1 - n;         /* Reflect over n = -1/2.    */
        else if (n >= N)
            n = 2 * N - 1 - n;  /* Reflect over n = N - 1/2. */
        else
            break;
        
    return n;
}

void recursive_filter_impulse(num *h, long N,
    const num *b, int p, const num *a, int q);

void init_recursive_filter(num *dest, const num *src, long N, long stride,
    const num *b, int p, const num *a, int q,
    num sum, num tol, long max_iter);

#endif /* _FILTER_UTIL_H_ */