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
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
|
/*
* 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-lib
* Created on: 31 мар. 2020 г.
*
* lsp-dsp-lib 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-lib 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-lib. If not, see <https://www.gnu.org/licenses/>.
*/
#ifndef LSP_PLUG_IN_DSP_COMMON_COMPLEX_H_
#define LSP_PLUG_IN_DSP_COMMON_COMPLEX_H_
#include <lsp-plug.in/dsp/common/types.h>
/** Calculate complex multiplication: { dst_re, dst_im } = { src1_re, src1_im } * { src2_re, src2_im }
*
* @param dst_re destination real part
* @param dst_im destination imaginary part
* @param src1_re source 1 real part
* @param src1_im source 1 imaginary part
* @param src2_re source 2 real part
* @param src2_im source 2 imaginary part
* @param count number of multiplications
*/
LSP_DSP_LIB_SYMBOL(void, complex_mul3,
float *dst_re, float *dst_im,
const float *src1_re, const float *src1_im,
const float *src2_re, const float *src2_im,
size_t count
);
/** Calculate complex multiplication: { dst_re, dst_im } = { dst_re, dst_im } * { src_re, src_im }
*
* @param dst_re destination and source real part
* @param dst_im destination and source imaginary part
* @param src_re source real part
* @param src_im source imaginary part
* @param count number of multiplications
*/
LSP_DSP_LIB_SYMBOL(void, complex_mul2,
float *dst_re, float *dst_im,
const float *src_re, const float *src_im,
size_t count
);
/**
* Divide complex numbers:
* { dst_re, dst_im } = { dst_re, dst_im } / { src_re, src_im }
*
* @param dst_re destination real part, source top real part
* @param dst_im destination imaginary part, source top imaginary part
* @param src_re bottom real part
* @param src_im bottom imaginary part
* @param count number of divisions
*/
LSP_DSP_LIB_SYMBOL(void, complex_div2, float *dst_re, float *dst_im, const float *src_re, const float *src_im, size_t count);
/**
* Divide complex numbers:
* { dst_re, dst_im } = { src_re, src_im } / { dst_re, dst_im }
*
* @param dst_re destination real part, source bottom real part
* @param dst_im destination imaginary part, source bottom imaginary part
* @param src_re top real part
* @param src_im top imaginary part
* @param count number of divisions
*/
LSP_DSP_LIB_SYMBOL(void, complex_rdiv2, float *dst_re, float *dst_im, const float *src_re, const float *src_im, size_t count);
/**
* Divide complex numbers:
* { dst_re, dst_im } = { t_re, t_im } / { b_re, b_im }
*
* @param dst_re destination real part, source bottom real part
* @param dst_im destination imaginary part, source bottom imaginary part
* @param t_re top real part
* @param t_im top imaginary part
* @param b_re bottom real part
* @param b_im bottom imaginary part
* @param count number of divisions
*/
LSP_DSP_LIB_SYMBOL(void, complex_div3, float *dst_re, float *dst_im, const float *t_re, const float *t_im, const float *b_re, const float *b_im, size_t count);
/** Calculate complex reciprocal: 1 / (re + j * im)
*
* @param dst_re source, destination real part
* @param dst_im source, destination imaginary part
* @param count number of multiplications
*/
LSP_DSP_LIB_SYMBOL(void, complex_rcp1,
float *dst_re, float *dst_im,
size_t count
);
/** Calculate complex reciprocal: 1 / (re + j * im)
*
* @param dst_re destination real part
* @param dst_im destination imaginary part
* @param src_re source real part
* @param src_im source imaginary part
* @param count number of multiplications
*/
LSP_DSP_LIB_SYMBOL(void, complex_rcp2,
float *dst_re, float *dst_im,
const float *src_re, const float *src_im,
size_t count
);
/** Convert real+imaginary complex number to polar form
*
* @param dst_mod module of the complex number
* @param dst_arg argument of the complex number in range of [-PI, +PI]
* @param src_re real part of complex number
* @param src_im imaginary part of complex number
* @param count number of elements to process
*/
LSP_DSP_LIB_SYMBOL(void, complex_cvt2modarg,
float *dst_mod, float *dst_arg,
const float *src_re, const float *src_im,
size_t count
);
/** Get module for complex numbers: mod = sqrt(re*re + im*im)
*
* @param dst_mod array to sore module
* @param src_re real part of complex number
* @param src_im imaginary part of complex number
* @param count number of elements to process
*/
LSP_DSP_LIB_SYMBOL(void, complex_mod,
float *dst_mod,
const float *src_re, const float *src_im,
size_t count
);
/** Get argument for complex numbers in range of [-PI; +PI]
*
* @param dst array to sore argument
* @param src packed complex number data
* @param count count number of elements to process
*/
LSP_DSP_LIB_SYMBOL(void, complex_arg, float *dst, const float *re, const float *im, size_t count);
/** Convert polar-form of complex number to real+imaginary
*
* @param dst_re real part of complex number
* @param dst_im imaginary part of complex number
* @param src_mod module of the complex number
* @param src_arg argument of the complex number
* @param count number of elements to process
*/
LSP_DSP_LIB_SYMBOL(void, complex_cvt2reim,
float *dst_re, float *dst_im,
const float *src_mod, const float *src_arg,
size_t count
);
#endif /* LSP_PLUG_IN_DSP_COMMON_COMPLEX_H_ */
|
