polynomials_bernstein.h

#include <deal.II/base/polynomial.h>
#include <fstream>
#include <iostream>

#ifndef dealii__polynomials_bernstein_h
#define dealii__polynomials_bernstein_h


DEAL_II_NAMESPACE_OPEN

template <typename number>
class PolynomialsBernstein : public Polynomials::Polynomial<number>
{
public:
  PolynomialsBernstein (
    const unsigned int index,
    const unsigned int degree);
};


template <typename number>
std::vector<Polynomials::Polynomial<number> >
generate_complete_bernstein_basis (
  const unsigned int degree)
{
  std::vector<Polynomials::Polynomial<number> > v;
  for (unsigned int i = 0; i < degree + 1; ++i)
    v.push_back(PolynomialsBernstein<number>(i, degree));
  return v;
}

DEAL_II_NAMESPACE_CLOSE

#endif