1 /*! ========================================================================
2 ** Extended Template and Library
3 ** Calculus Functional Classes Implementation
4 ** $Id: _calculus.h,v 1.1.1.1 2005/01/04 01:31:46 darco Exp $
6 ** Copyright (c) 2002 Robert B. Quattlebaum Jr.
8 ** This package is free software; you can redistribute it and/or
9 ** modify it under the terms of the GNU General Public License as
10 ** published by the Free Software Foundation; either version 2 of
11 ** the License, or (at your option) any later version.
13 ** This package is distributed in the hope that it will be useful,
14 ** but WITHOUT ANY WARRANTY; without even the implied warranty of
15 ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
16 ** General Public License for more details.
18 ** === N O T E S ===========================================================
20 ** ========================================================================= */
22 /* === S T A R T =========================================================== */
24 #ifndef __ETL__CALCULUS_H
25 #define __ETL__CALCULUS_H
27 /* === H E A D E R S ======================================================= */
31 /* === M A C R O S ========================================================= */
34 //#define _EPSILON 0.0000001
37 #define ETL_FIXED_DERIVATIVE 1
39 /* === T Y P E D E F S ===================================================== */
41 /* === C L A S S E S & S T R U C T S ======================================= */
46 class derivative : public std::unary_function<typename T::argument_type,typename T::result_type>
49 typename T::argument_type epsilon;
51 explicit derivative(const T &x, const typename T::argument_type &epsilon=0.000001):func(x),epsilon(epsilon) { }
53 typename T::result_type
54 operator()(const typename T::argument_type &x)const
56 #ifdef ETL_FIXED_DERIVATIVE
57 return (func(x+epsilon)-func(x))/epsilon;
59 return (func(x)-func(x+epsilon))/epsilon;
65 class integral : public std::binary_function<typename T::argument_type,typename T::argument_type,typename T::result_type>
70 explicit integral(const T &x, const int &samples=500):func(x),samples(samples) { }
72 typename T::result_type
73 operator()(typename T::argument_type x,typename T::argument_type y)const
75 typename T::result_type ret=0;
77 const typename T::argument_type increment=(y-x)/i;
79 for(;i;i--,x+=increment)
80 ret+=(func(x)+func(x+increment))*increment/2;
87 /* === E N D =============================================================== */