File indexing completed on 2024-05-12 15:27:08
0001 /*************************************************************************** 0002 File : nsl_int.h 0003 Project : LabPlot 0004 Description : NSL numerical integration functions 0005 -------------------------------------------------------------------- 0006 Copyright : (C) 2016 by Stefan Gerlach (stefan.gerlach@uni.kn) 0007 0008 ***************************************************************************/ 0009 0010 /*************************************************************************** 0011 * * 0012 * This program is free software; you can redistribute it and/or modify * 0013 * it under the terms of the GNU General Public License as published by * 0014 * the Free Software Foundation; either version 2 of the License, or * 0015 * (at your option) any later version. * 0016 * * 0017 * This program is distributed in the hope that it will be useful, * 0018 * but WITHOUT ANY WARRANTY; without even the implied warranty of * 0019 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * 0020 * GNU General Public License for more details. * 0021 * * 0022 * You should have received a copy of the GNU General Public License * 0023 * along with this program; if not, write to the Free Software * 0024 * Foundation, Inc., 51 Franklin Street, Fifth Floor, * 0025 * Boston, MA 02110-1301 USA * 0026 * * 0027 ***************************************************************************/ 0028 0029 #ifndef NSL_INT_H 0030 #define NSL_INT_H 0031 0032 #include <stdlib.h> 0033 0034 #define NSL_INT_NETHOD_COUNT 4 0035 typedef enum {nsl_int_method_rectangle, nsl_int_method_trapezoid, nsl_int_method_simpson, nsl_int_method_simpson_3_8} nsl_int_method_type; 0036 extern const char* nsl_int_method_name[]; 0037 0038 /* 0039 numerical integration for non-uniform samples 0040 rectangle rule (1-point) 0041 trapezoid rule (2-point) 0042 Simpson-1/3 rule (3-point) returns number of points, abs not supported yet 0043 Simpson-3/8 rule (4-point) returns number of points, abs not supported yet 0044 abs - 0:return mathem. area, 1: return absolute area 0045 */ 0046 int nsl_int_rectangle(const double *x, double *y, const size_t n, int abs); 0047 int nsl_int_trapezoid(const double *x, double *y, const size_t n, int abs); 0048 size_t nsl_int_simpson(double *x, double *y, const size_t n, int abs); 0049 size_t nsl_int_simpson_3_8(double *x, double *y, const size_t n, int abs); 0050 0051 #endif /* NSL_INT_H */