]> https://gitweb.dealii.org/ - dealii.git/commitdiff
When defining abbreviations for non-scalar data, use references,
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Sat, 14 Aug 2010 03:38:34 +0000 (03:38 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Sat, 14 Aug 2010 03:38:34 +0000 (03:38 +0000)
i.e. use
      const std::vector<Tensor<2, dim> > &
jacobians = fe_values.get_jacobians ();
instead of
      std::vector<Tensor<2, dim> >
        jacobians = fe_values.get_jacobians ();

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

deal.II/deal.II/include/numerics/vectors.templates.h

index 976202d974a479b8336b828e859424772f8dec72..1ecd8b47f9ae5bf811d4f7474a1d3fb62518cc8f 100644 (file)
@@ -2947,16 +2947,18 @@ namespace internals
 
                                       // Initialize the required
                                       // objects.
-      std::vector<Tensor<2, dim> > jacobians = fe_values.get_jacobians ();
+      const std::vector<Tensor<2, dim> > &
+       jacobians = fe_values.get_jacobians ();
+      const std::vector<Point<dim> > &
+       quadrature_points = quadrature.get_points ();
+
       std::vector<Point<dim> > tangentials (fe_values.n_quadrature_points);
-      std::vector<Point<dim> > quadrature_points = fe_values.get_quadrature_points ();
       std::vector<Vector<double> > values (fe_values.n_quadrature_points,
                                           Vector<double> (dim));
 
                                       // Get boundary function values
                                       // at quadrature points.
       boundary_function.vector_value_list (quadrature_points, values);
-      quadrature_points = quadrature.get_points ();
 
       const unsigned int superdegree = cell->get_fe ().degree;
       const unsigned int degree = superdegree - 1;
@@ -3167,8 +3169,11 @@ namespace internals
       fe_values.reinit (cell);
 
                                       // Initialize the required objects.
-      std::vector<Tensor<2, dim> > jacobians = fe_values.get_jacobians ();
-      std::vector<Point<dim> > quadrature_points = fe_values.get_quadrature_points ();
+      const std::vector<Tensor<2, dim> > &
+       jacobians = fe_values.get_jacobians ();
+      const std::vector<Point<dim> > &
+       quadrature_points = fe_values.get_quadrature_points ();
+
       std::vector<Vector<double> > values (fe_values.n_quadrature_points,
                                           Vector<double> (dim));
 
@@ -3416,16 +3421,18 @@ namespace internals
       fe_values.reinit (cell);
 
                                       // Initialize the required objects.
+      const std::vector<Tensor<2, dim> > &
+       jacobians = fe_values.get_jacobians ();
+      const std::vector<Point<dim> > &
+       quadrature_points = quadrature.get_points ();
+
       std::vector<Point<dim> > tangentials (fe_values.n_quadrature_points);
-      std::vector<Tensor<2, dim> > jacobians = fe_values.get_jacobians ();
-      std::vector<Point<dim> > quadrature_points = fe_values.get_quadrature_points ();
       std::vector<Vector<double> > values (fe_values.n_quadrature_points,
                                           Vector<double> (dim));
 
                                       // Get boundary function values
                                       // at quadrature points.
       boundary_function.vector_value_list (quadrature_points, values);
-      quadrature_points = quadrature.get_points ();
 
       const unsigned int degree = cell->get_fe ().degree - 1;
       Point<dim> shifted_reference_point_1;

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.