]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Use LagrangeEquidistant::generate_complete_basis instead of creating this basis by...
authorhartmann <hartmann@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 27 Jun 2005 07:49:00 +0000 (07:49 +0000)
committerhartmann <hartmann@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 27 Jun 2005 07:49:00 +0000 (07:49 +0000)
git-svn-id: https://svn.dealii.org/trunk@10949 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/fe/fe_q.cc

index 552316d2a99a294bfd71b80d3a97907cf9f6ee47..d9c52f7f49b83775935be6e848b7fa61f0492e5b 100644 (file)
@@ -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<Polynomials::LagrangeEquidistant> v;
-  for (unsigned int i=0;i<=this->degree;++i)
-    v.push_back(Polynomials::LagrangeEquidistant(this->degree,i));
-
-  const TensorProductPolynomials<dim-1> face_polynomials(v);
+  const std::vector<Polynomials::Polynomial<double> > polynomial_basis=
+    Polynomials::LagrangeEquidistant::generate_complete_basis(this->degree);
+  const TensorProductPolynomials<dim-1> face_polynomials(polynomial_basis);
 
   this->interface_constraints
     .TableBase<2,double>::reinit (this->interface_constraints_size());

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.