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
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
|
// -*- Mode: C++; tab-width: 2; -*-
// vi: set ts=2:
//
// $Id: piecewisePolynomial.C,v 1.9 2003/08/26 09:17:51 oliver Exp $
//
#include <BALL/MATHS/piecewisePolynomial.h>
#include <cmath>
using namespace std;
namespace BALL
{
PiecewisePolynomial::PiecewisePolynomial()
: PiecewiseFunction(),
degree_(0)
{
}
PiecewisePolynomial::PiecewisePolynomial(const PiecewisePolynomial& poly)
: PiecewiseFunction(poly),
degree_(poly.degree_)
{
}
PiecewisePolynomial::PiecewisePolynomial(Size degree,
const std::vector<Interval>& intervals,
const std::vector<Coefficients>& coefficients)
: PiecewiseFunction(intervals, coefficients),
degree_(degree)
{
}
PiecewisePolynomial::~PiecewisePolynomial()
{
clear();
valid_ = false;
}
void PiecewisePolynomial::set(Size degree,
const std::vector<Interval>& intervals,
const std::vector<Coefficients>& coeffs)
{
PiecewiseFunction::set(intervals, coeffs);
degree_ = degree;
}
void PiecewisePolynomial::clear()
{
degree_ = 0;
PiecewiseFunction::clear();
}
PiecewisePolynomial& PiecewisePolynomial::operator =
(const PiecewisePolynomial& poly)
{
PiecewiseFunction::operator = (poly);
degree_ = poly.degree_;
return *this;
}
void PiecewisePolynomial::setDegree(Size degree)
{
degree_ = degree;
}
Size PiecewisePolynomial::getDegree() const
{
return degree_;
}
double PiecewisePolynomial::operator () (double x) const
{
// ?????
if (!isInRange(x))
{
Log.error() << "PiecewisePolynomial::operator () (): x is not in range"
<< endl;
return 0;
}
else
{
Position k = getIntervalIndex(x);
double val = 0.0;
double sub = intervals_[k].first;
for (Size i = 0; i < degree_; ++i)
{
val += coefficients_[k][i] * pow((x - sub), (int) i);
}
return val;
}
}
bool PiecewisePolynomial::operator == (const PiecewisePolynomial& poly)
const
{
return (PiecewiseFunction::operator == (poly)
&& (degree_ == poly.degree_));
}
void PiecewisePolynomial::dump(ostream& stream, Size /* depth */) const
{
stream << "[PiecewisePolynomial: degree " << degree_ << "]" << endl;
PiecewiseFunction::dump();
}
}
|
