libdspl-2.0
Библиотека алгоритмов цифровой обработки сигналов
complex.c
1 /*
2 * Copyright (c) 2015-2020 Sergey Bakhurin
3 * Digital Signal Processing Library [http://dsplib.org]
4 *
5 * This file is part of DSPL.
6 *
7 * is free software: you can redistribute it and/or modify
8 * it under the terms of the GNU General Public License as published by
9 * the Free Software Foundation, either version 3 of the License, or
10 * (at your option) any later version.
11 *
12 * DSPL is distributed in the hope that it will be useful,
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15 * GNU General Public License for more details.
16 *
17 * You should have received a copy of the GNU General Public License
18 * along with Foobar. If not, see <http://www.gnu.org/licenses/>.
19 */
20 
21 
22 #include <stdio.h>
23 #include <stdlib.h>
24 #include "dspl.h"
25 
26 
27 
28 
29 /******************************************************************************
30 Separate complex vector to the real and image vectors
31 --------------------------------------------------------------------------------
32 Documented: RU, EN
33 *******************************************************************************/
34 int DSPL_API cmplx2re(complex_t* x, int n, double* re, double* im)
35 {
36  int k;
37  if(!x)
38  return ERROR_PTR;
39  if(n < 1)
40  return ERROR_SIZE;
41 
42  if(re)
43  {
44  for(k = 0; k < n; k++)
45  re[k] = RE(x[k]);
46  }
47  if(im)
48  {
49  for(k = 0; k < n; k++)
50  im[k] = IM(x[k]);
51  }
52  return RES_OK;
53 }
54 
55 
56 
57 
58 
59 /******************************************************************************
60 Convert real array to the complex array.
61 --------------------------------------------------------------------------------
62 Documented: RU, EN
63 *******************************************************************************/
64 int DSPL_API re2cmplx(double* x, int n, complex_t* y)
65 {
66  int k;
67  if(!x || !y)
68  return ERROR_PTR;
69  if(n < 1)
70  return ERROR_SIZE;
71 
72  for(k = 0; k < n; k++)
73  {
74  RE(y[k]) = x[k];
75  IM(y[k]) = 0.0;
76  }
77  return RES_OK;
78 }
79 
80 
#define RE(x)
Макрос определяющий реальную часть комплексного числа.
Definition: dspl.h:77
#define ERROR_PTR
Ошибка указателя. Данная ошибка означает, что один из обязательных указателей (память под который дол...
Definition: dspl.h:142
#define ERROR_SIZE
Ошибка при передаче размера массива. Данная ошибка возникает когда помимо указателя на массив входных...
Definition: dspl.h:150
int re2cmplx(double *x, int n, complex_t *y)
Преобразование массива вещественных данных в массив комплексных данных.
Definition: complex.c:64
int cmplx2re(complex_t *x, int n, double *re, double *im)
Преобразование массива комплексных данных в два массива вещественных данных, содержащих реальную и мн...
Definition: complex.c:34
double complex_t[2]
Описание комплексного типа данных.
Definition: dspl.h:41
#define RES_OK
Функция завершилась корректно. Ошибки отсутствуют.
Definition: dspl.h:94
#define IM(x)
Макрос определяющий мнимую часть комплексного числа.
Definition: dspl.h:78