00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032 #ifndef SACADO_PCE_STANDARDPOLY_HPP
00033 #define SACADO_PCE_STANDARDPOLY_HPP
00034
00035 #include <vector>
00036 #include <ostream>
00037
00038 namespace Sacado {
00039
00040 namespace PCE {
00041
00042 template <typename T>
00043 class StandardPoly {
00044 public:
00045
00047 StandardPoly(unsigned int deg);
00048
00050 StandardPoly(const std::vector<T>& coefficients);
00051
00053 StandardPoly(const StandardPoly& p);
00054
00056 ~StandardPoly();
00057
00059 StandardPoly& operator=(const StandardPoly& p);
00060
00062 unsigned int degree() const;
00063
00065 const T& coeff(unsigned int i) const;
00066
00068 T& coeff(unsigned int i);
00069
00071 const T& operator[](unsigned int i) const;
00072
00074 T& operator[](unsigned int i);
00075
00077
00080 void multiply(const T& alpha,
00081 const StandardPoly<T>& a,
00082 const StandardPoly<T>& b,
00083 const T& beta);
00084
00086
00089 void add(const T& alpha,
00090 const StandardPoly<T>& a,
00091 const T& gamma);
00092
00094
00097 void add(const T& alpha,
00098 const StandardPoly<T>& a,
00099 const T& beta,
00100 const StandardPoly<T>& b,
00101 const T& gamma);
00102
00103 void print(std::ostream& os) const;
00104
00105 protected:
00106
00108 std::vector<T> coeffs;
00109
00110 };
00111
00112 template <typename T>
00113 std::ostream& operator << (std::ostream& os, const StandardPoly<T>& p) {
00114 p.print(os);
00115 return os;
00116 }
00117
00118 }
00119
00120 }
00121
00122
00123 #include "Sacado_PCE_StandardPolyImp.hpp"
00124
00125 #endif // SACADO_PCE_STANDARDPOLY_HPP