/usr/include/ETL/_calculus.h is in etl-dev 0.04.19-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 | /*! ========================================================================
** Extended Template and Library
** Calculus Functional Classes Implementation
** $Id$
**
** Copyright (c) 2002 Robert B. Quattlebaum Jr.
** Copyright (c) 2008 Chris Moore
**
** This package is free software; you can redistribute it and/or
** modify it under the terms of the GNU General Public License as
** published by the Free Software Foundation; either version 2 of
** the License, or (at your option) any later version.
**
** This package is distributed in the hope that it will be useful,
** but WITHOUT ANY WARRANTY; without even the implied warranty of
** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
** General Public License for more details.
**
** === N O T E S ===========================================================
**
** ========================================================================= */
/* === S T A R T =========================================================== */
#ifndef __ETL__CALCULUS_H
#define __ETL__CALCULUS_H
/* === H E A D E R S ======================================================= */
#include <functional>
#include "hermite"
/* === M A C R O S ========================================================= */
//#ifndef _EPSILON
//#define _EPSILON 0.0000001
//#endif
/* === T Y P E D E F S ===================================================== */
/* === C L A S S E S & S T R U C T S ======================================= */
_ETL_BEGIN_NAMESPACE
template <typename T>
class derivative : public std::unary_function<typename T::argument_type,typename T::result_type>
{
T func;
typename T::argument_type epsilon;
public:
explicit derivative(const T &x, const typename T::argument_type &epsilon=0.000001):func(x),epsilon(epsilon) { }
typename T::result_type
operator()(const typename T::argument_type &x)const
{
return (func(x+epsilon)-func(x))/epsilon;
}
};
template <typename T>
class derivative<hermite<T> > : public std::unary_function<typename hermite<T>::argument_type,typename hermite<T>::result_type>
{
hermite<T> func;
public:
explicit derivative(const hermite<T> &x):func(x) { }
typename hermite<T>::result_type
operator()(const typename hermite<T>::argument_type &x)const
{
T a = func[0], b = func[1], c = func[2], d = func[3];
typename hermite<T>::argument_type y(1-x);
return ((b-a)*y*y + (c-b)*x*y*2 + (d-c)*x*x) * 3;
}
};
template <typename T>
class integral : public std::binary_function<typename T::argument_type,typename T::argument_type,typename T::result_type>
{
T func;
int samples;
public:
explicit integral(const T &x, const int &samples=500):func(x),samples(samples) { }
typename T::result_type
operator()(typename T::argument_type x,typename T::argument_type y)const
{
typename T::result_type ret=0;
int i=samples;
const typename T::argument_type increment=(y-x)/i;
for(;i;i--,x+=increment)
ret+=(func(x)+func(x+increment))*increment/2;
return ret;
}
};
_ETL_END_NAMESPACE
/* === E N D =============================================================== */
#endif
|