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
71
72
73
74
75
76
77
78
79
80
81
82
83
84
|
#define TYPEDEPARGS 0, 1
#define SINGLEARGS
#define COMPLEXARGS
#define OCTFILENAME comp_ifilterbank_fft // change to filename
#define OCTFILEHELP "This function calls the C-library \n\
F = comp_ifilterbank_fft(c,G,a) \n\
Yeah."
#include "ltfat_oct_template_helper.h"
// octave_idx_type 32 or 64 bit signed integer
static inline void
fwd_ifilterbank_fft(const Complex *c[], const Complex *filt[],
const ltfat_int L, const ltfat_int W,
ltfat_int a[], const ltfat_int M,
Complex *f)
{
ltfat_ifilterbank_fft_d(reinterpret_cast<const ltfat_complex_d **>(c),
reinterpret_cast<const ltfat_complex_d **>(filt),
L, W, a, M,
reinterpret_cast<ltfat_complex_d *>(f));
}
static inline void
fwd_ifilterbank_fft(const FloatComplex *c[], const FloatComplex *filt[],
const ltfat_int L, const ltfat_int W,
ltfat_int a[], const ltfat_int M,
FloatComplex *f)
{
ltfat_ifilterbank_fft_s(reinterpret_cast<const ltfat_complex_s **>(c),
reinterpret_cast<const ltfat_complex_s **>(filt),
L, W, a, M,
reinterpret_cast<ltfat_complex_s *>(f));
}
template <class LTFAT_TYPE, class LTFAT_REAL, class LTFAT_COMPLEX>
octave_value_list octFunction(const octave_value_list& args, int nargout)
{
// Input data
Cell c = args(0).cell_value();
// Cell aray containing impulse responses
Cell G = args(1).cell_value();
// Subsampling factors
Matrix aDouble = args(2).matrix_value();
// Output length
const octave_idx_type L = G.elem(0).rows();
// Number of channels
const octave_idx_type W = c.elem(0).columns();
// Number of filters
const octave_idx_type M = G.numel();
// Output signal
MArray<LTFAT_COMPLEX> F(dim_vector(L, W));
// Allocating temporary arrays
// Output subband lengths
// Impulse responses pointers
OCTAVE_LOCAL_BUFFER (const LTFAT_COMPLEX*, GPtrs, M);
// Output subbands pointers
OCTAVE_LOCAL_BUFFER (const LTFAT_COMPLEX*, cPtrs, M);
// Output cell elements array,
OCTAVE_LOCAL_BUFFER (MArray<LTFAT_COMPLEX>, cElems, M);
//
OCTAVE_LOCAL_BUFFER (MArray<LTFAT_COMPLEX>, GElems, M);
OCTAVE_LOCAL_BUFFER (ltfat_int, a, M);
for (octave_idx_type m = 0; m < M; m++)
{
a[m] = (ltfat_int) aDouble(m);
GElems[m] = ltfatOctArray<LTFAT_COMPLEX>(G.elem(m));
GPtrs[m] = GElems[m].data();
cElems[m] = ltfatOctArray<LTFAT_COMPLEX>(c.elem(m));
cPtrs[m] = cElems[m].data();
}
fwd_ifilterbank_fft(cPtrs, GPtrs, L, W, a, M, F.fortran_vec());
return octave_value(F);
}
|
