1 | /*
|
---|
2 | * vmg - a versatile multigrid solver
|
---|
3 | * Copyright (C) 2012 Institute for Numerical Simulation, University of Bonn
|
---|
4 | *
|
---|
5 | * vmg is free software: you can redistribute it and/or modify
|
---|
6 | * it under the terms of the GNU General Public License as published by
|
---|
7 | * the Free Software Foundation, either version 3 of the License, or
|
---|
8 | * (at your option) any later version.
|
---|
9 | *
|
---|
10 | * vmg is distributed in the hope that it will be useful,
|
---|
11 | * but WITHOUT ANY WARRANTY; without even the implied warranty of
|
---|
12 | * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
---|
13 | * GNU General Public License for more details.
|
---|
14 | *
|
---|
15 | * You should have received a copy of the GNU General Public License
|
---|
16 | * along with this program. If not, see <http://www.gnu.org/licenses/>.
|
---|
17 | */
|
---|
18 |
|
---|
19 | /**
|
---|
20 | * @file polynomial.hpp
|
---|
21 | * @author Julian Iseringhausen <isering@ins.uni-bonn.de>
|
---|
22 | * @date Mon Nov 21 13:26:56 2011
|
---|
23 | *
|
---|
24 | * @brief Class to represent and evaluate a polynomial.
|
---|
25 | *
|
---|
26 | */
|
---|
27 |
|
---|
28 | #ifndef POLYNOMIAL_HPP_
|
---|
29 | #define POLYNOMIAL_HPP_
|
---|
30 |
|
---|
31 | #include <cstdarg>
|
---|
32 | #include <vector>
|
---|
33 |
|
---|
34 | namespace VMG
|
---|
35 | {
|
---|
36 |
|
---|
37 | class Polynomial
|
---|
38 | {
|
---|
39 | public:
|
---|
40 | Polynomial()
|
---|
41 | {}
|
---|
42 |
|
---|
43 | Polynomial(const int& degree,...) :
|
---|
44 | coeff(degree+1)
|
---|
45 | {
|
---|
46 | va_list vl;
|
---|
47 | va_start(vl, degree);
|
---|
48 | for (int i=0; i<=degree; ++i)
|
---|
49 | coeff[i] = va_arg(vl, vmg_float);
|
---|
50 | va_end(vl);
|
---|
51 | }
|
---|
52 |
|
---|
53 | Polynomial(const Polynomial& rhs) :
|
---|
54 | coeff(rhs.coeff.size())
|
---|
55 | {
|
---|
56 | for (unsigned int i=0; i<coeff.size(); ++i)
|
---|
57 | coeff[i] = rhs.coeff[i];
|
---|
58 | }
|
---|
59 |
|
---|
60 | vmg_float operator()(const vmg_float& val) const
|
---|
61 | {
|
---|
62 | switch (coeff.size())
|
---|
63 | {
|
---|
64 | case 1:
|
---|
65 | return coeff[0];
|
---|
66 | case 2:
|
---|
67 | return coeff[1] * val + coeff[0];
|
---|
68 | case 3:
|
---|
69 | return (coeff[2] * val + coeff[1]) * val + coeff[0];
|
---|
70 | case 4:
|
---|
71 | return ((coeff[3] * val + coeff[2]) * val + coeff[1]) * val + coeff[0];
|
---|
72 | case 5:
|
---|
73 | return (((coeff[4] * val + coeff[3]) * val + coeff[2]) * val + coeff[1]) * val + coeff[0];
|
---|
74 | case 6:
|
---|
75 | return ((((coeff[5] * val + coeff[4]) * val + coeff[3]) * val + coeff[2]) * val + coeff[1]) * val + coeff[0];
|
---|
76 | case 7:
|
---|
77 | return (((((coeff[6]*val+coeff[5])*val+coeff[4])*val+coeff[3])*val+coeff[2])*val+coeff[1])*val+coeff[0];
|
---|
78 | case 8:
|
---|
79 | return ((((((coeff[7]*val+coeff[6])*val+coeff[5])*val+coeff[4])*val+coeff[3])*val+coeff[2])*val+coeff[1])*val+coeff[0];
|
---|
80 | default:
|
---|
81 | vmg_float result = coeff.back();
|
---|
82 | for (int i=coeff.size()-2; i>=0; --i)
|
---|
83 | result = coeff[i] + result * val;
|
---|
84 | return result;
|
---|
85 | }
|
---|
86 |
|
---|
87 | }
|
---|
88 |
|
---|
89 | Polynomial& operator=(const Polynomial& rhs)
|
---|
90 | {
|
---|
91 | coeff.resize(rhs.coeff.size());
|
---|
92 | for (unsigned int i=0; i<coeff.size(); ++i)
|
---|
93 | coeff[i] = rhs.coeff[i];
|
---|
94 |
|
---|
95 | return *this;
|
---|
96 | }
|
---|
97 |
|
---|
98 | private:
|
---|
99 | std::vector<vmg_float> coeff;
|
---|
100 | };
|
---|
101 |
|
---|
102 | }
|
---|
103 |
|
---|
104 | #endif /* POLYNOMIAL_HPP_ */
|
---|