]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Rename types::invalid_material_id to numbers::invalid_material_id and similar for...
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 29 Aug 2012 20:54:41 +0000 (20:54 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 29 Aug 2012 20:54:41 +0000 (20:54 +0000)
git-svn-id: https://svn.dealii.org/trunk@26165 0785d39b-7218-0410-832d-ea1e28bc413d

19 files changed:
deal.II/doc/news/changes.h
deal.II/include/deal.II/base/types.h
deal.II/include/deal.II/grid/grid_out.h
deal.II/include/deal.II/grid/tria.h
deal.II/include/deal.II/grid/tria_accessor.h
deal.II/include/deal.II/grid/tria_accessor.templates.h
deal.II/include/deal.II/grid/tria_objects.h
deal.II/include/deal.II/numerics/error_estimator.h
deal.II/include/deal.II/numerics/vector_tools.h
deal.II/include/deal.II/numerics/vector_tools.templates.h
deal.II/source/dofs/dof_handler.cc
deal.II/source/dofs/dof_tools.cc
deal.II/source/grid/grid_generator.cc
deal.II/source/grid/grid_in.cc
deal.II/source/grid/grid_out.cc
deal.II/source/grid/tria.cc
deal.II/source/grid/tria_accessor.cc
deal.II/source/hp/dof_handler.cc
deal.II/source/numerics/error_estimator.cc

index df12f77d71835e710c26bc52efcc8e70c492eec9..2ff709a695482acae85fa017f88a54f0978f0f9f 100644 (file)
@@ -38,7 +38,7 @@ and then copy it into a vector with ghost elements.
 
 <li> Changed: The argument material_id of the estimate-functions of
 the KellyErrorEstimator class is now of type types::material_id
-with default value types::invalid_material_id, instead of type
+with default value numbers::invalid_material_id, instead of type
 unsigned int with default value numbers::invalid_unsigned_int. This
 should not make a difference to most users unless you specified the
 argument's default value by hand.
@@ -63,12 +63,12 @@ used to store boundary indicators internally.
 
 <ol>
 <li>
-New: we now support parallel output using HDF5/xdmf.
+New: We now support parallel output using HDF5/xdmf.
 <br>
 (Eric Heien, Timo Heister, 2012/08/28)
 
 <li>
-New: we are now compatible with Trilinos 10.4.2, 10.8.5, and 10.12.2. See the
+New: We are now compatible with Trilinos 10.4.2, 10.8.5, and 10.12.2. See the
 readme for more information.
 <br>
 (Timo Heister, 2012/08/24)
@@ -210,10 +210,10 @@ type unsigned char. Throughout the library, we changed their datatype
 to <code>types::material_id</code>
 resp. <code>types::boundary_id</code>, both typedefs of unsigned
 char. Internal faces are now characterized by
-types::internal_face_boundary_id(<code>=static_cast@<types::boundary_id@>(-1)</code>)
+numbers::internal_face_boundary_id(<code>=static_cast@<types::boundary_id@>(-1)</code>)
 instead of 255, so we get rid of that mysterious number in the source
 code.  Material_ids are also assumed to lie in the range from 0 to
-types::invalid_material_id-1 (where <code>types::invalid_material_id =
+numbers::invalid_material_id-1 (where <code>numbers::invalid_material_id =
 static_cast<types::material_id>(-1)</code>). With this change, it is now
 much easier to extend the range of boundary or material ids, if
 needed.
index b94be7f37b340906ea9687b2d78ccf69c7c05778..a59a8f06e6e9e8f9cb3fe4b6bca3c44bab4c1d2b 100644 (file)
@@ -29,6 +29,11 @@ namespace types
                                     *
                                     * See the @ref GlossSubdomainId
                                     * "glossary" for more information.
+                                   *
+                                   * There is a special value,
+                                   * numbers::invalid_subdomain_id
+                                   * that is used to indicate an
+                                   * invalid value of this type.
                                     */
   typedef unsigned int subdomain_id;
 
@@ -62,6 +67,12 @@ namespace types
            * The type used to denote boundary indicators associated with every
            * piece of the boundary and, in the case of meshes that describe
            * manifolds in higher dimensions, associated with every cell.
+          *
+          * There is a special value, numbers::internal_face_boundary_id
+          * that is used to indicate an invalid value of this type and that
+          * is used as the boundary indicator for faces that are in the interior
+          * of the domain and therefore not part of any addressable boundary
+          * component.
            */
   typedef unsigned char boundary_id;
 
@@ -70,14 +81,12 @@ namespace types
                                    */
   typedef boundary_id boundary_id_t;
 
-          /**
-           * @deprecated Use numbers::internal_face_boundary_id
-           */
-  const boundary_id internal_face_boundary_id = static_cast<boundary_id>(-1);
-
           /**
            * The type used to denote material indicators associated with every
            * cell.
+          *
+          * There is a special value, numbers::invalid_material_id
+          * that is used to indicate an invalid value of this type.
            */
   typedef unsigned char material_id;
 
@@ -86,11 +95,6 @@ namespace types
                                    */
   typedef material_id material_id_t;
 
-          /**
-          * @deprecated Use numbers::invalid_material_id
-           */
-  const material_id invalid_material_id = static_cast<material_id>(-1);
-
 }
 
 
index 824d896e5fea04524ca11f2617cc741ba3412d31..4fdc09c9429942b0227cf12a08b2178214e4588b 100644 (file)
@@ -1545,7 +1545,7 @@ class GridOut
                                       * these faces are explicitly printed
                                       * in the <tt>write_*</tt> functions;
                                       * all faces with indicator
-                                      * types::internal_face_boundary_id are
+                                      * numbers::internal_face_boundary_id are
                                       * interior ones and an indicator with
                                       * value zero for faces at the boundary
                                       * are considered default.
@@ -1580,7 +1580,7 @@ class GridOut
                                       * these lines are explicitly printed
                                       * in the <tt>write_*</tt> functions;
                                       * all lines with indicator
-                                      * types::internal_face_boundary_id are
+                                      * numbers::internal_face_boundary_id are
                                       * interior ones and an indicator with
                                       * value zero for faces at the boundary
                                       * are considered default.
index 39a678440c79afc69a20b4f1db53329f9ed6cd12..366b45cdbd00d80a0ec85a3afdde605c888161d8 100644 (file)
@@ -136,8 +136,8 @@ struct CellData
  *  number describing the boundary condition to hold on this part of the
  *  boundary. The triangulation creation function gives lines not in this
  *  list either the boundary indicator zero (if on the boundary) or
- *  types::internal_face_boundary_id (if in the interior). Explicitely giving a
- *  line the indicator types::internal_face_boundary_id will result in an error, as well as giving
+ *  numbers::internal_face_boundary_id (if in the interior). Explicitely giving a
+ *  line the indicator numbers::internal_face_boundary_id will result in an error, as well as giving
  *  an interior line a boundary indicator.
  *
  * @ingroup grid
@@ -783,12 +783,12 @@ namespace internal
  *   use is like the material id of cells).
 
  *   Boundary indicators may be in the range from zero to
- *   types::internal_face_boundary_id-1. The value
- *   types::internal_face_boundary_id is reserved to denote interior lines (in 2D)
+ *   numbers::internal_face_boundary_id-1. The value
+ *   numbers::internal_face_boundary_id is reserved to denote interior lines (in 2D)
  *   and interior lines and quads (in 3D), which do not have a
  *   boundary indicator. This way, a program can easily
  *   determine, whether such an object is at the boundary or not.
- *   Material indicators may be in the range from zero to types::invalid_material_id-1.
+ *   Material indicators may be in the range from zero to numbers::invalid_material_id-1.
  *
  *   Lines in two dimensions and quads in three dimensions inherit their
  *   boundary indicator to their children upon refinement. You should therefore
index d9a1a5a4a585821d32a503915abf38f2aad221e7..b12710e3f56c4cb819b50f3d9304fd224bd0a14b 100644 (file)
@@ -1019,7 +1019,7 @@ class TriaAccessor : public TriaAccessorBase<structdim, dim, spacedim>
                                       * object.
                                       *
                                       * If the return value is the special
-                                      * value types::internal_face_boundary_id,
+                                      * value numbers::internal_face_boundary_id,
                                       * then this object is in the
                                       * interior of the domain.
                                       */
@@ -1052,14 +1052,14 @@ class TriaAccessor : public TriaAccessorBase<structdim, dim, spacedim>
                                       * (a face not at the boundary of the
                                       * domain), or set set the boundary
                                       * indicator of an exterior face to
-                                      * types::internal_face_boundary_id
+                                      * numbers::internal_face_boundary_id
                                       * (this value is reserved for another
                                       * purpose). Algorithms may not work or
                                       * produce very confusing results if
                                       * boundary cells have a boundary
-                                      * indicator of types::internal_face_boundary_id
+                                      * indicator of numbers::internal_face_boundary_id
                                       * or if interior cells have boundary
-                                      * indicators other than types::internal_face_boundary_id.
+                                      * indicators other than numbers::internal_face_boundary_id.
                                       * Unfortunately, the current object
                                       * has no means of finding out whether it
                                       * really is at the boundary of the
@@ -1930,7 +1930,7 @@ class TriaAccessor<0, 1, spacedim>
                                       * boundary indicator one.
                                       *
                                       * If the return value is the special
-                                      * value types::internal_face_boundary_id,
+                                      * value numbers::internal_face_boundary_id,
                                       * then this object is in the
                                       * interior of the domain.
                                       */
@@ -2054,14 +2054,14 @@ class TriaAccessor<0, 1, spacedim>
                                       * (a face not at the boundary of the
                                       * domain), or set set the boundary
                                       * indicator of an exterior face to
-                                      *  types::internal_face_boundary_id
+                                      *  numbers::internal_face_boundary_id
                                       * (this value is reserved for another
                                       * purpose). Algorithms may not work or
                                       * produce very confusing results if
                                       * boundary cells have a boundary
-                                      * indicator of types::internal_face_boundary_id
+                                      * indicator of numbers::internal_face_boundary_id
                                       * or if interior cells have boundary
-                                      * indicators other than types::internal_face_boundary_id.
+                                      * indicators other than numbers::internal_face_boundary_id.
                                       * Unfortunately, the current object
                                       * has no means of finding out whether it
                                       * really is at the boundary of the
index 47c48cb45b93e7424ea344f2a1b9471071d9c1c8..fd90e894a4e3377109a3f5fe2a35bce674e33abd 100644 (file)
@@ -1878,7 +1878,7 @@ TriaAccessor<structdim, dim, spacedim>::at_boundary () const
 {
                                    // error checking is done
                                    // in boundary_indicator()
-  return (boundary_indicator() != types::internal_face_boundary_id);
+  return (boundary_indicator() != numbers::internal_face_boundary_id);
 }
 
 
@@ -2252,7 +2252,7 @@ TriaAccessor<0, 1, spacedim>::boundary_indicator () const
       }
 
       default:
-            return types::internal_face_boundary_id;
+            return numbers::internal_face_boundary_id;
     }
 
 }
index 73d369e317f39e06c958e8ea96d0547355266550..f64e7fd589b05764daeb4ae9b2f7dd1874f053cd 100644 (file)
@@ -174,7 +174,7 @@ namespace internal
                                           * example, in one dimension, this field
                                           * stores the material id of a line, which
                                           * is a number between 0 and
-                                          * types::invalid_material_id-1. In more
+                                          * numbers::invalid_material_id-1. In more
                                           * than one dimension, lines have no
                                           * material id, but they may be at the
                                           * boundary; then, we store the
@@ -184,8 +184,8 @@ namespace internal
                                           * boundary conditions hold on this
                                           * part. The boundary indicator also
                                           * is a number between zero and
-                                          * types::internal_face_boundary_id-1;
-                                          * the id types::internal_face_boundary_id
+                                          * numbers::internal_face_boundary_id-1;
+                                          * the id numbers::internal_face_boundary_id
                                           * is reserved for lines
                                           * in the interior and may be used
                                           * to check whether a line is at the
@@ -712,7 +712,7 @@ namespace internal
     inline
     TriaObjects<G>::BoundaryOrMaterialId::BoundaryOrMaterialId ()
     {
-      material_id = types::invalid_material_id;
+      material_id = numbers::invalid_material_id;
     }
 
 
index adfec866bbdc588a6a8d96ccfb29e4b6d0143a2b..e856482be96af875f59d51bb361404bef4eb4502 100644 (file)
@@ -300,7 +300,7 @@ class KellyErrorEstimator
                                       *
                                       * The @p material_id parameter has a similar
                                       * meaning: if not set to its default value
-                                      * (which is types::invalid_material_id),
+                                      * (which is numbers::invalid_material_id),
                                       * it means that indicators will only be
                                       * computed for cells with this particular
                                       * material id.
@@ -349,7 +349,7 @@ class KellyErrorEstimator
                           const Function<spacedim>     *coefficients   = 0,
                           const unsigned int       n_threads = numbers::invalid_unsigned_int,
                           const types::subdomain_id subdomain_id = types::invalid_subdomain_id,
-                          const types::material_id       material_id = types::invalid_material_id);
+                          const types::material_id       material_id = numbers::invalid_material_id);
 
                                      /**
                                       * Calls the @p estimate
@@ -366,7 +366,7 @@ class KellyErrorEstimator
                           const Function<spacedim>     *coefficients   = 0,
                           const unsigned int       n_threads = multithread_info.n_default_threads,
                           const types::subdomain_id subdomain_id = types::invalid_subdomain_id,
-                          const types::material_id       material_id = types::invalid_material_id);
+                          const types::material_id       material_id = numbers::invalid_material_id);
 
                                      /**
                                       * Same function as above, but
@@ -406,7 +406,7 @@ class KellyErrorEstimator
                           const Function<spacedim>         *coefficients   = 0,
                           const unsigned int           n_threads = multithread_info.n_default_threads,
                           const types::subdomain_id subdomain_id = types::invalid_subdomain_id,
-                          const types::material_id           material_id = types::invalid_material_id);
+                          const types::material_id           material_id = numbers::invalid_material_id);
 
                                      /**
                                       * Calls the @p estimate
@@ -423,7 +423,7 @@ class KellyErrorEstimator
                           const Function<spacedim>         *coefficients   = 0,
                           const unsigned int           n_threads = multithread_info.n_default_threads,
                           const types::subdomain_id subdomain_id = types::invalid_subdomain_id,
-                          const types::material_id          material_id = types::invalid_material_id);
+                          const types::material_id          material_id = numbers::invalid_material_id);
 
 
                                      /**
@@ -443,7 +443,7 @@ class KellyErrorEstimator
                           const Function<spacedim>     *coefficients   = 0,
                           const unsigned int       n_threads = multithread_info.n_default_threads,
                           const types::subdomain_id subdomain_id = types::invalid_subdomain_id,
-                          const types::material_id       material_id = types::invalid_material_id);
+                          const types::material_id       material_id = numbers::invalid_material_id);
 
 
                                      /**
@@ -462,7 +462,7 @@ class KellyErrorEstimator
                           const Function<spacedim>     *coefficients   = 0,
                           const unsigned int       n_threads = multithread_info.n_default_threads,
                           const types::subdomain_id subdomain_id = types::invalid_subdomain_id,
-                          const types::material_id       material_id = types::invalid_material_id);
+                          const types::material_id       material_id = numbers::invalid_material_id);
 
 
                                      /**
@@ -482,7 +482,7 @@ class KellyErrorEstimator
                           const Function<spacedim>         *coefficients   = 0,
                           const unsigned int           n_threads = multithread_info.n_default_threads,
                           const types::subdomain_id subdomain_id = types::invalid_subdomain_id,
-                          const types::material_id           material_id = types::invalid_material_id);
+                          const types::material_id           material_id = numbers::invalid_material_id);
 
 
                                      /**
@@ -501,7 +501,7 @@ class KellyErrorEstimator
                           const Function<spacedim>*    coefficients   = 0,
                           const unsigned int           n_threads = multithread_info.n_default_threads,
                           const types::subdomain_id subdomain_id = types::invalid_subdomain_id,
-                          const types::material_id           material_id = types::invalid_material_id);
+                          const types::material_id           material_id = numbers::invalid_material_id);
 
 
                                      /**
@@ -612,7 +612,7 @@ class KellyErrorEstimator<1,spacedim>
                           const Function<spacedim>     *coefficients   = 0,
                           const unsigned int       n_threads = multithread_info.n_default_threads,
                           const types::subdomain_id subdomain_id = types::invalid_subdomain_id,
-                          const types::material_id       material_id = types::invalid_material_id);
+                          const types::material_id       material_id = numbers::invalid_material_id);
 
                                      /**
                                       * Calls the @p estimate
@@ -629,7 +629,7 @@ class KellyErrorEstimator<1,spacedim>
                           const Function<spacedim>     *coefficients   = 0,
                           const unsigned int       n_threads = multithread_info.n_default_threads,
                           const types::subdomain_id subdomain_id = types::invalid_subdomain_id,
-                          const types::material_id       material_id = types::invalid_material_id);
+                          const types::material_id       material_id = numbers::invalid_material_id);
 
                                      /**
                                       * Same function as above, but
@@ -669,7 +669,7 @@ class KellyErrorEstimator<1,spacedim>
                           const Function<spacedim>         *coefficients   = 0,
                           const unsigned int           n_threads = multithread_info.n_default_threads,
                           const types::subdomain_id subdomain_id = types::invalid_subdomain_id,
-                          const types::material_id           material_id = types::invalid_material_id);
+                          const types::material_id           material_id = numbers::invalid_material_id);
 
                                      /**
                                       * Calls the @p estimate
@@ -686,7 +686,7 @@ class KellyErrorEstimator<1,spacedim>
                           const Function<spacedim>         *coefficients   = 0,
                           const unsigned int           n_threads = multithread_info.n_default_threads,
                           const types::subdomain_id subdomain_id = types::invalid_subdomain_id,
-                          const types::material_id           material_id = types::invalid_material_id);
+                          const types::material_id           material_id = numbers::invalid_material_id);
 
 
                                      /**
@@ -706,7 +706,7 @@ class KellyErrorEstimator<1,spacedim>
                           const Function<spacedim>     *coefficients   = 0,
                           const unsigned int       n_threads = multithread_info.n_default_threads,
                           const types::subdomain_id subdomain_id = types::invalid_subdomain_id,
-                          const types::material_id       material_id = types::invalid_material_id);
+                          const types::material_id       material_id = numbers::invalid_material_id);
 
 
                                      /**
@@ -725,7 +725,7 @@ class KellyErrorEstimator<1,spacedim>
                           const Function<spacedim>     *coefficients   = 0,
                           const unsigned int       n_threads = multithread_info.n_default_threads,
                           const types::subdomain_id subdomain_id = types::invalid_subdomain_id,
-                          const types::material_id       material_id = types::invalid_material_id);
+                          const types::material_id       material_id = numbers::invalid_material_id);
 
 
                                      /**
@@ -745,7 +745,7 @@ class KellyErrorEstimator<1,spacedim>
                           const Function<spacedim>         *coefficients   = 0,
                           const unsigned int           n_threads = multithread_info.n_default_threads,
                           const types::subdomain_id subdomain_id = types::invalid_subdomain_id,
-                          const types::material_id           material_id = types::invalid_material_id);
+                          const types::material_id           material_id = numbers::invalid_material_id);
 
 
                                      /**
@@ -764,7 +764,7 @@ class KellyErrorEstimator<1,spacedim>
                           const Function<spacedim>         *coefficients   = 0,
                           const unsigned int           n_threads = multithread_info.n_default_threads,
                           const types::subdomain_id subdomain_id = types::invalid_subdomain_id,
-                          const types::material_id           material_id = types::invalid_material_id);
+                          const types::material_id           material_id = numbers::invalid_material_id);
 
                                      /**
                                       * Exception
index f6a1387560ce0487deb6e9a6a955923f61beed3e..157c56b0c59dc2afa3b275bdc2e8f437f0467f99 100644 (file)
@@ -121,8 +121,8 @@ class ConstraintMatrix;
  *   the trace of the function to the boundary is done with the
  *   VectorTools::project_boundary_values() (see below) function,
  *   which is called with a map of boundary functions FunctioMap in
- *   which all boundary indicators from zero to types::internal_face_boundary_id-1
- *   (types::internal_face_boundary_id is used for other purposes,
+ *   which all boundary indicators from zero to numbers::internal_face_boundary_id-1
+ *   (numbers::internal_face_boundary_id is used for other purposes,
  *   see the Triangulation class documentation) point
  *   to the function to be projected. The projection to the boundary
  *   takes place using a second quadrature formula on the boundary
@@ -694,7 +694,7 @@ namespace VectorTools
                                     * keys of this map
                                     * correspond to the number
                                     * @p boundary_indicator of the
-                                    * face.  types::internal_face_boundary_id
+                                    * face.  numbers::internal_face_boundary_id
                                     * is an illegal value for this key since
                                     * it is reserved for interior faces.
                                     *
@@ -841,7 +841,7 @@ namespace VectorTools
                                     * The parameter @p boundary_component
                                     * corresponds to the number @p
                                     * boundary_indicator of the face.
-                                    * types::internal_face_boundary_id
+                                    * numbers::internal_face_boundary_id
                                     * is an illegal value, since it is
                                     * reserved for interior faces.
                                     *
@@ -1147,7 +1147,7 @@ namespace VectorTools
                                     * The parameter @p boundary_component
                                     * corresponds to the number
                                     * @p boundary_indicator of the face.
-                                    * types::internal_face_boundary_id
+                                    * numbers::internal_face_boundary_id
                                     * is an illegal value, since it is
                                     * reserved for interior faces.
                                     *
@@ -1242,7 +1242,7 @@ namespace VectorTools
                                     * The parameter @p boundary_component
                                     * corresponds to the number
                                     * @p boundary_indicator of the face.
-                                    *  types::internal_face_boundary_id
+                                    *  numbers::internal_face_boundary_id
                                     * is an illegal value, since it is
                                     * reserved for interior faces.
                                     *
index 2df4c30f05ff8302998d7219d82f3be636275f1f..bea6de958ed0df221069784d7d7ea7e7127ac9be 100644 (file)
@@ -308,7 +308,7 @@ namespace VectorTools
 
                                              // count, how often we have
                                              // added to this dof
-            Assert (touch_count[local_dof_indices[j]] < types::internal_face_boundary_id,
+            Assert (touch_count[local_dof_indices[j]] < numbers::internal_face_boundary_id,
                     ExcInternalError());
             ++touch_count[local_dof_indices[j]];
           };
@@ -571,7 +571,7 @@ namespace VectorTools
                                            // function to hold on
                                            // all parts of the boundary
           typename FunctionMap<spacedim>::type boundary_functions;
-          for (types::boundary_id c=0; c<types::internal_face_boundary_id; ++c)
+          for (types::boundary_id c=0; c<numbers::internal_face_boundary_id; ++c)
             boundary_functions[c] = &function;
           project_boundary_values (dof, boundary_functions, q_boundary,
                                    boundary_values);
@@ -1588,7 +1588,7 @@ namespace VectorTools
       if (function_map.size() == 0)
         return;
 
-      Assert (function_map.find(types::internal_face_boundary_id) == function_map.end(),
+      Assert (function_map.find(numbers::internal_face_boundary_id) == function_map.end(),
               ExcInvalidBoundaryIndicator());
 
       const unsigned int        n_components = DoFTools::n_components(dof);
index 9e466cd2d0bb311e284a669b5f3f255daf9a91d7..7a793fd719597c1407af57e894d007ebe7e74c1d 100644 (file)
@@ -1579,7 +1579,7 @@ template<int dim, int spacedim>
 unsigned int
 DoFHandler<dim,spacedim>::n_boundary_dofs (const FunctionMap &boundary_indicators) const
 {
-  Assert (boundary_indicators.find(types::internal_face_boundary_id) == boundary_indicators.end(),
+  Assert (boundary_indicators.find(numbers::internal_face_boundary_id) == boundary_indicators.end(),
           ExcInvalidBoundaryIndicator());
 
   std::set<int> boundary_dofs;
@@ -1605,7 +1605,7 @@ template<int dim, int spacedim>
 unsigned int
 DoFHandler<dim,spacedim>::n_boundary_dofs (const std::set<types::boundary_id> &boundary_indicators) const
 {
-  Assert (boundary_indicators.find (types::internal_face_boundary_id) == boundary_indicators.end(),
+  Assert (boundary_indicators.find (numbers::internal_face_boundary_id) == boundary_indicators.end(),
           ExcInvalidBoundaryIndicator());
 
   std::set<int> boundary_dofs;
index e23dd691a8399107fb8ca80851567498f4c01fec..6b7d7276663c6390ef8b6851b96615c7ea7cee34 100644 (file)
@@ -472,7 +472,7 @@ namespace DoFTools
     const unsigned int n_dofs = dof.n_dofs();
 
     AssertDimension (dof_to_boundary_mapping.size(), n_dofs);
-    Assert (boundary_indicators.find(types::internal_face_boundary_id) == boundary_indicators.end(),
+    Assert (boundary_indicators.find(numbers::internal_face_boundary_id) == boundary_indicators.end(),
             typename DH::ExcInvalidBoundaryIndicator());
     Assert (sparsity.n_rows() == dof.n_boundary_dofs (boundary_indicators),
             ExcDimensionMismatch (sparsity.n_rows(), dof.n_boundary_dofs (boundary_indicators)));
@@ -3948,7 +3948,7 @@ namespace DoFTools
                          const std::set<types::boundary_id> &boundary_indicators)
   {
     AssertDimension (component_select.size(), n_components(dof_handler));
-    Assert (boundary_indicators.find (types::internal_face_boundary_id) == boundary_indicators.end(),
+    Assert (boundary_indicators.find (numbers::internal_face_boundary_id) == boundary_indicators.end(),
             ExcInvalidBoundaryIndicator());
     const unsigned int dim=DH::dimension;
 
@@ -4076,7 +4076,7 @@ namespace DoFTools
                                          const std::set<types::boundary_id> &boundary_indicators)
   {
     AssertDimension (component_select.size(), n_components(dof_handler));
-    Assert (boundary_indicators.find (types::internal_face_boundary_id) == boundary_indicators.end(),
+    Assert (boundary_indicators.find (numbers::internal_face_boundary_id) == boundary_indicators.end(),
             ExcInvalidBoundaryIndicator());
 
                                      // let's see whether we have to
@@ -5875,7 +5875,7 @@ namespace DoFTools
     std::vector<unsigned int>     &mapping)
   {
     Assert (&dof_handler.get_fe() != 0, ExcNoFESelected());
-    Assert (boundary_indicators.find (types::internal_face_boundary_id) == boundary_indicators.end(),
+    Assert (boundary_indicators.find (numbers::internal_face_boundary_id) == boundary_indicators.end(),
             ExcInvalidBoundaryIndicator());
 
     mapping.clear ();
index 24a3432374292785327db8df36a486145ae64958..20136921f4e936946bac62e8e87a171a08c4405c 100644 (file)
@@ -920,12 +920,12 @@ GridGenerator::subdivided_hyper_rectangle (
                                    // create the cells
   unsigned int n_val_cells = 0;
   for (unsigned int i=0; i<n_cells; i++)
-    if (material_id[i]!=types::invalid_material_id) n_val_cells++;
+    if (material_id[i]!=numbers::invalid_material_id) n_val_cells++;
 
   std::vector<CellData<1> > cells(n_val_cells);
   unsigned int id = 0;
   for (unsigned int x=0; x<n_cells; ++x)
-    if (material_id[x] != types::invalid_material_id)
+    if (material_id[x] != numbers::invalid_material_id)
       {
         cells[id].vertices[0] = x;
         cells[id].vertices[1] = x+1;
@@ -993,14 +993,14 @@ GridGenerator::subdivided_hyper_rectangle (
   unsigned int n_val_cells = 0;
   for (unsigned int i=0; i<material_id.size(0); i++)
     for (unsigned int j=0; j<material_id.size(1); j++)
-      if (material_id[i][j] != types::invalid_material_id)
+      if (material_id[i][j] != numbers::invalid_material_id)
         n_val_cells++;
 
   std::vector<CellData<2> > cells(n_val_cells);
   unsigned int id = 0;
   for (unsigned int y=0; y<repetitions[1]; ++y)
     for (unsigned int x=0; x<repetitions[0]; ++x)
-      if (material_id[x][y]!=types::invalid_material_id)
+      if (material_id[x][y]!=numbers::invalid_material_id)
         {
           cells[id].vertices[0] = y*(repetitions[0]+1)+x;
           cells[id].vertices[1] = y*(repetitions[0]+1)+x+1;
@@ -1100,7 +1100,7 @@ GridGenerator::subdivided_hyper_rectangle (
   for (unsigned int i=0; i<material_id.size(0); i++)
     for (unsigned int j=0; j<material_id.size(1); j++)
       for (unsigned int k=0; k<material_id.size(2); k++)
-        if (material_id[i][j][k]!=types::invalid_material_id)
+        if (material_id[i][j][k]!=numbers::invalid_material_id)
           n_val_cells++;
 
   std::vector<CellData<dim> > cells(n_val_cells);
@@ -1110,7 +1110,7 @@ GridGenerator::subdivided_hyper_rectangle (
   for (unsigned int z=0; z<repetitions[2]; ++z)
     for (unsigned int y=0; y<repetitions[1]; ++y)
       for (unsigned int x=0; x<repetitions[0]; ++x)
-        if (material_id[x][y][z]!=types::invalid_material_id)
+        if (material_id[x][y][z]!=numbers::invalid_material_id)
           {
             cells[id].vertices[0] = z*n_xy + y*n_x + x;
             cells[id].vertices[1] = z*n_xy + y*n_x + x+1;
@@ -1751,7 +1751,7 @@ GridGenerator::half_hyper_ball (Triangulation<2> &tria,
     {
         for (unsigned int i=0;i<GeometryInfo<2>::faces_per_cell;++i)
         {
-            if (cell->face(i)->boundary_indicator() == types::internal_face_boundary_id)
+            if (cell->face(i)->boundary_indicator() == numbers::internal_face_boundary_id)
                 continue;
 
             // If x is zero, then this is part of the plane
index a48f534ef33a39eb6047bf6e48cb811c6f55053f..25c769810166bd5e51661f3e2d1fa09c705ff05b 100644 (file)
@@ -357,9 +357,9 @@ void GridIn<dim, spacedim>::read_ucd (std::istream &in)
           // to make sure that the cast wont fail
     Assert(material_id<= std::numeric_limits<types::material_id>::max(),
         ExcIndexRange(material_id,0,std::numeric_limits<types::material_id>::max()));
-    // we use only material_ids in the range from 0 to types::invalid_material_id-1
-    Assert(material_id < types::invalid_material_id,
-        ExcIndexRange(material_id,0,types::invalid_material_id));
+    // we use only material_ids in the range from 0 to numbers::invalid_material_id-1
+    Assert(material_id < numbers::invalid_material_id,
+        ExcIndexRange(material_id,0,numbers::invalid_material_id));
 
           cells.back().material_id = static_cast<types::material_id>(material_id);
 
@@ -387,9 +387,9 @@ void GridIn<dim, spacedim>::read_ucd (std::istream &in)
             // to make sure that the cast wont fail
             Assert(material_id<= std::numeric_limits<types::boundary_id>::max(),
                 ExcIndexRange(material_id,0,std::numeric_limits<types::boundary_id>::max()));
-            // we use only boundary_ids in the range from 0 to types::internal_face_boundary_id-1
-            Assert(material_id < types::internal_face_boundary_id,
-                ExcIndexRange(material_id,0,types::internal_face_boundary_id));
+            // we use only boundary_ids in the range from 0 to numbers::internal_face_boundary_id-1
+            Assert(material_id < numbers::internal_face_boundary_id,
+                ExcIndexRange(material_id,0,numbers::internal_face_boundary_id));
 
             subcelldata.boundary_lines.back().boundary_id
                 = static_cast<types::boundary_id>(material_id);
@@ -425,9 +425,9 @@ void GridIn<dim, spacedim>::read_ucd (std::istream &in)
               // to make sure that the cast wont fail
               Assert(material_id<= std::numeric_limits<types::boundary_id>::max(),
                   ExcIndexRange(material_id,0,std::numeric_limits<types::boundary_id>::max()));
-              // we use only boundary_ids in the range from 0 to types::internal_face_boundary_id-1
-              Assert(material_id < types::internal_face_boundary_id,
-                  ExcIndexRange(material_id,0,types::internal_face_boundary_id));
+              // we use only boundary_ids in the range from 0 to numbers::internal_face_boundary_id-1
+              Assert(material_id < numbers::internal_face_boundary_id,
+                  ExcIndexRange(material_id,0,numbers::internal_face_boundary_id));
 
               subcelldata.boundary_quads.back().boundary_id
                   = static_cast<types::boundary_id>(material_id);
@@ -1000,9 +1000,9 @@ void GridIn<dim, spacedim>::read_msh (std::istream &in)
     // to make sure that the cast wont fail
     Assert(material_id<= std::numeric_limits<types::material_id>::max(),
         ExcIndexRange(material_id,0,std::numeric_limits<types::material_id>::max()));
-    // we use only material_ids in the range from 0 to types::invalid_material_id-1
-    Assert(material_id < types::invalid_material_id,
-        ExcIndexRange(material_id,0,types::invalid_material_id));
+    // we use only material_ids in the range from 0 to numbers::invalid_material_id-1
+    Assert(material_id < numbers::invalid_material_id,
+        ExcIndexRange(material_id,0,numbers::invalid_material_id));
 
           cells.back().material_id = static_cast<types::material_id>(material_id);
 
@@ -1029,9 +1029,9 @@ void GridIn<dim, spacedim>::read_msh (std::istream &in)
       // to make sure that the cast wont fail
       Assert(material_id<= std::numeric_limits<types::boundary_id>::max(),
           ExcIndexRange(material_id,0,std::numeric_limits<types::boundary_id>::max()));
-      // we use only boundary_ids in the range from 0 to types::internal_face_boundary_id-1
-      Assert(material_id < types::internal_face_boundary_id,
-          ExcIndexRange(material_id,0,types::internal_face_boundary_id));
+      // we use only boundary_ids in the range from 0 to numbers::internal_face_boundary_id-1
+      Assert(material_id < numbers::internal_face_boundary_id,
+          ExcIndexRange(material_id,0,numbers::internal_face_boundary_id));
 
             subcelldata.boundary_lines.back().boundary_id
                 = static_cast<types::boundary_id>(material_id);
@@ -1067,9 +1067,9 @@ void GridIn<dim, spacedim>::read_msh (std::istream &in)
               // to make sure that the cast wont fail
               Assert(material_id<= std::numeric_limits<types::boundary_id>::max(),
                   ExcIndexRange(material_id,0,std::numeric_limits<types::boundary_id>::max()));
-              // we use only boundary_ids in the range from 0 to types::internal_face_boundary_id-1
-              Assert(material_id < types::internal_face_boundary_id,
-                  ExcIndexRange(material_id,0,types::internal_face_boundary_id));
+              // we use only boundary_ids in the range from 0 to numbers::internal_face_boundary_id-1
+              Assert(material_id < numbers::internal_face_boundary_id,
+                  ExcIndexRange(material_id,0,numbers::internal_face_boundary_id));
 
               subcelldata.boundary_quads.back().boundary_id
                   = static_cast<types::boundary_id>(material_id);
@@ -1618,10 +1618,10 @@ void GridIn<3>::read_netcdf (const std::string &filename)
                                    // we only handle boundary
                                    // indicators that fit into an
                                    // types::boundary_id. Also, we don't
-                                   // take types::internal_face_boundary_id
+                                   // take numbers::internal_face_boundary_id
            // as it denotes an internal face
   for (unsigned int i=0; i<bmarker.size(); ++i)
-    Assert(0<=bmarker[i] && bmarker[i]<types::internal_face_boundary_id, ExcIO());
+    Assert(0<=bmarker[i] && bmarker[i]<numbers::internal_face_boundary_id, ExcIO());
 
                                    // finally we setup the boundary
                                    // information
index a58d2634a19dd7d7af6a0eac2989389172218658..56eae72ae698377bdf1b414eac353089e5d70691 100644 (file)
@@ -1163,7 +1163,7 @@ void GridOut::write_xfig (
             Triangulation<dim>::face_iterator
               face = cell->face(face_reorder[f]);
             const types::boundary_id bi = face->boundary_indicator();
-            if (bi != types::internal_face_boundary_id)
+            if (bi != numbers::internal_face_boundary_id)
               {
                                                  // Code for polyline
                 out << "2 1 "
index 31a0ef4830929427e90e54f0027cc93bb3a6b412..adef525f7f0787b165c91080844da12c0c1faa69 100644 (file)
@@ -1789,8 +1789,8 @@ namespace internal
                 if (n_adj_cells == 1)
                   line->set_boundary_indicator (0);
                 else
-                                                   // interior line -> types::internal_face_boundary_id
-                  line->set_boundary_indicator (types::internal_face_boundary_id);
+                                                   // interior line -> numbers::internal_face_boundary_id
+                  line->set_boundary_indicator (numbers::internal_face_boundary_id);
               }
 
                                              // set boundary indicators where
@@ -1824,13 +1824,13 @@ namespace internal
                                                  // assert that we only set
                                                  // boundary info once
                 AssertThrow (! (line->boundary_indicator() != 0 &&
-                                line->boundary_indicator() != types::internal_face_boundary_id),
+                                line->boundary_indicator() != numbers::internal_face_boundary_id),
                              ExcMultiplySetLineInfoOfLine(line_vertices.first,
                                                           line_vertices.second));
 
                                                  // Assert that only exterior lines
                                                  // are given a boundary indicator
-                AssertThrow (! (line->boundary_indicator() == types::internal_face_boundary_id),
+                AssertThrow (! (line->boundary_indicator() == numbers::internal_face_boundary_id),
                              ExcInteriorLineCantBeBoundary());
 
                 line->set_boundary_indicator (boundary_line->boundary_id);
@@ -2509,8 +2509,8 @@ namespace internal
                 if (n_adj_cells == 1)
                   quad->set_boundary_indicator (0);
                 else
-                                                   // interior quad -> types::internal_face_boundary_id
-                  quad->set_boundary_indicator (types::internal_face_boundary_id);
+                                                   // interior quad -> numbers::internal_face_boundary_id
+                  quad->set_boundary_indicator (numbers::internal_face_boundary_id);
               }
 
                                              /////////////////////////////////////////
@@ -2522,7 +2522,7 @@ namespace internal
                                              // as interior
             for (typename Triangulation<dim,spacedim>::line_iterator
                    line=triangulation.begin_line(); line!=triangulation.end_line(); ++line)
-              line->set_boundary_indicator (types::internal_face_boundary_id);
+              line->set_boundary_indicator (numbers::internal_face_boundary_id);
                                              // next reset all lines bounding
                                              // boundary quads as on the
                                              // boundary also. note that since
@@ -4117,7 +4117,7 @@ namespace internal
                 new_lines[l]->clear_user_data();
                 new_lines[l]->clear_children();
                                                  // interior line
-                new_lines[l]->set_boundary_indicator(types::internal_face_boundary_id);
+                new_lines[l]->set_boundary_indicator(numbers::internal_face_boundary_id);
               }
 
                                              // Now add the four (two)
@@ -6375,7 +6375,7 @@ namespace internal
                           new_lines[i]->clear_user_data();
                           new_lines[i]->clear_children();
                                                            // interior line
-                          new_lines[i]->set_boundary_indicator(types::internal_face_boundary_id);
+                          new_lines[i]->set_boundary_indicator(numbers::internal_face_boundary_id);
                         }
 
                                                        // find some space for the newly to
@@ -6394,7 +6394,7 @@ namespace internal
                           new_quads[i]->clear_user_data();
                           new_quads[i]->clear_children();
                                                            // interior quad
-                          new_quads[i]->set_boundary_indicator (types::internal_face_boundary_id);
+                          new_quads[i]->set_boundary_indicator (numbers::internal_face_boundary_id);
                                                            // set all line orientation
                                                            // flags to true by default,
                                                            // change this afterwards, if
@@ -9472,7 +9472,7 @@ void
 Triangulation<dim, spacedim>::set_boundary (const types::boundary_id number,
                                             const Boundary<dim, spacedim>& boundary_object)
 {
-  Assert(number < types::internal_face_boundary_id, ExcIndexRange(number,0,types::internal_face_boundary_id));
+  Assert(number < numbers::internal_face_boundary_id, ExcIndexRange(number,0,numbers::internal_face_boundary_id));
 
   boundary[number] = &boundary_object;
 }
@@ -9483,7 +9483,7 @@ template <int dim, int spacedim>
 void
 Triangulation<dim, spacedim>::set_boundary (const types::boundary_id number)
 {
-  Assert(number < types::internal_face_boundary_id, ExcIndexRange(number,0,types::internal_face_boundary_id));
+  Assert(number < numbers::internal_face_boundary_id, ExcIndexRange(number,0,numbers::internal_face_boundary_id));
 
   //delete the entry located at number.
   boundary.erase(number);
@@ -9495,7 +9495,7 @@ template <int dim, int spacedim>
 const Boundary<dim, spacedim> &
 Triangulation<dim, spacedim>::get_boundary (const types::boundary_id number) const
 {
-  Assert(number<types::internal_face_boundary_id, ExcIndexRange(number,0,types::internal_face_boundary_id));
+  Assert(number<numbers::internal_face_boundary_id, ExcIndexRange(number,0,numbers::internal_face_boundary_id));
 
   //look, if there is a boundary stored at
   //boundary_id number.
@@ -9539,7 +9539,7 @@ Triangulation<dim, spacedim>::get_boundary_indicators () const
     {
 //TODO: if boundary_id_t is a real integer type, this might become
 //a humongously large array...
-      std::vector<bool> bi_exists(types::internal_face_boundary_id, false);
+      std::vector<bool> bi_exists(numbers::internal_face_boundary_id, false);
       active_cell_iterator cell=begin_active();
       for (; cell!=end(); ++cell)
         for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
index 972b4696d744091700069115f687b8a956bb3b33..cb713aa258a712775000cdbde6abf9c341176414 100644 (file)
@@ -1285,7 +1285,7 @@ template <int dim, int spacedim>
 void CellAccessor<dim, spacedim>::set_material_id (const types::material_id mat_id) const
 {
   Assert (this->used(), TriaAccessorExceptions::ExcCellNotUsed());
-  Assert ( mat_id < types::invalid_material_id, ExcIndexRange(mat_id,0,types::invalid_material_id));
+  Assert ( mat_id < numbers::invalid_material_id, ExcIndexRange(mat_id,0,numbers::invalid_material_id));
   this->tria->levels[this->present_level]->cells.boundary_or_material_id[this->present_index].material_id = mat_id;
 }
 
index 75f26683925ed98b433c3c4fd83e6c72d86e04a5..35a74c96705f2a12f68a6d060fe08e7a769afe8e 100644 (file)
@@ -2869,7 +2869,7 @@ template <>
   DoFHandler<dim,spacedim>::n_boundary_dofs (const FunctionMap &boundary_indicators) const
   {
     Assert (finite_elements != 0, ExcNoFESelected());
-    Assert (boundary_indicators.find(types::internal_face_boundary_id) == boundary_indicators.end(),
+    Assert (boundary_indicators.find(numbers::internal_face_boundary_id) == boundary_indicators.end(),
             ExcInvalidBoundaryIndicator());
 
                                      // same as above, but with
@@ -2905,7 +2905,7 @@ template <>
   DoFHandler<dim,spacedim>::n_boundary_dofs (const std::set<types::boundary_id> &boundary_indicators) const
   {
     Assert (finite_elements != 0, ExcNoFESelected());
-    Assert (boundary_indicators.find (types::internal_face_boundary_id) == boundary_indicators.end(),
+    Assert (boundary_indicators.find (numbers::internal_face_boundary_id) == boundary_indicators.end(),
             ExcInvalidBoundaryIndicator());
 
                                      // same as above, but with
index c598979682aa4206b9b6c3d350f8426f838d49f2..95517074b6d7edc12bfc6f7354678a4ffc554497 100644 (file)
@@ -907,7 +907,7 @@ namespace internal
                     ||
                     (cell->subdomain_id() == subdomain_id))
                    &&
-                   ((material_id == types::invalid_material_id)
+                   ((material_id == numbers::invalid_material_id)
                     ||
                     (cell->material_id() == material_id))) )
             {
@@ -927,7 +927,7 @@ namespace internal
                                  &&
                                  ((cell->neighbor(face_no)->material_id()
                                    == material_id) ||
-                                  (material_id == types::invalid_material_id));
+                                  (material_id == numbers::invalid_material_id));
               else
                 {
                   for (unsigned int sf=0; sf<face->n_children(); ++sf)
@@ -935,7 +935,7 @@ namespace internal
                           ->subdomain_id() == subdomain_id)
                          &&
                          (material_id ==
-                          types::invalid_material_id))
+                          numbers::invalid_material_id))
                         ||
                         ((cell->neighbor_child_on_subface(face_no,sf)
                           ->material_id() == material_id)
@@ -1207,7 +1207,7 @@ estimate (const Mapping<1,spacedim>                    &mapping,
   const unsigned int n_solution_vectors = solutions.size();
 
                                    // sanity checks
-  Assert (neumann_bc.find(types::internal_face_boundary_id) == neumann_bc.end(),
+  Assert (neumann_bc.find(numbers::internal_face_boundary_id) == neumann_bc.end(),
           ExcInvalidBoundaryIndicator());
 
   for (typename FunctionMap<spacedim>::type::const_iterator i=neumann_bc.begin();
@@ -1305,7 +1305,7 @@ estimate (const Mapping<1,spacedim>                    &mapping,
          ||
          (cell->subdomain_id() == subdomain_id))
         &&
-        ((material_id == types::invalid_material_id)
+        ((material_id == numbers::invalid_material_id)
          ||
          (cell->material_id() == material_id)))
       {
@@ -1668,7 +1668,7 @@ estimate (const Mapping<dim, spacedim>                  &mapping,
           ||
           (cell->subdomain_id() == subdomain_id))
          &&
-         ((material_id == types::invalid_material_id)
+         ((material_id == numbers::invalid_material_id)
           ||
           (cell->material_id() == material_id)))
       {

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.