]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
revert more whitespace crap
authorheister <heister@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 21 Dec 2012 01:40:23 +0000 (01:40 +0000)
committerheister <heister@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 21 Dec 2012 01:40:23 +0000 (01:40 +0000)
git-svn-id: https://svn.dealii.org/branches/branch_merge_mg_into_dof_handler@27841 0785d39b-7218-0410-832d-ea1e28bc413d

34 files changed:
deal.II/source/base/data_out_base.cc
deal.II/source/base/function_lib.cc
deal.II/source/base/function_parser.cc
deal.II/source/base/parameter_handler.cc
deal.II/source/base/parsed_function.cc
deal.II/source/base/polynomial_space.cc
deal.II/source/fe/fe_dgp_nonparametric.cc
deal.II/source/fe/fe_nothing.cc
deal.II/source/fe/fe_poly_tensor.cc
deal.II/source/fe/fe_q.cc
deal.II/source/fe/fe_q_hierarchical.cc
deal.II/source/fe/fe_raviart_thomas_nodal.cc
deal.II/source/fe/fe_system.cc
deal.II/source/fe/mapping_cartesian.cc
deal.II/source/fe/mapping_q.cc
deal.II/source/fe/mapping_q1.cc
deal.II/source/grid/grid_generator.cc
deal.II/source/grid/grid_in.cc
deal.II/source/grid/grid_reordering.cc
deal.II/source/grid/tria.cc
deal.II/source/grid/tria_boundary.cc
deal.II/source/lac/constraint_matrix.cc
deal.II/source/lac/petsc_matrix_free.cc
deal.II/source/lac/petsc_parallel_vector.cc
deal.II/source/lac/petsc_solver.cc
deal.II/source/lac/petsc_vector_base.cc
deal.II/source/lac/slepc_solver.cc
deal.II/source/lac/sparse_direct.cc
deal.II/source/lac/trilinos_solver.cc
deal.II/source/numerics/data_postprocessor.cc
deal.II/source/numerics/derivative_approximation.cc
deal.II/source/numerics/error_estimator.cc
deal.II/source/numerics/histogram.cc
deal.II/source/numerics/matrix_tools.cc

index ba57a738b61490fbf08d88eed7fada4c0188c9e9..5a5d9c9dd740a0aa039366c606927db0792975f3 100644 (file)
@@ -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());
 
index 3d8b0e073624852020dda7dc485871828c7f8655..6f4f5cc5f8f23cb0ba2e16cff4e968570396f9a1 100644 (file)
@@ -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(),
index 7dc72e6e071d58afe10d506ce05e4efca36ae607..5b4b800a46e56ba3d5925d22d7fdf8c5ae75c602 100644 (file)
@@ -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());
index 8227166e21a6726ec85790c53d43215a74759471..1975331b27c59e69aa94711eaf04eeebe84c6284 100644 (file)
@@ -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)
     :
index 42540badf2967455f16ea16d8fb1ebee6922d492..da058af424076b2a813ca40f2270fec63ffd9f9c 100644 (file)
@@ -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());
index 2852f213af8d00b8057342e034cbe1d1794f0d99..f97d13f01fa30f514fafc4c68d9cbef367300ed4 100644 (file)
@@ -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);
index fccbec69f84824647d492e6b565e1d50f6957510..7f891fd4dcef2d3b9eef1971bf8406efe8a43df5 100644 (file)
@@ -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
index 71b32f25d719c96014a0dd7d4609cc7f15155dd6..527bc4f4cf7ba0c3e907cf7e73d6c0a94ebddad1 100644 (file)
@@ -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
index 8d5562849c6710ac953e6b185fa01f5522aaf9cd..00c281f83a12894f009de63d80fc309d2e881324 100644 (file)
@@ -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)
   {
index ede4249d59a565bdfdcf4fa6f2132ed75ed78428..93b8de41dc9570513f0479f8201313c88f9c3cd5 100644 (file)
@@ -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
 }
index ceef17fdac10a8aed9cf1310c7a23fadde49966e..0f712f1aec8081ab5c8690534ef7fbcf32186743 100644 (file)
@@ -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));
 }
index f8dfcd8bbb6fec88ec1799a32eece9e066145b27..23597aadf3d98d5534a384984f23d43522ad19b0 100644 (file)
@@ -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));
 }
index 88048074d2c1e84f6dfa8f89f85f50694b954915..08678aaa3cac4d9f350bfa8dfe783c08b12cbc15 100644 (file)
@@ -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;
         }
     }
index 48ec6cfaea29775f2c173f7c12e63d527d18a017..5529fc527b4eb3e59fee1b83df9eb1afae914b30 100644 (file)
@@ -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
index 213905d2a72be35a1610f8ac6de70f6a7d4af798..665748712bc0743ac998630a6b46edaffd8e664c 100644 (file)
@@ -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
index 953aea7fb390abb7a3ecd58ee84c35bdd69a7d51..211ad2c8b91a7a887c97ea07fc46a42da045849c 100644 (file)
@@ -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>();
index 8e0966cb7c3d272eab6a9d4525f53d11ae161bb4..a2b6d02c6386bb0e0e88620bb8b079a4eb486b7e 100644 (file)
@@ -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)
 {
index e61a7551aedbfc00ef9f5a28108d517d4fb03cfc..022757a363cfe8749debea93706f7e3c3c0ff83e 100644 (file)
@@ -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());
 }
index dd8aeee74d210faf8e9e5780023ebdc0c5c5e58a..f6b79371d04b000296b31ebfaaa070047b7b0f3d 100644 (file)
@@ -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)
     {
index 5f9f470826d9b0c90668390455f2a5fbfc1dc22b..6b0834b82b16cd9932af02f705923e177f4e7853 100644 (file)
@@ -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"));
index 247b780afced1cd986c46a289fb700576f0790f9..f79e2f4e03368e58fc10f39ee5ccbc473eaf5f93 100644 (file)
@@ -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;
index 8f69b8e04cc7495035de2a4c5099af319922e48e..550a0a393c23460986d8628aa7b333cbe0c964a1 100644 (file)
@@ -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());
index 16f51fe291e2374657a9898937d9826c29540b98..f688cf5517a8636fac56373149ae2185bfbf1ade 100755 (executable)
@@ -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
index fb25d77f615620ed4f31fe18568a2b37e35c16fd..225c1ebda18ef4016d881fdc45616d564113640f 100644 (file)
@@ -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)
index ad9f1dfd5069705fdd3477cc0c3d401401c6eb46..caaa9ea4b94d0ae054739af016ac612673f44be0 100644 (file)
@@ -42,7 +42,7 @@ namespace PETScWrappers
 
 
 
-  SolverBase::SolverBase (SolverControl  &cn,
+  SolverBase::SolverBase (SolverControl &cn,
                           const MPI_Comm &mpi_communicator)
     :
     solver_control (cn),
index 135319b0d90564724c1853a1e3be6eeeab17acbf..4b80d74837c2981c49e74085d05e1907e49d9086 100644 (file)
@@ -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"));
index 27c0895bcd147b43011c8484e8bb73961de3fc10..a0a5c6e865d4ab6b574237e9750e70b450a91bd8 100644 (file)
@@ -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),
index 5bdb4eed7a953aea4b7e75705a915eed2f31d821..200a2d09def86233c7b8bcf4118d42ed7bfba9d6 100644 (file)
@@ -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);
 
 
index 3f1e06cf197cec9145fdd249138a5855c9d627a7..061a203401119a7adac09bbfc9d71a5fb8e5d32e 100644 (file)
@@ -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),
index 99149e8bd2301627d5ace9443728110373018806..5a8a00b3f8bbd9cd0978b8ba75872d7f9bd2b337 100644 (file)
@@ -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);
index a0a72efe1a7d3309f54f9ffb1ed1e3ee3c89b14d..1f4c109a2d7d9ad5586e8c830efbc31d00b37e83 100644 (file)
@@ -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;
index 374b01091480912702be529f2f3bf3eb7c504304..786fbc34ec94962d10834828c04fa8e0d298f5bb 100644 (file)
@@ -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*/)
index b80e62aa3e8d5435556381e9d0d527769468fb03..ca6b9cf992c13fe266d18ede880102563ab6bc89 100644 (file)
@@ -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);
index 24bb447f4173e705a8dddca929862c9c5bfb8a21..9fe2a393d1ec7171ef113040336b71e0f1974f35 100644 (file)
@@ -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);

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.