summaryrefslogtreecommitdiffstats
path: root/src/3rdparty/libdepixelize/priv/integral.h
blob: fc1a49f273d93cc91e9cb23c5ea4e0d0410ba7d5 (plain)
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
/*  This file is part of the libdepixelize project
    Copyright (C) 2013 Vinícius dos Santos Oliveira <vini.ipsmaker@gmail.com>

    GNU Lesser General Public License Usage
    This library is free software; you can redistribute it and/or modify it
    under the terms of the GNU Lesser General Public License as published by the
    Free Software Foundation; either version 2.1 of the License, or (at your
    option) any later version.
    You should have received a copy of the GNU Lesser General Public License
    along with this library.  If not, see <http://www.gnu.org/licenses/>.

    GNU General Public License Usage
    Alternatively, this library may be used 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.
    You should have received a copy of the GNU General Public License along with
    this library.  If not, see <http://www.gnu.org/licenses/>.

    This library 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
    Lesser General Public License for more details.
*/

#ifndef LIBDEPIXELIZE_TRACER_INTEGRAL_H
#define LIBDEPIXELIZE_TRACER_INTEGRAL_H

#include <2geom/coord.h>

namespace Tracer {

/**
 * Compute the integral numerically using Gaussian Quadrature rule with
 * \p samples number of samples.
 */
template<class T, class F>
Geom::Coord integral(F f, T begin, T end, unsigned samples)
{
    T ret = 0;
    const T width = (end - begin) / samples;

    for ( unsigned i = 0 ; i != samples ; ++i )
        ret += width * f(begin + width * (i + .5));

    return ret;
}

} // namespace Tracer

#endif // LIBDEPIXELIZE_TRACER_INTEGRAL_H

/*
  Local Variables:
  mode:c++
  c-file-style:"stroustrup"
  c-file-offsets:((innamespace . 0)(inline-open . 0)(case-label . +))
  indent-tabs-mode:nil
  fill-column:99
  End:
*/
// vim: filetype=cpp:expandtab:shiftwidth=4:tabstop=8:softtabstop=4:encoding=utf-8:textwidth=99 :