Polynomial.h

Go to the documentation of this file.

#ifndef POLY_H
#define POLY_H

#include <cmath>
#include <sstream>

template <StatType T>
class NORMALIZED_POLYNOMIAL
{

//definition of member variables

private:

        string name;
        
        unsigned int dim;
        
        vector<Mdouble> coefficients; //starting with the highest order
        
        vector<Mdouble> averaged_coefficients;

//definition of member functions

public:

        NORMALIZED_POLYNOMIAL () {      
                set_name("Polynomial");
                coefficients.resize(0);
                dim = 0;
        }

        void set_polynomial(vector<Mdouble> new_coefficients, unsigned int new_dim);

        void set_polynomial(Mdouble* new_coefficients, unsigned int num_coeff, unsigned int new_dim);


private:
        void set_polynomial();

public:
        void set_name(const char* new_name) {name = new_name;}

        string get_name() {return name;}

        Mdouble get_volume();
        
        void set_average();

        void set_average_1D();

        void set_average_2D();

        Mdouble evaluate(Mdouble r);
        
        Mdouble evaluate_1D(Mdouble r);
        
        Mdouble evaluate_2D(Mdouble r);
        
        Mdouble evaluateIntegral(Mdouble a, Mdouble b, Mdouble t);
        
        Mdouble evaluateIntegral_1D(Mdouble a, Mdouble b, Mdouble t);
        
        Mdouble evaluateIntegral_2D(Mdouble a, Mdouble b, Mdouble t);

private:
        Mdouble operator[](int i) const {
          return coefficients[i];
        }

public:
        int order (void) {
          return coefficients.size()-1;
        }

        friend inline std::ostream& operator<<(std::ostream& os, const NORMALIZED_POLYNOMIAL &P) {
                unsigned int N = P.coefficients.size();
                for (unsigned int i=0; i<N; i++) {
                        if (!P[i]) continue;
                        if (P[i]>=0) os << "+";
                        os << setprecision(2) << P[i];
                        if (N-1-i>1) os << "r^" << N-1-i; 
                        else if (N-1-i==1) os << "r"; 
                }
                return os;
        }
        
};
#include "Polynomial.hcc"
#endif