1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
|
#include "iir.h"
#include <ippi.h>
#include <ipps.h>
#include <ippcc.h>
/* 32 bit floating point */
int nsutil_iir_Create_F32(const float *coefficients, int order, nsutil_iir_t *out_iir)
{
IppsIIRState_32f *state;
ippsIIRInitAlloc_32f(&state, coefficients, order, 0);
*out_iir = state;
return 0;
}
int nsutil_iir_Filter_F32_IP(nsutil_iir_t *iir, float *samples, int num_samples)
{
ippsIIR_32f_I(samples, num_samples, (IppsIIRState_32f *)iir);
return 0;
}
int nsutil_iir_Filter_F32(nsutil_iir_t *iir, const float *input, float *output, int num_samples)
{
ippsIIR_32f(input, output, num_samples, (IppsIIRState_32f *)iir);
return 0;
}
int nsutil_iir_Reset_F32(nsutil_iir_t *iir)
{
ippsIIRSetDlyLine_32f((IppsIIRState_32f *)iir, 0);
return 0;
}
int nsutil_iir_Destroy_F32(nsutil_iir_t *iir)
{
ippsIIRFree_32f((IppsIIRState_32f *)iir);
return 0;
}
/* 64 bit floating point */
int nsutil_iir_Create_F64(const double *coefficients, int order, nsutil_iir_t *out_iir)
{
IppsIIRState_64f *state;
ippsIIRInitAlloc_64f(&state, coefficients, order, 0);
*out_iir = state;
return 0;
}
int nsutil_iir_Filter_F64_IP(nsutil_iir_t *iir, double *samples, int num_samples)
{
ippsIIR_64f_I(samples, num_samples, (IppsIIRState_64f *)iir);
return 0;
}
int nsutil_iir_Filter_642(nsutil_iir_t *iir, const double *input, double *output, int num_samples)
{
ippsIIR_64f(input, output, num_samples, (IppsIIRState_64f *)iir);
return 0;
}
int nsutil_iir_Reset_F64(nsutil_iir_t *iir)
{
ippsIIRSetDlyLine_64f((IppsIIRState_64f *)iir, 0);
return 0;
}
int nsutil_iir_Destroy_F64(nsutil_iir_t *iir)
{
ippsIIRFree_64f((IppsIIRState_64f *)iir);
return 0;
}
|
