From: Guido Kanschat Date: Tue, 24 May 2005 10:34:01 +0000 (+0000) Subject: use formulae from Brezzi-Fortin for number of polynomials X-Git-Tag: v8.0.0~13835 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=235ba500828e39b0d9909b15dad54b5dbba26f02;p=dealii.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; }