]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
get_degree and minor clean-ups.
authorwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 6 Mar 2001 16:32:29 +0000 (16:32 +0000)
committerwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 6 Mar 2001 16:32:29 +0000 (16:32 +0000)
git-svn-id: https://svn.dealii.org/trunk@4131 0785d39b-7218-0410-832d-ea1e28bc413d

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

index a603519233840176d57336e1afe82154e35d8e13..1d65c0b48697b8f961691ad7b3ce7fa620cac5b5 100644 (file)
@@ -30,7 +30,8 @@ class FE_Q : public FiniteElement<dim>
                                      * Constructor for tensor product
                                      * polynomials of degree @p{k}.
                                      */
-    FE_Q (unsigned int k);
+    FE_Q (const unsigned int k);
+    
                                     /**
                                      * Destructor.
                                      */
@@ -49,18 +50,21 @@ class FE_Q : public FiniteElement<dim>
     virtual UpdateFlags update_each (UpdateFlags flags) const;
 
                                     /**
-                                     * Return the support points of the
-                                     * trial functions on the unit cell.
+                                     * Return the support points of
+                                     * the trial functions on the
+                                     * unit cell.
                                      *
-                                     * The order of points in
-                                     * the array matches that returned by
-                                     * the #cell->get_dof_indices# function, but:
+                                     * The order of points in the
+                                     * array matches that returned by
+                                     * the @p{cell->get_dof_indices}
+                                     * function, but:
                                      *
                                      * If the shape functions are not
-                                     * Lagrangian interpolants at some
-                                     * points, the size of the array
-                                     * will be zero. This is the standard behavior,
-                                     * if the function is not overloaded.
+                                     * Lagrangian interpolants at
+                                     * some points, the size of the
+                                     * array will be zero. This is
+                                     * the standard behavior, if the
+                                     * function is not overloaded.
                                      */
     virtual void get_unit_support_points (std::vector<Point<dim> > &) const;    
     
@@ -69,15 +73,17 @@ class FE_Q : public FiniteElement<dim>
                                      * the trial functions on the
                                      * first face of the unit cell.
                                      *
-                                     * The order of points in
-                                     * the array matches that returned by
-                                     * the #cell->get_dof_indices# function.
+                                     * The order of points in the
+                                     * array matches that returned by
+                                     * the @p{cell->get_dof_indices}
+                                     * function.
                                      *
                                      * If the shape functions are not
-                                     * Lagrangian interpolants at some
-                                     * points, the size of the array
-                                     * will be zero. This is the standard behavior,
-                                     * if the function is not overloaded.
+                                     * Lagrangian interpolants at
+                                     * some points, the size of the
+                                     * array will be zero. This is
+                                     * the standard behavior, if the
+                                     * function is not overloaded.
                                      */
     virtual void get_unit_face_support_points (std::vector<Point<dim-1> > &) const;
 
@@ -85,8 +91,15 @@ class FE_Q : public FiniteElement<dim>
                                      * Read-only access to the
                                      * renumber vector.
                                      */
-    const std::vector<unsigned int> &get_renumber() const;
+    const std::vector<unsigned int> & get_renumber() const;
 
+                                    /**
+                                     * Return the polynomial degree
+                                     * of this finite element,
+                                     * i.e. the value passed to the
+                                     * constructor.
+                                     */
+    unsigned int get_degree () const;
     
                                     /**
                                      * Determine an estimate for the
@@ -110,7 +123,7 @@ class FE_Q : public FiniteElement<dim>
                                      * This function is needed by the
                                      * constructors of @p{FESystem}.
                                      */
-    virtual FiniteElement<dim> *clone() const;
+    virtual FiniteElement<dim> * clone() const;
   
                                     /**
                                      * Prepare internal data
@@ -175,7 +188,7 @@ class FE_Q : public FiniteElement<dim>
                                      * be passed to the constructor of
                                      * @p{FiniteElementData}.
                                      */
-    static std::vector<unsigned int> get_dpo_vector(unsigned int degree);
+    static std::vector<unsigned int> get_dpo_vector(const unsigned int degree);
     
                                     /**
                                      * Map tensor product data to shape
@@ -212,8 +225,9 @@ class FE_Q : public FiniteElement<dim>
                                   std::vector<unsigned int>& numbering);
 
                                     /**
-                                     * Map tensor product data to shape function
-                                     * numbering on first face.
+                                     * Map tensor product data to
+                                     * shape function numbering on
+                                     * first face.
                                      *
                                      * This function does the same as
                                      * @p{build_renumbering}, only on
@@ -237,16 +251,19 @@ class FE_Q : public FiniteElement<dim>
                                      * Degree of the polynomials.
                                      */  
     const unsigned int degree;
+
                                     /**
                                      * Mapping from lexicographic to
                                      * shape function numbering.
                                      */
     std::vector<unsigned int> renumber;
+    
                                     /**
                                      * Mapping from lexicographic to
                                      * shape function numbering on first face.
                                      */
     std::vector<unsigned int> face_renumber;
+    
                                     /**
                                      * Vector of one-dimensional
                                      * polynomials used.
@@ -272,6 +289,7 @@ class FE_Q : public FiniteElement<dim>
                                          * values for each quadrature point.
                                          */
        std::vector<std::vector<double> > shape_values;
+
                                         /**
                                          * Array with shape function gradients
                                          * in quadrature points. There is one
@@ -288,13 +306,13 @@ class FE_Q : public FiniteElement<dim>
   
                                     /**
                                      * Allows @p{MappingQ} class to
-                                     * access to build_renumbering
+                                     * access to @p{build_renumbering}
                                      * function.
                                      */
     friend class MappingQ<dim>;
 };
 
-#endif
+
 
 
 template<int dim>
@@ -304,3 +322,8 @@ FE_Q<dim>::get_renumber() const
 {
   return renumber;
 }
+
+
+
+#endif
+

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.