From: guido Date: Tue, 24 May 2005 10:34:01 +0000 (+0000) Subject: use formulae from Brezzi-Fortin for number of polynomials X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=b98f5d46439455a282ad939d878ee542212034c3;p=dealii-svn.git use formulae from Brezzi-Fortin for number of polynomials git-svn-id: https://svn.dealii.org/trunk@10720 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/base/source/polynomials_bdm.cc b/deal.II/base/source/polynomials_bdm.cc index 4213d9ffea..83574b543d 100644 --- a/deal.II/base/source/polynomials_bdm.cc +++ b/deal.II/base/source/polynomials_bdm.cc @@ -196,12 +196,9 @@ template unsigned int PolynomialsBDM::compute_n_pols(unsigned int k) { - if (dim == 1) - return PolynomialSpace::compute_n_pols(k)+1; - if (dim == 2) - return 2*PolynomialSpace::compute_n_pols(k)+2; - if (dim == 3) - return 3*PolynomialSpace::compute_n_pols(k)+3*(k+1); + if (dim == 1) return k+1; + if (dim == 2) return (k+1)*(k+2)+2; + if (dim == 3) return ((k+1)*(k+2)*(k+3))/2+3*(k+1); Assert(false, ExcNotImplemented()); return 0; }