libdspl-2.0
Digital Signal Processing Algorithm Library
trapint.c
1/*
2* Copyright (c) 2015-2022 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#ifdef DOXYGEN_ENGLISH
28
29#endif
30#ifdef DOXYGEN_RUSSIAN
31
32#endif
33int DSPL_API trapint(double* x, double* y, int n, double* sum)
34{
35 int k;
36
37 if(!x || !y)
38 return ERROR_PTR;
39 if(n<2)
40 return ERROR_SIZE;
41 *sum = 0.0;
42
43 for(k = 1; k < n; k++)
44 *sum += 0.5 * (x[k] - x[k-1]) * (y[k] + y[k-1]);
45
46 return RES_OK;
47}
int DSPL_API sum(double *x, int n, double *s)
Definition: sum.c:43
#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_SIZE
Error array size. This error occurs when in addition to the pointer the wrong input is passed to the ...
Definition: dspl.h:618