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
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
|
/*
* Copyright (C) 2020 Linux Studio Plugins Project <https://lsp-plug.in/>
* (C) 2020 Vladimir Sadovnikov <sadko4u@gmail.com>
*
* This file is part of lsp-dsp-units
* Created on: 29 янв. 2016 г.
*
* lsp-dsp-units is free software: you can redistribute it and/or modify
* it under the terms of the GNU Lesser General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* any later version.
*
* lsp-dsp-units is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public License
* along with lsp-dsp-units. If not, see <https://www.gnu.org/licenses/>.
*/
#ifndef LSP_PLUG_IN_DSP_UNITS_UTIL_CONVOLVER_H_
#define LSP_PLUG_IN_DSP_UNITS_UTIL_CONVOLVER_H_
#include <lsp-plug.in/dsp-units/version.h>
#include <lsp-plug.in/dsp-units/iface/IStateDumper.h>
#define CONVOLVER_RANK_MIN 8 /* buffer of 256 samples (128 effective) */
#define CONVOLVER_RANK_MAX 16 /* buffer of 8192 samples (4096 effective) */
namespace lsp
{
namespace dspu
{
class LSP_DSP_UNITS_PUBLIC Convolver
{
private:
Convolver & operator = (const Convolver &);
Convolver(const Convolver &);
private:
float *vDataBuffer; // Buffer for storing convolution tail data
float *vFrame; // Pointer to the beginning of the input data frame
float *vConvBuffer; // Convolution buffer to perform convolution
float *vTaskData; // Task data for tail convolution
float *vConvData; // FFT convolution data
float *vDirectData; // Direct convolution data
size_t nDataBufferSize; // Size of data buffer
size_t nDirectSize; // Size of direct convolution data
size_t nFrameSize; // Size of input data frame
size_t nFrameOff; // Offset from the beginning of the input data frame
size_t nConvSize; // The actual convolution size in samples
size_t nLevels; // Number of raising convolution levels
size_t nBlocks; // Number of constant-size blocks
size_t nBlocksDone; // Number of applied constant-size blocks
size_t nRank; // The actual rank of the convolution
size_t nBlkInit; // Initial number of blocks to apply at step # 0
float fBlkCoef; // The actual coefficient to compute proper block number per formula
uint8_t *vData; // Non-aligned pointer to the whole allocated data
public:
explicit Convolver();
~Convolver();
/** Construct the convolver
*
*/
void construct();
/** Destroy convolver
*
*/
void destroy();
public:
/** Initialize convolver
*
* @param data convolution data
* @param count number of samples in convolution
* @param rank convolution rank
* @return true on success
*/
bool init(const float *data, size_t count, size_t rank, float phase);
/** Process samples
*
* @param dst destination buffer
* @param src source buffer
* @param count number of samples to process
*/
void process(float *dst, const float *src, size_t count);
/** Get the actual convolution size in samples
*
* @return actual convolution size in samples
*/
inline size_t data_size() const { return nConvSize; }
/**
* Get actual convolution rank of the convolver
* @return convolution rank
*/
inline size_t rank() const { return nRank; }
/**
* Dump internal state
* @param v state dumper
*/
void dump(IStateDumper *v) const;
};
}
} /* namespace lsp */
#endif /* LSP_PLUG_IN_DSP_UNITS_UTIL_CONVOLVER_H_ */
|
