]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
new heads
authorguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 22 Feb 2000 20:05:29 +0000 (20:05 +0000)
committerguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 22 Feb 2000 20:05:29 +0000 (20:05 +0000)
git-svn-id: https://svn.dealii.org/trunk@2474 0785d39b-7218-0410-832d-ea1e28bc413d

39 files changed:
deal.II/lac/include/lac/block_vector.h
deal.II/lac/include/lac/block_vector.templates.h
deal.II/lac/include/lac/full_matrix.h
deal.II/lac/include/lac/full_matrix.templates.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
deal.II/lac/include/lac/solver.h
deal.II/lac/include/lac/solver_bicgstab.h
deal.II/lac/include/lac/solver_cg.h
deal.II/lac/include/lac/solver_control.h
deal.II/lac/include/lac/solver_gmres.h
deal.II/lac/include/lac/solver_qmrs.h
deal.II/lac/include/lac/solver_richardson.h
deal.II/lac/include/lac/solver_selector.h
deal.II/lac/include/lac/sparse_ilu.h
deal.II/lac/include/lac/sparse_ilu.templates.h
deal.II/lac/include/lac/sparse_matrix.2.templates
deal.II/lac/include/lac/sparse_matrix.h
deal.II/lac/include/lac/sparse_matrix.templates.h
deal.II/lac/include/lac/sparse_vanka.h
deal.II/lac/include/lac/sparse_vanka.templates.h
deal.II/lac/include/lac/sparsity_pattern.h
deal.II/lac/include/lac/vector.h
deal.II/lac/include/lac/vector.templates.h
deal.II/lac/include/lac/vector_memory.h
deal.II/lac/source/full_matrix.double.cc
deal.II/lac/source/full_matrix.float.cc
deal.II/lac/source/precondition_block.cc
deal.II/lac/source/solver_control.cc
deal.II/lac/source/sparse_ilu.cc
deal.II/lac/source/sparse_matrix.double.cc
deal.II/lac/source/sparse_matrix.float.cc
deal.II/lac/source/sparse_matrix.long_double.cc
deal.II/lac/source/sparse_vanka.cc
deal.II/lac/source/sparsity_pattern.cc
deal.II/lac/source/vector.cc
deal.II/lac/source/vector.long_double.cc

index 3c5754e7bcd0da92fd4f2e33c87a4fe692a974ea..82a5fb649e207cb55d2fe6ddbee55f3cac78932b 100644 (file)
@@ -1,6 +1,18 @@
-/*      $Id$                 */
-#ifndef __lac__block_vector_H
-#define __lac__block_vector_H
+//----------------------------  block_vector.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  block_vector.h  ---------------------------
+#ifndef __deal2__block_vector_h
+#define __deal2__block_vector_h
+
 
 #include <lac/vector.h>
 #include <base/exceptions.h>
@@ -136,9 +148,9 @@ class BlockVector
                                      * $l_\infty$-norm of a vector.
                                      */
     Number linfty_norm () const;
-    
-    
-                                    /**
+
+
+/**
                                      * Change the dimension of the vector to
                                      * #N#. The reserved memory for this vector
                                      * remains unchanged if possible, to make
@@ -214,9 +226,8 @@ class BlockVector
 //                                   */
 //      const_iterator end () const;  
 
-    
-    
-                                    /**
+
+/**
                                      * @name 2: Data-Access
                                      */
                                     //@{
@@ -231,9 +242,9 @@ class BlockVector
                                      */
     Number& operator() (const unsigned int i);
                                     //@}
-    
-    
-                                    /**
+
+
+/**
                                      * @name 3: Modification of vectors
                                      */
                                     //@{
@@ -329,9 +340,9 @@ class BlockVector
              const Number b, const BlockVector<n_blocks,Number>& W);
 
                                     //@}
-    
-    
-                                    /**
+
+
+/**
                                      * @name 5: Mixed stuff
                                      */
                                     //@{
@@ -429,10 +440,6 @@ class BlockVector
 };
 
 
-
-
-
-
 /*----------------------- Inline functions ----------------------------------*/
 
 
@@ -444,7 +451,6 @@ 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 () {
@@ -452,7 +458,6 @@ 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 {
@@ -460,7 +465,6 @@ 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 () {
@@ -468,7 +472,6 @@ unsigned int BlockVector<n_blocks,Number>::size () const
 //  };
 
 
-
 //  template <typename Number>
 //  inline
 //  BlockVector<n_blocks,Number>::const_iterator BlockVector<n_blocks,Number>::end () const {
@@ -476,7 +479,6 @@ 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
@@ -490,7 +492,6 @@ 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)
@@ -524,5 +525,4 @@ BlockVector<n_blocks,Number>::block(unsigned int i) const
 }
 
 
-
 #endif
index eaa93a2c1321c548f9965fba94af74644c643ac4..6d54d753bfa7dcc282cc1caf38d9514d83fb0e77 100644 (file)
@@ -1,9 +1,18 @@
-// $Id$
+//----------------------------  block_vector.templates.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  block_vector.templates.h  ---------------------------
+#ifndef __deal2__block_vector_templates_h
+#define __deal2__block_vector_templates_h
 
-// This file was once part of the DEAL Library
-// DEAL is Copyright(1995) by
-// Roland Becker, Guido Kanschat, Franz-Theo Suttmeier
-// Revised and extended by Wolfgang Bangerth, 1998, 1999
 
 #include <lac/block_vector.h>
 #include <cmath>
@@ -17,7 +26,6 @@ static inline Number sqr (const Number x)
 };
 
 
-
 template <int n_blocks, typename Number>
 BlockVector<n_blocks,Number>::BlockVector ()
 {
@@ -26,7 +34,6 @@ BlockVector<n_blocks,Number>::BlockVector ()
 }
 
 
-
 template <int n_blocks, typename Number>
 BlockVector<n_blocks,Number>::BlockVector (const vector<unsigned int>& n)
 {
@@ -36,7 +43,6 @@ 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)
 {
@@ -49,7 +55,6 @@ BlockVector<n_blocks,Number>::BlockVector (const BlockVector<n_blocks,Number>& v
 }
 
 
-
 // see the .h file for why this function was disabled
 //
 // template <int n_blocks, typename Number>
@@ -68,7 +73,6 @@ BlockVector<n_blocks,Number>::BlockVector (const BlockVector<n_blocks,Number>& v
 // }
 
 
-
 template <int n_blocks, typename Number>
 void BlockVector<n_blocks,Number>::reinit (const vector<unsigned int>& n,
                                            const bool fast)
@@ -87,7 +91,6 @@ 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,
                                            const bool fast)
@@ -101,15 +104,12 @@ void BlockVector<n_blocks,Number>::reinit (const BlockVector<n_blocks,Number>& v
 }
 
 
-
-
 template <int n_blocks, typename Number>
 BlockVector<n_blocks,Number>::~BlockVector ()
 {
 }
 
 
-
 template <int n_blocks, typename Number>
 void BlockVector<n_blocks,Number>::clear ()
 {
@@ -120,7 +120,6 @@ void BlockVector<n_blocks,Number>::clear ()
 }
 
 
-
 template <int n_blocks, typename Number>
 bool BlockVector<n_blocks,Number>::all_zero () const
 {
@@ -133,7 +132,6 @@ 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
 {
@@ -146,7 +144,6 @@ Number BlockVector<n_blocks,Number>::operator * (const BlockVector<n_blocks,Numb
 }
 
 
-
 template <int n_blocks, typename Number>
 Number BlockVector<n_blocks,Number>::norm_sqr () const
 {
@@ -159,7 +156,6 @@ Number BlockVector<n_blocks,Number>::norm_sqr () const
 };
 
 
-
 template <int n_blocks, typename Number>
 Number BlockVector<n_blocks,Number>::mean_value () const
 {
@@ -172,7 +168,6 @@ Number BlockVector<n_blocks,Number>::mean_value () const
 }
 
 
-
 template <int n_blocks, typename Number>
 Number BlockVector<n_blocks,Number>::l1_norm () const
 {
@@ -185,7 +180,6 @@ Number BlockVector<n_blocks,Number>::l1_norm () const
 }
 
 
-
 template <int n_blocks, typename Number>
 Number BlockVector<n_blocks,Number>::l2_norm () const
 {
@@ -193,7 +187,6 @@ Number BlockVector<n_blocks,Number>::l2_norm () const
 }
 
 
-
 template <int n_blocks, typename Number>
 Number BlockVector<n_blocks,Number>::linfty_norm () const
 {
@@ -206,9 +199,6 @@ Number BlockVector<n_blocks,Number>::linfty_norm () const
     }
   return sum;
 }
-  
-
-
 
 
 template <int n_blocks, typename Number>
@@ -219,7 +209,6 @@ BlockVector<n_blocks,Number>& BlockVector<n_blocks,Number>::operator += (const B
 }
 
 
-
 template <int n_blocks, typename Number>
 BlockVector<n_blocks,Number>& BlockVector<n_blocks,Number>::operator -= (const BlockVector<n_blocks,Number>& v)
 {
@@ -231,7 +220,6 @@ BlockVector<n_blocks,Number>& BlockVector<n_blocks,Number>::operator -= (const B
 }
 
 
-
 template <int n_blocks, typename Number>
 void BlockVector<n_blocks,Number>::add (const Number v)
 {
@@ -242,7 +230,6 @@ 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)
 {
@@ -253,7 +240,6 @@ 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)
@@ -265,7 +251,6 @@ 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,
@@ -279,7 +264,6 @@ 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)
@@ -291,7 +275,6 @@ 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)
@@ -303,7 +286,6 @@ 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,
@@ -317,7 +299,6 @@ 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,
@@ -334,7 +315,6 @@ 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)
 {
@@ -345,7 +325,6 @@ 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,
@@ -359,7 +338,6 @@ 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)
@@ -382,7 +360,6 @@ BlockVector<n_blocks,Number>& BlockVector<n_blocks,Number>::operator = (const Nu
 }
 
 
-
 template <int n_blocks, typename Number>
 BlockVector<n_blocks,Number>&
 BlockVector<n_blocks,Number>::operator = (const BlockVector<n_blocks,Number>& v)
@@ -395,7 +372,6 @@ BlockVector<n_blocks,Number>::operator = (const BlockVector<n_blocks,Number>& v)
 }
 
 
-
 template <int n_blocks, typename Number>
 template<typename Number2>
 BlockVector<n_blocks,Number>&
@@ -426,7 +402,6 @@ void BlockVector<n_blocks,Number>::print (ostream &out,
 }
 
 
-
 template <int n_blocks, typename Number>
 void BlockVector<n_blocks,Number>::block_write (ostream &out) const
 {
@@ -445,3 +420,5 @@ void BlockVector<n_blocks,Number>::block_read (istream &in)
       components[i].block_read(in);
     }  
 }
+
+#endif
index 5eadbb8a256756901d337a89d8f122bedd4db4ad..a7c7aa1a5408d092c63e290a90d7b39f56c6aa0f 100644 (file)
@@ -1,14 +1,20 @@
-/*----------------------------   fullmatrix.h     ---------------------------*/
-//      $Id$
-#ifndef __lac_fullmatrix_H
-#define __lac_fullmatrix_H
-/*----------------------------   fullmatrix.h     ---------------------------*/
+//----------------------------  full_matrix.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  full_matrix.h  ---------------------------
+#ifndef __deal2__full_matrix_h
+#define __deal2__full_matrix_h
 
-// This file is part of the DEAL Library
-// DEAL is Copyright(1995) by
-// Roland Becker, Guido Kanschat, Franz-Theo Suttmeier
-// Revised by Wolfgang Bangerth
 
+/*----------------------------   fullmatrix.h     ---------------------------*/
 
 #include <base/exceptions.h>
 #include <base/subscriptor.h>
@@ -19,7 +25,6 @@
 template<typename number> class Vector;
 
 
-
 /**
  * Rectangular/quadratic full matrix.
  *
@@ -128,9 +133,9 @@ class FullMatrix : public Subscriptor
                                      */
     template<typename number2>
     FullMatrix<number>& operator = (const FullMatrix<number2>& src);
-    
-    
-                                    /**
+
+
+/**
                                      * Fill rectangular block.
                                      *
                                      * The matrix #src# is copied
@@ -655,9 +660,6 @@ class FullMatrix : public Subscriptor
 };
 
 
-
-
-
 /*-------------------------Inline functions -------------------------------*/
 
 template <typename number>
@@ -720,9 +722,7 @@ FullMatrix<number>::operator() (const unsigned int i, const unsigned int j)
 }
 
 
-
-
 /*----------------------------   fullmatrix.h     ---------------------------*/
-/* end of #ifndef __lac_fullmatrix_H */
+
 #endif
 /*----------------------------   fullmatrix.h     ---------------------------*/
index d441ee34507421438087434c5e97a4c4a0938e10..d73fa76fe43b5d895f5c22042b8e77196a96efb2 100644 (file)
@@ -1,4 +1,18 @@
-// $Id$
+//----------------------------  full_matrix.templates.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  full_matrix.templates.h  ---------------------------
+#ifndef __deal2__full_matrix_templates_h
+#define __deal2__full_matrix_templates_h
+
 
 #include <lac/vector.h>
 #include <lac/full_matrix.h>
@@ -57,7 +71,6 @@ FullMatrix<number>::~FullMatrix ()
 };
 
 
-
 template <typename number>
 void
 FullMatrix<number>::reinit (const unsigned int mm,
@@ -103,7 +116,6 @@ FullMatrix<number>::reinit (const unsigned int mm,
 };
 
 
-
 template <typename number>
 void
 FullMatrix<number>::reinit (const unsigned int n)
@@ -137,7 +149,6 @@ FullMatrix<number>::all_zero () const
 };
 
 
-
 template <typename number>
 template <typename number2>
 void
@@ -269,7 +280,6 @@ FullMatrix<number>::vmult (Vector<number2>& dst,
 }
 
 
-
 template <typename number>
 template <typename number2>
 void FullMatrix<number>::Tvmult (Vector<number2>       &dst,
@@ -296,7 +306,6 @@ void FullMatrix<number>::Tvmult (Vector<number2>       &dst,
 }
 
 
-
 template <typename number>
 template <typename number2, typename number3>
 double FullMatrix<number>::residual (Vector<number2>& dst,
@@ -324,7 +333,6 @@ double FullMatrix<number>::residual (Vector<number2>& dst,
 }
 
 
-
 template <typename number>
 template <typename number2>
 void FullMatrix<number>::forward (Vector<number2>& dst,
@@ -347,7 +355,6 @@ void FullMatrix<number>::forward (Vector<number2>& dst,
 }
 
 
-
 template <typename number>
 template <typename number2>
 void FullMatrix<number>::backward (Vector<number2>& dst,
@@ -367,7 +374,6 @@ void FullMatrix<number>::backward (Vector<number2>& dst,
 }
 
 
-
 template <typename number>
 FullMatrix<number>&
 FullMatrix<number>::operator = (const FullMatrix<number>& m) 
@@ -381,7 +387,6 @@ FullMatrix<number>::operator = (const FullMatrix<number>& m)
 }
 
 
-
 template <typename number>
 template <typename number2>
 FullMatrix<number>&
@@ -396,7 +401,6 @@ FullMatrix<number>::operator = (const FullMatrix<number2>& m)
 }
 
 
-
 template <typename number>
 template <typename number2>
 void FullMatrix<number>::fill (const FullMatrix<number2>& src,
@@ -412,7 +416,6 @@ void FullMatrix<number>::fill (const FullMatrix<number2>& src,
 }
 
 
-
 template <typename number>
 void FullMatrix<number>::add_row (const unsigned int i,
                                  const number s,
@@ -425,7 +428,6 @@ void FullMatrix<number>::add_row (const unsigned int i,
 }
 
 
-
 template <typename number>
 void FullMatrix<number>::add_row (const unsigned int i,
                                  const number s,
@@ -441,7 +443,6 @@ void FullMatrix<number>::add_row (const unsigned int i,
 }
 
 
-
 template <typename number>
 void FullMatrix<number>::add_col (const unsigned int i, const number s,
                        const unsigned int j)
@@ -453,7 +454,6 @@ void FullMatrix<number>::add_col (const unsigned int i, const number s,
 }
 
 
-
 template <typename number>
 void FullMatrix<number>::add_col (const unsigned int i, const number s,
                        const unsigned int j, const number t,
@@ -466,7 +466,6 @@ void FullMatrix<number>::add_col (const unsigned int i, const number s,
 }
 
 
-
 template <typename number>
 void FullMatrix<number>::swap_row (const unsigned int i, const unsigned int j)
 {
@@ -480,7 +479,6 @@ void FullMatrix<number>::swap_row (const unsigned int i, const unsigned int j)
 }
 
 
-
 template <typename number>
 void FullMatrix<number>::swap_col (const unsigned int i, const unsigned int j)
 {
@@ -494,7 +492,6 @@ void FullMatrix<number>::swap_col (const unsigned int i, const unsigned int j)
 }
 
 
-
 template <typename number>
 void FullMatrix<number>::diagadd (const number src)
 {
@@ -506,7 +503,6 @@ void FullMatrix<number>::diagadd (const number src)
 }
 
 
-
 template <typename number>
 template <typename number2>
 void FullMatrix<number>::mmult (FullMatrix<number2>& dst,
@@ -529,7 +525,6 @@ void FullMatrix<number>::mmult (FullMatrix<number2>& dst,
 }
 
 
-
 /*void FullMatrix<number>::mmult (FullMatrix& dst, const FullMatrix& src) const
 {
   Assert (val != 0, ExcEmptyMatrix());  
@@ -550,7 +545,6 @@ void FullMatrix<number>::mmult (FullMatrix<number2>& dst,
 }*/
 
 
-
 template <typename number>
 template <typename number2>
 void FullMatrix<number>::Tmmult (FullMatrix<number2>& dst, const FullMatrix<number2>& src) const
@@ -572,7 +566,6 @@ void FullMatrix<number>::Tmmult (FullMatrix<number2>& dst, const FullMatrix<numb
 }
 
 
-
 /*void FullMatrix<number>::Tmmult(FullMatrix& dst, const FullMatrix& src) const
 {
   Assert (val != 0, ExcEmptyMatrix());  
@@ -593,7 +586,6 @@ void FullMatrix<number>::Tmmult (FullMatrix<number2>& dst, const FullMatrix<numb
 }*/
 
 
-
 template <typename number>
 template <typename number2>
 double FullMatrix<number>::matrix_norm (const Vector<number2> &v) const
@@ -623,7 +615,6 @@ double FullMatrix<number>::matrix_norm (const Vector<number2> &v) const
 };
 
 
-
 template <typename number>
 template <typename number2>
 double FullMatrix<number>::matrix_scalar_product (const Vector<number2> &u,
@@ -655,7 +646,6 @@ double FullMatrix<number>::matrix_scalar_product (const Vector<number2> &u,
 };
 
 
-
 template <typename number>
 number FullMatrix<number>::l1_norm () const
 {
@@ -676,7 +666,6 @@ number FullMatrix<number>::l1_norm () const
 };
 
 
-
 template <typename number>
 number FullMatrix<number>::linfty_norm () const
 {
@@ -697,7 +686,6 @@ number FullMatrix<number>::linfty_norm () const
 };
 
 
-
 template <typename number>
 void
 FullMatrix<number>::print (ostream& s, int w, int p) const
@@ -713,7 +701,6 @@ FullMatrix<number>::print (ostream& s, int w, int p) const
 }
 
 
-
 template <typename number>
 template <typename number2>
 void
@@ -836,7 +823,6 @@ FullMatrix<number>::add (const number s,const FullMatrix<number2>& src)
 }
 
 
-
 template <typename number>
 template <typename number2>
 void
@@ -959,7 +945,6 @@ FullMatrix<number>::add_diag (const number s, const FullMatrix<number2>& src)
 }
 
 
-
 template <typename number>
 template <typename number2>
 void
@@ -1086,7 +1071,6 @@ FullMatrix<number>::Tadd (const number s, const FullMatrix<number2>& src)
 }
 
 
-
 template <typename number>
 bool
 FullMatrix<number>::operator == (const FullMatrix<number> &m) const
@@ -1101,7 +1085,6 @@ FullMatrix<number>::operator == (const FullMatrix<number> &m) const
 };
 
 
-
 template <typename number>
 double
 FullMatrix<number>::determinant () const
@@ -1131,7 +1114,6 @@ FullMatrix<number>::determinant () const
 };
 
 
-
 template <typename number>
 number
 FullMatrix<number>::norm2 () const
@@ -1145,7 +1127,6 @@ FullMatrix<number>::norm2 () const
 }
 
 
-
 template <typename number>
 void FullMatrix<number>::clear ()
 {
@@ -1154,7 +1135,6 @@ void FullMatrix<number>::clear ()
 };
 
 
-
 template <typename number>
 void
 FullMatrix<number>::invert (const FullMatrix<number> &M)
@@ -1311,7 +1291,6 @@ FullMatrix<number>::invert (const FullMatrix<number> &M)
            AssertThrow (false, ExcNotImplemented(dim_range));
     };    
 };
-  
 
 
 template <typename number>
@@ -1453,7 +1432,6 @@ FullMatrix<number>::householder(Vector<number2>& src)
 }
 
 
-
 template <typename number>
 template <typename number2>
 double
@@ -1470,3 +1448,6 @@ FullMatrix<number>::least_squares(Vector<number2>& dst, Vector<number2>& src)
   for (unsigned int i=n() ; i<m() ; ++i) sum += src(i) * src(i);
   return sqrt(sum);
 }
+
+
+#endif
index c59b6c51a15ea5bf62ee87b53673bf92e5b238dc..2d9390bf6bff5e8e84bfadeceacdca33079865bd 100644 (file)
@@ -1,10 +1,17 @@
-/*----------------------------   precondition.h     ---------------------------*/
-/*      $Id$                 */
-#ifndef __precondition_H
-#define __precondition_H
-/*----------------------------   precondition.h     ---------------------------*/
-
-
+//----------------------------  precondition.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  precondition.h  ---------------------------
+#ifndef __deal2__precondition_h
+#define __deal2__precondition_h
 
 
 /**
@@ -26,8 +33,6 @@ class PreconditionIdentity
 };
 
 
-
-
 /**
  * Preconditioner using a matrix-builtin function.
  * This class forms a preconditioner suitable for the LAC solver
@@ -115,7 +120,6 @@ class PreconditionUseMatrix
 };
 
 
-
 /**
  * Preconditioner for builtin relaxation methods.
  * Application of this preconditioner includes
@@ -250,7 +254,6 @@ class PreconditionLACSolver
 };
 
 
-
 /* ---------------------------------- Inline functions ------------------- */
 
 template<class VECTOR>
@@ -314,9 +317,4 @@ PreconditionLACSolver<SOLVER,MATRIX,PRECONDITION>::operator() (VECTOR& dst,
 }
 
 
-
-
-/*----------------------------   precondition.h     ---------------------------*/
-/* end of #ifndef __precondition_H */
 #endif
-/*----------------------------   precondition.h     ---------------------------*/
index beca0f1f7f3fd899c869e93ec60cc6e8a2f4329b..5c64cf7236feb7917b6a6d4b8c6573d81a72922c 100644 (file)
@@ -1,8 +1,18 @@
-/*----------------------------   precondition_block.h     ---------------------------*/
-/*      $Id$                 */
-#ifndef __precondition_block_H
-#define __precondition_block_H
-/*----------------------------   precondition_block.h     ---------------------------*/
+//----------------------------  precondition_block.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  precondition_block.h  ---------------------------
+#ifndef __deal2__precondition_block_h
+#define __deal2__precondition_block_h
+
 
 #include <lac/forward_declarations.h>
 #include <base/exceptions.h>
@@ -190,7 +200,6 @@ class PreconditionBlock: public Subscriptor
 };
 
 
-
 /**
  * Block SOR preconditioning.
  *
@@ -243,9 +252,5 @@ class PreconditionBlockSOR : public PreconditionBlock<number,inverse_type>
 };
 
 
-
-/*------------------------   precondition_block.h     ---------------------------*/
-/* end of #ifndef __precondition_block_H */
 #endif
-/*------------------------   precondition_block.h     ---------------------------*/
 
index 6f373177be998429a81c1bdb9c61b8fdc8476a06..3084c1fb19e934b9ecf3969954fc48915becbc5b 100644 (file)
@@ -1,6 +1,17 @@
-/*-----------------   precondition_block.templates.h     ----------------------*/
-/*      $Id$                 */
-/*           Ralf Hartmann, University of Heidelberg                           */
+//----------------------------  precondition_block.templates.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  precondition_block.templates.h  ---------------------------
+#ifndef __deal2__precondition_block_templates_h
+#define __deal2__precondition_block_templates_h
 
 
 #include <lac/precondition_block.h>
@@ -15,7 +26,6 @@ PreconditionBlock<number,inverse_type>::PreconditionBlock ():
                A(0)  {};
 
 
-
 template <typename number, typename inverse_type>
 PreconditionBlock<number,inverse_type>::~PreconditionBlock ()
 {
@@ -24,7 +34,6 @@ PreconditionBlock<number,inverse_type>::~PreconditionBlock ()
 }
 
 
-
 template <typename number, typename inverse_type>
 void PreconditionBlock<number,inverse_type>::clear ()
 {
@@ -34,7 +43,6 @@ void PreconditionBlock<number,inverse_type>::clear ()
 }
 
 
-
 template <typename number, typename inverse_type>
 void PreconditionBlock<number,inverse_type>::use_matrix(
   const SparseMatrix<number> &M)
@@ -43,21 +51,18 @@ void PreconditionBlock<number,inverse_type>::use_matrix(
 }
 
 
-
 template <typename number, typename inverse_type>
 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 {
   return blocksize;
 }
 
 
-
 template <typename number, typename inverse_type>
 void PreconditionBlock<number,inverse_type>::invert_diagblocks()
 {
@@ -119,7 +124,6 @@ void PreconditionBlock<number,inverse_type>::invert_diagblocks()
 }
 
 
-
 /*--------------------- PreconditionBlockSOR -----------------------*/
 
 template<typename number, typename inverse_type>
@@ -127,12 +131,10 @@ PreconditionBlockSOR<number,inverse_type>::PreconditionBlockSOR(const number ome
                omega(omega)  {}
 
 
-
 template<typename number, typename inverse_type>
 PreconditionBlockSOR<number,inverse_type>::~PreconditionBlockSOR(){}
 
 
-
 template <typename number, typename inverse_type>
 template <typename number2>
 void PreconditionBlockSOR<number,inverse_type>::operator() (Vector<number2>       &dst,
@@ -211,7 +213,4 @@ void PreconditionBlockSOR<number,inverse_type>::operator() (Vector<number2>
 }
 
 
-
-
-
-/*----------------------   precondition_block.templates.h     ----------------*/
+#endif
index 0b321c392c370ed6565138565da065184d8a7e2e..b7e8c7cde3c302c5a30063cfbaef8c0afbb81570 100644 (file)
@@ -1,9 +1,17 @@
-/*----------------------   precondition_selector.h     -------------------------*/
-/*      $Id$                 */
-/*                Ralf Hartmann, University of Heidelberg                       */
-#ifndef __precondition_selector_H
-#define __precondition_selector_H
-/*----------------------   precondition_selector.h     -------------------------*/
+//----------------------------  precondition_selector.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  precondition_selector.h  ---------------------------
+#ifndef __deal2__precondition_selector_h
+#define __deal2__precondition_selector_h
 
 
 #include <base/smartpointer.h>
@@ -116,7 +124,7 @@ class PreconditionSelector
     DeclException0 (ExcNoMatrixGivenToUse);
 
 
-  protected:
+protected:
     
                                     /**
                                      * Stores the name of the
@@ -202,8 +210,4 @@ string PreconditionSelector<Matrix,Vector>::get_precondition_names()
 }
 
 
-
-/*----------------------------   precondition_selector.h     ---------------------------*/
-/* end of #ifndef __precondition_selector_H */
 #endif
-/*----------------------------   precondition_selector.h     ---------------------------*/
index 530a8b2f5c9542d46606fc9a1ce2382bb50e51a1..98a1742317f1ddb5e45329d3368180e504475d50 100644 (file)
@@ -1,15 +1,22 @@
-/*----------------------------   solver.h     ---------------------------*/
-/*      $Id$                 */
-#ifndef __solver_H
-#define __solver_H
-/*----------------------------   solver.h     ---------------------------*/
-
+//----------------------------  solver.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  solver.h  ---------------------------
+#ifndef __deal2__solver_h
+#define __deal2__solver_h
 
 
 #include <lac/forward_declarations.h>
 
 
-
 /**
  * Base class for iterative solvers.
  *
@@ -134,10 +141,9 @@ class Solver
     enum ReturnState {
          success=0, exceeded, breakdown
     };
-    
-    
-    
-                                    /**
+
+
+/**
                                      * Constructor. Assign a control
                                      * object which stores the required
                                      * precision and an object to provide
@@ -181,7 +187,6 @@ SolverControl & Solver<Matrix,Vector>::control() const
 };
 
 
-
 template<class Matrix, class Vector>
 inline
 Solver<Matrix, Vector>::Solver(SolverControl &cn, VectorMemory<Vector> &mem)
@@ -190,7 +195,4 @@ Solver<Matrix, Vector>::Solver(SolverControl &cn, VectorMemory<Vector> &mem)
 {};
 
 
-/*----------------------------   solver.h     ---------------------------*/
-/* end of #ifndef __solver_H */
 #endif
-/*----------------------------   solver.h     ---------------------------*/
index 097b58b144e94998ab034735e3d4b2a0df0192c1..32e15a9fbf05cb7931373de3832d3e9f3ab8b6f0 100644 (file)
@@ -1,8 +1,18 @@
-/*----------------------------   solver_bicgstab.h     ---------------------------*/
-/*      $Id$                 */
-#ifndef __solver_bicgstab_H
-#define __solver_bicgstab_H
-/*----------------------------   solver_bicgstab.h     ---------------------------*/
+//----------------------------  solver_bicgstab.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  solver_bicgstab.h  ---------------------------
+#ifndef __deal2__solver_bicgstab_h
+#define __deal2__solver_bicgstab_h
+
 
 #include <base/logstream.h>
 #include <lac/solver.h>
@@ -185,7 +195,6 @@ SolverBicgstab<Matrix, Vector>::start()
 }
 
 
-
 template<class Matrix, class Vector>
 template<class Preconditioner>
 typename Solver<Matrix,Vector>::ReturnState
@@ -284,7 +293,4 @@ SolverBicgstab<Matrix, Vector>::solve(const Matrix &A,
 }
 
 
-/*----------------------------   solver_bicgstab.h     ---------------------------*/
-/* end of #ifndef __solver_bicgstab_H */
 #endif
-/*----------------------------   solver_bicgstab.h     ---------------------------*/
index e71d5026915b36ee4b2b2f5e3010c00b1c196165..f046c7752943ac5b39df0d2680fac10d9b130a79 100644 (file)
@@ -1,9 +1,17 @@
-/*----------------------------   solver_cg.h     ---------------------------*/
-/*      $Id$                 */
-#ifndef __solver_cg_H
-#define __solver_cg_H
-/*----------------------------   solver_cg.h     ---------------------------*/
-
+//----------------------------  solver_cg.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  solver_cg.h  ---------------------------
+#ifndef __deal2__solver_cg_h
+#define __deal2__solver_cg_h
 
 
 #include <lac/solver.h>
@@ -12,7 +20,6 @@
 #include <cmath>
 
 
-
 /**
  * Preconditioned cg method.
  *
@@ -98,10 +105,8 @@ class SolverCG : public Solver<Matrix,Vector>
 };
 
 
-
-
 /*------------------------- Implementation ----------------------------*/
+
 
 template<class Matrix, class Vector>
 SolverCG<Matrix,Vector>::SolverCG(SolverControl &cn,
@@ -118,7 +123,6 @@ SolverCG<Matrix,Vector>::criterion()
 };
 
 
-
 template<class Matrix, class Vector>
 template<class Preconditioner>
 typename Solver<Matrix,Vector>::ReturnState 
@@ -196,8 +200,8 @@ SolverCG<Matrix,Vector>::solve (const Matrix &A,
       it++;
     };
 
-    
-  // Deallocate Memory
+
+// Deallocate Memory
  
   memory.free(Vr);
   memory.free(Vp);
@@ -215,9 +219,4 @@ SolverCG<Matrix,Vector>::solve (const Matrix &A,
 };
 
 
-
-
-/*----------------------------   solver_cg.h     ---------------------------*/
-/* end of #ifndef __solver_cg_H */
 #endif
-/*----------------------------   solver_cg.h     ---------------------------*/
index 4e24bcd645e9d9bcf6a7670e99c58c887b701597..343463c643a65ac0be968b7d6dd792e2cf652bb3 100644 (file)
@@ -1,8 +1,17 @@
-/*----------------------------   solver_control.h     ---------------------------*/
-/*      $Id$                 */
-#ifndef __solver_control_H
-#define __solver_control_H
-/*----------------------------   solver_control.h     ---------------------------*/
+//----------------------------  solver_control.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  solver_control.h  ---------------------------
+#ifndef __deal2__solver_control_h
+#define __deal2__solver_control_h
 
 
 #include <base/subscriptor.h>
@@ -161,10 +170,6 @@ class SolverControl : public Subscriptor
 };
 
 
-
-
-
-
 /**
  * Specialization of #SolverControl# which returns #success# if either
  * the specified tolerance is achieved or if the initial residual (or
@@ -216,9 +221,9 @@ class ReductionControl : public SolverControl
                                      * criterion.
                                      */
     double initial_value() const;
-    
 
-  protected:
+
+protected:
                                     /**
                                      * Desired reduction factor.
                                      */
@@ -239,9 +244,4 @@ class ReductionControl : public SolverControl
 };
 
 
-
-
-/*----------------------------   solver_control.h     ---------------------------*/
-/* end of #ifndef __solver_control_H */
 #endif
-/*----------------------------   solver_control.h     ---------------------------*/
index 1dbd82e1c6b1027e037a6b0ba44b2e2d6032aeb1..d20e0bdb5f0e8f7a59a3638037aab6187d8e3a2c 100644 (file)
@@ -1,7 +1,19 @@
-/*----------------------------   solver_pgmres.h     ---------------------------*/
-/*      $Id$                 */
-#ifndef __solver_pgmres_H
-#define __solver_pgmres_H
+//----------------------------  solver_gmres.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  solver_gmres.h  ---------------------------
+#ifndef __deal2__solver_gmres_h
+#define __deal2__solver_gmres_h
+
+
 /*----------------------------   solver_pgmres.h     ---------------------------*/
 
 
@@ -12,7 +24,6 @@
 #include <vector>
 
 
-
 /**
  * Implementation of the Restarted Preconditioned Direct Generalized
  * Minimal Residual Method. The stopping criterion is the norm of the
@@ -135,8 +146,6 @@ class SolverGMRES : public Solver<MATRIX, VECTOR>
 };
 
 
-
-
 /* --------------------- Inline and template functions ------------------- */
 
 
@@ -152,7 +161,6 @@ SolverGMRES<MATRIX,VECTOR>::SolverGMRES (SolverControl        &cn,
 };
 
 
-
 template <class MATRIX, class VECTOR>
 inline
 void
@@ -180,8 +188,6 @@ SolverGMRES<MATRIX,VECTOR>::givens_rotation (Vector<double> &h,
 }
 
 
-
-
 template<class MATRIX, class VECTOR>
 template<class Preconditioner>
 typename Solver<MATRIX,VECTOR>::ReturnState
@@ -236,7 +242,7 @@ SolverGMRES<MATRIX,VECTOR>::solve (const MATRIX& A,
                    h    (n_tmp_vectors-1);
 
 
-  unsigned int dim = 0;
+unsigned int dim = 0;
 
   SolverControl::State iteration_state = SolverControl::iterate;
   
@@ -250,8 +256,8 @@ SolverGMRES<MATRIX,VECTOR>::solve (const MATRIX& A,
   VECTOR &v = *tmp_vectors[0];
   VECTOR &p = *tmp_vectors[n_tmp_vectors-1];
 
-  
-                                  ///////////////////////////////////
+
+///////////////////////////////////
                                   // outer iteration: loop until we
                                   // either reach convergence or the
                                   // maximum number of iterations is
@@ -394,9 +400,6 @@ SolverGMRES<MATRIX,VECTOR>::solve (const MATRIX& A,
 };
 
 
-  
-
-
 template<class MATRIX, class VECTOR>
 double
 SolverGMRES<MATRIX,VECTOR>::criterion () 
@@ -409,8 +412,7 @@ SolverGMRES<MATRIX,VECTOR>::criterion ()
 };
 
 
-
 /*----------------------------   solver_pgmres.h     ---------------------------*/
-/* end of #ifndef __solver_pgmres_H */
+
 #endif
 /*----------------------------   solver_pgmres.h     ---------------------------*/
index 1dcda2651e0a8450bca56724f4d25f26c0887041..b9e7a04ca3f651ba64eb217a7122fcf405af871a 100644 (file)
@@ -1,8 +1,18 @@
-/*----------------------------   solver_qmrs.h     ---------------------------*/
-/*      $Id$                 */
-#ifndef __lac__solver_qmrs_H
-#define __lac__solver_qmrs_H
-/*----------------------------   solver_qmrs.h     ---------------------------*/
+//----------------------------  solver_qmrs.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  solver_qmrs.h  ---------------------------
+#ifndef __deal2__solver_qmrs_h
+#define __deal2__solver_qmrs_h
+
 
 //TODO: Check for programming errors!!!
 
@@ -163,10 +173,8 @@ class SolverQMRS : public Solver<Matrix,Vector>
 };
 
 
-
-
 /*------------------------- Implementation ----------------------------*/
+
 
 template<class Matrix, class Vector>
 SolverQMRS<Matrix,Vector>::SolverQMRS(SolverControl &cn,
@@ -185,7 +193,6 @@ SolverQMRS<Matrix,Vector>::criterion()
 };
 
 
-
 template<class Matrix, class Vector>
 template<class Preconditioner>
 typename Solver<Matrix,Vector>::ReturnState 
@@ -289,8 +296,8 @@ SolverQMRS<Matrix,Vector>::iterate(const Preconditioner& precondition)
   rho = q*v;
   //deallog << "rho:" << rho << endl;
 
-  
-  while (state == SolverControl::iterate)
+
+while (state == SolverControl::iterate)
     {
       step++; it++;
                                       // Step 1
@@ -345,8 +352,4 @@ SolverQMRS<Matrix,Vector>::iterate(const Preconditioner& precondition)
 }
 
 
-
-/*----------------------------   solver_qmrs.h     ---------------------------*/
-/* end of #ifndef __solver_qmrs_H */
 #endif
-/*----------------------------   solver_qmrs.h     ---------------------------*/
index c1db38f4bf67098848f2f2deba6d4aa9eb98db28..cc327c6860e145e20344a8f9d88bed25ce9e5b43 100644 (file)
@@ -1,17 +1,23 @@
-/*----------------------------   solver_richardson.h     ---------------------------*/
-/*   $Id$              */
-/*            Ralf Hartmann, University of Heidelberg                               */
-#ifndef __solver_richardson_H
-#define __solver_richardson_H
-/*----------------------------   solver_richardson.h     ---------------------------*/
-
+//----------------------------  solver_richardson.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  solver_richardson.h  ---------------------------
+#ifndef __deal2__solver_richardson_h
+#define __deal2__solver_richardson_h
 
 
 #include <lac/solver.h>
 #include <lac/solver_control.h>
 
 
-
 /**
  * Implementation of the richardson iteration method. The stopping criterion
  * is the norm of the residual.
@@ -115,8 +121,6 @@ class SolverRichardson : public Solver<Matrix, Vector>
 };
 
 
-
-
 /*----------------- Implementation of the Richardson Method ------------------*/
 
 
@@ -186,7 +190,4 @@ SolverRichardson<Matrix,Vector>::set_omega(double om)
 }
 
 
-/*------------------   solver_richardson.h     ----------------------*/
-/* end of #ifndef __solver_richardson_H */
 #endif
-/*------------------   solver_richardson.h     ----------------------*/
index c1b4c3c48f827be893d81ff49187ea3e94b62f2c..aa47cfbfc65beeffa631aef7a4fd0af4dcfe8113 100644 (file)
@@ -1,9 +1,18 @@
-/*----------------------------   solver_selector.h     ---------------------------*/
-/*      $Id$                 */
-/*                Ralf Hartmann, University of Heidelberg                         */
-#ifndef __solver_selector_H
-#define __solver_selector_H
-/*----------------------------   solver_selector.h     ---------------------------*/
+//----------------------------  solver_selector.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  solver_selector.h  ---------------------------
+#ifndef __deal2__solver_selector_h
+#define __deal2__solver_selector_h
+
 
 #include <base/smartpointer.h>
 #include <lac/solver.h>
@@ -193,7 +202,6 @@ class SolverSelector
 };
 
 
-
 /* --------------------- Inline and template functions ------------------- */
 
 
@@ -207,13 +215,11 @@ SolverSelector<Matrix, Vector>::SolverSelector(string solver_name,
 {};
 
 
-
 template <class Matrix, class Vector>
 SolverSelector<Matrix, Vector>::~SolverSelector()
 {};
 
 
-
 template <class Matrix, class Vector>
 template<class Preconditioner>
 typename Solver<Matrix,Vector>::ReturnState 
@@ -253,7 +259,6 @@ SolverSelector<Matrix, Vector>::solve(const Matrix &A,
 };
 
 
-
 template <class Matrix, class Vector>
 string SolverSelector<Matrix, Vector>::get_solver_names()
 {
@@ -261,7 +266,6 @@ string SolverSelector<Matrix, Vector>::get_solver_names()
 };
 
 
-
 template <class Matrix, class Vector>
 void SolverSelector<Matrix, Vector>::set_data(
   const typename SolverGMRES<Matrix,Vector>::AdditionalData &data)
@@ -294,8 +298,4 @@ void SolverSelector<Matrix, Vector>::set_data(
 };
 
 
-
-/*----------------------------   solver_selector.h     ---------------------------*/
-/* end of #ifndef __solver_selector_H */
 #endif
-/*----------------------------   solver_selector.h     ---------------------------*/
index 51b9fac3c7b3fe6e13d641e6f1156a4d2c0aeb4b..0a51693e6641cc953b9d0372f81d8793fc3deb88 100644 (file)
@@ -1,15 +1,22 @@
-/*----------------------------   matrix_decompositions.h     ---------------------------*/
-/*      $Id$                 */
-#ifndef __sparse_ilu_H
-#define __sparse_ilu_H
-/*----------------------------   sparse_ilu.h     ---------------------------*/
+//----------------------------  sparse_ilu.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  sparse_ilu.h  ---------------------------
+#ifndef __deal2__sparse_ilu_h
+#define __deal2__sparse_ilu_h
 
 
 #include <lac/sparse_matrix.h>
 
 
-
-
 /**
  * Incomplete LU decomposition of a sparse matrix into another sparse matrix.
  * A given matrix is decomposed into a incomplete LU factorization, where
@@ -207,9 +214,4 @@ SparseILU<number>::operator() (Vector<somenumber>       &dst,
 }
 
 
-
-
-/*----------------------------   sparse_ilu.h     ---------------------------*/
-/* end of #ifndef __sparse_ilu_H */
 #endif
-/*----------------------------   sparse_ilu.h     ---------------------------*/
index e3569dbdf2c6c7bd0c7dda6c060cd973083028e4..d930662c08f54f586551cbfbb9c9c6a69e4c185f 100644 (file)
@@ -1,7 +1,19 @@
-/*----------------------------   matrix_decompositions.templates.h     ---------------------------*/
-/*      $Id$                 */
-#ifndef __matrix_decompositions_templates_H
-#define __matrix_decompositions_templates_H
+//----------------------------  sparse_ilu.templates.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  sparse_ilu.templates.h  ---------------------------
+#ifndef __deal2__sparse_ilu_templates_h
+#define __deal2__sparse_ilu_templates_h
+
+
 /*----------------------------   matrix_decompositions.templates.h     ---------------------------*/
 
 
 #include <cmath>
 
 
-
 template <typename number>
 SparseILU<number>::SparseILU () 
 {};
 
 
-
 template <typename number>
 SparseILU<number>::SparseILU (const SparsityPattern &sparsity) :
                SparseMatrix<number> (sparsity)
 {};
 
 
-
 template <typename number>
 void SparseILU<number>::reinit ()
 {
@@ -33,7 +42,6 @@ void SparseILU<number>::reinit ()
 };
 
 
-
 template <typename number>
 void SparseILU<number>::reinit (const SparsityPattern &sparsity)
 {
@@ -41,8 +49,6 @@ void SparseILU<number>::reinit (const SparsityPattern &sparsity)
 };
 
 
-
-
 template <typename number>
 template <typename somenumber>
 void SparseILU<number>::decompose (const SparseMatrix<somenumber> &matrix,
@@ -53,9 +59,9 @@ void SparseILU<number>::decompose (const SparseMatrix<somenumber> &matrix,
   Assert (matrix.m()==m(),        ExcSizeMismatch(matrix.m(), m()));
   
   Assert (strengthen_diagonal>=0, ExcInvalidStrengthening (strengthen_diagonal));
-  
-  
-                                  // first thing: copy over all elements
+
+
+// first thing: copy over all elements
                                   // of #matrix# to the present object
                                   //
                                   // note that some elements in this
@@ -104,9 +110,9 @@ void SparseILU<number>::decompose (const SparseMatrix<somenumber> &matrix,
 
        *diagonal_element += strengthen_diagonal * rowsum;
       };
-  
-         
-                                  // now work only on this
+
+
+// now work only on this
                                   // matrix
   const SparsityPattern             &sparsity = get_sparsity_pattern();
   const unsigned int * const rowstart_indices = sparsity.get_rowstart_indices();
@@ -127,8 +133,8 @@ void SparseILU<number>::decompose (const SparseMatrix<somenumber> &matrix,
           a[i,j] -= a[i,k] * a[k,j]
 */
 
-  
-                                  // i := row
+
+// i := row
   for (unsigned int row=1; row<m(); ++row)
     {
                                       // invert diagonal element of the
@@ -260,8 +266,6 @@ void SparseILU<number>::decompose (const SparseMatrix<somenumber> &matrix,
 };
 
 
-
-
 template <typename number>
 template <typename somenumber>
 void SparseILU<number>::apply_decomposition (Vector<somenumber>       &dst,
@@ -332,9 +336,7 @@ void SparseILU<number>::apply_decomposition (Vector<somenumber>       &dst,
 };
 
 
-
-
 /*----------------------------   matrix_decompositions.templates.h     ---------------------------*/
-/* end of #ifndef __matrix_decompositions_templates_H */
+
 #endif
 /*----------------------------   matrix_decompositions.templates.h     ---------------------------*/
index 865b835b97cf6b90f299b1b0842665e406a561cc..0755a0dcce48f71d1561795dc0c1714b3c3d28d0 100644 (file)
@@ -1,4 +1,18 @@
-// $Id$
+//----------------------------  sparse_matrix.2.templates  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  sparse_matrix.2.templates  ---------------------------
+
+
+
 // Driver file for SparseMatrix functions with two types.
 
 // TYPEMAT and TYPE2 are defined in sparsematrix?.cc
index 911b599239e2a1e0aa30ebc06dac6df1fcb236c1..3aa504e119a1446451a21294c016bd3cf3001e7d 100644 (file)
@@ -1,14 +1,20 @@
-/*----------------------------   sparsematrix.h     ---------------------------*/
-/*      $Id$                 */
-#ifndef __sparsematrix_H
-#define __sparsematrix_H
-/*----------------------------   sparsematrix.h     ---------------------------*/
+//----------------------------  sparse_matrix.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  sparse_matrix.h  ---------------------------
+#ifndef __deal2__sparse_matrix_h
+#define __deal2__sparse_matrix_h
 
 
-// This file was once part of the DEAL Library
-// DEAL is Copyright(1995) by
-// Roland Becker, Guido Kanschat, Franz-Theo Suttmeier
-// Revised, modified and extended by Wolfgang Bangerth
+/*----------------------------   sparsematrix.h     ---------------------------*/
 
 
 #include <base/exceptions.h>
@@ -17,7 +23,6 @@
 #include <lac/sparsity_pattern.h>
 
 
-
 /**
  * Sparse matrix.
  *
@@ -63,9 +68,9 @@ class SparseMatrix : public Subscriptor
                                      * #copy_from# function.
                                      */
     SparseMatrix (const SparseMatrix &);
-    
-    
-                                    /**
+
+
+/**
                                      * Constructor. Takes the given matrix
                                      * sparsity structure to represent the
                                      * sparsity pattern of this matrix. You
@@ -91,9 +96,8 @@ class SparseMatrix : public Subscriptor
                                      */
     SparseMatrix<number>& operator = (const SparseMatrix<number> &);
 
-    
 
-                                    /**
+/**
                                      * Reinitialize the object but keep to
                                      * the sparsity pattern previously used.
                                      * This may be necessary if you #reinit#'d
@@ -645,18 +649,14 @@ class SparseMatrix : public Subscriptor
                            somenumber               *partial_norm) const;
 
 
-                                    // make all other sparse matrices
+// make all other sparse matrices
                                     // friends
     template <typename somenumber> friend class SparseMatrix;
 };
 
 
-
-
-
 /*---------------------- Inline functions -----------------------------------*/
 
 
 template <typename number>
 inline
@@ -667,7 +667,6 @@ unsigned int SparseMatrix<number>::m () const
 };
 
 
-
 template <typename number>
 inline
 unsigned int SparseMatrix<number>::n () const
@@ -677,7 +676,6 @@ unsigned int SparseMatrix<number>::n () const
 };
 
 
-
 template <typename number>
 inline
 void SparseMatrix<number>::set (const unsigned int i,
@@ -694,7 +692,6 @@ void SparseMatrix<number>::set (const unsigned int i,
 };
 
 
-
 template <typename number>
 inline
 void SparseMatrix<number>::add (const unsigned int i, const unsigned int j,
@@ -709,9 +706,6 @@ void SparseMatrix<number>::add (const unsigned int i, const unsigned int j,
 };
 
 
-
-
-
 template <typename number>
 inline
 number SparseMatrix<number>::operator () (const unsigned int i,
@@ -724,7 +718,6 @@ number SparseMatrix<number>::operator () (const unsigned int i,
 };
 
 
-
 template <typename number>
 inline
 number SparseMatrix<number>::diag_element (const unsigned int i) const
@@ -740,7 +733,6 @@ number SparseMatrix<number>::diag_element (const unsigned int i) const
 };
 
 
-
 template <typename number>
 inline
 number & SparseMatrix<number>::diag_element (const unsigned int i)
@@ -756,7 +748,6 @@ number & SparseMatrix<number>::diag_element (const unsigned int i)
 };
 
 
-
 template <typename number>
 inline
 number
@@ -770,7 +761,6 @@ SparseMatrix<number>::raw_entry (const unsigned int row,
 };
 
 
-
 template <typename number>
 inline
 number SparseMatrix<number>::global_entry (const unsigned int j) const
@@ -780,7 +770,6 @@ number SparseMatrix<number>::global_entry (const unsigned int j) const
 };
 
 
-
 template <typename number>
 inline
 number & SparseMatrix<number>::global_entry (const unsigned int j)
@@ -790,9 +779,8 @@ number & SparseMatrix<number>::global_entry (const unsigned int j)
 };
 
 
-
 /*----------------------------   sparsematrix.h     ---------------------------*/
-/* end of #ifndef __sparsematrix_H */
+
 #endif
 /*----------------------------   sparsematrix.h     ---------------------------*/
 
index 55a63344b43b3359b4d51da9dad2092bfc974193..c9e7d1a7ef68b578600f9dbcc2c00fd61b0682be 100644 (file)
@@ -1,9 +1,17 @@
-// $Id$
-
-// This file was once part of the DEAL Library
-// DEAL is Copyright(1995) by
-// Roland Becker, Guido Kanschat, Franz-Theo Suttmeier
-// Revised, modified and extended by Wolfgang Bangerth, 1998, 1999
+//----------------------------  sparse_matrix.templates.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  sparse_matrix.templates.h  ---------------------------
+#ifndef __deal2__sparse_matrix_templates_h
+#define __deal2__sparse_matrix_templates_h
 
 
 #include <lac/sparse_matrix.h>
@@ -24,9 +32,6 @@
 #endif
 
 
-
-
-
 template <typename number>
 SparseMatrix<number>::SparseMatrix () :
                cols(0),
@@ -48,7 +53,6 @@ SparseMatrix<number>::SparseMatrix (const SparseMatrix &m) :
 };
 
 
-
 template <typename number>
 SparseMatrix<number>&
 SparseMatrix<number>::operator = (const SparseMatrix<number> &m)
@@ -61,7 +65,6 @@ SparseMatrix<number>::operator = (const SparseMatrix<number> &m)
 };
 
 
-
 template <typename number>
 SparseMatrix<number>::SparseMatrix (const SparsityPattern &c) :
                cols(&c),
@@ -73,7 +76,6 @@ SparseMatrix<number>::SparseMatrix (const SparsityPattern &c) :
 };
 
 
-
 template <typename number>
 SparseMatrix<number>::~SparseMatrix ()
 {
@@ -86,7 +88,6 @@ SparseMatrix<number>::~SparseMatrix ()
 };
 
 
-
 template <typename number>
 void
 SparseMatrix<number>::reinit ()
@@ -115,7 +116,6 @@ SparseMatrix<number>::reinit ()
 }
 
 
-
 template <typename number>
 void
 SparseMatrix<number>::reinit (const SparsityPattern &sparsity)
@@ -130,7 +130,6 @@ SparseMatrix<number>::reinit (const SparsityPattern &sparsity)
 };
 
 
-
 template <typename number>
 void
 SparseMatrix<number>::clear ()
@@ -144,7 +143,6 @@ SparseMatrix<number>::clear ()
 };
 
 
-
 template <typename number>
 bool
 SparseMatrix<number>::empty () const
@@ -156,7 +154,6 @@ SparseMatrix<number>::empty () const
 };
 
 
-
 template <typename number>
 unsigned int
 SparseMatrix<number>::n_nonzero_elements () const
@@ -166,7 +163,6 @@ SparseMatrix<number>::n_nonzero_elements () const
 };
 
 
-
 template <typename number>
 template <typename somenumber>
 SparseMatrix<number> &
@@ -183,7 +179,6 @@ SparseMatrix<number>::copy_from (const SparseMatrix<somenumber> &matrix)
 };
 
 
-
 template <typename number>
 template <typename somenumber>
 void
@@ -203,7 +198,6 @@ SparseMatrix<number>::add_scaled (const number factor,
 };
 
 
-
 template <typename number>
 template <typename somenumber>
 void
@@ -274,7 +268,6 @@ SparseMatrix<number>::vmult (Vector<somenumber>& dst, const Vector<somenumber>&
 };
 
 
-
 template <typename number>
 template <typename somenumber>
 void
@@ -303,7 +296,6 @@ SparseMatrix<number>::threaded_vmult (Vector<somenumber>       &dst,
 };
 
 
-
 template <typename number>
 template <typename somenumber>
 void
@@ -327,7 +319,6 @@ SparseMatrix<number>::Tvmult (Vector<somenumber>& dst, const Vector<somenumber>&
 }
 
 
-
 template <typename number>
 template <typename somenumber>
 void
@@ -374,7 +365,6 @@ SparseMatrix<number>::Tvmult_add (Vector<somenumber>& dst, const Vector<somenumb
 }
 
 
-
 template <typename number>
 template <typename somenumber>
 somenumber
@@ -453,7 +443,6 @@ SparseMatrix<number>::matrix_norm (const Vector<somenumber>& v) const
 };
 
 
-
 template <typename number>
 template <typename somenumber>
 void
@@ -639,7 +628,6 @@ SparseMatrix<number>::threaded_residual (Vector<somenumber>       &dst,
 };
 
 
-
 template <typename number>
 template <typename somenumber>
 void
@@ -665,7 +653,6 @@ SparseMatrix<number>::precondition_Jacobi (Vector<somenumber>& dst,
 };
 
 
-
 template <typename number>
 template <typename somenumber>
 void
@@ -733,7 +720,6 @@ SparseMatrix<number>::precondition_SSOR (Vector<somenumber>& dst,
 }
 
 
-
 template <typename number>
 template <typename somenumber>
 void
@@ -749,7 +735,6 @@ SparseMatrix<number>::precondition_SOR (Vector<somenumber>& dst, const Vector<so
 };
 
 
-
 template <typename number>
 template <typename somenumber>
 void
@@ -877,7 +862,6 @@ SparseMatrix<number>::SSOR (Vector<somenumber>& dst, const number om) const
 }
 
 
-
 template <typename number>
 const SparsityPattern &
 SparseMatrix<number>::get_sparsity_pattern () const
@@ -887,7 +871,6 @@ SparseMatrix<number>::get_sparsity_pattern () const
 };
 
 
-
 template <typename number>
 void SparseMatrix<number>::print (ostream &out) const {
   Assert (cols != 0, ExcMatrixNotInitialized());
@@ -901,7 +884,6 @@ void SparseMatrix<number>::print (ostream &out) const {
 };
 
 
-
 template <typename number>
 void SparseMatrix<number>::print_formatted (ostream &out,
                                            const unsigned int precision,
@@ -941,3 +923,4 @@ void SparseMatrix<number>::print_formatted (ostream &out,
 //  out.setf (0, ios::floatfield);                 // reset output format
 };
 
+#endif
index 76b0d1ad813ebd162a3f2cd314150309f3606b40..328d2b43e0e7e76d96e84d2805c6803ebc1954b8 100644 (file)
@@ -1,10 +1,20 @@
-/*----------------------------   sparse_vanka.h     ---------------------------*/
-/*      $Id$                 */
-#ifndef __sparse_vanka_H
-#define __sparse_vanka_H
-/* Copyright Guido Kanschat, 1999 */
-/*----------------------------   sparse_vanka.h     ---------------------------*/
+//----------------------------  sparse_vanka.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  sparse_vanka.h  ---------------------------
+#ifndef __deal2__sparse_vanka_h
+#define __deal2__sparse_vanka_h
+
 
+/* Copyright Guido Kanschat, 1999 */
 
 
 #include <base/smartpointer.h>
@@ -310,7 +320,6 @@ class SparseVanka
 };
 
 
-
 /**
  * Block version of the sparse Vanka preconditioner. This class
  * divides the matrix into blocks and works on the diagonal blocks
@@ -525,10 +534,4 @@ class SparseBlockVanka : public SparseVanka<number>
 };
 
 
-
-
-
-/*----------------------------   sparse_vanka.h     ---------------------------*/
-/* end of #ifndef __sparse_vanka_H */
 #endif
-/*----------------------------   sparse_vanka.h     ---------------------------*/
index fa2be95f8b3a1fd286212cfb24c5dfe2f647b234..013a9e9e5249e29c10e2ecd452d0c308b9bec20a 100644 (file)
@@ -1,5 +1,18 @@
-// $Id$
-// Copyright Guido Kanschat, 1999
+//----------------------------  sparse_vanka.templates.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  sparse_vanka.templates.h  ---------------------------
+#ifndef __deal2__sparse_vanka_templates_h
+#define __deal2__sparse_vanka_templates_h
+
 
 #include <lac/sparse_vanka.h>
 #include <lac/full_matrix.h>
@@ -14,7 +27,6 @@
 #include <map>
 
 
-
 template<typename number>
 SparseVanka<number>::SparseVanka(const SparseMatrix<number> &M,
                                 const vector<bool>         &selected,
@@ -48,7 +60,6 @@ SparseVanka<number>::~SparseVanka()
 }
 
 
-
 template <typename number>
 void
 SparseVanka<number>::compute_inverses () 
@@ -113,7 +124,6 @@ SparseVanka<number>::compute_inverses ()
 };
 
 
-
 template <typename number>
 void
 SparseVanka<number>::compute_inverses (const unsigned int begin,
@@ -132,7 +142,6 @@ SparseVanka<number>::compute_inverses (const unsigned int begin,
 };
 
 
-
 template <typename number>
 void
 SparseVanka<number>::compute_inverse (const unsigned int               row,
@@ -205,8 +214,6 @@ SparseVanka<number>::compute_inverse (const unsigned int               row,
 };
 
 
-
-
 template<typename number>
 template<typename number2>
 void
@@ -221,8 +228,6 @@ SparseVanka<number>::operator ()(Vector<number2>       &dst,
 };
 
 
-
-
 template<typename number>
 template<typename number2>
 void
@@ -242,7 +247,7 @@ SparseVanka<number>::apply_preconditioner (Vector<number2>       &dst,
     = matrix->get_sparsity_pattern();
 
 
-                                  // store whether we shall work on
+// store whether we shall work on
                                   // the whole matrix, or only on
                                   // blocks. this variable is used to
                                   // optimize access to vectors a
@@ -383,8 +388,6 @@ SparseVanka<number>::apply_preconditioner (Vector<number2>       &dst,
 };
 
 
-
-
 template <typename number>
 SparseBlockVanka<number>::SparseBlockVanka (const SparseMatrix<number> &M,
                                            const vector<bool>         &selected,
@@ -402,7 +405,6 @@ SparseBlockVanka<number>::SparseBlockVanka (const SparseMatrix<number> &M,
 };
 
 
-
 template <typename number>
 void
 SparseBlockVanka<number>::compute_dof_masks (const SparseMatrix<number> &M,
@@ -567,8 +569,6 @@ SparseBlockVanka<number>::compute_dof_masks (const SparseMatrix<number> &M,
 };
 
 
-
-
 template <typename number>
 template <typename number2>
 void SparseBlockVanka<number>::operator() (Vector<number2>       &dst,
@@ -608,5 +608,7 @@ void SparseBlockVanka<number>::operator() (Vector<number2>       &dst,
        apply_preconditioner (dst, src,
                              &dof_masks[block]);
 #endif
-    };
-};
+    }
+}
+
+#endif
index 992a1e48c4f159e0a5ae7d534754e06370b2cd7c..29bb0bdfc1f31ff4c73acdcfc98879c46083f634 100644 (file)
@@ -1,9 +1,17 @@
-/*----------------------------   sparsity_pattern.h     ---------------------------*/
-/*      $Id$                 */
-#ifndef __sparsity_pattern_H
-#define __sparsity_pattern_H
-/*----------------------------   sparsity_pattern.h     ---------------------------*/
-
+//----------------------------  sparsity_pattern.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  sparsity_pattern.h  ---------------------------
+#ifndef __deal2__sparsity_pattern_h
+#define __deal2__sparsity_pattern_h
 
 
 #include <base/exceptions.h>
@@ -13,8 +21,6 @@
 #include <vector>
 
 
-
-
 /**
  * Structure representing the sparsity pattern of a sparse matrix.
  * 
@@ -491,9 +497,9 @@ class SparsityPattern : public Subscriptor
                                      * information!
                                      */
     const unsigned int * get_column_numbers () const;
-    
-    
-                                    /**
+
+
+/**
                                      * Exception
                                      */
     DeclException1 (ExcInvalidNumber,
@@ -670,17 +676,14 @@ class SparsityPattern : public Subscriptor
                                      */
     bool compressed;
 
-    
-    template <typename number> friend class SparseMatrix;
-};
-
 
+template <typename number> friend class SparseMatrix;
+};
 
 
 /*---------------------- Inline functions -----------------------------------*/
 
 
-
 inline
 unsigned int
 SparsityPattern::n_rows () const
@@ -689,7 +692,6 @@ SparsityPattern::n_rows () const
 };
 
 
-
 inline
 unsigned int
 SparsityPattern::n_cols () const
@@ -698,7 +700,6 @@ SparsityPattern::n_cols () const
 };
 
 
-
 inline
 bool
 SparsityPattern::is_compressed () const
@@ -707,7 +708,6 @@ SparsityPattern::is_compressed () const
 };
 
 
-
 inline
 const unsigned int *
 SparsityPattern::get_rowstart_indices () const
@@ -716,7 +716,6 @@ SparsityPattern::get_rowstart_indices () const
 };
 
 
-
 inline
 const unsigned int *
 SparsityPattern::get_column_numbers () const
@@ -725,7 +724,6 @@ SparsityPattern::get_column_numbers () const
 };
 
 
-
 inline
 unsigned int
 SparsityPattern::row_length (const unsigned int row) const
@@ -735,7 +733,6 @@ SparsityPattern::row_length (const unsigned int row) const
 }
 
 
-
 inline
 unsigned int
 SparsityPattern::column_number (const unsigned int row,
@@ -748,7 +745,6 @@ SparsityPattern::column_number (const unsigned int row,
 }
 
 
-
 inline
 unsigned int
 SparsityPattern::n_nonzero_elements () const
@@ -759,8 +755,4 @@ SparsityPattern::n_nonzero_elements () const
 };
 
 
-
-/*----------------------------   sparsity_pattern.h     ---------------------------*/
-/* end of #ifndef __sparsity_pattern_H */
 #endif
-/*----------------------------   sparsity_pattern.h     ---------------------------*/
index b20d240a086fbfd8856b3d0d9fa35bcd2d060dfc..b5017424a15b6659c02eb62c83f5401a9bff1a0e 100644 (file)
@@ -1,19 +1,23 @@
-/*----------------------------   vector.h     ---------------------------*/
-/*      $Id$                 */
-#ifndef __vector_H
-#define __vector_H
-/*----------------------------   vector.h     ---------------------------*/
+//----------------------------  vector.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  vector.h  ---------------------------
+#ifndef __deal2__vector_h
+#define __deal2__vector_h
 
-// This file was once part of the DEAL Library
-// DEAL is Copyright(1995) by
-// Roland Becker, Guido Kanschat, Franz-Theo Suttmeier
-// Revised and extended by Wolfgang Bangerth
 
 #include <base/exceptions.h>
 #include <cstdio>
 
 
-
 /**
  *  Numerical Vector of data.  For this class there are different
  * types of functions available. The first type of function mesures
@@ -43,9 +47,9 @@ class Vector {
     typedef value_type& reference;
     typedef const value_type& const_reference;
     typedef size_t size_type;
-    
 
-                                    /**
+
+/**
                                      * @name 1: Basic Object-handling 
                                      */
                                     //@{
@@ -153,9 +157,9 @@ class Vector {
                                      * $l_\infty$-norm of a vector.
                                      */
     Number linfty_norm () const;
-    
-    
-                                    /**
+
+
+/**
                                      * Change the dimension of the vector to
                                      * #N#. The reserved memory for this vector
                                      * remains unchanged if possible, to make
@@ -233,9 +237,9 @@ class Vector {
                                      */
     const_iterator end () const;  
                                     //@}
-    
-    
-                                    /**
+
+
+/**
                                      * @name 2: Data-Access
                                      */
                                     //@{
@@ -250,9 +254,9 @@ class Vector {
                                      */
     Number& operator() (const unsigned int i);
                                     //@}
-    
-    
-                                    /**
+
+
+/**
                                      * @name 3: Modification of vectors
                                      */
                                     //@{
@@ -358,9 +362,9 @@ class Vector {
     void ratio (const Vector<Number> &a,
                const Vector<Number> &b);
                                     //@}
-    
-    
-                                    /**
+
+
+/**
                                      * @name 5: Mixed stuff
                                      */
                                     //@{
@@ -463,10 +467,6 @@ class Vector {
 };
 
 
-
-
-
-
 /*----------------------- Inline functions ----------------------------------*/
 
 
@@ -478,7 +478,6 @@ unsigned int Vector<Number>::size () const
 }
 
 
-
 template <typename Number>
 inline
 Vector<Number>::iterator Vector<Number>::begin () {
@@ -486,7 +485,6 @@ Vector<Number>::iterator Vector<Number>::begin () {
 };
 
 
-
 template <typename Number>
 inline
 Vector<Number>::const_iterator Vector<Number>::begin () const {
@@ -494,7 +492,6 @@ Vector<Number>::const_iterator Vector<Number>::begin () const {
 };
 
 
-
 template <typename Number>
 inline
 Vector<Number>::iterator Vector<Number>::end () {
@@ -502,7 +499,6 @@ Vector<Number>::iterator Vector<Number>::end () {
 };
 
 
-
 template <typename Number>
 inline
 Vector<Number>::const_iterator Vector<Number>::end () const {
@@ -510,7 +506,6 @@ Vector<Number>::const_iterator Vector<Number>::end () const {
 };
 
 
-
 template <typename Number>
 inline
 Number Vector<Number>::operator() (const unsigned int i) const
@@ -520,7 +515,6 @@ Number Vector<Number>::operator() (const unsigned int i) const
 }
 
 
-
 template <typename Number>
 inline
 Number& Vector<Number>::operator() (const unsigned int i)
@@ -530,10 +524,4 @@ Number& Vector<Number>::operator() (const unsigned int i)
 }
 
 
-
-
-
-/*----------------------------   vector.h     ---------------------------*/
-/* end of #ifndef __vector_H */
 #endif
-/*----------------------------   vector.h     ---------------------------*/
index d1855941f34f2a78330db4f3cd411192fcbb919b..5d79eb4257131de1cd02742d39a544e05c70920c 100644 (file)
@@ -1,9 +1,18 @@
-// $Id$
+//----------------------------  vector.templates.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  vector.templates.h  ---------------------------
+#ifndef __deal2__vector_templates_h
+#define __deal2__vector_templates_h
 
-// This file was once part of the DEAL Library
-// DEAL is Copyright(1995) by
-// Roland Becker, Guido Kanschat, Franz-Theo Suttmeier
-// Revised and extended by Wolfgang Bangerth, 1998, 1999
 
 #include <lac/vector.h>
 #include <cmath>
@@ -16,7 +25,6 @@ static inline Number sqr (const Number x) {
 };
 
 
-
 template <typename Number>
 Vector<Number>::Vector () :
                dim(0),
@@ -25,7 +33,6 @@ Vector<Number>::Vector () :
 {}
 
 
-
 template <typename Number>
 Vector<Number>::Vector (const unsigned int n) :
                dim(0),
@@ -36,7 +43,6 @@ Vector<Number>::Vector (const unsigned int n) :
 }
 
 
-
 template <typename Number>
 Vector<Number>::Vector (const Vector<Number>& v) :
                dim(v.size()),
@@ -52,7 +58,6 @@ Vector<Number>::Vector (const Vector<Number>& v) :
 }
 
 
-
 // see the .h file for why this function was disabled
 //
 // template <typename Number>
@@ -71,7 +76,6 @@ Vector<Number>::Vector (const Vector<Number>& v) :
 // }
 
 
-
 template <typename Number>
 void Vector<Number>::reinit (const unsigned int n, const bool fast) {
   if (n==0) 
@@ -95,15 +99,12 @@ void Vector<Number>::reinit (const unsigned int n, const bool fast) {
 }
 
 
-
 template <typename Number>
 void Vector<Number>::reinit (const Vector<Number>& v, const bool fast) {
   reinit (v.size(), fast);
 };
 
 
-
-
 template <typename Number>
 Vector<Number>::~Vector ()
 {
@@ -111,7 +112,6 @@ Vector<Number>::~Vector ()
 }
 
 
-
 template <typename Number>
 void Vector<Number>::clear () {
   if (dim>0)
@@ -119,7 +119,6 @@ void Vector<Number>::clear () {
 }
 
 
-
 template <typename Number>
 bool Vector<Number>::all_zero () const {
   Assert (dim!=0, ExcEmptyVector());
@@ -133,7 +132,6 @@ bool Vector<Number>::all_zero () const {
 };
 
 
-
 template <typename Number>
 Number Vector<Number>::operator * (const Vector<Number>& v) const
 {
@@ -170,7 +168,6 @@ Number Vector<Number>::operator * (const Vector<Number>& v) const
 }
 
 
-
 template <typename Number>
 Number Vector<Number>::norm_sqr () const
 {
@@ -201,7 +198,6 @@ Number Vector<Number>::norm_sqr () const
 };
 
 
-
 template <typename Number>
 Number Vector<Number>::mean_value () const
 {
@@ -232,7 +228,6 @@ Number Vector<Number>::mean_value () const
 };
 
 
-
 template <typename Number>
 Number Vector<Number>::l1_norm () const
 {
@@ -263,7 +258,6 @@ Number Vector<Number>::l1_norm () const
 };
 
 
-
 template <typename Number>
 Number Vector<Number>::l2_norm () const
 {
@@ -271,7 +265,6 @@ Number Vector<Number>::l2_norm () const
 };
 
 
-
 template <typename Number>
 Number Vector<Number>::linfty_norm () const {
   Assert (dim!=0, ExcEmptyVector());
@@ -295,9 +288,6 @@ Number Vector<Number>::linfty_norm () const {
   return max (max(max0, max1),
              max(max2, max3));
 };
-  
-
-
 
 
 template <typename Number>
@@ -310,7 +300,6 @@ Vector<Number>& Vector<Number>::operator += (const Vector<Number>& v)
 }
 
 
-
 template <typename Number>
 Vector<Number>& Vector<Number>::operator -= (const Vector<Number>& v)
 {
@@ -327,7 +316,6 @@ Vector<Number>& Vector<Number>::operator -= (const Vector<Number>& v)
 }
 
 
-
 template <typename Number>
 void Vector<Number>::add (const Number v)
 {
@@ -340,7 +328,6 @@ void Vector<Number>::add (const Number v)
 }
 
 
-
 template <typename Number>
 void Vector<Number>::add (const Vector<Number>& v)
 {
@@ -355,7 +342,6 @@ void Vector<Number>::add (const Vector<Number>& v)
 }
 
 
-
 template <typename Number>
 void Vector<Number>::add (const Number a, const Vector<Number>& v)
 {
@@ -370,7 +356,6 @@ void Vector<Number>::add (const Number a, const Vector<Number>& v)
 }
 
 
-
 template <typename Number>
 void Vector<Number>::add (const Number a, const Vector<Number>& v,
                          const Number b, const Vector<Number>& w)
@@ -387,7 +372,6 @@ void Vector<Number>::add (const Number a, const Vector<Number>& v,
 }
 
 
-
 template <typename Number>
 void Vector<Number>::sadd (const Number x, const Vector<Number>& v)
 {
@@ -401,7 +385,6 @@ void Vector<Number>::sadd (const Number x, const Vector<Number>& v)
 }
 
 
-
 template <typename Number>
 void Vector<Number>::sadd (const Number x, const Number a, const Vector<Number>& v)
 {
@@ -415,7 +398,6 @@ void Vector<Number>::sadd (const Number x, const Number a, const Vector<Number>&
 }
 
 
-
 template <typename Number>
 void Vector<Number>::sadd (const Number x, const Number a,
                           const Vector<Number>& v, const Number b, const Vector<Number>& w)
@@ -432,7 +414,6 @@ void Vector<Number>::sadd (const Number x, const Number a,
 }
 
 
-
 template <typename Number>
 void Vector<Number>::sadd (const Number x, const Number a,
                           const Vector<Number>& v, const Number b,
@@ -453,7 +434,6 @@ void Vector<Number>::sadd (const Number x, const Number a,
 }
 
 
-
 template <typename Number>
 void Vector<Number>::scale (const Number factor)
 {
@@ -465,7 +445,6 @@ void Vector<Number>::scale (const Number factor)
 }
 
 
-
 template <typename Number>
 void Vector<Number>::equ (const Number a, const Vector<Number>& u,
                          const Number b, const Vector<Number>& v)
@@ -482,7 +461,6 @@ void Vector<Number>::equ (const Number a, const Vector<Number>& u,
 }
 
 
-
 template <typename Number>
 void Vector<Number>::equ (const Number a, const Vector<Number>& u)
 {
@@ -496,7 +474,6 @@ void Vector<Number>::equ (const Number a, const Vector<Number>& u)
 }
 
 
-
 template <typename Number>
 void Vector<Number>::ratio (const Vector<Number> &a, const Vector<Number> &b) {
   Assert (dim!=0, ExcEmptyVector());
@@ -514,7 +491,6 @@ void Vector<Number>::ratio (const Vector<Number> &a, const Vector<Number> &b) {
 };
 
 
-
 template <typename Number>
 Vector<Number>& Vector<Number>::operator = (const Number s)
 {
@@ -524,7 +500,6 @@ Vector<Number>& Vector<Number>::operator = (const Number s)
 }
 
 
-
 template <typename Number>
 Vector<Number>&
 Vector<Number>::operator = (const Vector<Number>& v)
@@ -538,7 +513,6 @@ Vector<Number>::operator = (const Vector<Number>& v)
 }
 
 
-
 template <typename Number>
 template<typename Number2>
 Vector<Number>&
@@ -553,7 +527,6 @@ Vector<Number>::operator = (const Vector<Number2>& v)
 }
 
 
-
 template <typename Number>
 void Vector<Number>::print (FILE* f, const char* format) const
 {
@@ -565,7 +538,6 @@ void Vector<Number>::print (FILE* f, const char* format) const
 }
 
 
-
 template <typename Number>
 void Vector<Number>::print (const char* format) const
 {
@@ -577,7 +549,6 @@ void Vector<Number>::print (const char* format) const
 }
 
 
-
 template <typename Number>
 void Vector<Number>::print (ostream &out,
                            unsigned int precision,
@@ -606,7 +577,6 @@ void Vector<Number>::print (ostream &out,
 };
 
 
-
 template <typename Number>
 void Vector<Number>::block_write (ostream &out) const {
   AssertThrow (out, ExcIO());
@@ -621,7 +591,6 @@ void Vector<Number>::block_write (ostream &out) const {
 };
 
 
-
 template <typename Number>
 void Vector<Number>::block_read (istream &in) {
   AssertThrow (in, ExcIO());
@@ -643,5 +612,6 @@ void Vector<Number>::block_read (istream &in) {
   in >> c;
   AssertThrow (c==']', ExcIO());
   AssertThrow (in, ExcIO());
-};
+}
 
+#endif
index aab9e728c5a9c777665b3464bbea251312e81e9c..c0960b7ca20e1a24c3c6939ac4474ce98cb266c5 100644 (file)
@@ -1,8 +1,17 @@
-/*----------------------------   vector_memory.h     ---------------------------*/
-/*      $Id$                 */
-#ifndef __lac__vector_memory_H
-#define __lac__vector_memory_H
-/*----------------------------   vector_memory.h     ---------------------------*/
+//----------------------------  vector_memory.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  vector_memory.h  ---------------------------
+#ifndef __deal2__vector_memory_h
+#define __deal2__vector_memory_h
 
 
 #include <base/subscriptor.h>
@@ -11,7 +20,6 @@
 #include <vector>
 
 
-
 /**
  * Memory management for vectors. This class is used by all
  * iterative methods to allocate space for auxilliary
@@ -58,9 +66,6 @@ class VectorMemory : public Subscriptor
 };
 
 
-
-
-
 /**
  * Simple memory management.  This memory class is just made for
  * tests. It just allocates and deletes
@@ -160,7 +165,6 @@ class GrowingVectorMemory : public VectorMemory<Vector>
 };
 
 
-
 /* --------------------- inline functions ---------------------- */
 
 
@@ -178,7 +182,6 @@ GrowingVectorMemory<Vector>::GrowingVectorMemory(const unsigned int initial_size
 }
 
 
-
 template<typename Vector>
 GrowingVectorMemory<Vector>::~GrowingVectorMemory()
 {
@@ -199,7 +202,6 @@ GrowingVectorMemory<Vector>::~GrowingVectorMemory()
 }
 
 
-
 template<typename Vector>
 Vector*
 GrowingVectorMemory<Vector>::alloc()
@@ -220,8 +222,7 @@ GrowingVectorMemory<Vector>::alloc()
   return t.second;
 }
 
-  
-  
+
 template<typename Vector>
 void
 GrowingVectorMemory<Vector>::free(const Vector* const v)
@@ -238,8 +239,5 @@ GrowingVectorMemory<Vector>::free(const Vector* const v)
 }
 
 
-/*----------------------------   vector_memory.h     ---------------------------*/
-/* end of #ifndef __vector_memory_H */
 #endif
-/*----------------------------   vector_memory.h     ---------------------------*/
 
index 1818ba60d0da90af999766c6f240e0c52f0af05f..4af35357f1ec7fd3fca3d867cb1b8e06aa98bd0c 100644 (file)
@@ -1,14 +1,16 @@
-// $Id$
+//----------------------------  full_matrix.double.cc  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  full_matrix.double.cc  ---------------------------
 
-// Driver for FullMatrix template instantiation.
-
-/* Instantiation is controlled by preprocessor symbols:
- *
- * 1. TYPEMAT : numerical type used in the matrix
- * 2. TYPEVEC : numerical type for vector entries
- * 3. TYPERES : numerical type for entries in the right hand side vector
- * 4. TYPEMAT2: numerical type for the second matrix
- */
 
 #include <lac/full_matrix.templates.h>
 
index 1871cb4d528d83468d414a7685acc7adb3e8fae1..ec174e7f2684ac21fb034b5ec060de590142b14f 100644 (file)
@@ -1,14 +1,16 @@
-// $Id$
+//----------------------------  full_matrix.float.cc  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  full_matrix.float.cc  ---------------------------
 
-// Driver for FullMatrix template instantiation.
-
-/* Instantiation is controlled by preprocessor symbols:
- *
- * 1. TYPEMAT : numerical type used in the matrix
- * 2. TYPEVEC : numerical type for vector entries
- * 3. TYPERES : numerical type for entries in the right hand side vector
- * 4. TYPEMAT2: numerical type for the second matrix
- */
 
 #include <lac/full_matrix.templates.h>
 
index 528dd90c024b80aaaa90acd750a04edd73f7390a..823b11b63827f221775256547e8edbd25f689628 100644 (file)
@@ -1,14 +1,20 @@
-/*----------------------------   precondition_block.cc     ---------------------------*/
-/*      $Id$                 */
-/*                Ralf Hartmann, University of Heidelberg                            */
-
+//----------------------------  precondition_block.cc  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  precondition_block.cc  ---------------------------
 
 
 #include <lac/precondition_block.templates.h>
 
 
-
-
 // explicit instantiations for "float" PreconditionBlock
 template class PreconditionBlock<float, float>;
 
@@ -40,7 +46,6 @@ template void PreconditionBlockSOR<float, float>::operator() (
 // higher precision than the matrix itself
 
 
-
 // explicit instantiations for "double" PreconditionBlockSOR
 template class PreconditionBlockSOR<double, float>;
 
@@ -58,8 +63,3 @@ template void PreconditionBlockSOR<double, double>::operator() (
   Vector<double> &, const Vector<double> &) const;
 
 
-
-
-
-
-/*----------------------------   precondition_block.cc     ---------------------------*/
index 1ef8150cba9d1106db151fbdc7849aba096b1851..d2577d40423d10fb16a336fa1120cb23150d0913 100644 (file)
@@ -1,4 +1,16 @@
-// $Id$
+//----------------------------  solver_control.cc  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  solver_control.cc  ---------------------------
+
 
 #include <base/logstream.h>
 #include <lac/solver_control.h>
@@ -39,9 +51,9 @@ SolverControl::check (const unsigned int step,
     {
       deallog << "Starting value " << check_value << endl;
     }
-  
 
-  if (step>=maxsteps)
+
+if (step>=maxsteps)
     {
       if (log_result)
        deallog << "Failure step " << step
@@ -59,7 +71,6 @@ SolverControl::check (const unsigned int step,
 };
 
 
-
 double
 SolverControl::last_value() const
 {
@@ -67,7 +78,6 @@ SolverControl::last_value() const
 };
 
 
-
 unsigned int
 SolverControl::last_step() const
 {
@@ -75,12 +85,9 @@ SolverControl::last_step() const
 };
 
 
-
-
 /*----------------------- ReductionControl ---------------------------------*/
 
 
-
 ReductionControl::ReductionControl(const unsigned int n,
                                   const double tol,
                                   const double red,
@@ -100,7 +107,6 @@ double
 ReductionControl::initial_value() const {
   return initial_val;
 };
-    
 
 
 SolverControl::State
index 7d48252d6b65f9e40c48c29b407f97ab30281988..3aa33b6997ce8962f9b80d3f363404a7297e749f 100644 (file)
@@ -1,4 +1,16 @@
-/* $Id$ */
+//----------------------------  sparse_ilu.cc  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  sparse_ilu.cc  ---------------------------
+
 
 #include <lac/sparse_ilu.templates.h>
 
index 5768bb4f3786f0dc927c9a571cbe75ff0cf9793c..69ca4421c0534bd93db7f77ed1bddb7287260759 100644 (file)
@@ -1,13 +1,16 @@
-// $Id$
+//----------------------------  sparse_matrix.double.cc  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  sparse_matrix.double.cc  ---------------------------
 
-// SparseMatrix template instantiation
-
-
-/* Instantiation is controlled by preprocessor symbols:
- *
- * 1. TYPEMAT : numerical type used in the matrix
- * 2. TYPE2 : numerical type for function arguments
- */
 
 #include <cmath>
 #include <lac/sparse_matrix.templates.h>
index 7a8a327aeb60d79704145731f6d0d7892cdf6c65..b2395269120c4ac9b2e623ccad51e4228280887a 100644 (file)
@@ -1,13 +1,16 @@
-// $Id$
+//----------------------------  sparse_matrix.float.cc  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  sparse_matrix.float.cc  ---------------------------
 
-// SparseMatrix template instantiation
-
-
-/* Instantiation is controlled by preprocessor symbols:
- *
- * 1. TYPEMAT : numerical type used in the matrix
- * 2. TYPE2 : numerical type for function arguments
- */
 
 #include <cmath>
 #include <lac/sparse_matrix.templates.h>
index 62f326ca38b3b2de9bb5b204aa513a4b13ed0933..50a8cadfb38cdde8e4933e334878caa6599c00dc 100644 (file)
@@ -1,13 +1,16 @@
-// $Id$
+//----------------------------  sparse_matrix.long_double.cc  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  sparse_matrix.long_double.cc  ---------------------------
 
-// SparseMatrix template instantiation
-
-
-/* Instantiation is controlled by preprocessor symbols:
- *
- * 1. TYPEMAT : numerical type used in the matrix
- * 2. TYPE2 : numerical type for function arguments
- */
 
 #include <cmath>
 #include <lac/sparse_matrix.templates.h>
index 687577c6f12227b9dc8f49a3eecf65b088b2b06e..b23b455b4dece7ecb68f309c4905858a7ca9ca76 100644 (file)
@@ -1,4 +1,16 @@
-/* $Id$ */
+//----------------------------  sparse_vanka.cc  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  sparse_vanka.cc  ---------------------------
+
 
 #include <lac/sparse_vanka.templates.h>
 
index 6c6bf1689ff1b5b96715619a9baa6ac1d501cccf..e010645ef536c15f551ea7432eb3b8c8f4424ad3 100644 (file)
@@ -1,9 +1,15 @@
-// $Id$
-
-// This file was once part of the DEAL Library
-// DEAL is Copyright(1995) by
-// Roland Becker, Guido Kanschat, Franz-Theo Suttmeier
-// Revised, modified and extended by Wolfgang Bangerth, 1998, 1999
+//----------------------------  sparsity_pattern.cc  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  sparsity_pattern.cc  ---------------------------
 
 
 #include <lac/sparse_matrix.h>
@@ -25,7 +31,6 @@ SparsityPattern::SparsityPattern () :
 };
 
 
-
 SparsityPattern::SparsityPattern (const SparsityPattern &s) :
                Subscriptor(),
                max_dim(0),
@@ -42,7 +47,6 @@ SparsityPattern::SparsityPattern (const SparsityPattern &s) :
 };
 
 
-
 SparsityPattern::SparsityPattern (const unsigned int m,
                                  const unsigned int n,
                                  const unsigned int max_per_row) 
@@ -55,7 +59,6 @@ SparsityPattern::SparsityPattern (const unsigned int m,
 };
 
 
-
 SparsityPattern::SparsityPattern (const unsigned int          m,
                                  const unsigned int          n,
                                  const vector<unsigned int> &row_lengths) 
@@ -68,7 +71,6 @@ SparsityPattern::SparsityPattern (const unsigned int          m,
 };
 
 
-
 SparsityPattern::SparsityPattern (const unsigned int n,
                                  const unsigned int max_per_row)
                : max_dim(0),
@@ -80,7 +82,6 @@ SparsityPattern::SparsityPattern (const unsigned int n,
 };
 
 
-
 SparsityPattern::SparsityPattern (const unsigned int          m,
                                  const vector<unsigned int> &row_lengths) 
                : max_dim(0),
@@ -92,7 +93,6 @@ SparsityPattern::SparsityPattern (const unsigned int          m,
 };
 
 
-
 SparsityPattern::SparsityPattern (const SparsityPattern &original,
                                  const unsigned int        max_per_row,
                                  const unsigned int        extra_off_diagonals)
@@ -176,7 +176,6 @@ SparsityPattern::SparsityPattern (const SparsityPattern &original,
 };
 
 
-
 SparsityPattern::~SparsityPattern ()
 {
   if (rowstart != 0)  delete[] rowstart;
@@ -184,7 +183,6 @@ SparsityPattern::~SparsityPattern ()
 }
 
 
-
 SparsityPattern &
 SparsityPattern::operator = (const SparsityPattern &s)
 {
@@ -202,8 +200,6 @@ SparsityPattern::operator = (const SparsityPattern &s)
 };
 
 
-
-
 void
 SparsityPattern::reinit (const unsigned int m,
                         const unsigned int n,
@@ -216,7 +212,6 @@ SparsityPattern::reinit (const unsigned int m,
 };
 
 
-
 void
 SparsityPattern::reinit (const unsigned int m,
                         const unsigned int n,
@@ -256,7 +251,7 @@ SparsityPattern::reinit (const unsigned int m,
                    *max_element(row_lengths.begin(), row_lengths.end()));
 
 
-                                  // allocate memory for the rowstart
+// allocate memory for the rowstart
                                   // values, if necessary
   if (rows > max_dim)
     {
@@ -320,8 +315,8 @@ SparsityPattern::compress ()
                                   // now allocate the respective memory
   unsigned int *new_colnums = new unsigned int[nonzero_elements];
 
-  
-                                  // reserve temporary storage to
+
+// reserve temporary storage to
                                   // store the entries of one row
   vector<unsigned int> tmp_entries (max_row_length);
   
@@ -403,7 +398,6 @@ SparsityPattern::compress ()
 };
 
 
-
 bool
 SparsityPattern::empty () const {
                                   // let's try to be on the safe side of
@@ -429,7 +423,6 @@ SparsityPattern::empty () const {
 };
 
 
-
 unsigned int
 SparsityPattern::max_entries_per_row () const 
 {
@@ -451,7 +444,6 @@ SparsityPattern::max_entries_per_row () const
 };
 
 
-
 unsigned int
 SparsityPattern::operator () (const unsigned int i,
                              const unsigned int j) const
@@ -496,7 +488,6 @@ SparsityPattern::operator () (const unsigned int i,
 }
 
 
-
 void
 SparsityPattern::add (const unsigned int i,
                      const unsigned int j)
@@ -526,7 +517,6 @@ SparsityPattern::add (const unsigned int i,
 }
 
 
-
 void
 SparsityPattern::print_gnuplot (ostream &out) const
 {
@@ -540,7 +530,6 @@ SparsityPattern::print_gnuplot (ostream &out) const
 }
 
 
-
 unsigned int
 SparsityPattern::bandwidth () const
 {
index ce8f4baeb597d6d2dda620c4fb192eac9edf2661..79bc6a33bb20665139128db2b5c9998025bb4ea5 100644 (file)
@@ -1,4 +1,16 @@
-// $Id$
+//----------------------------  vector.cc  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  vector.cc  ---------------------------
+
 
 #include <lac/vector.templates.h>
 
index 2fc8e16f28721690ae744636f1f5886cfb9d6d4f..7e49c8ce2676f1cb225d814c4cf0ceb7cc68068a 100644 (file)
@@ -1,4 +1,16 @@
-// $Id$
+//----------------------------  vector.long_double.cc  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  vector.long_double.cc  ---------------------------
+
 
 #include <lac/vector.templates.h>
 

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.