This file is indexed.

/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