libdspl-2.0
Библиотека алгоритмов цифровой обработки сигналов
verif_cmplx.c
1 /*
2 * Copyright (c) 2015-2024 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 <string.h>
25 #include "dspl.h"
26 
27 
28 
29 
30 
31 #ifdef DOXYGEN_ENGLISH
32 
85 #endif
86 #ifdef DOXYGEN_RUSSIAN
87 
138 #endif
139 int DSPL_API verif_cmplx(complex_t* x, complex_t* y, size_t n,
140  double eps, double* err)
141 {
142 
143  complex_t d;
144  double mx, md, maxd;
145  size_t k;
146  int res;
147  if(!x || !y)
148  return ERROR_PTR;
149  if(n < 1)
150  return ERROR_SIZE;
151  if(eps <= 0.0 )
152  return ERROR_NEGATIVE;
153 
154  maxd = -100.0;
155 
156  for(k = 0; k < n; k++)
157  {
158  RE(d) = RE(x[k]) - RE(y[k]);
159  IM(d) = IM(x[k]) - IM(y[k]);
160  md = ABS(d);
161  mx = ABS(x[k]);
162  if(mx > 0.0)
163  {
164  md = md / mx;
165  if(md > maxd)
166  maxd = md;
167  }
168  }
169  if(err)
170  *err = maxd;
171 
172  if(maxd > eps)
173  res = DSPL_VERIF_FAILED;
174  else
175  res = DSPL_VERIF_SUCCESS;
176 
177  return res;
178 }
179 
#define RE(x)
Макрос определяющий реальную часть комплексного числа.
Definition: dspl.h:420
#define ERROR_PTR
Ошибка указателя. Данная ошибка означает, что один из обязательных указателей (память под который дол...
Definition: dspl.h:610
#define ERROR_SIZE
Ошибка при передаче размера массива. Данная ошибка возникает когда помимо указателя на массив входных...
Definition: dspl.h:618
double complex_t[2]
Описание комплексного типа данных.
Definition: dspl.h:86
#define IM(x)
Макрос определяющий мнимую часть комплексного числа.
Definition: dspl.h:478
#define ERROR_NEGATIVE
Отрицательный параметр. Функция возвращает данный код ошибки, когда принимает отрицательный параметр ...
Definition: dspl.h:604
int verif_cmplx(complex_t *x, complex_t *y, size_t n, double eps, double *err)
Верификация комплексных массивов
Definition: verif_cmplx.c:139