]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Doc update and use the GeometryInfo<dim> class in two places to enhance readability.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Mon, 25 May 1998 22:31:20 +0000 (22:31 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Mon, 25 May 1998 22:31:20 +0000 (22:31 +0000)
git-svn-id: https://svn.dealii.org/trunk@363 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/include/fe/fe.h

index 523e5fc65f4a5c3f22154b393d816aa447c9f46b..ad4212407d626278a56838a71feb4c739d9302ae 100644 (file)
@@ -7,6 +7,7 @@
 #include <base/exceptions.h>
 #include <grid/point.h>
 #include <grid/dof.h>
+#include <grid/geometry_info.h>
 #include <lac/dfmatrix.h>
 
 
@@ -290,7 +291,7 @@ struct FiniteElementBase : public FiniteElementData<dim> {
                                      * matrix are discarded and will not fill
                                      * up the transfer matrix.
                                      */
-    dFMatrix restriction[(1<<dim)];
+    dFMatrix restriction[GeometryInfo<dim>::children_per_cell];
 
                                     /**
                                      * Have #N=2^dim# matrices keeping the
@@ -303,7 +304,12 @@ struct FiniteElementBase : public FiniteElementData<dim> {
                                      * to the destination cell, i.e. the
                                      * refined one, while the column indices
                                      * are for the unrefined cell's degrees of
-                                     * freedom.
+                                     * freedom. Thus, if #u0# is the vector
+                                     * of values of degrees of freedom on the
+                                     * coarse cell, * #prolongation[i]*u0#
+                                     * yields the vector of values of the
+                                     * degrees of freedom on the #i#th child
+                                     * cell.
                                      *
                                      * Upon assembling the transfer matrix
                                      * between cells using this matrix array,
@@ -311,7 +317,7 @@ struct FiniteElementBase : public FiniteElementData<dim> {
                                      * matrix are discarded and will not fill
                                      * up the transfer matrix.
                                      */
-    dFMatrix prolongation[(1<<dim)];
+    dFMatrix prolongation[GeometryInfo<dim>::children_per_cell];
 
                                     /**
                                      * Specify the constraints which the

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.