49 const double a = 1.0 / k;
50 const double b = 2.0 * k - 3.0;
51 const double c = k - 3.0;
55 for (
unsigned int i = 1; i <= k - 2; ++i)
57 coefficients[i] = b * coefficients_km1[i - 1] - c * coefficients_km2[i];
74std::vector<Polynomials::Polynomial<double>>
77 std::vector<Polynomials::Polynomial<double>> v;
79 for (
unsigned int i = 0; i <=
degree; ++i)
static std::vector< Polynomials::Polynomial< double > > generate_complete_basis(const unsigned int degree)
IntegratedLegendreSZ(const unsigned int p)
static std::vector< double > get_coefficients(const unsigned int k)
Polynomial(const std::vector< double > &coefficients)
std::vector< double > coefficients
unsigned int degree() const
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE