]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Update documentation of FE_Q constructor. 12661/head
authorWolfgang Bangerth <bangerth@colostate.edu>
Mon, 16 Aug 2021 23:53:09 +0000 (17:53 -0600)
committerWolfgang Bangerth <bangerth@colostate.edu>
Tue, 17 Aug 2021 19:18:31 +0000 (13:18 -0600)
include/deal.II/fe/fe_q.h

index b77dcef2f6bde9a4e5c6fd3a04f692596fe7b46e..f451b7df8d63e8827cca8bf31bc04688ac8d1b90 100644 (file)
@@ -559,11 +559,16 @@ public:
    * Constructor for tensor product polynomials with support points @p points
    * based on a one-dimensional quadrature formula. The degree of the finite
    * element is <tt>points.size()-1</tt>. Note that the first point has to be
-   * 0 and the last one 1. Constructing
+   * 0 and the last one 1, to ensure that the resulting element is continuous.
+   * As a consequence, only *certain* quadrature formulas make sense as
+   * argument to this constructor. Specifically, constructing
    * <tt>FE_Q<dim>(QGaussLobatto<1>(fe_degree+1))</tt> is equivalent to the
    * constructor that specifies the polynomial degree only. For selecting
-   * equidistant nodes at <tt>fe_degree > 2</tt>, construct
-   * <tt>FE_Q<dim>(QIterated<1>(QTrapezoid<1>(),fe_degree))</tt>.
+   * equidistant nodes at <tt>fe_degree > 2</tt>, use
+   * <tt>FE_Q<dim>(QIterated<1>(QTrapezoid<1>(),fe_degree))</tt>, though the
+   * use of equidistant node points is known to lead to poorly conditioned
+   * linear systems with poor interpolation properties for high polynomial
+   * degrees (which is why the default uses the QGaussLobatto points).
    *
    * The *space* created by this constructor is the
    * same as if you had called `FE_Q<dim>(points.size()-1)`, but what *is*

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.