libdspl-2.0
Digital Signal Processing Algorithm Library
polyval.c
1/*
2* Copyright (c) 2015-2022 Sergey Bakhurin
3* Digital Signal Processing Library [http://dsplib.org]
4*
5* This file is part of libdspl-2.0.
6*
7* is free software: you can redistribute it and/or modify
8* it under the terms of the GNU Lesser 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 Lesser General Public License
18* along with Foobar. If not, see <http://www.gnu.org/licenses/>.
19*/
20
21#include <stdlib.h>
22#include <string.h>
23#include <math.h>
24#include "dspl.h"
25
26
27#ifdef DOXYGEN_ENGLISH
28
29#endif
30#ifdef DOXYGEN_RUSSIAN
76#endif
77int DSPL_API polyval(double* a, int ord, double* x, int n, double* y)
78{
79 int k, m;
80
81 if(!a || !x || !y)
82 return ERROR_PTR;
83 if(ord<0)
84 return ERROR_POLY_ORD;
85 if(n<1)
86 return ERROR_SIZE;
87
88 for(k = 0; k < n; k++)
89 {
90 y[k] = a[ord];
91 for(m = ord-1; m>-1; m--)
92 y[k] = y[k]*x[k] + a[m];
93 }
94 return RES_OK;
95}
96
#define RES_OK
The function completed correctly. No errors.
Definition: dspl.h:558
#define ERROR_PTR
Pointer error. This error means that one of the required pointers (memory to be allocated for) is tra...
Definition: dspl.h:610
#define ERROR_POLY_ORD
The polynomial order is incorrect. The order of the polynomial must be a positive integer.
Definition: dspl.h:609
#define ERROR_SIZE
Error array size. This error occurs when in addition to the pointer the wrong input is passed to the ...
Definition: dspl.h:618