Only formatting changes.

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

diff --git a/deal.II/deal.II/source/fe/fe_lib.quadratic.cc b/deal.II/deal.II/source/fe/fe_lib.quadratic.cc
index 5755c6752f28090e3fee3c3505a2d8d858f088dc..229af02ed4bad82bf7e2a3cfd7affbf6bb7e96d6 100644 (file)
--- a/deal.II/deal.II/source/fe/fe_lib.quadratic.cc
+++ b/deal.II/deal.II/source/fe/fe_lib.quadratic.cc
@@ -34,6 +34,7 @@ FEQ2<1>::FEQ2 () :
 };
 
 
+
 template <>
 FEQ2<1>::FEQ2 (const int) :
                FEQ1Mapping<1> (0, 3, 0, 0, 1,
@@ -43,6 +44,7 @@ FEQ2<1>::FEQ2 (const int) :
 };
 
 
+
 template <>
 void FEQ2<1>::initialize_matrices ()
 {
@@ -127,6 +129,7 @@ void FEQ2<1>::initialize_matrices ()
 };
 
 
+
 template <>
 double
 FEQ2<1>::shape_value(const unsigned int i,
@@ -144,6 +147,7 @@ FEQ2<1>::shape_value(const unsigned int i,
 };
 
 
+
 template <>
 Tensor<1,1>
 FEQ2<1>::shape_grad(const unsigned int i,
@@ -166,6 +170,7 @@ FEQ2<1>::shape_grad(const unsigned int i,
 };
 
 
+
 template <>
 Tensor<2,1>
 FEQ2<1>::shape_grad_grad (const unsigned int i,
@@ -190,29 +195,37 @@ FEQ2<1>::shape_grad_grad (const unsigned int i,
 };
 
 
+
 template <>
-void FEQ2<1>::get_unit_support_points (vector<Point<1> > &unit_points) const {
+void FEQ2<1>::get_unit_support_points (vector<Point<1> > &unit_points) const
+{
   FiniteElement<1>::get_unit_support_points (unit_points);
 };
 
 
+
 template <>
 void FEQ2<1>::get_support_points (const DoFHandler<1>::cell_iterator &cell,
-                                 vector<Point<1> >  &support_points) const {
+                                 vector<Point<1> >  &support_points) const
+{
   FiniteElement<1>::get_support_points (cell, support_points);
 };
 
 
+
 template <>
 void FEQ2<1>::get_face_support_points (const DoFHandler<1>::face_iterator &,
-                                      vector<Point<1> >  &) const {
+                                      vector<Point<1> >  &) const
+{
   Assert (false, ExcInternalError());
 };
 
 
+
 template <>
 void FEQ2<1>::get_local_mass_matrix (const DoFHandler<1>::cell_iterator &cell,
-                                    FullMatrix<double> &local_mass_matrix) const {
+                                    FullMatrix<double> &local_mass_matrix) const
+{
   Assert (local_mass_matrix.n() == dofs_per_cell,
          ExcWrongFieldDimension(local_mass_matrix.n(),dofs_per_cell));
   Assert (local_mass_matrix.m() == dofs_per_cell,
@@ -250,6 +263,7 @@ FEQ2<2>::FEQ2 () :
 };
 
 
+
 template <>
 FEQ2<2>::FEQ2 (const int) :
                FEQ1Mapping<2> (0, 0, 9, 0, 1,
@@ -259,6 +273,7 @@ FEQ2<2>::FEQ2 (const int) :
 };
 
 
+
 template <>
 void FEQ2<2>::initialize_matrices ()
 {
@@ -645,10 +660,11 @@ void FEQ2<2>::initialize_matrices ()
 };
 
 
+
 template <>
 double
 FEQ2<2>::shape_value (const unsigned int i,
-                               const Point<2>    &p) const
+                     const Point<2>    &p) const
 {
   Assert (i<dofs_per_cell, ExcIndexRange(i, 0, dofs_per_cell));
 
@@ -670,10 +686,11 @@ FEQ2<2>::shape_value (const unsigned int i,
 };
 
 
+
 template <>
 Tensor<1,2>
 FEQ2<2>::shape_grad (const unsigned int i,
-                              const Point<2>    &p) const
+                    const Point<2>    &p) const
 {
   Assert (i<dofs_per_cell, ExcIndexRange(i, 0, dofs_per_cell));
 
@@ -709,10 +726,11 @@ FEQ2<2>::shape_grad (const unsigned int i,
 };
 
 
+
 template <>
 Tensor<2,2>
 FEQ2<2>::shape_grad_grad (const unsigned int i,
-                                   const Point<2>    &p) const
+                         const Point<2>    &p) const
 {
   Assert (i<dofs_per_cell, ExcIndexRange(i, 0, dofs_per_cell));
 
@@ -781,9 +799,11 @@ FEQ2<2>::shape_grad_grad (const unsigned int i,
 };
 
 
+
 template <>
 void FEQ2<2>::get_local_mass_matrix (const DoFHandler<2>::cell_iterator &cell,
-                                              FullMatrix<double> &local_mass_matrix) const {
+                                    FullMatrix<double> &local_mass_matrix) const
+{
   Assert (local_mass_matrix.n() == dofs_per_cell,
          ExcWrongFieldDimension(local_mass_matrix.n(),dofs_per_cell));
   Assert (local_mass_matrix.m() == dofs_per_cell,
@@ -1055,8 +1075,10 @@ void FEQ2<2>::get_local_mass_matrix (const DoFHandler<2>::cell_iterator &cell,
 };
 
 
+
 template <>
-void FEQ2<2>::get_unit_support_points (vector<Point<2> > &unit_points) const {
+void FEQ2<2>::get_unit_support_points (vector<Point<2> > &unit_points) const
+{
   Assert (unit_points.size() == dofs_per_cell,
          ExcWrongFieldDimension (unit_points.size(), dofs_per_cell));
   
@@ -1072,9 +1094,11 @@ void FEQ2<2>::get_unit_support_points (vector<Point<2> > &unit_points) const {
 };
 
 
+
 template <>
 void FEQ2<2>::get_support_points (const DoFHandler<2>::cell_iterator &cell,
-                                 vector<Point<2> >  &support_points) const {
+                                 vector<Point<2> >  &support_points) const
+{
   Assert (support_points.size() == dofs_per_cell,
          ExcWrongFieldDimension (support_points.size(), dofs_per_cell));
   
@@ -1099,9 +1123,11 @@ void FEQ2<2>::get_support_points (const DoFHandler<2>::cell_iterator &cell,
 };
 
 
+
 template <>
 void FEQ2<2>::get_face_support_points (const DoFHandler<2>::face_iterator &face,
-                                      vector<Point<2> >  &support_points) const {
+                                      vector<Point<2> >  &support_points) const
+{
   Assert (support_points.size() == dofs_per_face,
          ExcWrongFieldDimension (support_points.size(), dofs_per_face));
 
@@ -1202,6 +1228,7 @@ FEQ2<3>::FEQ2 () :
 };
 
 
+
 template <>
 FEQ2<3>::FEQ2 (const int) :
                FEQ1Mapping<3> (0, 0, 0, 27, 1,
@@ -1211,6 +1238,7 @@ FEQ2<3>::FEQ2 (const int) :
 };
 
 
+
 template <>
 void FEQ2<3>::initialize_matrices ()
 {
@@ -2282,10 +2310,11 @@ void FEQ2<3>::initialize_matrices ()
 };
 
 
+
 template <>
 double
 FEQ2<3>::shape_value (const unsigned int i,
-                               const Point<3>    &p) const
+                     const Point<3>    &p) const
 {
   Assert (i<dofs_per_cell, ExcIndexRange(i, 0, dofs_per_cell));
 
@@ -2375,10 +2404,11 @@ xi)*eta*eta)*zeta+((-16.0+48.0*xi-32.0*xi*xi)*eta+(16.0-48.0*xi+32.0*xi*xi)*eta
 };
 
 
+
 template <>
 Tensor<1,3>
 FEQ2<3>::shape_grad (const unsigned int i,
-                              const Point<3>    &p) const
+                    const Point<3>    &p) const
 {
   Assert (i<dofs_per_cell, ExcIndexRange(i, 0, dofs_per_cell));
 
@@ -2478,10 +2508,11 @@ FEQ2<3>::shape_grad (const unsigned int i,
 };
 
 
+
 template <>
 Tensor<2,3>
 FEQ2<3>::shape_grad_grad (const unsigned int i,
-                                   const Point<3>    &p) const
+                         const Point<3>    &p) const
 {
   Assert (i<dofs_per_cell, ExcIndexRange(i, 0, dofs_per_cell));
 
@@ -2761,9 +2792,11 @@ return_value[2][2] = 2.0*(-64.0*xi+64.0*xi*xi)*eta+2.0*(64.0*xi-64.0*xi*xi)*eta*
 };
 
 
+
 template <>
 void FEQ2<3>::get_local_mass_matrix (const DoFHandler<3>::cell_iterator &,
-                                              FullMatrix<double> &local_mass_matrix) const {
+                                    FullMatrix<double> &local_mass_matrix) const
+{
   Assert (local_mass_matrix.n() == dofs_per_cell,
          ExcWrongFieldDimension(local_mass_matrix.n(),dofs_per_cell));
   Assert (local_mass_matrix.m() == dofs_per_cell,
@@ -2773,8 +2806,10 @@ void FEQ2<3>::get_local_mass_matrix (const DoFHandler<3>::cell_iterator &,
 };
 
 
+
 template <>
-void FEQ2<3>::get_unit_support_points (vector<Point<3> > &unit_points) const {
+void FEQ2<3>::get_unit_support_points (vector<Point<3> > &unit_points) const
+{
   Assert (unit_points.size() == dofs_per_cell,
          ExcWrongFieldDimension (unit_points.size(), dofs_per_cell));
   
@@ -2808,9 +2843,11 @@ void FEQ2<3>::get_unit_support_points (vector<Point<3> > &unit_points) const {
 };
 
 
+
 template <>
 void FEQ2<3>::get_support_points (const DoFHandler<3>::cell_iterator &cell,
-                                 vector<Point<3> >  &support_points) const {
+                                 vector<Point<3> >  &support_points) const
+{
   Assert (support_points.size() == dofs_per_cell,
          ExcWrongFieldDimension (support_points.size(), dofs_per_cell));
 
@@ -2907,9 +2944,11 @@ void FEQ2<3>::get_support_points (const DoFHandler<3>::cell_iterator &cell,
 };
 
 
+
 template <>
 void FEQ2<3>::get_face_support_points (const DoFHandler<3>::face_iterator &face,
-                                      vector<Point<3> >  &support_points) const {
+                                      vector<Point<3> >  &support_points) const
+{
   Assert (support_points.size() == dofs_per_face,
          ExcWrongFieldDimension (support_points.size(), dofs_per_face));