From: heister <heister@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Fri, 21 Dec 2012 01:40:23 +0000 (+0000)
Subject: revert more whitespace crap
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=99fd5c32b522bf0b205f99c2e14f33374b05b964;p=dealii-svn.git

revert more whitespace crap

git-svn-id: https://svn.dealii.org/branches/branch_merge_mg_into_dof_handler@27841 0785d39b-7218-0410-832d-ea1e28bc413d
---

diff --git a/deal.II/source/base/data_out_base.cc b/deal.II/source/base/data_out_base.cc
index ba57a738b6..5a5d9c9dd7 100644
--- a/deal.II/source/base/data_out_base.cc
+++ b/deal.II/source/base/data_out_base.cc
@@ -1819,12 +1819,12 @@ DataOutBase::GnuplotFlags::GnuplotFlags ()
 
 
 
-void DataOutBase::GnuplotFlags::declare_parameters (ParameterHandler &/*prm*/)
+void DataOutBase::GnuplotFlags::declare_parameters (ParameterHandler & /*prm*/)
 {}
 
 
 
-void DataOutBase::GnuplotFlags::parse_parameters (const ParameterHandler &/*prm*/) const
+void DataOutBase::GnuplotFlags::parse_parameters (const ParameterHandler & /*prm*/) const
 {}
 
 
@@ -2128,12 +2128,12 @@ DataOutBase::GmvFlags::GmvFlags ()
 
 
 
-void DataOutBase::GmvFlags::declare_parameters (ParameterHandler &/*prm*/)
+void DataOutBase::GmvFlags::declare_parameters (ParameterHandler & /*prm*/)
 {}
 
 
 
-void DataOutBase::GmvFlags::parse_parameters (const ParameterHandler &/*prm*/) const
+void DataOutBase::GmvFlags::parse_parameters (const ParameterHandler & /*prm*/) const
 {}
 
 
@@ -2157,12 +2157,12 @@ TecplotFlags (const char *tecplot_binary_file_name,
 
 
 
-void DataOutBase::TecplotFlags::declare_parameters (ParameterHandler &/*prm*/)
+void DataOutBase::TecplotFlags::declare_parameters (ParameterHandler & /*prm*/)
 {}
 
 
 
-void DataOutBase::TecplotFlags::parse_parameters (const ParameterHandler &/*prm*/) const
+void DataOutBase::TecplotFlags::parse_parameters (const ParameterHandler & /*prm*/) const
 {}
 
 
@@ -2181,12 +2181,12 @@ DataOutBase::VtkFlags::VtkFlags ()
 
 
 
-void DataOutBase::VtkFlags::declare_parameters (ParameterHandler &/*prm*/)
+void DataOutBase::VtkFlags::declare_parameters (ParameterHandler & /*prm*/)
 {}
 
 
 
-void DataOutBase::VtkFlags::parse_parameters (const ParameterHandler &/*prm*/) const
+void DataOutBase::VtkFlags::parse_parameters (const ParameterHandler & /*prm*/) const
 {}
 
 
@@ -2208,12 +2208,12 @@ DataOutBase::Deal_II_IntermediateFlags::Deal_II_IntermediateFlags ()
 
 
 
-void DataOutBase::Deal_II_IntermediateFlags::declare_parameters (ParameterHandler &/*prm*/)
+void DataOutBase::Deal_II_IntermediateFlags::declare_parameters (ParameterHandler & /*prm*/)
 {}
 
 
 
-void DataOutBase::Deal_II_IntermediateFlags::parse_parameters (const ParameterHandler &/*prm*/) const
+void DataOutBase::Deal_II_IntermediateFlags::parse_parameters (const ParameterHandler & /*prm*/) const
 {}
 
 
@@ -2873,7 +2873,7 @@ template <int dim, int spacedim>
 void DataOutBase::write_gnuplot (const std::vector<Patch<dim,spacedim> > &patches,
                                  const std::vector<std::string>          &data_names,
                                  const std::vector<std_cxx1x::tuple<unsigned int, unsigned int, std::string> > &,
-                                 const GnuplotFlags                      &/*flags*/,
+                                 const GnuplotFlags &                      /*flags*/,
                                  std::ostream                            &out)
 {
   AssertThrow (out, ExcIO());
@@ -3473,7 +3473,7 @@ void DataOutBase::write_povray (const std::vector<Patch<dim,spacedim> > &patches
 
 template <int dim, int spacedim>
 void DataOutBase::write_eps (const std::vector<Patch<dim,spacedim> > &patches,
-                             const std::vector<std::string>          &/*data_names*/,
+                             const std::vector<std::string> &          /*data_names*/,
                              const std::vector<std_cxx1x::tuple<unsigned int, unsigned int, std::string> > &,
                              const EpsFlags                          &flags,
                              std::ostream                            &out)
@@ -5253,7 +5253,7 @@ DataOutBase::
 write_deal_II_intermediate (const std::vector<Patch<dim,spacedim> > &patches,
                             const std::vector<std::string>          &data_names,
                             const std::vector<std_cxx1x::tuple<unsigned int, unsigned int, std::string> > &vector_data_ranges,
-                            const Deal_II_IntermediateFlags         &/*flags*/,
+                            const Deal_II_IntermediateFlags &         /*flags*/,
                             std::ostream                            &out)
 {
   AssertThrow (out, ExcIO());
@@ -5535,7 +5535,7 @@ template <int dim, int spacedim>
 void
 DataOutInterface<dim,spacedim>::
 write_pvd_record (std::ostream &out,
-                  const std::vector<std::pair<double,std::string> >  &times_and_names) const
+                  const std::vector<std::pair<double,std::string> > &times_and_names) const
 {
   AssertThrow (out, ExcIO());
 
diff --git a/deal.II/source/base/function_lib.cc b/deal.II/source/base/function_lib.cc
index 3d8b0e0736..6f4f5cc5f8 100644
--- a/deal.II/source/base/function_lib.cc
+++ b/deal.II/source/base/function_lib.cc
@@ -1363,7 +1363,7 @@ namespace Functions
 
   Tensor<1,2>
   LSingularityGradFunction::gradient (
-    const Point<2>   &/*p*/,
+    const Point<2> &   /*p*/,
     const unsigned int /*component*/) const
   {
     Assert(false, ExcNotImplemented());
@@ -1802,7 +1802,7 @@ namespace Functions
   template<int dim>
   void
   JumpFunction<dim>::gradient_list (const std::vector<Point<dim> > &p,
-                                    std::vector<Tensor<1,dim> >  &gradients,
+                                    std::vector<Tensor<1,dim> > &gradients,
                                     const unsigned int) const
   {
     Assert (gradients.size() == p.size(),
diff --git a/deal.II/source/base/function_parser.cc b/deal.II/source/base/function_parser.cc
index 7dc72e6e07..5b4b800a46 100644
--- a/deal.II/source/base/function_parser.cc
+++ b/deal.II/source/base/function_parser.cc
@@ -197,7 +197,7 @@ void FunctionParser<dim>::initialize (const std::string &variables,
 
 
 template <int dim>
-double FunctionParser<dim>::value (const Point<dim>  &p,
+double FunctionParser<dim>::value (const Point<dim> &p,
                                    const unsigned int component) const
 {
   Assert (initialized==true, ExcNotInitialized());
diff --git a/deal.II/source/base/parameter_handler.cc b/deal.II/source/base/parameter_handler.cc
index 8227166e21..1975331b27 100644
--- a/deal.II/source/base/parameter_handler.cc
+++ b/deal.II/source/base/parameter_handler.cc
@@ -436,7 +436,7 @@ namespace Patterns
   const char *List::description_init = "[List";
 
 
-  List::List (const PatternBase  &p,
+  List::List (const PatternBase &p,
               const unsigned int  min_elements,
               const unsigned int  max_elements)
     :
@@ -571,8 +571,8 @@ namespace Patterns
   const char *Map::description_init = "[Map";
 
 
-  Map::Map (const PatternBase  &p_key,
-            const PatternBase  &p_value,
+  Map::Map (const PatternBase &p_key,
+            const PatternBase &p_value,
             const unsigned int  min_elements,
             const unsigned int  max_elements)
     :
diff --git a/deal.II/source/base/parsed_function.cc b/deal.II/source/base/parsed_function.cc
index 42540badf2..da058af424 100644
--- a/deal.II/source/base/parsed_function.cc
+++ b/deal.II/source/base/parsed_function.cc
@@ -30,7 +30,7 @@ namespace Functions
 
   template <int dim>
   void
-  ParsedFunction<dim>::declare_parameters(ParameterHandler  &prm,
+  ParsedFunction<dim>::declare_parameters(ParameterHandler &prm,
                                           const unsigned int n_components)
   {
     Assert(n_components > 0, ExcZero());
diff --git a/deal.II/source/base/polynomial_space.cc b/deal.II/source/base/polynomial_space.cc
index 2852f213af..f97d13f01f 100644
--- a/deal.II/source/base/polynomial_space.cc
+++ b/deal.II/source/base/polynomial_space.cc
@@ -125,7 +125,7 @@ PolynomialSpace<dim>::set_numbering(
 template <int dim>
 double
 PolynomialSpace<dim>::compute_value (const unsigned int i,
-                                     const Point<dim>  &p) const
+                                     const Point<dim> &p) const
 {
   unsigned int ix[dim];
   compute_index(i,ix);
@@ -143,7 +143,7 @@ PolynomialSpace<dim>::compute_value (const unsigned int i,
 template <int dim>
 Tensor<1,dim>
 PolynomialSpace<dim>::compute_grad (const unsigned int i,
-                                    const Point<dim>  &p) const
+                                    const Point<dim> &p) const
 {
   unsigned int ix[dim];
   compute_index(i,ix);
@@ -169,7 +169,7 @@ PolynomialSpace<dim>::compute_grad (const unsigned int i,
 template <int dim>
 Tensor<2,dim>
 PolynomialSpace<dim>::compute_grad_grad (const unsigned int i,
-                                         const Point<dim>  &p) const
+                                         const Point<dim> &p) const
 {
   unsigned int ix[dim];
   compute_index(i,ix);
diff --git a/deal.II/source/fe/fe_dgp_nonparametric.cc b/deal.II/source/fe/fe_dgp_nonparametric.cc
index fccbec69f8..7f891fd4dc 100644
--- a/deal.II/source/fe/fe_dgp_nonparametric.cc
+++ b/deal.II/source/fe/fe_dgp_nonparametric.cc
@@ -292,7 +292,7 @@ FE_DGPNonparametric<dim,spacedim>::fill_fe_values (
   typename Mapping<dim,spacedim>::InternalDataBase &,
   typename Mapping<dim,spacedim>::InternalDataBase &fedata,
   FEValuesData<dim,spacedim> &data,
-  CellSimilarity::Similarity &/*cell_similarity*/) const
+  CellSimilarity::Similarity & /*cell_similarity*/) const
 {
   // convert data object to internal
   // data for this class. fails with
diff --git a/deal.II/source/fe/fe_nothing.cc b/deal.II/source/fe/fe_nothing.cc
index 71b32f25d7..527bc4f4cf 100644
--- a/deal.II/source/fe/fe_nothing.cc
+++ b/deal.II/source/fe/fe_nothing.cc
@@ -168,7 +168,7 @@ compare_for_face_domination (const FiniteElement<dim> &) const
 template <int dim>
 std::vector<std::pair<unsigned int, unsigned int> >
 FE_Nothing<dim> ::
-hp_vertex_dof_identities (const FiniteElement<dim> &/*fe_other*/) const
+hp_vertex_dof_identities (const FiniteElement<dim> & /*fe_other*/) const
 {
   // the FE_Nothing has no
   // degrees of freedom, so there
@@ -181,7 +181,7 @@ hp_vertex_dof_identities (const FiniteElement<dim> &/*fe_other*/) const
 template <int dim>
 std::vector<std::pair<unsigned int, unsigned int> >
 FE_Nothing<dim> ::
-hp_line_dof_identities (const FiniteElement<dim> &/*fe_other*/) const
+hp_line_dof_identities (const FiniteElement<dim> & /*fe_other*/) const
 {
   // the FE_Nothing has no
   // degrees of freedom, so there
@@ -194,7 +194,7 @@ hp_line_dof_identities (const FiniteElement<dim> &/*fe_other*/) const
 template <int dim>
 std::vector<std::pair<unsigned int, unsigned int> >
 FE_Nothing<dim> ::
-hp_quad_dof_identities (const FiniteElement<dim> &/*fe_other*/) const
+hp_quad_dof_identities (const FiniteElement<dim> & /*fe_other*/) const
 {
   // the FE_Nothing has no
   // degrees of freedom, so there
@@ -216,7 +216,7 @@ hp_constraints_are_implemented () const
 template <int dim>
 void
 FE_Nothing<dim>::
-get_face_interpolation_matrix (const FiniteElement<dim> &/*source_fe*/,
+get_face_interpolation_matrix (const FiniteElement<dim> & /*source_fe*/,
                                FullMatrix<double>       &interpolation_matrix) const
 {
   // since this element has no face dofs, the
@@ -236,7 +236,7 @@ void
 FE_Nothing<dim>::
 get_subface_interpolation_matrix (const FiniteElement<dim> & /*source_fe*/,
                                   const unsigned int /*index*/,
-                                  FullMatrix<double>  &interpolation_matrix) const
+                                  FullMatrix<double> &interpolation_matrix) const
 {
   // since this element has no face dofs, the
   // interpolation matrix is necessarily empty
diff --git a/deal.II/source/fe/fe_poly_tensor.cc b/deal.II/source/fe/fe_poly_tensor.cc
index 8d5562849c..00c281f83a 100644
--- a/deal.II/source/fe/fe_poly_tensor.cc
+++ b/deal.II/source/fe/fe_poly_tensor.cc
@@ -83,7 +83,7 @@ namespace
 
 
   void
-  get_face_sign_change (const Triangulation<3>::cell_iterator &/*cell*/,
+  get_face_sign_change (const Triangulation<3>::cell_iterator & /*cell*/,
                         const unsigned int                     /*dofs_per_face*/,
                         std::vector<double>                   &face_sign)
   {
diff --git a/deal.II/source/fe/fe_q.cc b/deal.II/source/fe/fe_q.cc
index ede4249d59..93b8de41dc 100644
--- a/deal.II/source/fe/fe_q.cc
+++ b/deal.II/source/fe/fe_q.cc
@@ -116,7 +116,7 @@ namespace FE_Q_Helper
     template <int dim>
     inline
     std::vector<Polynomials::Polynomial<double> >
-    generate_poly_space1d (const std::vector<Point<dim> >  &unit_support_points,
+    generate_poly_space1d (const std::vector<Point<dim> > &unit_support_points,
                            const std::vector<unsigned int> &index_map_inverse,
                            const unsigned int               dofs1d)
     {
@@ -767,8 +767,8 @@ get_interpolation_matrix (const FiniteElement<dim,spacedim> &x_source_fe,
 template <>
 void
 FE_Q<1>::
-get_face_interpolation_matrix (const FiniteElement<1,1> &/*x_source_fe*/,
-                               FullMatrix<double>     &/*interpolation_matrix*/) const
+get_face_interpolation_matrix (const FiniteElement<1,1> & /*x_source_fe*/,
+                               FullMatrix<double> &     /*interpolation_matrix*/) const
 {
   Assert (false, ExcImpossibleInDim(1));
 }
@@ -778,9 +778,9 @@ get_face_interpolation_matrix (const FiniteElement<1,1> &/*x_source_fe*/,
 template <>
 void
 FE_Q<1>::
-get_subface_interpolation_matrix (const FiniteElement<1,1> &/*x_source_fe*/,
+get_subface_interpolation_matrix (const FiniteElement<1,1> & /*x_source_fe*/,
                                   const unsigned int      /*subface*/,
-                                  FullMatrix<double>     &/*interpolation_matrix*/) const
+                                  FullMatrix<double> &     /*interpolation_matrix*/) const
 {
   Assert (false, ExcImpossibleInDim(1));
 }
@@ -789,8 +789,8 @@ get_subface_interpolation_matrix (const FiniteElement<1,1> &/*x_source_fe*/,
 template <>
 void
 FE_Q<1,2>::
-get_face_interpolation_matrix (const FiniteElement<1,2> &/*x_source_fe*/,
-                               FullMatrix<double>     &/*interpolation_matrix*/) const
+get_face_interpolation_matrix (const FiniteElement<1,2> & /*x_source_fe*/,
+                               FullMatrix<double> &     /*interpolation_matrix*/) const
 {
   typedef FiniteElement<1,2> FEL;
   Assert (false,
@@ -802,9 +802,9 @@ get_face_interpolation_matrix (const FiniteElement<1,2> &/*x_source_fe*/,
 template <>
 void
 FE_Q<1,2>::
-get_subface_interpolation_matrix (const FiniteElement<1,2> &/*x_source_fe*/,
+get_subface_interpolation_matrix (const FiniteElement<1,2> & /*x_source_fe*/,
                                   const unsigned int      /*subface*/,
-                                  FullMatrix<double>     &/*interpolation_matrix*/) const
+                                  FullMatrix<double> &     /*interpolation_matrix*/) const
 {
   typedef FiniteElement<1,2> FEL;
   Assert (false,
@@ -816,8 +816,8 @@ get_subface_interpolation_matrix (const FiniteElement<1,2> &/*x_source_fe*/,
 template <>
 void
 FE_Q<1,3>::
-get_face_interpolation_matrix (const FiniteElement<1,3> &/*x_source_fe*/,
-                               FullMatrix<double>     &/*interpolation_matrix*/) const
+get_face_interpolation_matrix (const FiniteElement<1,3> & /*x_source_fe*/,
+                               FullMatrix<double> &     /*interpolation_matrix*/) const
 {
   typedef FiniteElement<1,3> FEL;
   Assert (false,
@@ -829,9 +829,9 @@ get_face_interpolation_matrix (const FiniteElement<1,3> &/*x_source_fe*/,
 template <>
 void
 FE_Q<1,3>::
-get_subface_interpolation_matrix (const FiniteElement<1,3> &/*x_source_fe*/,
+get_subface_interpolation_matrix (const FiniteElement<1,3> & /*x_source_fe*/,
                                   const unsigned int      /*subface*/,
-                                  FullMatrix<double>     &/*interpolation_matrix*/) const
+                                  FullMatrix<double> &     /*interpolation_matrix*/) const
 {
   typedef FiniteElement<1,3> FEL;
   Assert (false,
@@ -1335,7 +1335,7 @@ void FE_Q<1>::initialize_unit_face_support_points ()
 }
 
 template <>
-void FE_Q<1>::initialize_unit_face_support_points (const Quadrature<1> &/*points*/)
+void FE_Q<1>::initialize_unit_face_support_points (const Quadrature<1> & /*points*/)
 {
   // no faces in 1d, so nothing to do
 }
@@ -1347,7 +1347,7 @@ void FE_Q<1,2>::initialize_unit_face_support_points ()
 }
 
 template <>
-void FE_Q<1,2>::initialize_unit_face_support_points (const Quadrature<1> &/*points*/)
+void FE_Q<1,2>::initialize_unit_face_support_points (const Quadrature<1> & /*points*/)
 {
   // no faces in 1d, so nothing to do
 }
@@ -1359,7 +1359,7 @@ void FE_Q<1,3>::initialize_unit_face_support_points ()
 }
 
 template <>
-void FE_Q<1,3>::initialize_unit_face_support_points (const Quadrature<1> &/*points*/)
+void FE_Q<1,3>::initialize_unit_face_support_points (const Quadrature<1> & /*points*/)
 {
   // no faces in 1d, so nothing to do
 }
diff --git a/deal.II/source/fe/fe_q_hierarchical.cc b/deal.II/source/fe/fe_q_hierarchical.cc
index ceef17fdac..0f712f1aec 100644
--- a/deal.II/source/fe/fe_q_hierarchical.cc
+++ b/deal.II/source/fe/fe_q_hierarchical.cc
@@ -587,8 +587,8 @@ void FE_Q_Hierarchical<1>::initialize_unit_face_support_points ()
 
 template <>
 void FE_Q_Hierarchical<1>::
-get_face_interpolation_matrix (const FiniteElement<1,1> &/*x_source_fe*/,
-                               FullMatrix<double>     &/*interpolation_matrix*/) const
+get_face_interpolation_matrix (const FiniteElement<1,1> & /*x_source_fe*/,
+                               FullMatrix<double> &     /*interpolation_matrix*/) const
 {
   Assert (false, ExcImpossibleInDim(1));
 }
@@ -597,9 +597,9 @@ get_face_interpolation_matrix (const FiniteElement<1,1> &/*x_source_fe*/,
 template <>
 void
 FE_Q_Hierarchical<1>::
-get_subface_interpolation_matrix (const FiniteElement<1,1> &/*x_source_fe*/,
+get_subface_interpolation_matrix (const FiniteElement<1,1> & /*x_source_fe*/,
                                   const unsigned int      /*subface*/,
-                                  FullMatrix<double>     &/*interpolation_matrix*/) const
+                                  FullMatrix<double> &     /*interpolation_matrix*/) const
 {
   Assert (false, ExcImpossibleInDim(1));
 }
diff --git a/deal.II/source/fe/fe_raviart_thomas_nodal.cc b/deal.II/source/fe/fe_raviart_thomas_nodal.cc
index f8dfcd8bbb..23597aadf3 100644
--- a/deal.II/source/fe/fe_raviart_thomas_nodal.cc
+++ b/deal.II/source/fe/fe_raviart_thomas_nodal.cc
@@ -540,8 +540,8 @@ FE_RaviartThomasNodal<dim>::compare_for_face_domination (
 template <>
 void
 FE_RaviartThomasNodal<1>::get_face_interpolation_matrix (
-  const FiniteElement<1,1> &/*x_source_fe*/,
-  FullMatrix<double>     &/*interpolation_matrix*/) const
+  const FiniteElement<1,1> & /*x_source_fe*/,
+  FullMatrix<double> &     /*interpolation_matrix*/) const
 {
   Assert (false, ExcImpossibleInDim(1));
 }
@@ -550,9 +550,9 @@ FE_RaviartThomasNodal<1>::get_face_interpolation_matrix (
 template <>
 void
 FE_RaviartThomasNodal<1>::get_subface_interpolation_matrix (
-  const FiniteElement<1,1> &/*x_source_fe*/,
+  const FiniteElement<1,1> & /*x_source_fe*/,
   const unsigned int      /*subface*/,
-  FullMatrix<double>     &/*interpolation_matrix*/) const
+  FullMatrix<double> &     /*interpolation_matrix*/) const
 {
   Assert (false, ExcImpossibleInDim(1));
 }
diff --git a/deal.II/source/fe/fe_system.cc b/deal.II/source/fe/fe_system.cc
index 88048074d2..08678aaa3c 100644
--- a/deal.II/source/fe/fe_system.cc
+++ b/deal.II/source/fe/fe_system.cc
@@ -310,7 +310,7 @@ FESystem<dim,spacedim>::FESystem (const FiniteElement<dim,spacedim> &fe1,
 
 template <int dim, int spacedim>
 FESystem<dim,spacedim>::FESystem (
-  const std::vector<const FiniteElement<dim,spacedim>*>  &fes,
+  const std::vector<const FiniteElement<dim,spacedim>*> &fes,
   const std::vector<unsigned int>                  &multiplicities)
   :
   FiniteElement<dim,spacedim> (multiply_dof_numbers(fes, multiplicities),
@@ -392,7 +392,7 @@ FESystem<dim,spacedim>::shape_value (const unsigned int i,
 template <int dim, int spacedim>
 double
 FESystem<dim,spacedim>::shape_value_component (const unsigned int i,
-                                               const Point<dim>  &p,
+                                               const Point<dim> &p,
                                                const unsigned int component) const
 {
   Assert (i<this->dofs_per_cell, ExcIndexRange(i, 0, this->dofs_per_cell));
@@ -444,7 +444,7 @@ FESystem<dim,spacedim>::shape_grad (const unsigned int i,
 template <int dim, int spacedim>
 Tensor<1,dim>
 FESystem<dim,spacedim>::shape_grad_component (const unsigned int i,
-                                              const Point<dim>  &p,
+                                              const Point<dim> &p,
                                               const unsigned int component) const
 {
   Assert (i<this->dofs_per_cell, ExcIndexRange(i, 0, this->dofs_per_cell));
@@ -496,7 +496,7 @@ FESystem<dim,spacedim>::shape_grad_grad (const unsigned int i,
 template <int dim, int spacedim>
 Tensor<2,dim>
 FESystem<dim,spacedim>::shape_grad_grad_component (const unsigned int i,
-                                                   const Point<dim>  &p,
+                                                   const Point<dim> &p,
                                                    const unsigned int component) const
 {
   Assert (i<this->dofs_per_cell, ExcIndexRange(i, 0, this->dofs_per_cell));
@@ -1855,9 +1855,9 @@ void FESystem<dim,spacedim>::initialize (const std::vector<const FiniteElement<d
       if (multiplicities[i]>0)
         {
           base_elements[ind] =
-            std::make_pair (std_cxx1x::shared_ptr<const FiniteElement<dim,spacedim> >
-                            (fes[i]->clone()),
-                            multiplicities[i]);
+	    std::make_pair (std_cxx1x::shared_ptr<const FiniteElement<dim,spacedim> >
+			    (fes[i]->clone()),
+			    multiplicities[i]);
           ++ind;
         }
     }
diff --git a/deal.II/source/fe/mapping_cartesian.cc b/deal.II/source/fe/mapping_cartesian.cc
index 48ec6cfaea..5529fc527b 100644
--- a/deal.II/source/fe/mapping_cartesian.cc
+++ b/deal.II/source/fe/mapping_cartesian.cc
@@ -467,7 +467,7 @@ MappingCartesian<dim, spacedim>::fill_fe_face_values (
   typename Mapping<dim, spacedim>::InternalDataBase &mapping_data,
   std::vector<Point<dim> >     &quadrature_points,
   std::vector<double>          &JxW_values,
-  std::vector<Tensor<1,dim> >  &boundary_forms,
+  std::vector<Tensor<1,dim> > &boundary_forms,
   std::vector<Point<spacedim> >     &normal_vectors) const
 {
   // convert data object to internal
@@ -522,7 +522,7 @@ MappingCartesian<dim, spacedim>::fill_fe_subface_values (
   typename Mapping<dim, spacedim>::InternalDataBase &mapping_data,
   std::vector<Point<dim> >     &quadrature_points,
   std::vector<double>          &JxW_values,
-  std::vector<Tensor<1,dim> >  &boundary_forms,
+  std::vector<Tensor<1,dim> > &boundary_forms,
   std::vector<Point<spacedim> >     &normal_vectors) const
 {
   // convert data object to internal
diff --git a/deal.II/source/fe/mapping_q.cc b/deal.II/source/fe/mapping_q.cc
index 213905d2a7..665748712b 100644
--- a/deal.II/source/fe/mapping_q.cc
+++ b/deal.II/source/fe/mapping_q.cc
@@ -367,7 +367,7 @@ MappingQ<dim,spacedim>::fill_fe_face_values (
   typename Mapping<dim,spacedim>::InternalDataBase &mapping_data,
   std::vector<Point<spacedim> >     &quadrature_points,
   std::vector<double>          &JxW_values,
-  std::vector<Tensor<1,spacedim> >  &exterior_forms,
+  std::vector<Tensor<1,spacedim> > &exterior_forms,
   std::vector<Point<spacedim> >     &normal_vectors) const
 {
   // convert data object to internal
@@ -429,7 +429,7 @@ MappingQ<dim,spacedim>::fill_fe_subface_values (const typename Triangulation<dim
                                                 typename Mapping<dim,spacedim>::InternalDataBase &mapping_data,
                                                 std::vector<Point<spacedim> >     &quadrature_points,
                                                 std::vector<double>          &JxW_values,
-                                                std::vector<Tensor<1,spacedim> >  &exterior_forms,
+                                                std::vector<Tensor<1,spacedim> > &exterior_forms,
                                                 std::vector<Point<spacedim> >     &normal_vectors) const
 {
   // convert data object to internal
diff --git a/deal.II/source/fe/mapping_q1.cc b/deal.II/source/fe/mapping_q1.cc
index 953aea7fb3..211ad2c8b9 100644
--- a/deal.II/source/fe/mapping_q1.cc
+++ b/deal.II/source/fe/mapping_q1.cc
@@ -602,7 +602,7 @@ MappingQ1<dim,spacedim>::compute_fill (const typename Triangulation<dim,spacedim
                                        const unsigned int  n_q_points,
                                        const DataSetDescriptor  data_set,
                                        const CellSimilarity::Similarity cell_similarity,
-                                       InternalData  &data,
+                                       InternalData &data,
                                        std::vector<Point<spacedim> > &quadrature_points) const
 {
   const UpdateFlags update_flags(data.current_update_flags());
@@ -1139,7 +1139,7 @@ fill_fe_subface_values (const typename Triangulation<dim,spacedim>::cell_iterato
                         typename Mapping<dim,spacedim>::InternalDataBase &mapping_data,
                         std::vector<Point<spacedim> >     &quadrature_points,
                         std::vector<double>          &JxW_values,
-                        std::vector<Tensor<1,spacedim> >  &boundary_forms,
+                        std::vector<Tensor<1,spacedim> > &boundary_forms,
                         std::vector<Point<spacedim> >     &normal_vectors) const
 {
   // ensure that the following cast
@@ -1173,7 +1173,7 @@ void
 MappingQ1<dim,spacedim>::transform (
   const VectorSlice<const std::vector<Tensor<1, dim> > >  input,
   VectorSlice<std::vector<Tensor<1, spacedim> > >         output,
-  const typename Mapping<dim,spacedim>::InternalDataBase  &mapping_data,
+  const typename Mapping<dim,spacedim>::InternalDataBase &mapping_data,
   const MappingType                                       mapping_type) const
 {
   transform_fields(input, output, mapping_data, mapping_type);
@@ -1858,10 +1858,10 @@ transform_real_to_unit_cell_internal (const Triangulation<1,2>::cell_iterator &c
 template<>
 Point<1>
 MappingQ1<1,3>::
-transform_real_to_unit_cell_internal (const Triangulation<1,3>::cell_iterator &/*cell*/,
-                                      const Point<3> &/*p*/,
-                                      const Point<1> &/*initial_p_unit*/,
-                                      InternalData   &/*mdata*/) const
+transform_real_to_unit_cell_internal (const Triangulation<1,3>::cell_iterator & /*cell*/,
+                                      const Point<3> & /*p*/,
+                                      const Point<1> & /*initial_p_unit*/,
+                                      InternalData &   /*mdata*/) const
 {
   Assert(false, ExcNotImplemented());
   return Point<1>();
diff --git a/deal.II/source/grid/grid_generator.cc b/deal.II/source/grid/grid_generator.cc
index 8e0966cb7c..a2b6d02c63 100644
--- a/deal.II/source/grid/grid_generator.cc
+++ b/deal.II/source/grid/grid_generator.cc
@@ -180,7 +180,7 @@ void GridGenerator::hyper_cube (Triangulation<dim,spacedim> &tria,
 
 void
 GridGenerator::moebius (
-  Triangulation<3>  &tria,
+  Triangulation<3> &tria,
   const unsigned int   n_cells,
   const unsigned int   n_rotations,
   const double         R,
@@ -233,7 +233,7 @@ GridGenerator::moebius (
 
 
 void
-GridGenerator::torus (Triangulation<2,3>  &tria,
+GridGenerator::torus (Triangulation<2,3> &tria,
                       const double         R,
                       const double         r)
 {
@@ -372,7 +372,7 @@ GridGenerator::torus (Triangulation<2,3>  &tria,
 template<>
 void
 GridGenerator::parallelogram (
-  Triangulation<2>  &tria,
+  Triangulation<2> &tria,
   const Tensor<2,2> &corners,
   const bool      colorize)
 {
diff --git a/deal.II/source/grid/grid_in.cc b/deal.II/source/grid/grid_in.cc
index e61a7551ae..022757a363 100644
--- a/deal.II/source/grid/grid_in.cc
+++ b/deal.II/source/grid/grid_in.cc
@@ -2149,9 +2149,9 @@ void GridIn<dim, spacedim>::skip_comment_lines (std::istream &in,
 
 
 template <int dim, int spacedim>
-void GridIn<dim, spacedim>::debug_output_grid (const std::vector<CellData<dim> > &/*cells*/,
-                                               const std::vector<Point<spacedim> >    &/*vertices*/,
-                                               std::ostream                      &/*out*/)
+void GridIn<dim, spacedim>::debug_output_grid (const std::vector<CellData<dim> > & /*cells*/,
+                                               const std::vector<Point<spacedim> > &    /*vertices*/,
+                                               std::ostream &                      /*out*/)
 {
   Assert (false, ExcNotImplemented());
 }
diff --git a/deal.II/source/grid/grid_reordering.cc b/deal.II/source/grid/grid_reordering.cc
index dd8aeee74d..f6b79371d0 100644
--- a/deal.II/source/grid/grid_reordering.cc
+++ b/deal.II/source/grid/grid_reordering.cc
@@ -225,7 +225,7 @@ namespace internal
                   const unsigned int s1,
                   const unsigned int s2,
                   const unsigned int s3,
-                  const CellData<2>  &cd)
+                  const CellData<2> &cd)
       :
       original_cell_data (cd)
     {
diff --git a/deal.II/source/grid/tria.cc b/deal.II/source/grid/tria.cc
index 5f9f470826..6b0834b82b 100644
--- a/deal.II/source/grid/tria.cc
+++ b/deal.II/source/grid/tria.cc
@@ -1440,7 +1440,7 @@ namespace internal
       void
       create_triangulation (const std::vector<Point<spacedim> > &v,
                             const std::vector<CellData<1> >     &cells,
-                            const SubCellData                   &/*subcelldata*/,
+                            const SubCellData &                   /*subcelldata*/,
                             Triangulation<1,spacedim>           &triangulation)
       {
         AssertThrow (v.size() > 0, ExcMessage ("No vertices given"));
diff --git a/deal.II/source/grid/tria_boundary.cc b/deal.II/source/grid/tria_boundary.cc
index 247b780afc..f79e2f4e03 100644
--- a/deal.II/source/grid/tria_boundary.cc
+++ b/deal.II/source/grid/tria_boundary.cc
@@ -920,7 +920,7 @@ template <>
 Point<3>
 StraightBoundary<1,3>::
 project_to_surface (const Triangulation<1, 3>::quad_iterator &,
-                    const Point<3>  &y) const
+                    const Point<3> &y) const
 {
   return y;
 }
@@ -930,7 +930,7 @@ template <int dim, int spacedim>
 Point<spacedim>
 StraightBoundary<dim, spacedim>::
 project_to_surface (const typename Triangulation<dim, spacedim>::quad_iterator &quad,
-                    const Point<spacedim>  &y) const
+                    const Point<spacedim> &y) const
 {
   if (spacedim <= 2)
     return y;
diff --git a/deal.II/source/lac/constraint_matrix.cc b/deal.II/source/lac/constraint_matrix.cc
index 8f69b8e04c..550a0a393c 100644
--- a/deal.II/source/lac/constraint_matrix.cc
+++ b/deal.II/source/lac/constraint_matrix.cc
@@ -2131,7 +2131,7 @@ ConstraintMatrix::distribute (PETScWrappers::MPI::Vector &vec) const
 
 template<>
 void
-ConstraintMatrix::distribute (PETScWrappers::MPI::BlockVector &/*vec*/) const
+ConstraintMatrix::distribute (PETScWrappers::MPI::BlockVector & /*vec*/) const
 {
   Assert (sorted==true, ExcMatrixIsClosed());
   AssertThrow (false, ExcNotImplemented());
diff --git a/deal.II/source/lac/petsc_matrix_free.cc b/deal.II/source/lac/petsc_matrix_free.cc
index 16f51fe291..f688cf5517 100755
--- a/deal.II/source/lac/petsc_matrix_free.cc
+++ b/deal.II/source/lac/petsc_matrix_free.cc
@@ -179,7 +179,7 @@ namespace PETScWrappers
 
 
 
-  void MatrixFree::vmult (Vec  &dst, const Vec  &src) const
+  void MatrixFree::vmult (Vec &dst, const Vec &src) const
   {
 
 //TODO: Translate the given PETSc Vec* vector into a deal.II
diff --git a/deal.II/source/lac/petsc_parallel_vector.cc b/deal.II/source/lac/petsc_parallel_vector.cc
index fb25d77f61..225c1ebda1 100644
--- a/deal.II/source/lac/petsc_parallel_vector.cc
+++ b/deal.II/source/lac/petsc_parallel_vector.cc
@@ -53,7 +53,7 @@ namespace PETScWrappers
 
 
     Vector::Vector (const MPI_Comm    &communicator,
-                    const VectorBase  &v,
+                    const VectorBase &v,
                     const unsigned int local_size)
       :
       communicator (communicator)
diff --git a/deal.II/source/lac/petsc_solver.cc b/deal.II/source/lac/petsc_solver.cc
index ad9f1dfd50..caaa9ea4b9 100644
--- a/deal.II/source/lac/petsc_solver.cc
+++ b/deal.II/source/lac/petsc_solver.cc
@@ -42,7 +42,7 @@ namespace PETScWrappers
 
 
 
-  SolverBase::SolverBase (SolverControl  &cn,
+  SolverBase::SolverBase (SolverControl &cn,
                           const MPI_Comm &mpi_communicator)
     :
     solver_control (cn),
diff --git a/deal.II/source/lac/petsc_vector_base.cc b/deal.II/source/lac/petsc_vector_base.cc
index 135319b0d9..4b80d74837 100644
--- a/deal.II/source/lac/petsc_vector_base.cc
+++ b/deal.II/source/lac/petsc_vector_base.cc
@@ -321,7 +321,7 @@ namespace PETScWrappers
 
   void
   VectorBase::set (const std::vector<unsigned int> &indices,
-                   const std::vector<PetscScalar>  &values)
+                   const std::vector<PetscScalar> &values)
   {
     Assert (indices.size() == values.size(),
             ExcMessage ("Function called with arguments of different sizes"));
@@ -332,7 +332,7 @@ namespace PETScWrappers
 
   void
   VectorBase::add (const std::vector<unsigned int> &indices,
-                   const std::vector<PetscScalar>  &values)
+                   const std::vector<PetscScalar> &values)
   {
     Assert (indices.size() == values.size(),
             ExcMessage ("Function called with arguments of different sizes"));
diff --git a/deal.II/source/lac/slepc_solver.cc b/deal.II/source/lac/slepc_solver.cc
index 27c0895bcd..a0a5c6e865 100644
--- a/deal.II/source/lac/slepc_solver.cc
+++ b/deal.II/source/lac/slepc_solver.cc
@@ -44,7 +44,7 @@ namespace SLEPcWrappers
     AssertThrow (ierr == 0, ExcSLEPcError(ierr));
   }
 
-  SolverBase::SolverBase (SolverControl  &cn,
+  SolverBase::SolverBase (SolverControl &cn,
                           const MPI_Comm &mpi_communicator)
     :
     solver_control (cn),
diff --git a/deal.II/source/lac/sparse_direct.cc b/deal.II/source/lac/sparse_direct.cc
index 5bdb4eed7a..200a2d09de 100644
--- a/deal.II/source/lac/sparse_direct.cc
+++ b/deal.II/source/lac/sparse_direct.cc
@@ -2106,7 +2106,7 @@ void SparseDirectMA27::solve (const SparseMatrix<double> &matrix,
                               Vector<double>             &rhs_and_solution);
 
 template
-void SparseDirectMA27::solve (const SparseMatrix<float>  &matrix,
+void SparseDirectMA27::solve (const SparseMatrix<float> &matrix,
                               Vector<double>             &rhs_and_solution);
 
 
diff --git a/deal.II/source/lac/trilinos_solver.cc b/deal.II/source/lac/trilinos_solver.cc
index 3f1e06cf19..061a203401 100644
--- a/deal.II/source/lac/trilinos_solver.cc
+++ b/deal.II/source/lac/trilinos_solver.cc
@@ -37,7 +37,7 @@ namespace TrilinosWrappers
 
 
 
-  SolverBase::SolverBase (SolverControl  &cn)
+  SolverBase::SolverBase (SolverControl &cn)
     :
     solver_name    (gmres),
     solver_control (cn)
@@ -384,7 +384,7 @@ namespace TrilinosWrappers
 
 
 
-  SolverDirect::SolverDirect (SolverControl  &cn,
+  SolverDirect::SolverDirect (SolverControl &cn,
                               const AdditionalData &data)
     :
     solver_control (cn),
diff --git a/deal.II/source/numerics/data_postprocessor.cc b/deal.II/source/numerics/data_postprocessor.cc
index 99149e8bd2..5a8a00b3f8 100644
--- a/deal.II/source/numerics/data_postprocessor.cc
+++ b/deal.II/source/numerics/data_postprocessor.cc
@@ -27,10 +27,10 @@ DataPostprocessor<dim>::~DataPostprocessor()
 template <int dim>
 void
 DataPostprocessor<dim>::
-compute_derived_quantities_scalar (const std::vector<double>         &/*uh*/,
-                                   const std::vector<Tensor<1,dim> > &/*duh*/,
-                                   const std::vector<Tensor<2,dim> > &/*dduh*/,
-                                   const std::vector<Point<dim> >    &/*normals*/,
+compute_derived_quantities_scalar (const std::vector<double> &         /*uh*/,
+                                   const std::vector<Tensor<1,dim> > & /*duh*/,
+                                   const std::vector<Tensor<2,dim> > & /*dduh*/,
+                                   const std::vector<Point<dim> > &    /*normals*/,
                                    std::vector<Vector<double> >      &computed_quantities) const
 {
   computed_quantities.clear();
@@ -45,7 +45,7 @@ compute_derived_quantities_scalar (const std::vector<double>         &uh,
                                    const std::vector<Tensor<1,dim> > &duh,
                                    const std::vector<Tensor<2,dim> > &dduh,
                                    const std::vector<Point<dim> >    &normals,
-                                   const std::vector<Point<dim> >    &/*evaluation_points*/,
+                                   const std::vector<Point<dim> > &    /*evaluation_points*/,
                                    std::vector<Vector<double> >      &computed_quantities) const
 {
   compute_derived_quantities_scalar(uh, duh, dduh, normals, computed_quantities);
@@ -56,10 +56,10 @@ compute_derived_quantities_scalar (const std::vector<double>         &uh,
 template <int dim>
 void
 DataPostprocessor<dim>::
-compute_derived_quantities_vector (const std::vector<Vector<double> > &/*uh*/,
-                                   const std::vector<std::vector<Tensor<1,dim> > > &/*duh*/,
-                                   const std::vector<std::vector<Tensor<2,dim> > > &/*dduh*/,
-                                   const std::vector<Point<dim> >                  &/*normals*/,
+compute_derived_quantities_vector (const std::vector<Vector<double> > & /*uh*/,
+                                   const std::vector<std::vector<Tensor<1,dim> > > & /*duh*/,
+                                   const std::vector<std::vector<Tensor<2,dim> > > & /*dduh*/,
+                                   const std::vector<Point<dim> > &                  /*normals*/,
                                    std::vector<Vector<double> >                    &computed_quantities) const
 {
   computed_quantities.clear();
@@ -75,7 +75,7 @@ compute_derived_quantities_vector (const std::vector<Vector<double> > &uh,
                                    const std::vector<std::vector<Tensor<1,dim> > > &duh,
                                    const std::vector<std::vector<Tensor<2,dim> > > &dduh,
                                    const std::vector<Point<dim> >                  &normals,
-                                   const std::vector<Point<dim> >                  &/*evaluation_points*/,
+                                   const std::vector<Point<dim> > &                  /*evaluation_points*/,
                                    std::vector<Vector<double> >                    &computed_quantities) const
 {
   compute_derived_quantities_vector(uh, duh, dduh, normals, computed_quantities);
diff --git a/deal.II/source/numerics/derivative_approximation.cc b/deal.II/source/numerics/derivative_approximation.cc
index a0a72efe1a..1f4c109a2d 100644
--- a/deal.II/source/numerics/derivative_approximation.cc
+++ b/deal.II/source/numerics/derivative_approximation.cc
@@ -65,7 +65,7 @@ template <class InputVector, int spacedim>
 inline
 typename DerivativeApproximation::Gradient<dim>::ProjectedDerivative
 DerivativeApproximation::Gradient<dim>::
-get_projected_derivative (const FEValues<dim,spacedim>  &fe_values,
+get_projected_derivative (const FEValues<dim,spacedim> &fe_values,
                           const InputVector    &solution,
                           const unsigned int    component)
 {
@@ -114,7 +114,7 @@ template <class InputVector, int spacedim>
 inline
 typename DerivativeApproximation::SecondDerivative<dim>::ProjectedDerivative
 DerivativeApproximation::SecondDerivative<dim>::
-get_projected_derivative (const FEValues<dim,spacedim>  &fe_values,
+get_projected_derivative (const FEValues<dim,spacedim> &fe_values,
                           const InputVector    &solution,
                           const unsigned int    component)
 {
@@ -396,7 +396,7 @@ template <class InputVector, int spacedim>
 inline
 typename DerivativeApproximation::ThirdDerivative<dim>::ProjectedDerivative
 DerivativeApproximation::ThirdDerivative<dim>::
-get_projected_derivative (const FEValues<dim,spacedim>  &fe_values,
+get_projected_derivative (const FEValues<dim,spacedim> &fe_values,
                           const InputVector    &solution,
                           const unsigned int    component)
 {
@@ -703,7 +703,7 @@ approximate_cell (const Mapping<dim,spacedim>                   &mapping,
                   const DH<dim,spacedim>                        &dof_handler,
                   const InputVector                             &solution,
                   const unsigned int                             component,
-                  const typename DH<dim,spacedim>::active_cell_iterator  &cell,
+                  const typename DH<dim,spacedim>::active_cell_iterator &cell,
                   typename DerivativeDescription::Derivative    &derivative)
 {
   QMidpoint<dim> midpoint_rule;
diff --git a/deal.II/source/numerics/error_estimator.cc b/deal.II/source/numerics/error_estimator.cc
index 374b010914..786fbc34ec 100644
--- a/deal.II/source/numerics/error_estimator.cc
+++ b/deal.II/source/numerics/error_estimator.cc
@@ -1140,14 +1140,14 @@ estimate (const DH   &dof_handler,
 template <int spacedim>
 template <typename InputVector, class DH>
 void KellyErrorEstimator<1,spacedim>::
-estimate (const Mapping<1,spacedim>                    &/*mapping*/,
-          const DH                            &/*dof_handler*/,
+estimate (const Mapping<1,spacedim> &                    /*mapping*/,
+          const DH &                            /*dof_handler*/,
           const hp::QCollection<0> &,
-          const typename FunctionMap<spacedim>::type          &/*neumann_bc*/,
-          const std::vector<const InputVector *> &/*solutions*/,
-          std::vector<Vector<float>*>            &/*errors*/,
-          const ComponentMask                &/*component_mask_*/,
-          const Function<spacedim>                   */*coefficient*/,
+          const typename FunctionMap<spacedim>::type &          /*neumann_bc*/,
+          const std::vector<const InputVector *> & /*solutions*/,
+          std::vector<Vector<float>*> &            /*errors*/,
+          const ComponentMask &                /*component_mask_*/,
+          const Function<spacedim> *                   /*coefficient*/,
           const unsigned int,
           const types::subdomain_id          /*subdomain_id*/,
           const types::material_id                   /*material_id*/)
diff --git a/deal.II/source/numerics/histogram.cc b/deal.II/source/numerics/histogram.cc
index b80e62aa3e..ca6b9cf992 100644
--- a/deal.II/source/numerics/histogram.cc
+++ b/deal.II/source/numerics/histogram.cc
@@ -332,7 +332,7 @@ void Histogram::evaluate<float> (const Vector<float>   &values,
 
 // explicit instantiations for double
 template
-void Histogram::evaluate<double> (const std::vector<Vector<double> >  &values,
+void Histogram::evaluate<double> (const std::vector<Vector<double> > &values,
                                   const std::vector<double>                    &y_values,
                                   const unsigned int                            n_intervals,
                                   const IntervalSpacing                         interval_spacing);
diff --git a/deal.II/source/numerics/matrix_tools.cc b/deal.II/source/numerics/matrix_tools.cc
index 24bb447f41..9fe2a393d1 100644
--- a/deal.II/source/numerics/matrix_tools.cc
+++ b/deal.II/source/numerics/matrix_tools.cc
@@ -1123,11 +1123,11 @@ namespace MatrixCreator
 
   template <int dim, int spacedim>
   void
-  create_boundary_mass_matrix (const Mapping<dim, spacedim>  &mapping,
+  create_boundary_mass_matrix (const Mapping<dim, spacedim> &mapping,
                                const DoFHandler<dim,spacedim> &dof,
                                const Quadrature<dim-1>  &q,
-                               SparseMatrix<double>  &matrix,
-                               const typename FunctionMap<spacedim>::type  &boundary_functions,
+                               SparseMatrix<double> &matrix,
+                               const typename FunctionMap<spacedim>::type &boundary_functions,
                                Vector<double>            &rhs_vector,
                                std::vector<unsigned int> &dof_to_boundary_mapping,
                                const Function<spacedim> *const coefficient,
@@ -1934,7 +1934,7 @@ namespace MatrixTools
   template <typename number>
   void
   apply_boundary_values (const std::map<unsigned int,double> &boundary_values,
-                         SparseMatrix<number>  &matrix,
+                         SparseMatrix<number> &matrix,
                          Vector<number>   &solution,
                          Vector<number>   &right_hand_side,
                          const bool        eliminate_columns)
@@ -2100,7 +2100,7 @@ namespace MatrixTools
   template <typename number>
   void
   apply_boundary_values (const std::map<unsigned int,double> &boundary_values,
-                         BlockSparseMatrix<number>  &matrix,
+                         BlockSparseMatrix<number> &matrix,
                          BlockVector<number>   &solution,
                          BlockVector<number>   &right_hand_side,
                          const bool             eliminate_columns)
@@ -2575,8 +2575,8 @@ namespace MatrixTools
 
 
   void
-  apply_boundary_values (const std::map<unsigned int,double>  &boundary_values,
-                         PETScWrappers::MPI::BlockSparseMatrix  &matrix,
+  apply_boundary_values (const std::map<unsigned int,double> &boundary_values,
+                         PETScWrappers::MPI::BlockSparseMatrix &matrix,
                          PETScWrappers::MPI::BlockVector        &solution,
                          PETScWrappers::MPI::BlockVector        &right_hand_side,
                          const bool                            eliminate_columns)
@@ -2887,8 +2887,8 @@ namespace MatrixTools
 
 
   void
-  apply_boundary_values (const std::map<unsigned int,double>  &boundary_values,
-                         TrilinosWrappers::BlockSparseMatrix  &matrix,
+  apply_boundary_values (const std::map<unsigned int,double> &boundary_values,
+                         TrilinosWrappers::BlockSparseMatrix &matrix,
                          TrilinosWrappers::BlockVector        &solution,
                          TrilinosWrappers::BlockVector        &right_hand_side,
                          const bool                            eliminate_columns)
@@ -2901,8 +2901,8 @@ namespace MatrixTools
 
 
   void
-  apply_boundary_values (const std::map<unsigned int,double>  &boundary_values,
-                         TrilinosWrappers::BlockSparseMatrix  &matrix,
+  apply_boundary_values (const std::map<unsigned int,double> &boundary_values,
+                         TrilinosWrappers::BlockSparseMatrix &matrix,
                          TrilinosWrappers::MPI::BlockVector   &solution,
                          TrilinosWrappers::MPI::BlockVector   &right_hand_side,
                          const bool                            eliminate_columns)
@@ -3040,14 +3040,14 @@ namespace MatrixTools
   template
   void
   apply_boundary_values<double> (const std::map<unsigned int,double> &boundary_values,
-                                 SparseMatrix<double>  &matrix,
+                                 SparseMatrix<double> &matrix,
                                  Vector<double>   &solution,
                                  Vector<double>   &right_hand_side,
                                  const bool        eliminate_columns);
   template
   void
   apply_boundary_values<float> (const std::map<unsigned int,double> &boundary_values,
-                                SparseMatrix<float>  &matrix,
+                                SparseMatrix<float> &matrix,
                                 Vector<float>   &solution,
                                 Vector<float>   &right_hand_side,
                                 const bool        eliminate_columns);
@@ -3055,14 +3055,14 @@ namespace MatrixTools
   template
   void
   apply_boundary_values<double> (const std::map<unsigned int,double> &boundary_values,
-                                 BlockSparseMatrix<double>  &matrix,
+                                 BlockSparseMatrix<double> &matrix,
                                  BlockVector<double>   &solution,
                                  BlockVector<double>   &right_hand_side,
                                  const bool        eliminate_columns);
   template
   void
   apply_boundary_values<float> (const std::map<unsigned int,double> &boundary_values,
-                                BlockSparseMatrix<float>  &matrix,
+                                BlockSparseMatrix<float> &matrix,
                                 BlockVector<float>   &solution,
                                 BlockVector<float>   &right_hand_side,
                                 const bool        eliminate_columns);