X-Git-Url: https://git.pterodactylus.net/?a=blobdiff_plain;f=ETL%2FETL%2F_calculus.h;fp=ETL%2FETL%2F_calculus.h;h=45c8ed48960378df1a83a173a9d958fc61f76b47;hb=a095981e18cc37a8ecc7cd237cc22b9c10329264;hp=0000000000000000000000000000000000000000;hpb=9459638ad6797b8139f1e9f0715c96076dbf0890;p=synfig.git

diff --git a/ETL/ETL/_calculus.h b/ETL/ETL/_calculus.h
new file mode 100644
index 0000000..45c8ed4
--- /dev/null
+++ b/ETL/ETL/_calculus.h
@@ -0,0 +1,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