]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Make things simpler for out documentation tools.
authorwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 13 Dec 1999 15:37:11 +0000 (15:37 +0000)
committerwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 13 Dec 1999 15:37:11 +0000 (15:37 +0000)
git-svn-id: https://svn.dealii.org/trunk@2052 0785d39b-7218-0410-832d-ea1e28bc413d

12 files changed:
deal.II/base/include/base/logstream.h
deal.II/deal.II/include/dofs/dof_accessor.templates.h
deal.II/deal.II/include/grid/tria_accessor.h
deal.II/deal.II/include/grid/tria_accessor.templates.h
deal.II/deal.II/include/multigrid/mg_base.h
deal.II/lac/include/lac/block_vector.h
deal.II/lac/include/lac/block_vector.templates.h
deal.II/lac/include/lac/mgbase.h
deal.II/lac/include/lac/precondition.h
deal.II/lac/include/lac/precondition_block.h
deal.II/lac/include/lac/precondition_block.templates.h
deal.II/lac/include/lac/precondition_selector.h

index 3eeb490e18fc3c549f7e664571f03309854b0272..14cc58019f579728ecb158e923869c086da0a70f 100644 (file)
@@ -163,7 +163,7 @@ class LogStream
                                      */
     template <typename T>
     LogStream & operator << (const T &t);
-    
+
                                     /**
                                      * Output a function. This really is not
                                      * to output the function, but calls the
index 0ecbae29e59ab22893880770e481b831a7752c31..161b79d68356cff45da0d044d7af1bb1ad7dad46 100644 (file)
@@ -23,7 +23,7 @@
 
 template <int dim>
 inline
-int DoFObjectAccessor<1, dim>::dof_index (const unsigned int i) const {
+int DoFObjectAccessor<1,dim>::dof_index (const unsigned int i) const {
   Assert (dof_handler != 0, ExcInvalidObject());
                                   // make sure a FE has been selected
                                   // and enough room was reserved
@@ -40,7 +40,7 @@ int DoFObjectAccessor<1, dim>::dof_index (const unsigned int i) const {
 
 template <int dim>
 inline
-int DoFObjectAccessor<1, dim>::vertex_dof_index (const unsigned int vertex,
+int DoFObjectAccessor<1,dim>::vertex_dof_index (const unsigned int vertex,
                                                      const unsigned int i) const {
   Assert (dof_handler != 0, ExcInvalidObject());
   Assert (dof_handler->selected_fe != 0, ExcInvalidObject());
@@ -59,7 +59,7 @@ int DoFObjectAccessor<1, dim>::vertex_dof_index (const unsigned int vertex,
 template <int dim>
 inline
 void
-DoFObjectAccessor<1, dim>::get_dof_indices (vector<int> &dof_indices) const {
+DoFObjectAccessor<1,dim>::get_dof_indices (vector<int> &dof_indices) const {
   Assert (dof_handler != 0, ExcInvalidObject());
   Assert (dof_handler->selected_fe != 0, ExcInvalidObject());
   Assert (dof_indices.size() == (2*dof_handler->get_fe().dofs_per_vertex +
@@ -80,9 +80,9 @@ DoFObjectAccessor<1, dim>::get_dof_indices (vector<int> &dof_indices) const {
 
 template <int dim>
 inline
-TriaIterator<dim,DoFObjectAccessor<1, dim> >
-DoFObjectAccessor<1, dim>::child (const unsigned int i) const {
-  TriaIterator<dim,DoFObjectAccessor<1, dim> > q (tria,
+TriaIterator<dim,DoFObjectAccessor<1,dim> >
+DoFObjectAccessor<1,dim>::child (const unsigned int i) const {
+  TriaIterator<dim,DoFObjectAccessor<1,dim> > q (tria,
                                                       present_level+1,
                                                       child_index (i),
                                                       dof_handler);
@@ -99,7 +99,7 @@ DoFObjectAccessor<1, dim>::child (const unsigned int i) const {
 template <int dim>
 inline
 void
-DoFObjectAccessor<1, dim>::copy_from (const DoFObjectAccessor<1, dim> &a) {
+DoFObjectAccessor<1,dim>::copy_from (const DoFObjectAccessor<1,dim> &a) {
   BaseClass::copy_from (a);
   set_dof_handler (a.dof_handler);
 };
@@ -107,11 +107,11 @@ DoFObjectAccessor<1, dim>::copy_from (const DoFObjectAccessor<1, dim> &a) {
 
 
 
-/*------------------------- Functions: DoFObjectAccessor<2, dim> -----------------------*/
+/*------------------------- Functions: DoFObjectAccessor<2,dim> -----------------------*/
 
 template <int dim>
 inline
-int DoFObjectAccessor<2, dim>::dof_index (const unsigned int i) const {
+int DoFObjectAccessor<2,dim>::dof_index (const unsigned int i) const {
   Assert (dof_handler != 0, ExcInvalidObject());
                                   // make sure a FE has been selected
                                   // and enough room was reserved
@@ -127,7 +127,7 @@ int DoFObjectAccessor<2, dim>::dof_index (const unsigned int i) const {
 
 template <int dim>
 inline
-int DoFObjectAccessor<2, dim>::vertex_dof_index (const unsigned int vertex,
+int DoFObjectAccessor<2,dim>::vertex_dof_index (const unsigned int vertex,
                                                      const unsigned int i) const {
   Assert (dof_handler != 0, ExcInvalidObject());
   Assert (dof_handler->selected_fe != 0, ExcInvalidObject());
@@ -146,7 +146,7 @@ int DoFObjectAccessor<2, dim>::vertex_dof_index (const unsigned int vertex,
 template <int dim>
 inline
 void
-DoFObjectAccessor<2, dim>::get_dof_indices (vector<int> &dof_indices) const {
+DoFObjectAccessor<2,dim>::get_dof_indices (vector<int> &dof_indices) const {
   Assert (dof_handler != 0, ExcInvalidObject());
   Assert (dof_handler->selected_fe != 0, ExcInvalidObject());
   Assert (dof_indices.size() == (4*dof_handler->get_fe().dofs_per_vertex +
@@ -172,11 +172,11 @@ DoFObjectAccessor<2, dim>::get_dof_indices (vector<int> &dof_indices) const {
 
 template <int dim>
 inline
-TriaIterator<dim,DoFObjectAccessor<1, dim> >
-DoFObjectAccessor<2, dim>::line (const unsigned int i) const {
+TriaIterator<dim,DoFObjectAccessor<1,dim> >
+DoFObjectAccessor<2,dim>::line (const unsigned int i) const {
   Assert (i<4, ExcIndexRange (i, 0, 4));
 
-  return TriaIterator<dim,DoFObjectAccessor<1, dim> >
+  return TriaIterator<dim,DoFObjectAccessor<1,dim> >
     (
       tria,
       present_level,
@@ -189,9 +189,9 @@ DoFObjectAccessor<2, dim>::line (const unsigned int i) const {
 
 template <int dim>
 inline
-TriaIterator<dim,DoFObjectAccessor<2, dim> >
-DoFObjectAccessor<2, dim>::child (const unsigned int i) const {
-  TriaIterator<dim,DoFObjectAccessor<2, dim> > q (tria,
+TriaIterator<dim,DoFObjectAccessor<2,dim> >
+DoFObjectAccessor<2,dim>::child (const unsigned int i) const {
+  TriaIterator<dim,DoFObjectAccessor<2,dim> > q (tria,
                                                  present_level+1,
                                                  child_index (i),
                                                  dof_handler);
@@ -208,7 +208,7 @@ DoFObjectAccessor<2, dim>::child (const unsigned int i) const {
 template <int dim>
 inline
 void
-DoFObjectAccessor<2, dim>::copy_from (const DoFObjectAccessor<2, dim> &a) {
+DoFObjectAccessor<2,dim>::copy_from (const DoFObjectAccessor<2,dim> &a) {
   BaseClass::copy_from (a);
   set_dof_handler (a.dof_handler);
 };
@@ -221,7 +221,7 @@ DoFObjectAccessor<2, dim>::copy_from (const DoFObjectAccessor<2, dim> &a) {
 
 template <int dim>
 inline
-int DoFObjectAccessor<3, dim>::dof_index (const unsigned int i) const {
+int DoFObjectAccessor<3,dim>::dof_index (const unsigned int i) const {
   Assert (dof_handler != 0, ExcInvalidObject());
                                   // make sure a FE has been selected
                                   // and enough room was reserved
@@ -237,7 +237,7 @@ int DoFObjectAccessor<3, dim>::dof_index (const unsigned int i) const {
 
 template <int dim>
 inline
-int DoFObjectAccessor<3, dim>::vertex_dof_index (const unsigned int vertex,
+int DoFObjectAccessor<3,dim>::vertex_dof_index (const unsigned int vertex,
                                                      const unsigned int i) const {
   Assert (dof_handler != 0, ExcInvalidObject());
   Assert (dof_handler->selected_fe != 0, ExcInvalidObject());
@@ -256,7 +256,7 @@ int DoFObjectAccessor<3, dim>::vertex_dof_index (const unsigned int vertex,
 template <int dim>
 inline
 void
-DoFObjectAccessor<3, dim>::get_dof_indices (vector<int> &dof_indices) const {
+DoFObjectAccessor<3,dim>::get_dof_indices (vector<int> &dof_indices) const {
   Assert (dof_handler != 0, ExcInvalidObject());
   Assert (dof_handler->selected_fe != 0, ExcInvalidObject());
   Assert (dof_indices.size() == (8*dof_handler->get_fe().dofs_per_vertex +
@@ -287,10 +287,10 @@ DoFObjectAccessor<3, dim>::get_dof_indices (vector<int> &dof_indices) const {
 
 template <int dim>
 inline
-TriaIterator<dim,DoFObjectAccessor<1, dim> >
-DoFObjectAccessor<3, dim>::line (const unsigned int i) const {
-  TriaIterator<dim,TriaObjectAccessor<1, dim> > l = BaseClass::line(i);
-  return TriaIterator<dim,DoFObjectAccessor<1, dim> >
+TriaIterator<dim,DoFObjectAccessor<1,dim> >
+DoFObjectAccessor<3,dim>::line (const unsigned int i) const {
+  TriaIterator<dim,TriaObjectAccessor<1,dim> > l = BaseClass::line(i);
+  return TriaIterator<dim,DoFObjectAccessor<1,dim> >
     (
       tria,
       present_level,
@@ -303,11 +303,11 @@ DoFObjectAccessor<3, dim>::line (const unsigned int i) const {
 
 template <int dim>
 inline
-TriaIterator<dim,DoFObjectAccessor<2, dim> >
-DoFObjectAccessor<3, dim>::quad (const unsigned int i) const {
+TriaIterator<dim,DoFObjectAccessor<2,dim> >
+DoFObjectAccessor<3,dim>::quad (const unsigned int i) const {
   Assert (i<6, ExcIndexRange (i, 0, 6));
 
-  return TriaIterator<dim,DoFObjectAccessor<2, dim> >
+  return TriaIterator<dim,DoFObjectAccessor<2,dim> >
     (
       tria,
       present_level,
@@ -320,9 +320,9 @@ DoFObjectAccessor<3, dim>::quad (const unsigned int i) const {
 
 template <int dim>
 inline
-TriaIterator<dim,DoFObjectAccessor<3, dim> >
-DoFObjectAccessor<3, dim>::child (const unsigned int i) const {
-  TriaIterator<dim,DoFObjectAccessor<3, dim> > q (tria,
+TriaIterator<dim,DoFObjectAccessor<3,dim> >
+DoFObjectAccessor<3,dim>::child (const unsigned int i) const {
+  TriaIterator<dim,DoFObjectAccessor<3,dim> > q (tria,
                                                  present_level+1,
                                                  child_index (i),
                                                  dof_handler);
@@ -337,7 +337,7 @@ DoFObjectAccessor<3, dim>::child (const unsigned int i) const {
 
 
 template <int dim>
-void DoFObjectAccessor<3, dim>::copy_from (const DoFObjectAccessor<3, dim> &a) {
+void DoFObjectAccessor<3,dim>::copy_from (const DoFObjectAccessor<3,dim> &a) {
   BaseClass::copy_from (a);
   set_dof_handler (a.dof_handler);
 };
index e3bf65f933b96775c6c275c345c0b8e3f2f20685..6d87bb6143b7b9d06411a68a1da51376cf5c38b5 100644 (file)
@@ -408,7 +408,7 @@ class TriaObjectAccessor :  public TriaAccessor<dim>
                                      *  Pointer to the #i#th
                                      *  child.
                                      */
-    TriaIterator<dim,TriaObjectAccessor<celldim, dim> >
+    TriaIterator<dim,TriaObjectAccessor<celldim,dim> >
     child (const unsigned int i) const;
 
                                     /**
index f1f5dcec9eb5514ba0acc18b615f8910acd49b24..d36f8cea2db81d258390c7c8894bd51d9f691a98 100644 (file)
@@ -102,7 +102,7 @@ TriaAccessor<dim>::get_triangulation () const
 template <int dim>
 inline
 bool
-TriaObjectAccessor<1, dim>::used () const {
+TriaObjectAccessor<1,dim>::used () const {
   Assert (state() == valid, ExcDereferenceInvalidObject());
   return tria->levels[present_level]->lines.used[present_index];
 };
@@ -113,7 +113,7 @@ TriaObjectAccessor<1, dim>::used () const {
 template <int dim>
 inline
 bool
-TriaObjectAccessor<1, dim>::user_flag_set () const {
+TriaObjectAccessor<1,dim>::user_flag_set () const {
   Assert (used(), ExcCellNotUsed());
   return tria->levels[present_level]->lines.user_flags[present_index];
 };
@@ -123,7 +123,7 @@ TriaObjectAccessor<1, dim>::user_flag_set () const {
 template <int dim>
 inline
 void
-TriaObjectAccessor<1, dim>::set_user_flag () const {
+TriaObjectAccessor<1,dim>::set_user_flag () const {
   Assert (used(), ExcCellNotUsed());
   tria->levels[present_level]->lines.user_flags[present_index] = true;
 };
@@ -133,7 +133,7 @@ TriaObjectAccessor<1, dim>::set_user_flag () const {
 template <int dim>
 inline
 void
-TriaObjectAccessor<1, dim>::clear_user_flag () const {
+TriaObjectAccessor<1,dim>::clear_user_flag () const {
   Assert (used(), ExcCellNotUsed());
   tria->levels[present_level]->lines.user_flags[present_index] = false;
 };
@@ -143,11 +143,11 @@ TriaObjectAccessor<1, dim>::clear_user_flag () const {
 
 template <int dim>
 inline
-TriaIterator<dim,TriaObjectAccessor<1, dim> >
-TriaObjectAccessor<1, dim>::child (const unsigned int i) const {
+TriaIterator<dim,TriaObjectAccessor<1,dim> >
+TriaObjectAccessor<1,dim>::child (const unsigned int i) const {
   Assert (i<2, ExcIndexRange(i,0,2));
   
-  TriaIterator<dim,TriaObjectAccessor<1, dim> > q (tria, present_level+1, child_index (i));
+  TriaIterator<dim,TriaObjectAccessor<1,dim> > q (tria, present_level+1, child_index (i));
   
 #ifdef DEBUG
   if (q.state() != past_the_end)
@@ -161,7 +161,7 @@ TriaObjectAccessor<1, dim>::child (const unsigned int i) const {
 template <int dim>
 inline
 int
-TriaObjectAccessor<1, dim>::child_index (unsigned int i) const {
+TriaObjectAccessor<1,dim>::child_index (unsigned int i) const {
   Assert (i<2, ExcIndexRange(i,0,2));
   return tria->levels[present_level]->lines.children[present_index]+i;
 };
@@ -172,7 +172,7 @@ TriaObjectAccessor<1, dim>::child_index (unsigned int i) const {
 template <int dim>
 inline
 bool
-TriaObjectAccessor<1, dim>::has_children () const {
+TriaObjectAccessor<1,dim>::has_children () const {
   Assert (state() == valid, ExcDereferenceInvalidObject());
   return (tria->levels[present_level]->lines.children[present_index] != -1);
 }
@@ -182,7 +182,7 @@ TriaObjectAccessor<1, dim>::has_children () const {
 template <int dim>
 inline
 unsigned int
-TriaObjectAccessor<1, dim>::max_refinement_depth () const
+TriaObjectAccessor<1,dim>::max_refinement_depth () const
 {
   if (!has_children())
     return 0;
@@ -198,7 +198,7 @@ TriaObjectAccessor<1, dim>::max_refinement_depth () const
 template <int dim>
 inline
 void
-TriaObjectAccessor<1, dim>::operator ++ () {
+TriaObjectAccessor<1,dim>::operator ++ () {
   ++present_index;
                                   // is index still in the range of
                                   // the vector?
@@ -225,7 +225,7 @@ TriaObjectAccessor<1, dim>::operator ++ () {
 template <int dim>
 inline
 void
-TriaObjectAccessor<1, dim>::operator -- () {
+TriaObjectAccessor<1,dim>::operator -- () {
   --present_index;
                                   // is index still in the range of
                                   // the vector?
@@ -255,7 +255,7 @@ TriaObjectAccessor<1, dim>::operator -- () {
 template <int dim>
 inline
 bool
-TriaObjectAccessor<2, dim>::used () const {
+TriaObjectAccessor<2,dim>::used () const {
   Assert (state() == valid, ExcDereferenceInvalidObject());
   return tria->levels[present_level]->quads.used[present_index];
 };
@@ -266,7 +266,7 @@ TriaObjectAccessor<2, dim>::used () const {
 template <int dim>
 inline
 bool
-TriaObjectAccessor<2, dim>::user_flag_set () const {
+TriaObjectAccessor<2,dim>::user_flag_set () const {
   Assert (used(), ExcCellNotUsed());
   return tria->levels[present_level]->quads.user_flags[present_index];
 };
@@ -276,7 +276,7 @@ TriaObjectAccessor<2, dim>::user_flag_set () const {
 template <int dim>
 inline
 void
-TriaObjectAccessor<2, dim>::set_user_flag () const {
+TriaObjectAccessor<2,dim>::set_user_flag () const {
   Assert (used(), ExcCellNotUsed());
   tria->levels[present_level]->quads.user_flags[present_index] = true;
 };
@@ -286,7 +286,7 @@ TriaObjectAccessor<2, dim>::set_user_flag () const {
 template <int dim>
 inline
 void
-TriaObjectAccessor<2, dim>::clear_user_flag () const {
+TriaObjectAccessor<2,dim>::clear_user_flag () const {
   Assert (used(), ExcCellNotUsed());
   tria->levels[present_level]->quads.user_flags[present_index] = false;
 };
@@ -296,10 +296,10 @@ TriaObjectAccessor<2, dim>::clear_user_flag () const {
 
 template <int dim>
 inline
-TriaIterator<dim,TriaObjectAccessor<1, dim> >
-TriaObjectAccessor<2, dim>::line (const unsigned int i) const {
+TriaIterator<dim,TriaObjectAccessor<1,dim> >
+TriaObjectAccessor<2,dim>::line (const unsigned int i) const {
   return
-    TriaIterator<dim,TriaObjectAccessor<1, dim> >
+    TriaIterator<dim,TriaObjectAccessor<1,dim> >
     (
       tria,
       present_level,
@@ -312,7 +312,7 @@ TriaObjectAccessor<2, dim>::line (const unsigned int i) const {
 template <int dim>
 inline
 unsigned int
-TriaObjectAccessor<2, dim>::line_index (unsigned int i) const {
+TriaObjectAccessor<2,dim>::line_index (unsigned int i) const {
   Assert (i<4, ExcIndexRange(i,0,4));
 
   return tria->levels[present_level]->quads.quads[present_index].line(i);
@@ -322,11 +322,11 @@ TriaObjectAccessor<2, dim>::line_index (unsigned int i) const {
 
 template <int dim>
 inline
-TriaIterator<dim,TriaObjectAccessor<2, dim> >
-TriaObjectAccessor<2, dim>::child (const unsigned int i) const {
+TriaIterator<dim,TriaObjectAccessor<2,dim> >
+TriaObjectAccessor<2,dim>::child (const unsigned int i) const {
   Assert (i<4, ExcIndexRange(i,0,4));
   
-  TriaIterator<dim,TriaObjectAccessor<2, dim> > q (tria, present_level+1, child_index (i));
+  TriaIterator<dim,TriaObjectAccessor<2,dim> > q (tria, present_level+1, child_index (i));
   
 #ifdef DEBUG
   if (q.state() != past_the_end)
@@ -339,7 +339,7 @@ TriaObjectAccessor<2, dim>::child (const unsigned int i) const {
 
 template <int dim>
 inline
-int TriaObjectAccessor<2, dim>::child_index (unsigned int i) const {
+int TriaObjectAccessor<2,dim>::child_index (unsigned int i) const {
   Assert (i<4, ExcIndexRange(i,0,4));
   return tria->levels[present_level]->quads.children[present_index]+i;
 };
@@ -350,7 +350,7 @@ int TriaObjectAccessor<2, dim>::child_index (unsigned int i) const {
 template <int dim>
 inline
 bool
-TriaObjectAccessor<2, dim>::has_children () const {
+TriaObjectAccessor<2,dim>::has_children () const {
   Assert (state() == valid, ExcDereferenceInvalidObject());
   return (tria->levels[present_level]->quads.children[present_index] != -1);
 };
@@ -360,7 +360,7 @@ TriaObjectAccessor<2, dim>::has_children () const {
 template <int dim>
 inline
 unsigned int
-TriaObjectAccessor<2, dim>::max_refinement_depth () const
+TriaObjectAccessor<2,dim>::max_refinement_depth () const
 {
   if (!has_children())
     return 0;
@@ -378,7 +378,7 @@ TriaObjectAccessor<2, dim>::max_refinement_depth () const
 template <int dim>
 inline
 void
-TriaObjectAccessor<2, dim>::operator ++ () {
+TriaObjectAccessor<2,dim>::operator ++ () {
   ++present_index;
                                   // is index still in the range of
                                   // the vector?
@@ -404,7 +404,7 @@ TriaObjectAccessor<2, dim>::operator ++ () {
 template <int dim>
 inline
 void
-TriaObjectAccessor<2, dim>::operator -- () {
+TriaObjectAccessor<2,dim>::operator -- () {
   --present_index;
                                   // is index still in the range of
                                   // the vector?
@@ -434,7 +434,7 @@ TriaObjectAccessor<2, dim>::operator -- () {
 template <int dim>
 inline
 bool
-TriaObjectAccessor<3, dim>::used () const {
+TriaObjectAccessor<3,dim>::used () const {
   Assert (state() == valid, ExcDereferenceInvalidObject());
   return tria->levels[present_level]->hexes.used[present_index];
 };
@@ -444,7 +444,7 @@ TriaObjectAccessor<3, dim>::used () const {
 template <int dim>
 inline
 bool
-TriaObjectAccessor<3, dim>::user_flag_set () const {
+TriaObjectAccessor<3,dim>::user_flag_set () const {
   Assert (used(), ExcCellNotUsed());
   return tria->levels[present_level]->hexes.user_flags[present_index];
 };
@@ -454,7 +454,7 @@ TriaObjectAccessor<3, dim>::user_flag_set () const {
 template <int dim>
 inline
 void
-TriaObjectAccessor<3, dim>::set_user_flag () const {
+TriaObjectAccessor<3,dim>::set_user_flag () const {
   Assert (used(), ExcCellNotUsed());
   tria->levels[present_level]->hexes.user_flags[present_index] = true;
 };
@@ -463,7 +463,7 @@ TriaObjectAccessor<3, dim>::set_user_flag () const {
 
 template <int dim>
 inline
-void TriaObjectAccessor<3, dim>::clear_user_flag () const {
+void TriaObjectAccessor<3,dim>::clear_user_flag () const {
   Assert (used(), ExcCellNotUsed());
   tria->levels[present_level]->hexes.user_flags[present_index] = false;
 };
@@ -473,8 +473,8 @@ void TriaObjectAccessor<3, dim>::clear_user_flag () const {
 
 template <int dim>
 inline
-TriaIterator<dim,TriaObjectAccessor<1, dim> >
-TriaObjectAccessor<3, dim>::line (const unsigned int i) const {
+TriaIterator<dim,TriaObjectAccessor<1,dim> >
+TriaObjectAccessor<3,dim>::line (const unsigned int i) const {
   Assert (used(), ExcCellNotUsed());
   Assert (i<12, ExcIndexRange (i,0,12));
 
@@ -504,18 +504,18 @@ TriaObjectAccessor<3, dim>::line (const unsigned int i) const {
                return this->quad(4)->line(3);
        };
   Assert (false, ExcIndexRange(i,0,12));
-  return TriaIterator<dim,TriaObjectAccessor<1, dim> >(tria, -1, -1, 0);
+  return TriaIterator<dim,TriaObjectAccessor<1,dim> >(tria, -1, -1, 0);
 };
 
 
 
 template <int dim>
 inline
-TriaIterator<dim,TriaObjectAccessor<2, dim> >
-TriaObjectAccessor<3, dim>::quad (const unsigned int i) const {
+TriaIterator<dim,TriaObjectAccessor<2,dim> >
+TriaObjectAccessor<3,dim>::quad (const unsigned int i) const {
   Assert (used(), ExcCellNotUsed());
   return
-    TriaIterator<dim,TriaObjectAccessor<2, dim> >
+    TriaIterator<dim,TriaObjectAccessor<2,dim> >
     (
       tria,
       present_level,
@@ -528,7 +528,7 @@ TriaObjectAccessor<3, dim>::quad (const unsigned int i) const {
 template <int dim>
 inline
 unsigned int
-TriaObjectAccessor<3, dim>::line_index (unsigned int i) const {
+TriaObjectAccessor<3,dim>::line_index (unsigned int i) const {
   Assert (i<12, ExcIndexRange(i,0,12));
 
   if (i<4)
@@ -557,7 +557,7 @@ TriaObjectAccessor<3, dim>::line_index (unsigned int i) const {
 template <int dim>
 inline
 unsigned int
-TriaObjectAccessor<3, dim>::quad_index (unsigned int i) const {
+TriaObjectAccessor<3,dim>::quad_index (unsigned int i) const {
   Assert (i<6, ExcIndexRange(i,0,6));
 
   return tria->levels[present_level]->hexes.hexes[present_index].quad(i);
@@ -567,11 +567,11 @@ TriaObjectAccessor<3, dim>::quad_index (unsigned int i) const {
 
 template <int dim>
 inline
-TriaIterator<dim,TriaObjectAccessor<3, dim> >
-TriaObjectAccessor<3, dim>::child (const unsigned int i) const {
+TriaIterator<dim,TriaObjectAccessor<3,dim> >
+TriaObjectAccessor<3,dim>::child (const unsigned int i) const {
   Assert (i<6, ExcIndexRange(i,0,6));
   
-  TriaIterator<dim,TriaObjectAccessor<3, dim> > q (tria, present_level+1, child_index (i));
+  TriaIterator<dim,TriaObjectAccessor<3,dim> > q (tria, present_level+1, child_index (i));
   
 #ifdef DEBUG
   if (q.state() != past_the_end)
@@ -584,7 +584,7 @@ TriaObjectAccessor<3, dim>::child (const unsigned int i) const {
 
 template <int dim>
 inline
-int TriaObjectAccessor<3, dim>::child_index (unsigned int i) const {
+int TriaObjectAccessor<3,dim>::child_index (unsigned int i) const {
   Assert (i<8, ExcIndexRange(i,0,8));
   return tria->levels[present_level]->hexes.children[present_index]+i;
 };
@@ -592,7 +592,7 @@ int TriaObjectAccessor<3, dim>::child_index (unsigned int i) const {
 
 
 template <int dim>
-bool TriaObjectAccessor<3, dim>::has_children () const {
+bool TriaObjectAccessor<3,dim>::has_children () const {
   Assert (state() == valid, ExcDereferenceInvalidObject());
   return (tria->levels[present_level]->hexes.children[present_index] != -1);
 };
@@ -602,7 +602,7 @@ bool TriaObjectAccessor<3, dim>::has_children () const {
 template <int dim>
 inline
 unsigned int
-TriaObjectAccessor<3, dim>::max_refinement_depth () const
+TriaObjectAccessor<3,dim>::max_refinement_depth () const
 {
   if (!has_children())
     return 0;
@@ -626,7 +626,7 @@ TriaObjectAccessor<3, dim>::max_refinement_depth () const
 template <int dim>
 inline
 void
-TriaObjectAccessor<3, dim>::operator ++ () {
+TriaObjectAccessor<3,dim>::operator ++ () {
   ++present_index;
                                   // is index still in the range of
                                   // the vector?
@@ -652,7 +652,7 @@ TriaObjectAccessor<3, dim>::operator ++ () {
 template <int dim>
 inline
 void
-TriaObjectAccessor<3, dim>::operator -- () {
+TriaObjectAccessor<3,dim>::operator -- () {
   --present_index;
                                   // is index still in the range of
                                   // the vector?
@@ -682,7 +682,7 @@ TriaObjectAccessor<3, dim>::operator -- () {
 template <int celldim, int dim>
 inline
 bool
-TriaObjectAccessor<celldim, dim>::used () const
+TriaObjectAccessor<celldim,dim>::used () const
 {
   Assert (state() == valid, ExcDereferenceInvalidObject());
   return tria->levels[present_level]->hexes.used[present_index];
@@ -693,7 +693,7 @@ TriaObjectAccessor<celldim, dim>::used () const
 template<int celldim, int dim>
 inline
 bool
-TriaObjectAccessor<celldim, dim>::user_flag_set () const
+TriaObjectAccessor<celldim,dim>::user_flag_set () const
 {
   Assert (used(), ExcCellNotUsed());
   return tria->levels[present_level]->hexes.user_flags[present_index];
@@ -704,7 +704,7 @@ TriaObjectAccessor<celldim, dim>::user_flag_set () const
 template<int celldim, int dim>
 inline
 void
-TriaObjectAccessor<celldim, dim>::set_user_flag () const
+TriaObjectAccessor<celldim,dim>::set_user_flag () const
 {
   Assert (used(), ExcCellNotUsed());
   tria->levels[present_level]->hexes.user_flags[present_index] = true;
@@ -714,7 +714,7 @@ TriaObjectAccessor<celldim, dim>::set_user_flag () const
 
 template<int celldim, int dim>
 inline
-void TriaObjectAccessor<celldim, dim>::clear_user_flag () const
+void TriaObjectAccessor<celldim,dim>::clear_user_flag () const
 {
   Assert (used(), ExcCellNotUsed());
   tria->levels[present_level]->hexes.user_flags[present_index] = false;
@@ -725,8 +725,8 @@ void TriaObjectAccessor<celldim, dim>::clear_user_flag () const
 
 template<int celldim, int dim>
 inline
-TriaIterator<dim,TriaObjectAccessor<1, dim> >
-TriaObjectAccessor<celldim, dim>::line (const unsigned int i) const
+TriaIterator<dim,TriaObjectAccessor<1,dim> >
+TriaObjectAccessor<celldim,dim>::line (const unsigned int i) const
 {
   Assert (used(), ExcCellNotUsed());
   Assert (i < GeometryInfo<celldim>::lines_per_cell,
@@ -736,7 +736,7 @@ TriaObjectAccessor<celldim, dim>::line (const unsigned int i) const
     {
       case 2:
            return
-             TriaIterator<dim,TriaObjectAccessor<1, dim> >
+             TriaIterator<dim,TriaObjectAccessor<1,dim> >
              (
                tria,
                present_level,
@@ -765,21 +765,21 @@ TriaObjectAccessor<celldim, dim>::line (const unsigned int i) const
            Assert(false, ExcNotImplemented());
     }
   
-  return TriaIterator<dim,TriaObjectAccessor<1, dim> >(tria, -1, -1, 0);
+  return TriaIterator<dim,TriaObjectAccessor<1,dim> >(tria, -1, -1, 0);
 };
 
 
 
 template<int celldim, int dim>
 inline
-TriaIterator<dim,TriaObjectAccessor<2, dim> >
-TriaObjectAccessor<celldim, dim>::quad (const unsigned int i) const
+TriaIterator<dim,TriaObjectAccessor<2,dim> >
+TriaObjectAccessor<celldim,dim>::quad (const unsigned int i) const
 {
   Assert (used(), ExcCellNotUsed());
   Assert (i < GeometryInfo<celldim>::quads_per_cell,
          ExcIndexRange(i,0,GeometryInfo<celldim>::quads_per_cell));
   return
-    TriaIterator<dim,TriaObjectAccessor<2, dim> >
+    TriaIterator<dim,TriaObjectAccessor<2,dim> >
     (
       tria,
       present_level,
@@ -792,7 +792,7 @@ TriaObjectAccessor<celldim, dim>::quad (const unsigned int i) const
 template<int celldim, int dim>
 inline
 unsigned int
-TriaObjectAccessor<celldim, dim>::line_index (unsigned int i) const
+TriaObjectAccessor<celldim,dim>::line_index (unsigned int i) const
 {
   Assert (i < GeometryInfo<celldim>::lines_per_cell,
          ExcIndexRange(i,0,GeometryInfo<celldim>::lines_per_cell));
@@ -832,7 +832,7 @@ TriaObjectAccessor<celldim, dim>::line_index (unsigned int i) const
 template<int celldim, int dim>
 inline
 unsigned int
-TriaObjectAccessor<celldim, dim>::quad_index (unsigned int i) const
+TriaObjectAccessor<celldim,dim>::quad_index (unsigned int i) const
 {
   Assert (i < GeometryInfo<celldim>::quads_per_cell,
          ExcIndexRange(i,0,GeometryInfo<celldim>::quads_per_cell));
@@ -844,13 +844,13 @@ TriaObjectAccessor<celldim, dim>::quad_index (unsigned int i) const
 
 template<int celldim, int dim>
 inline
-TriaIterator<dim,TriaObjectAccessor<celldim, dim> >
-TriaObjectAccessor<celldim, dim>::child (const unsigned int i) const
+TriaIterator<dim,TriaObjectAccessor<celldim,dim> >
+TriaObjectAccessor<celldim,dim>::child (const unsigned int i) const
 {
   Assert (i < GeometryInfo<celldim>::children_per_cell,
          ExcIndexRange(i,0,GeometryInfo<celldim>::children_per_cell));
   
-  TriaIterator<dim,TriaObjectAccessor<celldim, dim> > q (tria, present_level+1, child_index (i));
+  TriaIterator<dim,TriaObjectAccessor<celldim,dim> > q (tria, present_level+1, child_index (i));
   
 #ifdef DEBUG
   if (q.state() != past_the_end)
@@ -863,7 +863,7 @@ TriaObjectAccessor<celldim, dim>::child (const unsigned int i) const
 
 template<int celldim, int dim>
 inline
-int TriaObjectAccessor<celldim, dim>::child_index (unsigned int i) const
+int TriaObjectAccessor<celldim,dim>::child_index (unsigned int i) const
 {
   Assert (i < GeometryInfo<celldim>::children_per_cell,
          ExcIndexRange(i,0,GeometryInfo<celldim>::children_per_cell));
@@ -873,7 +873,7 @@ int TriaObjectAccessor<celldim, dim>::child_index (unsigned int i) const
 
 
 template<int celldim, int dim>
-bool TriaObjectAccessor<celldim, dim>::has_children () const
+bool TriaObjectAccessor<celldim,dim>::has_children () const
 {
   Assert (state() == valid, ExcDereferenceInvalidObject());
   return (tria->levels[present_level]->hexes.children[present_index] != -1);
@@ -884,7 +884,7 @@ bool TriaObjectAccessor<celldim, dim>::has_children () const
 template<int celldim, int dim>
 inline
 unsigned int
-TriaObjectAccessor<celldim, dim>::max_refinement_depth () const
+TriaObjectAccessor<celldim,dim>::max_refinement_depth () const
 {
   if (!has_children())
     return 0;
@@ -908,7 +908,7 @@ TriaObjectAccessor<celldim, dim>::max_refinement_depth () const
 template<int celldim, int dim>
 inline
 void
-TriaObjectAccessor<celldim, dim>::operator ++ () {
+TriaObjectAccessor<celldim,dim>::operator ++ () {
   ++present_index;
                                   // is index still in the range of
                                   // the vector?
@@ -934,7 +934,7 @@ TriaObjectAccessor<celldim, dim>::operator ++ () {
 template<int celldim, int dim>
 inline
 void
-TriaObjectAccessor<celldim, dim>::operator -- ()
+TriaObjectAccessor<celldim,dim>::operator -- ()
 {
   --present_index;
                                   // is index still in the range of
index baea5b2fecf700a8a8d09a35cb0926215cdd91e4..3e92c8a903098b5376e8aba73018a5cd939bf7af 100644 (file)
@@ -625,7 +625,7 @@ MGMatrix<MATRIX>::operator[](unsigned int i) const
 
 
 template<class SOLVER, class MATRIX, class PRECOND>
-MGCoarseGridLACIteration<SOLVER, MATRIX, PRECOND>::
+MGCoarseGridLACIteration<SOLVER,MATRIX,PRECOND>::
 MGCoarseGridLACIteration(SOLVER        &s,
                         const MATRIX  &m,
                         const PRECOND &p)
@@ -639,7 +639,7 @@ MGCoarseGridLACIteration(SOLVER        &s,
 
 template<class SOLVER, class MATRIX, class PRECOND>
 void
-MGCoarseGridLACIteration<SOLVER, MATRIX, PRECOND>::
+MGCoarseGridLACIteration<SOLVER,MATRIX,PRECOND>::
 operator() (const unsigned int    /* level */,
            Vector<double>       &dst,
            const Vector<double> &src) const
index ae184b777ebd44d608e253e4619efbe712b4340f..3c5754e7bcd0da92fd4f2e33c87a4fe692a974ea 100644 (file)
@@ -33,7 +33,7 @@ class BlockVector
                                      * that of V, all components are copied
                                      * from V
                                      */
-    BlockVector (const BlockVector<n_blocks, Number>& V);
+    BlockVector (const BlockVector<n_blocks,Number>& V);
 
 
 // note: I disabled this function for the time being, since egcs1.1.2
@@ -86,25 +86,25 @@ class BlockVector
                                     /**
                                      * $U(0-N) = s$: fill all components.
                                      */
-    BlockVector<n_blocks, Number> & operator= (const Number s);
+    BlockVector<n_blocks,Number> & operator= (const Number s);
     
                                     /**
                                      *  $U = V$: copy all components.
                                      */
-    BlockVector<n_blocks, Number> &
-    operator= (const BlockVector<n_blocks, Number>& V);
+    BlockVector<n_blocks,Number> &
+    operator= (const BlockVector<n_blocks,Number>& V);
 
                                     /**
                                      * $U = V$ for different types.
                                      */
     template<typename Number2>
-    BlockVector<n_blocks, Number> &
+    BlockVector<n_blocks,Number> &
     operator= (const BlockVector<n_blocks, Number2>& V);
     
                                     /**
                                      * $U = U * V$: scalar product.
                                      */
-    Number operator* (const BlockVector<n_blocks, Number>& V) const;
+    Number operator* (const BlockVector<n_blocks,Number>& V) const;
 
                                     /**
                                      * Return square of the $l_2$-norm.
@@ -167,7 +167,7 @@ class BlockVector
                                      * this function is the same as calling
                                      * #reinit (V.size(), fast)#.
                                      */
-    void reinit (const BlockVector<n_blocks, Number> &V,
+    void reinit (const BlockVector<n_blocks,Number> &V,
                 const bool            fast=false);
     
                                     /**
@@ -241,15 +241,15 @@ class BlockVector
                                      * Addition operator.
                                      * Fast equivalent to #U.add(1, V)#.
                                      */
-    BlockVector<n_blocks, Number> &
-    operator += (const BlockVector<n_blocks, Number> &V);
+    BlockVector<n_blocks,Number> &
+    operator += (const BlockVector<n_blocks,Number> &V);
 
                                     /**
                                      * Subtraction operator.
                                      * Fast equivalent to #U.add(-1, V)#.
                                      */
-    BlockVector<n_blocks, Number> &
-    operator -= (const BlockVector<n_blocks, Number> &V);
+    BlockVector<n_blocks,Number> &
+    operator -= (const BlockVector<n_blocks,Number> &V);
 
                                     /**
                                      * $U(0-DIM)+=s$.
@@ -263,49 +263,49 @@ class BlockVector
                                      * Simple vector addition, equal to the
                                      * #operator +=#.
                                      */
-    void add (const BlockVector<n_blocks, Number>& V);
+    void add (const BlockVector<n_blocks,Number>& V);
     
                                     /**
                                      * U+=a*V.
                                      * Simple addition of a scaled vector.
                                      */
-    void add (const Number a, const BlockVector<n_blocks, Number>& V);
+    void add (const Number a, const BlockVector<n_blocks,Number>& V);
     
                                     /**
                                      * U+=a*V+b*W.
                                      * Multiple addition of scaled vectors.
                                      */
-    void add (const Number a, const BlockVector<n_blocks, Number>& V,
-             const Number b, const BlockVector<n_blocks, Number>& W);
+    void add (const Number a, const BlockVector<n_blocks,Number>& V,
+             const Number b, const BlockVector<n_blocks,Number>& W);
     
                                     /**
                                      * U=s*U+V.
                                      * Scaling and simple vector addition.
                                      */
-    void sadd (const Number s, const BlockVector<n_blocks, Number>& V);
+    void sadd (const Number s, const BlockVector<n_blocks,Number>& V);
     
                                     /**
                                      * U=s*U+a*V.
                                      * Scaling and simple addition.
                                      */
-    void sadd (const Number s, const Number a, const BlockVector<n_blocks, Number>& V);
+    void sadd (const Number s, const Number a, const BlockVector<n_blocks,Number>& V);
     
                                     /**
                                      * U=s*U+a*V+b*W.
                                      * Scaling and multiple addition.
                                      */
     void sadd (const Number s, const Number a,
-              const BlockVector<n_blocks, Number>& V,
-              const Number b, const BlockVector<n_blocks, Number>& W);
+              const BlockVector<n_blocks,Number>& V,
+              const Number b, const BlockVector<n_blocks,Number>& W);
     
                                     /**
                                      * U=s*U+a*V+b*W+c*X.
                                      * Scaling and multiple addition.
                                      */
     void sadd (const Number s, const Number a,
-              const BlockVector<n_blocks, Number>& V,
-              const Number b, const BlockVector<n_blocks, Number>& W, 
-              const Number c, const BlockVector<n_blocks, Number>& X);
+              const BlockVector<n_blocks,Number>& V,
+              const Number b, const BlockVector<n_blocks,Number>& W, 
+              const Number c, const BlockVector<n_blocks,Number>& X);
     
                                     /**
                                      * Scale each element of the vector by the
@@ -319,14 +319,14 @@ class BlockVector
                                     /**
                                      *  U=a*V. Replacing.
                                      */
-    void equ (const Number a, const BlockVector<n_blocks, Number>& V);
+    void equ (const Number a, const BlockVector<n_blocks,Number>& V);
     
                                     /**
                                      * U=a*V+b*W.
                                      * Replacing by sum.
                                      */
-    void equ (const Number a, const BlockVector<n_blocks, Number>& V,
-             const Number b, const BlockVector<n_blocks, Number>& W);
+    void equ (const Number a, const BlockVector<n_blocks,Number>& V,
+             const Number b, const BlockVector<n_blocks,Number>& W);
 
                                     //@}
     
@@ -438,7 +438,7 @@ class BlockVector
 
 template <int n_blocks, typename Number>
 inline
-unsigned int BlockVector<n_blocks, Number>::size () const
+unsigned int BlockVector<n_blocks,Number>::size () const
 {
   return start[n_blocks];
 }
@@ -447,7 +447,7 @@ unsigned int BlockVector<n_blocks, Number>::size () const
 
 //  template <int n_blocks, typename Number>
 //  inline
-//  BlockVector<n_blocks, Number>::iterator BlockVector<n_blocks, Number>::begin () {
+//  BlockVector<n_blocks,Number>::iterator BlockVector<n_blocks,Number>::begin () {
 //    return &val[0];
 //  };
 
@@ -455,7 +455,7 @@ unsigned int BlockVector<n_blocks, Number>::size () const
 
 //  template <int n_blocks, typename Number>
 //  inline
-//  BlockVector<n_blocks, Number>::const_iterator BlockVector<n_blocks, Number>::begin () const {
+//  BlockVector<n_blocks,Number>::const_iterator BlockVector<n_blocks,Number>::begin () const {
 //    return &val[0];
 //  };
 
@@ -463,7 +463,7 @@ unsigned int BlockVector<n_blocks, Number>::size () const
 
 //  template <int n_blocks, typename Number>
 //  inline
-//  BlockVector<n_blocks, Number>::iterator BlockVector<n_blocks, Number>::end () {
+//  BlockVector<n_blocks,Number>::iterator BlockVector<n_blocks,Number>::end () {
 //    return &val[dim];
 //  };
 
@@ -471,7 +471,7 @@ unsigned int BlockVector<n_blocks, Number>::size () const
 
 //  template <typename Number>
 //  inline
-//  BlockVector<n_blocks, Number>::const_iterator BlockVector<n_blocks, Number>::end () const {
+//  BlockVector<n_blocks,Number>::const_iterator BlockVector<n_blocks,Number>::end () const {
 //    return &val[dim];
 //  };
 
@@ -479,7 +479,7 @@ unsigned int BlockVector<n_blocks, Number>::size () const
 
 template <int n_blocks, typename Number>
 inline
-Number BlockVector<n_blocks, Number>::operator() (const unsigned int i) const
+Number BlockVector<n_blocks,Number>::operator() (const unsigned int i) const
 {
   //  Assert (i<start[n_blocks], ExcInvalidIndex(i,dim));
 
@@ -493,7 +493,7 @@ Number BlockVector<n_blocks, Number>::operator() (const unsigned int i) const
 
 template <int n_blocks, typename Number>
 inline
-Number& BlockVector<n_blocks, Number>::operator() (const unsigned int i)
+Number& BlockVector<n_blocks,Number>::operator() (const unsigned int i)
 {
   //  Assert (i<dim, ExcInvalidIndex(i,dim));
 
@@ -506,7 +506,7 @@ Number& BlockVector<n_blocks, Number>::operator() (const unsigned int i)
 template <int n_blocks, typename Number>
 inline
 Vector<Number>&
-BlockVector<n_blocks, Number>::block(unsigned int i)
+BlockVector<n_blocks,Number>::block(unsigned int i)
 {
   Assert(i<n_blocks, ExcIndexRange(i,0,n_blocks));
 
@@ -516,7 +516,7 @@ BlockVector<n_blocks, Number>::block(unsigned int i)
 template <int n_blocks, typename Number>
 inline
 const Vector<Number>&
-BlockVector<n_blocks, Number>::block(unsigned int i) const
+BlockVector<n_blocks,Number>::block(unsigned int i) const
 {
   Assert(i<n_blocks, ExcIndexRange(i,0,n_blocks));
 
index 1e4e8f1a393bc9de8bb898a863c5227b5a0c440a..eaa93a2c1321c548f9965fba94af74644c643ac4 100644 (file)
@@ -19,7 +19,7 @@ static inline Number sqr (const Number x)
 
 
 template <int n_blocks, typename Number>
-BlockVector<n_blocks, Number>::BlockVector ()
+BlockVector<n_blocks,Number>::BlockVector ()
 {
   for (unsigned int i=0;i<=n_blocks;++i)
     start[n_blocks] = 0;
@@ -28,7 +28,7 @@ BlockVector<n_blocks, Number>::BlockVector ()
 
 
 template <int n_blocks, typename Number>
-BlockVector<n_blocks, Number>::BlockVector (const vector<unsigned int>& n)
+BlockVector<n_blocks,Number>::BlockVector (const vector<unsigned int>& n)
 {
   Assert(n.size()==n_blocks, ExcDimensionMismatch(n.size(), n_blocks));
 
@@ -38,7 +38,7 @@ BlockVector<n_blocks, Number>::BlockVector (const vector<unsigned int>& n)
 
 
 template <int n_blocks, typename Number>
-BlockVector<n_blocks, Number>::BlockVector (const BlockVector<n_blocks, Number>& v)
+BlockVector<n_blocks,Number>::BlockVector (const BlockVector<n_blocks,Number>& v)
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
@@ -54,7 +54,7 @@ BlockVector<n_blocks, Number>::BlockVector (const BlockVector<n_blocks, Number>&
 //
 // template <int n_blocks, typename Number>
 // template <typename OtherNumber>
-// BlockVector<n_blocks, Number>::Vector (const BlockVector<n_blocks, OtherNumber>& v) :
+// BlockVector<n_blocks,Number>::Vector (const BlockVector<n_blocks, OtherNumber>& v) :
 //             dim(v.size()),
 //             maxdim(v.size()),
 //             val(0)
@@ -70,7 +70,7 @@ BlockVector<n_blocks, Number>::BlockVector (const BlockVector<n_blocks, Number>&
 
 
 template <int n_blocks, typename Number>
-void BlockVector<n_blocks, Number>::reinit (const vector<unsigned int>& n,
+void BlockVector<n_blocks,Number>::reinit (const vector<unsigned int>& n,
                                            const bool fast)
 {
   Assert (n.size() == n_blocks,
@@ -89,7 +89,7 @@ void BlockVector<n_blocks, Number>::reinit (const vector<unsigned int>& n,
 
 
 template <int n_blocks, typename Number>
-void BlockVector<n_blocks, Number>::reinit (const BlockVector<n_blocks, Number>& v,
+void BlockVector<n_blocks,Number>::reinit (const BlockVector<n_blocks,Number>& v,
                                            const bool fast)
 {
   for (unsigned int i=0;i<n_blocks;++i)
@@ -104,14 +104,14 @@ void BlockVector<n_blocks, Number>::reinit (const BlockVector<n_blocks, Number>&
 
 
 template <int n_blocks, typename Number>
-BlockVector<n_blocks, Number>::~BlockVector ()
+BlockVector<n_blocks,Number>::~BlockVector ()
 {
 }
 
 
 
 template <int n_blocks, typename Number>
-void BlockVector<n_blocks, Number>::clear ()
+void BlockVector<n_blocks,Number>::clear ()
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
@@ -122,7 +122,7 @@ void BlockVector<n_blocks, Number>::clear ()
 
 
 template <int n_blocks, typename Number>
-bool BlockVector<n_blocks, Number>::all_zero () const
+bool BlockVector<n_blocks,Number>::all_zero () const
 {
   bool result = true;
   for (unsigned int i=0;i<n_blocks;++i)
@@ -135,7 +135,7 @@ bool BlockVector<n_blocks, Number>::all_zero () const
 
 
 template <int n_blocks, typename Number>
-Number BlockVector<n_blocks, Number>::operator * (const BlockVector<n_blocks, Number>& v) const
+Number BlockVector<n_blocks,Number>::operator * (const BlockVector<n_blocks,Number>& v) const
 {
   Number sum = 0.;
   for (unsigned int i=0;i<n_blocks;++i)
@@ -148,7 +148,7 @@ Number BlockVector<n_blocks, Number>::operator * (const BlockVector<n_blocks, Nu
 
 
 template <int n_blocks, typename Number>
-Number BlockVector<n_blocks, Number>::norm_sqr () const
+Number BlockVector<n_blocks,Number>::norm_sqr () const
 {
   Number sum = 0.;
   for (unsigned int i=0;i<n_blocks;++i)
@@ -161,7 +161,7 @@ Number BlockVector<n_blocks, Number>::norm_sqr () const
 
 
 template <int n_blocks, typename Number>
-Number BlockVector<n_blocks, Number>::mean_value () const
+Number BlockVector<n_blocks,Number>::mean_value () const
 {
   Number sum = 0.;
   for (unsigned int i=0;i<n_blocks;++i)
@@ -174,7 +174,7 @@ Number BlockVector<n_blocks, Number>::mean_value () const
 
 
 template <int n_blocks, typename Number>
-Number BlockVector<n_blocks, Number>::l1_norm () const
+Number BlockVector<n_blocks,Number>::l1_norm () const
 {
   Number sum = 0.;
   for (unsigned int i=0;i<n_blocks;++i)
@@ -187,7 +187,7 @@ Number BlockVector<n_blocks, Number>::l1_norm () const
 
 
 template <int n_blocks, typename Number>
-Number BlockVector<n_blocks, Number>::l2_norm () const
+Number BlockVector<n_blocks,Number>::l2_norm () const
 {
   return sqrt(norm_sqr());
 }
@@ -195,7 +195,7 @@ Number BlockVector<n_blocks, Number>::l2_norm () const
 
 
 template <int n_blocks, typename Number>
-Number BlockVector<n_blocks, Number>::linfty_norm () const
+Number BlockVector<n_blocks,Number>::linfty_norm () const
 {
   Number sum = 0.;
   for (unsigned int i=0;i<n_blocks;++i)
@@ -212,7 +212,7 @@ Number BlockVector<n_blocks, Number>::linfty_norm () const
 
 
 template <int n_blocks, typename Number>
-BlockVector<n_blocks, Number>& BlockVector<n_blocks, Number>::operator += (const BlockVector<n_blocks, Number>& v)
+BlockVector<n_blocks,Number>& BlockVector<n_blocks,Number>::operator += (const BlockVector<n_blocks,Number>& v)
 {
   add (v);
   return *this;
@@ -221,7 +221,7 @@ BlockVector<n_blocks, Number>& BlockVector<n_blocks, Number>::operator += (const
 
 
 template <int n_blocks, typename Number>
-BlockVector<n_blocks, Number>& BlockVector<n_blocks, Number>::operator -= (const BlockVector<n_blocks, Number>& v)
+BlockVector<n_blocks,Number>& BlockVector<n_blocks,Number>::operator -= (const BlockVector<n_blocks,Number>& v)
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
@@ -233,7 +233,7 @@ BlockVector<n_blocks, Number>& BlockVector<n_blocks, Number>::operator -= (const
 
 
 template <int n_blocks, typename Number>
-void BlockVector<n_blocks, Number>::add (const Number v)
+void BlockVector<n_blocks,Number>::add (const Number v)
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
@@ -244,7 +244,7 @@ void BlockVector<n_blocks, Number>::add (const Number v)
 
 
 template <int n_blocks, typename Number>
-void BlockVector<n_blocks, Number>::add (const BlockVector<n_blocks, Number>& v)
+void BlockVector<n_blocks,Number>::add (const BlockVector<n_blocks,Number>& v)
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
@@ -255,8 +255,8 @@ void BlockVector<n_blocks, Number>::add (const BlockVector<n_blocks, Number>& v)
 
 
 template <int n_blocks, typename Number>
-void BlockVector<n_blocks, Number>::add (const Number a,
-                                        const BlockVector<n_blocks, Number>& v)
+void BlockVector<n_blocks,Number>::add (const Number a,
+                                        const BlockVector<n_blocks,Number>& v)
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
@@ -267,10 +267,10 @@ void BlockVector<n_blocks, Number>::add (const Number a,
 
 
 template <int n_blocks, typename Number>
-void BlockVector<n_blocks, Number>::add (const Number a,
-                                        const BlockVector<n_blocks, Number>& v,
+void BlockVector<n_blocks,Number>::add (const Number a,
+                                        const BlockVector<n_blocks,Number>& v,
                                         const Number b,
-                                        const BlockVector<n_blocks, Number>& w)
+                                        const BlockVector<n_blocks,Number>& w)
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
@@ -281,8 +281,8 @@ void BlockVector<n_blocks, Number>::add (const Number a,
 
 
 template <int n_blocks, typename Number>
-void BlockVector<n_blocks, Number>::sadd (const Number x,
-                                         const BlockVector<n_blocks, Number>& v)
+void BlockVector<n_blocks,Number>::sadd (const Number x,
+                                         const BlockVector<n_blocks,Number>& v)
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
@@ -293,8 +293,8 @@ void BlockVector<n_blocks, Number>::sadd (const Number x,
 
 
 template <int n_blocks, typename Number>
-void BlockVector<n_blocks, Number>::sadd (const Number x, const Number a,
-                                         const BlockVector<n_blocks, Number>& v)
+void BlockVector<n_blocks,Number>::sadd (const Number x, const Number a,
+                                         const BlockVector<n_blocks,Number>& v)
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
@@ -305,10 +305,10 @@ void BlockVector<n_blocks, Number>::sadd (const Number x, const Number a,
 
 
 template <int n_blocks, typename Number>
-void BlockVector<n_blocks, Number>::sadd (const Number x, const Number a,
-                                         const BlockVector<n_blocks, Number>& v,
+void BlockVector<n_blocks,Number>::sadd (const Number x, const Number a,
+                                         const BlockVector<n_blocks,Number>& v,
                                          const Number b,
-                                         const BlockVector<n_blocks, Number>& w)
+                                         const BlockVector<n_blocks,Number>& w)
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
@@ -319,12 +319,12 @@ void BlockVector<n_blocks, Number>::sadd (const Number x, const Number a,
 
 
 template <int n_blocks, typename Number>
-void BlockVector<n_blocks, Number>::sadd (const Number x, const Number a,
-                                         const BlockVector<n_blocks, Number>& v,
+void BlockVector<n_blocks,Number>::sadd (const Number x, const Number a,
+                                         const BlockVector<n_blocks,Number>& v,
                                          const Number b,
-                                         const BlockVector<n_blocks, Number>& w,
+                                         const BlockVector<n_blocks,Number>& w,
                                          const Number c,
-                                         const BlockVector<n_blocks, Number>& y)
+                                         const BlockVector<n_blocks,Number>& y)
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
@@ -336,7 +336,7 @@ void BlockVector<n_blocks, Number>::sadd (const Number x, const Number a,
 
 
 template <int n_blocks, typename Number>
-void BlockVector<n_blocks, Number>::scale (const Number factor)
+void BlockVector<n_blocks,Number>::scale (const Number factor)
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
@@ -347,10 +347,10 @@ void BlockVector<n_blocks, Number>::scale (const Number factor)
 
 
 template <int n_blocks, typename Number>
-void BlockVector<n_blocks, Number>::equ (const Number a,
-                                        const BlockVector<n_blocks, Number>& v,
+void BlockVector<n_blocks,Number>::equ (const Number a,
+                                        const BlockVector<n_blocks,Number>& v,
                                         const Number b,
-                                        const BlockVector<n_blocks, Number>& w)
+                                        const BlockVector<n_blocks,Number>& w)
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
@@ -361,8 +361,8 @@ void BlockVector<n_blocks, Number>::equ (const Number a,
 
 
 template <int n_blocks, typename Number>
-void BlockVector<n_blocks, Number>::equ (const Number a,
-                                        const BlockVector<n_blocks, Number>& v)
+void BlockVector<n_blocks,Number>::equ (const Number a,
+                                        const BlockVector<n_blocks,Number>& v)
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
@@ -372,7 +372,7 @@ void BlockVector<n_blocks, Number>::equ (const Number a,
 
 
 template <int n_blocks, typename Number>
-BlockVector<n_blocks, Number>& BlockVector<n_blocks, Number>::operator = (const Number s)
+BlockVector<n_blocks,Number>& BlockVector<n_blocks,Number>::operator = (const Number s)
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
@@ -384,8 +384,8 @@ BlockVector<n_blocks, Number>& BlockVector<n_blocks, Number>::operator = (const
 
 
 template <int n_blocks, typename Number>
-BlockVector<n_blocks, Number>&
-BlockVector<n_blocks, Number>::operator = (const BlockVector<n_blocks, Number>& v)
+BlockVector<n_blocks,Number>&
+BlockVector<n_blocks,Number>::operator = (const BlockVector<n_blocks,Number>& v)
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
@@ -398,8 +398,8 @@ BlockVector<n_blocks, Number>::operator = (const BlockVector<n_blocks, Number>&
 
 template <int n_blocks, typename Number>
 template<typename Number2>
-BlockVector<n_blocks, Number>&
-BlockVector<n_blocks, Number>::operator = (const BlockVector<n_blocks, Number2>& v)
+BlockVector<n_blocks,Number>&
+BlockVector<n_blocks,Number>::operator = (const BlockVector<n_blocks, Number2>& v)
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
@@ -410,7 +410,7 @@ BlockVector<n_blocks, Number>::operator = (const BlockVector<n_blocks, Number2>&
 
 
 template <int n_blocks, typename Number>
-void BlockVector<n_blocks, Number>::print (ostream &out,
+void BlockVector<n_blocks,Number>::print (ostream &out,
                            unsigned int precision,
                            bool scientific,
                            bool across) const
@@ -428,7 +428,7 @@ void BlockVector<n_blocks, Number>::print (ostream &out,
 
 
 template <int n_blocks, typename Number>
-void BlockVector<n_blocks, Number>::block_write (ostream &out) const
+void BlockVector<n_blocks,Number>::block_write (ostream &out) const
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
@@ -438,7 +438,7 @@ void BlockVector<n_blocks, Number>::block_write (ostream &out) const
 
 
 template <int n_blocks, typename Number>
-void BlockVector<n_blocks, Number>::block_read (istream &in)
+void BlockVector<n_blocks,Number>::block_read (istream &in)
 {
   for (unsigned int i=0;i<n_blocks;++i)
     {
index 92fec377f0511fe58d11c74921ad2d3c003094b6..6012ba2b64717aa53334d27d63f40172bf2e2650 100644 (file)
@@ -625,7 +625,7 @@ MGMatrix<MATRIX>::operator[](unsigned int i) const
 
 
 template<class SOLVER, class MATRIX, class PRECOND>
-MGCoarseGridLACIteration<SOLVER, MATRIX, PRECOND>::
+MGCoarseGridLACIteration<SOLVER,MATRIX,PRECOND>::
 MGCoarseGridLACIteration(SOLVER        &s,
                         const MATRIX  &m,
                         const PRECOND &p)
@@ -639,7 +639,7 @@ MGCoarseGridLACIteration(SOLVER        &s,
 
 template<class SOLVER, class MATRIX, class PRECOND>
 void
-MGCoarseGridLACIteration<SOLVER, MATRIX, PRECOND>::
+MGCoarseGridLACIteration<SOLVER,MATRIX,PRECOND>::
 operator() (const unsigned int    /* level */,
            Vector<double>       &dst,
            const Vector<double> &src) const
index 6897d521cf37ea4924d0f0852600e4c7f3173161..4f9410e246d67639d3525eb5698af52879754235 100644 (file)
@@ -1,9 +1,11 @@
-/*-------------------- precondition.h --------------------*/
-//$Id$
-// Guido Kanschat, University of Heidelberg, 1999
+/*----------------------------   precondition.h     ---------------------------*/
+/*      $Id$                 */
+#ifndef __precondition_H
+#define __precondition_H
+/*----------------------------   precondition.h     ---------------------------*/
+
+
 
-#ifndef __lac_precondition_h
-#define __lac_precondition_h
 
 /**
  * No preconditioning.
@@ -243,7 +245,7 @@ PreconditionIdentity::operator() (VECTOR& dst, const VECTOR& src) const
 }
 
 template<class MATRIX, class VECTOR>
-PreconditionUseMatrix<MATRIX, VECTOR>::PreconditionUseMatrix(const MATRIX& M,
+PreconditionUseMatrix<MATRIX,VECTOR>::PreconditionUseMatrix(const MATRIX& M,
                                                             function_ptr method)
                :
                matrix(M), precondition(method)
@@ -252,14 +254,14 @@ PreconditionUseMatrix<MATRIX, VECTOR>::PreconditionUseMatrix(const MATRIX& M,
 
 template<class MATRIX, class VECTOR>
 void
-PreconditionUseMatrix<MATRIX, VECTOR>::operator() (VECTOR& dst,
+PreconditionUseMatrix<MATRIX,VECTOR>::operator() (VECTOR& dst,
                                                   const VECTOR& src) const
 {
   (matrix.*precondition)(dst, src);
 }
 
 template<class MATRIX, class VECTOR>
-PreconditionRelaxation<MATRIX, VECTOR>::PreconditionRelaxation(const MATRIX& M,
+PreconditionRelaxation<MATRIX,VECTOR>::PreconditionRelaxation(const MATRIX& M,
                                                             function_ptr method,
                                                               double omega)
                :
@@ -269,7 +271,7 @@ PreconditionRelaxation<MATRIX, VECTOR>::PreconditionRelaxation(const MATRIX& M,
 
 template<class MATRIX, class VECTOR>
 void
-PreconditionRelaxation<MATRIX, VECTOR>::operator() (VECTOR& dst,
+PreconditionRelaxation<MATRIX,VECTOR>::operator() (VECTOR& dst,
                                                   const VECTOR& src) const
 {
   (matrix.*precondition)(dst, src, omega);
@@ -278,7 +280,7 @@ PreconditionRelaxation<MATRIX, VECTOR>::operator() (VECTOR& dst,
 //////////////////////////////////////////////////////////////////////
 
 template<class SOLVER, class MATRIX, class PRECONDITION>
-PreconditionLACSolver<SOLVER, MATRIX, PRECONDITION>
+PreconditionLACSolver<SOLVER,MATRIX,PRECONDITION>
 ::PreconditionLACSolver(SOLVER& solver,
                        const MATRIX& matrix,
                        const PRECONDITION& precondition)
@@ -289,7 +291,7 @@ PreconditionLACSolver<SOLVER, MATRIX, PRECONDITION>
 template<class SOLVER, class MATRIX, class PRECONDITION>
 template<class VECTOR>
 void
-PreconditionLACSolver<SOLVER, MATRIX, PRECONDITION>::operator() (VECTOR& dst,
+PreconditionLACSolver<SOLVER,MATRIX,PRECONDITION>::operator() (VECTOR& dst,
                                                                 const VECTOR& src) const
 {
   solver.solve(matrix, dst, src, precondition);
@@ -297,4 +299,8 @@ PreconditionLACSolver<SOLVER, MATRIX, PRECONDITION>::operator() (VECTOR& dst,
 
 
 
+
+/*----------------------------   precondition.h     ---------------------------*/
+/* end of #ifndef __precondition_H */
 #endif
+/*----------------------------   precondition.h     ---------------------------*/
index 75d33a9a939803b6bfcaf55c023c9f14ed42ce99..beca0f1f7f3fd899c869e93ec60cc6e8a2f4329b 100644 (file)
@@ -200,7 +200,7 @@ class PreconditionBlock: public Subscriptor
  */
 template<typename number,
          typename inverse_type = number>
-class PreconditionBlockSOR : public PreconditionBlock<number, inverse_type>
+class PreconditionBlockSOR : public PreconditionBlock<number,inverse_type>
 {
   public:
                                     /**
index 4aad18d3a9ef77ff43090bdb7ba502c1c9e3a483..a7c147d60ec4ed11b45f51b6aff0fd871d1ccf04 100644 (file)
 
 
 template <typename number, typename inverse_type>
-PreconditionBlock<number, inverse_type>::PreconditionBlock ():
+PreconditionBlock<number,inverse_type>::PreconditionBlock ():
                blocksize(0),
                A(0)  {};
 
 
 template <typename number, typename inverse_type>
-PreconditionBlock<number, inverse_type>::~PreconditionBlock ()
+PreconditionBlock<number,inverse_type>::~PreconditionBlock ()
 {
   if (inverse.size()!=0)
     inverse.erase(inverse.begin(), inverse.end());
 }
 
 template <typename number, typename inverse_type>
-void PreconditionBlock<number, inverse_type>::clear ()
+void PreconditionBlock<number,inverse_type>::clear ()
 {
   if (inverse.size()!=0)
     inverse.erase(inverse.begin(), inverse.end());
@@ -31,7 +31,7 @@ void PreconditionBlock<number, inverse_type>::clear ()
 }
 
 template <typename number, typename inverse_type>
-void PreconditionBlock<number, inverse_type>::use_matrix(
+void PreconditionBlock<number,inverse_type>::use_matrix(
   const SparseMatrix<number> &M)
 {
   A = &M;
@@ -39,19 +39,19 @@ void PreconditionBlock<number, inverse_type>::use_matrix(
 
 
 template <typename number, typename inverse_type>
-void PreconditionBlock<number, inverse_type>::set_block_size(unsigned int bsize) {
+void PreconditionBlock<number,inverse_type>::set_block_size(unsigned int bsize) {
   blocksize=bsize;
 }
 
 
 template <typename number, typename inverse_type>
-unsigned int PreconditionBlock<number, inverse_type>::block_size() const {
+unsigned int PreconditionBlock<number,inverse_type>::block_size() const {
   return blocksize;
 }
 
 
 template <typename number, typename inverse_type>
-void PreconditionBlock<number, inverse_type>::invert_diagblocks()
+void PreconditionBlock<number,inverse_type>::invert_diagblocks()
 {
   Assert (A!=0, ExcNoMatrixGivenToUse());
   const SparseMatrix<number> &M=*A;
@@ -115,18 +115,18 @@ void PreconditionBlock<number, inverse_type>::invert_diagblocks()
 /*--------------------- PreconditionBlockSOR -----------------------*/
 
 template<typename number, typename inverse_type>
-PreconditionBlockSOR<number, inverse_type>::PreconditionBlockSOR(const number omega):
+PreconditionBlockSOR<number,inverse_type>::PreconditionBlockSOR(const number omega):
                omega(omega)  {}
 
 
 template<typename number, typename inverse_type>
-PreconditionBlockSOR<number, inverse_type>::~PreconditionBlockSOR(){}
+PreconditionBlockSOR<number,inverse_type>::~PreconditionBlockSOR(){}
 
 
 template <typename number, typename inverse_type>
 template <typename number2>
-void PreconditionBlockSOR<number, inverse_type>::operator() (Vector<number2> &dst,
-                                                            const Vector<number2> &src) const
+void PreconditionBlockSOR<number,inverse_type>::operator() (Vector<number2> &dst,
+                                                           const Vector<number2> &src) const
 {
   Assert(A!=0, ExcNoMatrixGivenToUse());
   const SparseMatrix<number> &M=*A;
index 72c08e56b3758847cae71977d798539d6ad31f55..0b321c392c370ed6565138565da065184d8a7e2e 100644 (file)
@@ -144,7 +144,7 @@ class PreconditionSelector
 
 
 template <class Matrix, class Vector>
-PreconditionSelector<Matrix, Vector>
+PreconditionSelector<Matrix,Vector>
 ::PreconditionSelector(string preconditioning,
                       const typename Vector::value_type &omega) :
                preconditioning(preconditioning),
@@ -152,7 +152,7 @@ PreconditionSelector<Matrix, Vector>
 
 
 template <class Matrix, class Vector>
-PreconditionSelector<Matrix, Vector>::~PreconditionSelector()
+PreconditionSelector<Matrix,Vector>::~PreconditionSelector()
 {
                                   // release the matrix A
   A=0;
@@ -160,14 +160,14 @@ PreconditionSelector<Matrix, Vector>::~PreconditionSelector()
 
 
 template <class Matrix, class Vector>
-void PreconditionSelector<Matrix, Vector>::use_matrix(const Matrix &M)
+void PreconditionSelector<Matrix,Vector>::use_matrix(const Matrix &M)
 {
   A=&M;
 }
 
 template <class Matrix, class Vector>
-void PreconditionSelector<Matrix, Vector>::operator() (Vector &dst,
-                                                      const Vector &src) const
+void PreconditionSelector<Matrix,Vector>::operator() (Vector &dst,
+                                                     const Vector &src) const
 {
   if (preconditioning=="none")
     {
@@ -196,7 +196,7 @@ void PreconditionSelector<Matrix, Vector>::operator() (Vector &dst,
 
 
 template <class Matrix, class Vector>
-string PreconditionSelector<Matrix, Vector>::get_precondition_names()
+string PreconditionSelector<Matrix,Vector>::get_precondition_names()
 {
   return "none|jacobi|sor|ssor";
 }

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.