Make the cell argument to transform_*_to_*_cell a reference, rather than a value.

git-svn-id: https://svn.dealii.org/trunk@7366 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/deal.II/include/fe/mapping.h b/deal.II/deal.II/include/fe/mapping.h
index 5f9806e5796b659859e4189301cb4fbd36df6c2d..db8e033abf8948d78b6290e3127690a049f0157e 100644 (file)
--- a/deal.II/deal.II/include/fe/mapping.h
+++ b/deal.II/deal.II/include/fe/mapping.h
@@ -67,9 +67,10 @@ class Mapping : public Subscriptor
                                      * @p{p_real} on the real cell
                                      * @p{cell} and returns @p{p_real}.
                                      */
-    virtual Point<dim> transform_unit_to_real_cell (
-      const typename Triangulation<dim>::cell_iterator cell,
-      const Point<dim> &p) const=0;
+    virtual Point<dim>
+    transform_unit_to_real_cell (
+      const typename Triangulation<dim>::cell_iterator &cell,
+      const Point<dim>                                 &p) const = 0;
     
                                     /**
                                      * Transforms the point @p{p} on
@@ -77,9 +78,10 @@ class Mapping : public Subscriptor
                                      * @p{p_unit} on the unit cell
                                      * @p{cell} and returns @p{p_unit}.
                                      */
-    virtual Point<dim> transform_real_to_unit_cell (
-      const typename Triangulation<dim>::cell_iterator cell,
-      const Point<dim> &p) const=0;
+    virtual Point<dim>
+    transform_real_to_unit_cell (
+      const typename Triangulation<dim>::cell_iterator &cell,
+      const Point<dim>                                 &p) const = 0;
     
                                     /**
                                      * Base class for internal data

diff --git a/deal.II/deal.II/include/fe/mapping_cartesian.h b/deal.II/deal.II/include/fe/mapping_cartesian.h
index f46808eb9d60f91612db7ec20b29e27d58733d8b..2161a1e25eaaa461b2faa3520d05cf0f42795f31 100644 (file)
--- a/deal.II/deal.II/include/fe/mapping_cartesian.h
+++ b/deal.II/deal.II/include/fe/mapping_cartesian.h
@@ -148,9 +148,10 @@ class MappingCartesian : public Mapping<dim>
                                      * @p{p_real} on the real cell
                                      * @p{cell} and returns @p{p_real}.
                                      */
-    virtual Point<dim> transform_unit_to_real_cell (
-      const typename Triangulation<dim>::cell_iterator cell,
-      const Point<dim> &p) const;
+    virtual Point<dim>
+    transform_unit_to_real_cell (
+      const typename Triangulation<dim>::cell_iterator &cell,
+      const Point<dim>                                 &p) const;
 
                                     /**
                                      * Transforms the point @p{p} on
@@ -162,9 +163,10 @@ class MappingCartesian : public Mapping<dim>
                                      * @p{transform_unit_to_real_cell}
                                      * function.
                                      */
-    virtual Point<dim> transform_real_to_unit_cell (
-      const typename Triangulation<dim>::cell_iterator cell,
-      const Point<dim> &p) const;
+    virtual Point<dim>
+    transform_real_to_unit_cell (
+      const typename Triangulation<dim>::cell_iterator &cell,
+      const Point<dim>                                 &p) const;
     
                                     /**
                                      * Implementation of the interface in

diff --git a/deal.II/deal.II/include/fe/mapping_q.h b/deal.II/deal.II/include/fe/mapping_q.h
index 2b1edbc99da3dbe2e07c29c8caf63235e211a671..9071c0cab2b17da844aa9dbf8cfa91fb886d5666 100644 (file)
--- a/deal.II/deal.II/include/fe/mapping_q.h
+++ b/deal.II/deal.II/include/fe/mapping_q.h
@@ -55,9 +55,10 @@ class MappingQ : public MappingQ1<dim>
                                      * @p{p_real} on the real cell
                                      * @p{cell} and returns @p{p_real}.
                                      */
-    virtual Point<dim> transform_unit_to_real_cell (
-      const typename Triangulation<dim>::cell_iterator cell,
-      const Point<dim> &p) const;
+    virtual Point<dim>
+    transform_unit_to_real_cell (
+      const typename Triangulation<dim>::cell_iterator &cell,
+      const Point<dim>                                 &p) const;
     
                                     /**
                                      * Transforms the point @p{p} on
@@ -69,9 +70,10 @@ class MappingQ : public MappingQ1<dim>
                                      * @p{transform_unit_to_real_cell}
                                      * function.
                                      */
-    virtual Point<dim> transform_real_to_unit_cell (
-      const typename Triangulation<dim>::cell_iterator cell,
-      const Point<dim> &p) const;
+    virtual Point<dim>
+    transform_real_to_unit_cell (
+      const typename Triangulation<dim>::cell_iterator &cell,
+      const Point<dim>                                 &p) const;
 
                                     /**
                                      * Implementation of the

diff --git a/deal.II/deal.II/include/fe/mapping_q1.h b/deal.II/deal.II/include/fe/mapping_q1.h
index f941af5cf1c88e0716d69385ac24a6a8fadb9a76..e53e5a4195b56e549d1cfd9efbacce853bee5f9c 100644 (file)
--- a/deal.II/deal.II/include/fe/mapping_q1.h
+++ b/deal.II/deal.II/include/fe/mapping_q1.h
@@ -52,9 +52,10 @@ class MappingQ1 : public Mapping<dim>
                                      * @p{p_real} on the real cell
                                      * @p{cell} and returns @p{p_real}.
                                      */
-    virtual Point<dim> transform_unit_to_real_cell (
-      const typename Triangulation<dim>::cell_iterator cell,
-      const Point<dim> &p) const;
+    virtual Point<dim>
+    transform_unit_to_real_cell (
+      const typename Triangulation<dim>::cell_iterator &cell,
+      const Point<dim>                                 &p) const;
     
                                     /**
                                      * Transforms the point @p{p} on
@@ -66,9 +67,10 @@ class MappingQ1 : public Mapping<dim>
                                      * @p{transform_unit_to_real_cell}
                                      * function.
                                      */
-    virtual Point<dim> transform_real_to_unit_cell (
-      const typename Triangulation<dim>::cell_iterator cell,
-      const Point<dim> &p) const;
+    virtual Point<dim>
+    transform_real_to_unit_cell (
+      const typename Triangulation<dim>::cell_iterator &cell,
+      const Point<dim>                                 &p) const;
     
                                     /**
                                      * Implementation of the interface in
@@ -527,7 +529,7 @@ class MappingQ1 : public Mapping<dim>
                                      * @p{mdata} will be changed by
                                      * this function.
                                      */
-    void transform_real_to_unit_cell_internal (const typename Triangulation<dim>::cell_iterator cell,
+    void transform_real_to_unit_cell_internal (const typename Triangulation<dim>::cell_iterator &cell,
                                               const Point<dim> &p,
                                               InternalData &mdata,
                                               Point<dim> &p_unit) const;

diff --git a/deal.II/deal.II/source/fe/mapping_cartesian.cc b/deal.II/deal.II/source/fe/mapping_cartesian.cc
index c948e405341ccfbf6fb1e45f697c9cf31afee66f..5214b53ae3cb04b90e9de0e89c3832de9e72717c 100644 (file)
--- a/deal.II/deal.II/source/fe/mapping_cartesian.cc
+++ b/deal.II/deal.II/source/fe/mapping_cartesian.cc
@@ -545,9 +545,10 @@ MappingCartesian<dim>::transform_contravariant (Tensor<2,dim>       *begin,
 
 
 template <int dim>
-Point<dim> MappingCartesian<dim>::transform_unit_to_real_cell (
-  const typename Triangulation<dim>::cell_iterator cell,
-  const Point<dim> &p) const
+Point<dim>
+MappingCartesian<dim>::transform_unit_to_real_cell (
+  const typename Triangulation<dim>::cell_iterator &cell,
+  const Point<dim>                                 &p) const
 {
   Tensor<1,dim> length;
   const Point<dim> start = cell->vertex(0);
@@ -579,9 +580,10 @@ Point<dim> MappingCartesian<dim>::transform_unit_to_real_cell (
 
 
 template <int dim>
-Point<dim> MappingCartesian<dim>::transform_real_to_unit_cell (
-  const typename Triangulation<dim>::cell_iterator cell,
-  const Point<dim> &p) const
+Point<dim>
+MappingCartesian<dim>::transform_real_to_unit_cell (
+  const typename Triangulation<dim>::cell_iterator &cell,
+  const Point<dim>                                 &p) const
 {
   const Point<dim> &start = cell->vertex(0);
   Point<dim> real = p;

diff --git a/deal.II/deal.II/source/fe/mapping_q.cc b/deal.II/deal.II/source/fe/mapping_q.cc
index 0524e035ddf0def48eee920d733bea90ed1b95c7..47543ddd12ee7d5b8cd3d3a60096c86a0d86edf7 100644 (file)
--- a/deal.II/deal.II/source/fe/mapping_q.cc
+++ b/deal.II/deal.II/source/fe/mapping_q.cc
@@ -1063,9 +1063,10 @@ MappingQ<dim>::transform_contravariant (Tensor<2,dim>       *begin,
 
 
 template <int dim>
-Point<dim> MappingQ<dim>::transform_unit_to_real_cell (
-  const typename Triangulation<dim>::cell_iterator cell,
-  const Point<dim> &p) const
+Point<dim>
+MappingQ<dim>::transform_unit_to_real_cell (
+  const typename Triangulation<dim>::cell_iterator &cell,
+  const Point<dim>                                 &p) const
 {
                                   // Use the get_data function to
                                   // create an InternalData with data
@@ -1096,9 +1097,10 @@ Point<dim> MappingQ<dim>::transform_unit_to_real_cell (
 
 
 template <int dim>
-Point<dim> MappingQ<dim>::transform_real_to_unit_cell (
-  const typename Triangulation<dim>::cell_iterator cell,
-  const Point<dim> &p) const
+Point<dim>
+MappingQ<dim>::transform_real_to_unit_cell (
+  const typename Triangulation<dim>::cell_iterator &cell,
+  const Point<dim>                                 &p) const
 {
                                   // first a Newton iteration based
                                   // on a Q1 mapping

diff --git a/deal.II/deal.II/source/fe/mapping_q1.cc b/deal.II/deal.II/source/fe/mapping_q1.cc
index b389455a2d44e4fcf88dafc69703df45113f6730..9c71143df16adcb149df950c612f439e0f124e82 100644 (file)
--- a/deal.II/deal.II/source/fe/mapping_q1.cc
+++ b/deal.II/deal.II/source/fe/mapping_q1.cc
@@ -884,9 +884,10 @@ MappingQ1<dim>::transform_contravariant (Tensor<2,dim>       *begin,
 
 
 template <int dim>
-Point<dim> MappingQ1<dim>::transform_unit_to_real_cell (
-  const typename Triangulation<dim>::cell_iterator cell,
-  const Point<dim> &p) const
+Point<dim>
+MappingQ1<dim>::transform_unit_to_real_cell (
+  const typename Triangulation<dim>::cell_iterator &cell,
+  const Point<dim>                                 &p) const
 {
                                   // Use the get_data function to
                                   // create an InternalData with data
@@ -910,7 +911,8 @@ Point<dim> MappingQ1<dim>::transform_unit_to_real_cell (
 
 
 template <int dim>
-Point<dim> MappingQ1<dim>::transform_unit_to_real_cell_internal (
+Point<dim>
+MappingQ1<dim>::transform_unit_to_real_cell_internal (
   const InternalData &data) const
 {
   const unsigned int n_mapping_points=data.mapping_support_points.size();
@@ -927,9 +929,10 @@ Point<dim> MappingQ1<dim>::transform_unit_to_real_cell_internal (
 
 
 template <int dim>
-Point<dim> MappingQ1<dim>::transform_real_to_unit_cell (
-  const typename Triangulation<dim>::cell_iterator cell,
-  const Point<dim> &p) const
+Point<dim>
+MappingQ1<dim>::transform_real_to_unit_cell (
+  const typename Triangulation<dim>::cell_iterator &cell,
+  const Point<dim>                                 &p) const
 {
                                   // Let the start value of the
                                   // newton iteration be the center
@@ -965,11 +968,12 @@ Point<dim> MappingQ1<dim>::transform_real_to_unit_cell (
 
 
 template <int dim>
-void MappingQ1<dim>::transform_real_to_unit_cell_internal (
-  const typename Triangulation<dim>::cell_iterator cell,
-  const Point<dim> &p,
-  InternalData &mdata,
-  Point<dim> &p_unit) const
+void
+MappingQ1<dim>::transform_real_to_unit_cell_internal (
+  const typename Triangulation<dim>::cell_iterator &cell,
+  const Point<dim>                                 &p,
+  InternalData                                     &mdata,
+  Point<dim>                                       &p_unit) const
 {
   const unsigned int n_shapes=mdata.shape_values.size();
   Assert(n_shapes!=0, ExcInternalError());

diff --git a/deal.II/deal.II/source/fe/mapping_q1_eulerian.cc b/deal.II/deal.II/source/fe/mapping_q1_eulerian.cc
index c001b4d4609e897607379be47ff96248c70acc82..5d1e6527ae15e0b561b4c1bd53f3a51422d8fef5 100644 (file)
--- a/deal.II/deal.II/source/fe/mapping_q1_eulerian.cc
+++ b/deal.II/deal.II/source/fe/mapping_q1_eulerian.cc
@@ -22,8 +22,8 @@
 
 
 template <int dim>
-MappingQ1Eulerian<dim>::MappingQ1Eulerian ( const Vector<double>  &euler_transform_vectors,
-                                           const DoFHandler<dim> &shiftmap_dof_handler)
+MappingQ1Eulerian<dim>::MappingQ1Eulerian (const Vector<double>  &euler_transform_vectors,
+                                           const DoFHandler<dim> &shiftmap_dof_handler)
                      :
                     euler_transform_vectors(euler_transform_vectors),
                     shiftmap_dof_handler(&shiftmap_dof_handler)