From ba01c42e376755d30d6ce24be34080d217887a00 Mon Sep 17 00:00:00 2001 From: hartmann Date: Mon, 27 Jun 2005 07:49:00 +0000 Subject: [PATCH] Use LagrangeEquidistant::generate_complete_basis instead of creating this basis by hand. git-svn-id: https://svn.dealii.org/trunk@10949 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/fe/fe_q.cc | 9 ++++----- 1 file changed, 4 insertions(+), 5 deletions(-) diff --git a/deal.II/deal.II/source/fe/fe_q.cc b/deal.II/deal.II/source/fe/fe_q.cc index 552316d2a9..d9c52f7f49 100644 --- a/deal.II/deal.II/source/fe/fe_q.cc +++ b/deal.II/deal.II/source/fe/fe_q.cc @@ -689,11 +689,10 @@ FE_Q<3>::initialize_constraints () // destination (child) and source (mother) // dofs. const unsigned int pnts=(this->degree+1)*(this->degree+1); - std::vector v; - for (unsigned int i=0;i<=this->degree;++i) - v.push_back(Polynomials::LagrangeEquidistant(this->degree,i)); - - const TensorProductPolynomials face_polynomials(v); + const std::vector > polynomial_basis= + Polynomials::LagrangeEquidistant::generate_complete_basis(this->degree); + + const TensorProductPolynomials face_polynomials(polynomial_basis); this->interface_constraints .TableBase<2,double>::reinit (this->interface_constraints_size()); -- 2.39.5