From: guido Date: Thu, 14 Apr 2005 09:39:23 +0000 (+0000) Subject: list numbering in documentation and allow output to dellog X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=feee71a7cc96ccbfa4aefdd932654ae68b16a03e;p=dealii-svn.git list numbering in documentation and allow output to dellog git-svn-id: https://svn.dealii.org/trunk@10483 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/base/include/base/polynomial_space.h b/deal.II/base/include/base/polynomial_space.h index 078775d5cc..08a7a43bf9 100644 --- a/deal.II/base/include/base/polynomial_space.h +++ b/deal.II/base/include/base/polynomial_space.h @@ -43,7 +43,40 @@ * dim-dimensional polynomials can be changed by using the * set_numbering() function. * - * @author Guido Kanschat, 2002, Wolfgang Bangerth, 2003, Ralf Hartmann 2004 + * The standard ordering of polynomials is that indices for the first + * space dimension vasry fastest and the last space dimension is + * slowest. In particular, if we take for simplicity the vector of + * monomials x0, x1, x2,..., + * xn, we get + * + *
+ *
1D
x0, x1,...,xn + *
2D:
x0y0, + * x1y0,..., + * xny0,
+ * x0y1, + * x1y1,..., + * xn-1y1,
+ * x0y2,... + * xn-2y2,
...
+ * x0yn-1, + * x1yn-1,
+ * x0yn + *
+ *
3D:
x0y0z0,..., + * xny0z0,
+ * x0y1z0,..., + * xn-1y1z0,
...
+ * x0ynz0,
+ * x0y0z1,... + * xn-1y0z1,
...
+ * x0yn-1z1,
+ * x0y0z2,... + * xn-2y0z2,
...
+ * x0y0zn + *
+ *
+ * @author Guido Kanschat, Wolfgang Bangerth, Ralf Hartmann 2002, 2003, 2004, 2005 */ template class PolynomialSpace @@ -70,7 +103,8 @@ class PolynomialSpace * Prints the list of the indices * to out. */ - void output_indices(std::ostream &out) const; + template + void output_indices(STREAM &out) const; /** * Sets the ordering of the @@ -274,6 +308,25 @@ PolynomialSpace::degree() const return polynomials.size(); } + +template +template +void +PolynomialSpace::output_indices(STREAM &out) const +{ + unsigned int ix[dim]; + for (unsigned int i=0; i -void -PolynomialSpace::output_indices(std::ostream &out) const -{ - unsigned int ix[dim]; - for (unsigned int i=0; i void PolynomialSpace::set_numbering(