]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Lots, lots, lots of small reindentations, etc (several thousands of them...).
authorwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 9 Jun 2003 17:01:06 +0000 (17:01 +0000)
committerwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 9 Jun 2003 17:01:06 +0000 (17:01 +0000)
git-svn-id: https://svn.dealii.org/trunk@7776 0785d39b-7218-0410-832d-ea1e28bc413d

66 files changed:
deal.II/base/include/base/auto_derivative_function.h
deal.II/base/include/base/convergence_table.h
deal.II/base/include/base/data_out_base.h
deal.II/base/include/base/exceptions.h
deal.II/base/include/base/function.h
deal.II/base/include/base/function_derivative.h
deal.II/base/include/base/function_lib.h
deal.II/base/include/base/function_time.h
deal.II/base/include/base/job_identifier.h
deal.II/base/include/base/logstream.h
deal.II/base/include/base/memory_consumption.h
deal.II/base/include/base/multithread_info.h
deal.II/base/include/base/parameter_handler.h
deal.II/base/include/base/point.h
deal.II/base/include/base/quadrature.h
deal.II/base/include/base/smartpointer.h
deal.II/base/include/base/table.h
deal.II/base/include/base/table_handler.h
deal.II/base/include/base/tensor.h
deal.II/base/include/base/tensor_base.h
deal.II/base/include/base/tensor_function.h
deal.II/base/include/base/thread_management.h
deal.II/base/include/base/timer.h
deal.II/base/source/auto_derivative_function.cc
deal.II/base/source/convergence_table.cc
deal.II/base/source/data_out_base.cc
deal.II/base/source/exceptions.cc
deal.II/base/source/function.cc
deal.II/base/source/function_lib.cc
deal.II/base/source/function_lib_cutoff.cc
deal.II/base/source/function_time.cc
deal.II/base/source/log.cc
deal.II/base/source/multithread_info.cc
deal.II/base/source/parameter_handler.cc
deal.II/base/source/polynomial.cc
deal.II/base/source/polynomial_space.cc
deal.II/base/source/quadrature.cc
deal.II/base/source/quadrature_lib.cc
deal.II/base/source/subscriptor.cc
deal.II/base/source/table_handler.cc
deal.II/base/source/tensor.cc
deal.II/base/source/tensor_function.cc
deal.II/base/source/thread_management.cc
deal.II/base/source/timer.cc
deal.II/contrib/boost/include/boost_local/config/compiler/gcc.hpp
deal.II/deal.II/include/grid/tria_accessor.h
deal.II/deal.II/source/grid/tria_accessor.cc
deal.II/examples/step-14/Makefile
deal.II/lac/include/lac/block_indices.h
deal.II/lac/include/lac/block_matrix.h
deal.II/lac/include/lac/block_sparse_matrix.h
deal.II/lac/include/lac/block_sparse_matrix_ez.h
deal.II/lac/include/lac/block_sparsity_pattern.h
deal.II/lac/include/lac/block_vector.h
deal.II/lac/include/lac/compressed_sparsity_pattern.h
deal.II/lac/include/lac/eigen.h
deal.II/lac/include/lac/filtered_matrix.templates.h
deal.II/lac/include/lac/full_matrix.h
deal.II/lac/include/lac/full_matrix.templates.h
deal.II/lac/include/lac/matrix_out.h
deal.II/lac/include/lac/sparse_matrix.h
deal.II/lac/include/lac/sparse_matrix.templates.h
deal.II/lac/include/lac/vector.h
deal.II/lac/include/lac/vector.templates.h
deal.II/lac/source/compressed_sparsity_pattern.cc
deal.II/lac/source/sparse_matrix_vector.in.h

index 526ebb7975fb6defa30792f98aff66041d812898..c0d772ab867728c3877bcd30c2662f7fe15b153b 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 2001, 2002 by the deal authors
+//    Copyright (C) 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -223,7 +223,9 @@ class AutoDerivativeFunction : public Function<dim>
                                      * @p{DifferenceFormula} of the
                                      * order @p{ord} at minimum.
                                      */
-    static DifferenceFormula get_formula_of_order(const unsigned int ord);
+    static
+    DifferenceFormula
+    get_formula_of_order (const unsigned int ord);
 
                                     /**
                                      * Exception.
index b237379433b283317e12fc00f104ee3375092135..de038f81f88c5fa8c72f3e058689e77ddd63ce0e 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002 by the deal authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -104,9 +104,10 @@ class ConvergenceTable: public TableHandler
                                      *
                                      * Still not implemented.
                                      */
-    void evaluate_convergence_rates (const std::string &data_column_key,
-                                    const std::string &reference_column_key,
-                                    const RateMode     rate_mode);
+    void
+    evaluate_convergence_rates (const std::string &data_column_key,
+                                const std::string &reference_column_key,
+                                const RateMode     rate_mode);
 
 
                                     /**
@@ -126,8 +127,9 @@ class ConvergenceTable: public TableHandler
                                      * the @p{set_tex_supercaption (...)} function
                                      * of the base class @p{TableHandler}.
                                      */
-    void evaluate_convergence_rates (const std::string &data_column_key, 
-                                    const RateMode     rate_mode);
+    void
+    evaluate_convergence_rates (const std::string &data_column_key, 
+                                const RateMode     rate_mode);
 
                                     /**
                                      * Omit this column @p{key} 
@@ -140,7 +142,8 @@ class ConvergenceTable: public TableHandler
                                      * omitting the column from convergence
                                      * rate evalution.
                                      */
-    void omit_column_from_convergence_rate_evaluation(const std::string &key);
+    void
+    omit_column_from_convergence_rate_evaluation(const std::string &key);
 
                                     /**
                                      * Evaluates convergence rates
@@ -168,8 +171,9 @@ class ConvergenceTable: public TableHandler
                                      * should omitted by calling the
                                      * @p{omit_column_from_convergence_rate_evaluation(..)}.
                                      */
-    void evaluate_all_convergence_rates(const std::string &reference_column_key,
-                                       const RateMode     rate_mode);
+    void
+    evaluate_all_convergence_rates(const std::string &reference_column_key,
+                                   const RateMode     rate_mode);
 
                                     /**
                                      * Evaluates convergence rates
@@ -195,7 +199,8 @@ class ConvergenceTable: public TableHandler
                                      * should omitted by calling the
                                      * @p{omit_column_from_convergence_rate_evaluation(..)}.
                                      */
-    void evaluate_all_convergence_rates(const RateMode rate_mode);
+    void
+    evaluate_all_convergence_rates(const RateMode rate_mode);
 
                                     /**
                                      * Exception
index e9d1c64dfb0aa0bd5a48792b29b1884c77d52c44..aaab631ad4c53fb8d48da62cbe76c674f289002e 100644 (file)
@@ -554,8 +554,7 @@ class DataOutBase
                                          */
        DXFlags (const bool write_multigrid = false,
                 const bool write_neighbors = false);
-
-      public:
+        
                                         /**
                                          * Declare all flags with name
                                          * and type as offered by this
@@ -1011,9 +1010,10 @@ class DataOutBase
                                          * This function was originally written
                                          * by Stefan Nauber.
                                          */
-       static RgbValues default_color_function (const double value,
-                                                const double min_value,
-                                                const double max_value);
+       static RgbValues
+        default_color_function (const double value,
+                                const double min_value,
+                                const double max_value);
 
                                         /**
                                          * This is an alternative color
@@ -1024,9 +1024,10 @@ class DataOutBase
                                          * @p{color_function} variable to the
                                          * address of this function.
                                          */
-       static RgbValues grey_scale_color_function (const double value,
-                                                   const double min_value,
-                                                   const double max_value);
+       static RgbValues
+        grey_scale_color_function (const double value,
+                                   const double min_value,
+                                   const double max_value);
        
                                         /**
                                          * This is one more
@@ -1042,9 +1043,10 @@ class DataOutBase
                                          * variable to the address of
                                          * this function.
                                          */
-       static RgbValues reverse_grey_scale_color_function (const double value,
-                                                           const double min_value,
-                                                           const double max_value);
+       static RgbValues
+        reverse_grey_scale_color_function (const double value,
+                                           const double min_value,
+                                           const double max_value);
        
                                         /**
                                          * Constructor.
@@ -2059,7 +2061,7 @@ class DataOutInterface : private DataOutBase
 
 inline
 bool
-DataOutBase::EpsFlags::RgbValues::is_grey () const 
+DataOutBase::EpsFlags::RgbValues::is_grey () const
 {
   return (red == green) && (red == blue);
 }
index 3c6a20806830f1d77e16797692fa6ae035f7b744..3cc0c21b4b7b1e997761bf360672bb7f13ebaaa2 100644 (file)
@@ -656,7 +656,7 @@ class Exception1 : public ExceptionBase {                             \
 #define DeclException2(Exception2, type1, type2, outsequence)         \
 class Exception2 : public ExceptionBase {                             \
   public:                                                             \
-      Exception2 (const type1 a1, const type2 a2) :                   \
+      Exception2 (const type1 a1, const type2 a2) :          \
              arg1 (a1), arg2(a2) {};                                 \
       virtual ~Exception2 () throw () {};                             \
       virtual void PrintInfo (std::ostream &out) const {              \
@@ -678,7 +678,7 @@ class Exception2 : public ExceptionBase {                             \
 #define DeclException3(Exception3, type1, type2, type3, outsequence)  \
 class Exception3 : public ExceptionBase {                             \
   public:                                                             \
-      Exception3 (const type1 a1, const type2 a2, const type3 a3) :   \
+      Exception3 (const type1 a1, const type2 a2, const type3 a3) : \
              arg1 (a1), arg2(a2), arg3(a3) {};                       \
       virtual ~Exception3 () throw () {};                             \
       virtual void PrintInfo (std::ostream &out) const {              \
@@ -702,7 +702,7 @@ class Exception3 : public ExceptionBase {                             \
 class Exception4 : public ExceptionBase {                             \
   public:                                                             \
       Exception4 (const type1 a1, const type2 a2,                     \
-           const type3 a3, const type4 a4) :                         \
+           const type3 a3, const type4 a4) :                \
              arg1 (a1), arg2(a2), arg3(a3), arg4(a4) {};             \
       virtual ~Exception4 () throw () {};                             \
       virtual void PrintInfo (std::ostream &out) const {              \
@@ -727,7 +727,7 @@ class Exception4 : public ExceptionBase {                             \
 class Exception5 : public ExceptionBase {                             \
   public:                                                             \
       Exception5 (const type1 a1, const type2 a2, const type3 a3,     \
-           const type4 a4, const type5 a5) :                         \
+           const type4 a4, const type5 a5) :                \
              arg1 (a1), arg2(a2), arg3(a3), arg4(a4), arg5(a5) {};   \
       virtual ~Exception5 () throw () {};                             \
       virtual void PrintInfo (std::ostream &out) const {              \
index 90349aaf96a9080edae774a17a65abea6f5b6c81..b0d925b8d1b13eedd94cf587b34c3d254a59be73 100644 (file)
@@ -23,7 +23,7 @@
 template <int dim> class Point;
 template <int rank_, int dim> class Tensor;
 template <int dim> class Tensor<1,dim>;
-template<typename number> class Vector;
+template <typename number> class Vector;
 
 /**
  * This class is a model for a general function. It serves the purpose
@@ -202,8 +202,8 @@ class Function : public FunctionTime,
                                      * size beforehand,
                                      * i.e. @p{n_components}.
                                      */
-    virtual void   vector_value (const Point<dim>   &p,
-                                Vector<double>     &values) const;
+    virtual void vector_value (const Point<dim>   &p,
+                               Vector<double>     &values) const;
     
                                     /**
                                      * Set @p{values} to the point
@@ -245,8 +245,8 @@ class Function : public FunctionTime,
                                      * components of the
                                      * function at the given point.
                                      */
-    virtual void          vector_gradient (const Point<dim>            &p,
-                                          std::vector<Tensor<1,dim> > &gradients) const;
+    virtual void vector_gradient (const Point<dim>            &p,
+                                  std::vector<Tensor<1,dim> > &gradients) const;
     
                                     /**
                                      * Set @p{gradients} to the
@@ -379,8 +379,8 @@ class ZeroFunction : public Function<dim>
                                      * at the given point for all
                                      * components.
                                      */
-    virtual void   vector_value (const Point<dim> &p,
-                                Vector<double>   &return_value) const;
+    virtual void vector_value (const Point<dim> &p,
+                               Vector<double>   &return_value) const;
 
                                     /**
                                      * Set @p{values} to the point values
@@ -421,8 +421,8 @@ class ZeroFunction : public Function<dim>
                                      * function at the given point,
                                      * for all components.
                                      */
-    virtual void          vector_gradient (const Point<dim>            &p,
-                                          std::vector<Tensor<1,dim> > &gradients) const;
+    virtual void vector_gradient (const Point<dim>            &p,
+                                  std::vector<Tensor<1,dim> > &gradients) const;
     
                                     /**
                                      * Set @p{gradients} to the gradients of
@@ -562,6 +562,8 @@ class ConstantFunction : public ZeroFunction<dim>
     const double function_value;
 };
 
+
+
 /**
  * This is a constant vector-valued function, that is different from
  * zero only in one component.
index d5a717f74e1de2ed0e86c83583bf95bd3a395882..a3cdfcaac39adf55799d45ceaa135fe2ddcb31ca 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002 by the deal authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -56,9 +56,9 @@ class FunctionDerivative : public Function<dim>
                                      * and the step size @p{h} of the
                                      * difference formula.
                                      */
-    FunctionDerivative (const Function<dim>f,
-                       const Point<dim>direction,
-                       const double h = 1.e-6);
+    FunctionDerivative (const Function<dim> &f,
+                       const Point<dim>    &direction,
+                       const double         h = 1.e-6);
     
                                     /**
                                      * Constructor. Provided are the
@@ -110,8 +110,8 @@ class FunctionDerivative : public Function<dim>
                                      * points.
                                      */
     virtual void value_list (const std::vector<Point<dim> > &points,
-                            std::vector<double>                     &values,
-                            const unsigned int                       component = 0) const;    
+                            std::vector<double>            &values,
+                            const unsigned int              component = 0) const;    
 
                                     /**
                                      * Determine an estimate for
@@ -147,10 +147,12 @@ class FunctionDerivative : public Function<dim>
                                      * formula.
                                      */
     double h;
+    
                                     /**
                                      * Difference formula.
                                      */
     DifferenceFormula formula;
+    
                                     /**
                                      * Helper object. Contains the
                                      * increment vector for the
index a43c781bbb8050f177aa5ad4b491ca1989c86ee9..376f1654713cc53721408ff62ccc8da905f54586 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002 by the deal authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -256,13 +256,15 @@ namespace Functions
                                   const unsigned int              component = 0) const;
       
                                       /**
-                                       * Gradient at a single point.
+                                       * Second derivatives at a
+                                       * single point.
                                        */
       virtual Tensor<2,dim> hessian (const Point<dim>   &p,
                                     const unsigned int  component = 0) const;
       
                                       /**
-                                       * Gradients at multiple points.
+                                       * Second derivatives at
+                                       * multiple points.
                                        */
       virtual void hessian_list (const std::vector<Point<dim> > &points,
                                 std::vector<Tensor<2,dim> >    &hessians,
@@ -502,7 +504,8 @@ namespace Functions
                            const unsigned int  component = 0) const;
       
                                       /**
-                                       * Function values at multiple points.
+                                       * Function values at multiple
+                                       * points.
                                        */
       virtual void value_list (const std::vector<Point<dim> > &points,
                               std::vector<double>            &values,
index a747284bb916f1391314ff4c33f5b4aef088f3f5..434c21bb7b9213f4d00a0b68b339339609ccc0dd 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002 by the deal authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
index d295606e6a2ec28ade4cd65feab690ae4c3f4314..c8ee8639e501ac1a06a0aa145ecd8d26a0969627 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002 by the deal authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -55,7 +55,7 @@ class JobIdentifier
                                     /**
                                      * Return the value of @p{id}.
                                      */
-    const std::string operator () ()const;
+    const std::string operator () () const;
     
   private:
                                     /**
index 09a9eb4ded02ea8fc678e2cc2e2f0c33fec0c7c6..9acea8a88601be2815899d838a8b6b4a28dcc6f5 100644 (file)
@@ -117,7 +117,7 @@ class LogStream
                                      * this function with @p{n=0}, no
                                      * console output will be written.
                                      */
-    void depth_console (unsigned int n);
+    void depth_console (const unsigned int n);
      
                                     /**
                                      * Maximum number of levels to be
@@ -128,7 +128,7 @@ class LogStream
                                      * with care, since it may spoile
                                      * the value of a log file.
                                      */
-    void depth_file (unsigned int n);
+    void depth_file (const unsigned int n);
 
                                      /**
                                      * Set time printing flag. If this flag
@@ -136,7 +136,7 @@ class LogStream
                                      * be prepended by the user time used
                                      * by the running program so far.
                                      */
-    void log_execution_time (bool flag);
+    void log_execution_time (const bool flag);
 
                                     /**
                                      * Output time differences
@@ -155,7 +155,7 @@ class LogStream
                                      * not changed by this function,
                                      * just the output.
                                      */
-    void log_time_differences (bool flag);
+    void log_time_differences (const bool flag);
     
                                     /**
                                      * Output a constant something
@@ -280,7 +280,7 @@ class LogStream
                                      * optional time information and
                                      * the contents of the prefix stack.
                                      */
-    void print_line_head();
+    void print_line_head ();
 
                                     /**
                                      * Actually do the work of
index 026b7e8ec981946d431fdfc6a83df905e985d2ef..e779d940da1df95f9d0376fbf1324c9aa78264a4 100644 (file)
@@ -390,7 +390,7 @@ namespace MemoryConsumption
 namespace MemoryConsumption
 {
   inline
-  unsigned int memory_consumption (const bool) 
+  unsigned int memory_consumption (const bool)
   {
     return sizeof(bool);
   }
index e1ad34062d2ce5beba41e2b48dea2fe4617662a1..e99b25bd9878b3cee50081270376b40abfeb2e9a 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 2000, 2001, 2002 by the deal authors
+//    Copyright (C) 2000, 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -46,7 +46,7 @@ class MultithreadInfo
                                      * detection failed or if
                                      * detection is not supported.
                                      */
-    MultithreadInfo();
+    MultithreadInfo ();
 
                                     /**
                                      * The number of CPUs in the
index 270ea4d554014c3807bf93a04e69f84a91a2795b..ced2614aadcf76718c62d51dd165763f71c70668 100644 (file)
@@ -1184,13 +1184,13 @@ class ParameterHandler
                                      * Inhibit automatic
                                      * CopyConstructor.
                                      */
-    ParameterHandler(const ParameterHandler&);
+    ParameterHandler (const ParameterHandler&);
 
                                     /**
                                      * Inhibit automatic
                                      * assignment operator.
                                      */
-    ParameterHandler& operator=(const ParameterHandler&);
+    ParameterHandler& operator= (const ParameterHandler&);
     
   public:
                                     /**
index 4c224d9983ec4c2a6fe1bc057532e5940e0c12d6..ad8f3933e6d1ab1935761a8c72a0f3bda26bfa42 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002 by the deal authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -47,12 +47,15 @@ class Point : public Tensor<1,dim>
 {
   public:
                                     /**
-                                     * Standard constructor. Creates an origin.
+                                     * Standard constructor. Creates
+                                     * an origin.
                                      */
     Point ();
+    
                                     /**
-                                     * Constructor. Initialize all entries
-                                     * to zero if @p{initialize==true}.
+                                     * Constructor. Initialize all
+                                     * entries to zero if
+                                     * @p{initialize==true}.
                                      */
     explicit Point (const bool initialize);
 
@@ -182,7 +185,8 @@ Point<dim>::Point ()
 
 template <int dim>
 inline
-Point<dim>::Point (const bool initialize) :
+Point<dim>::Point (const bool initialize)
+                :
                Tensor<1,dim>(initialize) 
 {}
 
@@ -190,7 +194,8 @@ Point<dim>::Point (const bool initialize) :
 
 template <int dim>
 inline
-Point<dim>::Point (const Tensor<1,dim> &t) :
+Point<dim>::Point (const Tensor<1,dim> &t)
+                :
                Tensor<1,dim>(t) 
 {}
 
@@ -241,7 +246,7 @@ double Point<dim>::operator () (const unsigned int index) const
 
 template <int dim>
 inline
-double & Point<dim>::operator () (const unsigned int index) 
+double & Point<dim>::operator () (const unsigned int index)
 {
   Assert (index<dim, ExcIndexRange (index, 0, dim));
   return this->values[index];
@@ -251,7 +256,7 @@ double & Point<dim>::operator () (const unsigned int index)
 
 template <int dim>
 inline
-Point<dim> Point<dim>::operator + (const Tensor<1,dim> &p) const 
+Point<dim> Point<dim>::operator + (const Tensor<1,dim> &p) const
 {
   return (Point<dim>(*this) += p);
 }
@@ -260,7 +265,7 @@ Point<dim> Point<dim>::operator + (const Tensor<1,dim> &p) const
 
 template <int dim>
 inline
-Point<dim> Point<dim>::operator - (const Tensor<1,dim> &p) const 
+Point<dim> Point<dim>::operator - (const Tensor<1,dim> &p) const
 {
   return (Point<dim>(*this) -= p);
 }
@@ -269,7 +274,7 @@ Point<dim> Point<dim>::operator - (const Tensor<1,dim> &p) const
 
 template <int dim>
 inline
-Point<dim> Point<dim>::operator - () const 
+Point<dim> Point<dim>::operator - () const
 {
   Point<dim> result;
   for (unsigned int i=0; i<dim; ++i)
@@ -281,7 +286,7 @@ Point<dim> Point<dim>::operator - () const
 
 template <int dim>
 inline
-Point<dim> Point<dim>::operator * (const double factor) const 
+Point<dim> Point<dim>::operator * (const double factor) const
 {
   return (Point<dim>(*this) *= factor);
 }
@@ -290,7 +295,7 @@ Point<dim> Point<dim>::operator * (const double factor) const
 
 template <int dim>
 inline
-double Point<dim>::operator * (const Tensor<1,dim> &p) const 
+double Point<dim>::operator * (const Tensor<1,dim> &p) const
 {
                                   // simply pass down
   return Tensor<1,dim>::operator * (p);
@@ -299,7 +304,7 @@ double Point<dim>::operator * (const Tensor<1,dim> &p) const
 
 template <int dim>
 inline
-double Point<dim>::square () const 
+double Point<dim>::square () const
 {
   double q=0;
   for (unsigned int i=0; i<dim; ++i)
@@ -325,7 +330,7 @@ double Point<dim>::distance (const Point<dim> &p) const
 
 template <int dim>
 inline
-Point<dim> Point<dim>::operator / (const double factor) const 
+Point<dim> Point<dim>::operator / (const double factor) const
 {
   return (Point<dim>(*this) /= factor);
 }
@@ -340,7 +345,7 @@ Point<dim> Point<dim>::operator / (const double factor) const
  */
 template <int dim>
 inline
-Point<dim> operator * (const double factor, const Point<dim> &p) 
+Point<dim> operator * (const double factor, const Point<dim> &p)
 {
   return p*factor;
 }
index 4a2682425e5d5bb247b4376b0a034b2074452886..2a94844ecac698b077b6b0470e767d7992599d45 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002 by the deal authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -104,7 +104,7 @@ class Quadrature : public Subscriptor
                                      * @p{Quadrature<dim-1>}.
                                      */
     Quadrature (const SubQuadrature &,
-               const Quadrature<1>                     &);
+               const Quadrature<1> &);
     
                                     /**
                                      * Construct a quadrature formula
@@ -149,13 +149,15 @@ class Quadrature : public Subscriptor
     virtual ~Quadrature ();
     
                                     /**
-                                     * Return the @p{i}th quadrature point.
+                                     * Return the @p{i}th quadrature
+                                     * point.
                                      */
     const Point<dim> & point (const unsigned int i) const;
 
                                     /**
-                                     * Return a reference to the whole array of
-                                     * quadrature points.
+                                     * Return a reference to the
+                                     * whole array of quadrature
+                                     * points.
                                      */
     const std::vector<Point<dim> > & get_points () const;
     
@@ -195,6 +197,7 @@ class Quadrature : public Subscriptor
 };
 
 
+
 /**
  * Quadrature formula constructed by iteration of another quadrature formula in
  * each direction. In more than one space dimension, the resulting quadrature
@@ -255,10 +258,12 @@ class QIterated : public Quadrature<dim>
                                      * points at the left and right end points
                                      * of the interval.
                                      */
-    static bool uses_both_endpoints (const Quadrature<1> &base_quadrature);
+    static bool
+    uses_both_endpoints (const Quadrature<1> &base_quadrature);
 };
 
 
+
 /**
  *  This class is a helper class to facilitate the usage of quadrature formulae
  *  on faces or subfaces of cells. It computes the locations of quadrature
@@ -417,45 +422,54 @@ template <>
 const std::vector<double> & Quadrature<0>::get_weights () const;
 
 template <>
-void QProjector<1>::project_to_face (const Quadrature<0> &,
-                                    const unsigned int,
-                                    std::vector<Point<1> > &);
+void
+QProjector<1>::project_to_face (const Quadrature<0> &,
+                                const unsigned int,
+                                std::vector<Point<1> > &);
 template <>
-void QProjector<2>::project_to_face (const Quadrature<1>      &quadrature,
-                                    const unsigned int        face_no,
-                                    std::vector<Point<2> >   &q_points);
+void
+QProjector<2>::project_to_face (const Quadrature<1>      &quadrature,
+                                const unsigned int        face_no,
+                                std::vector<Point<2> >   &q_points);
 template <>
-void QProjector<3>::project_to_face (const Quadrature<2>    &quadrature,
-                                    const unsigned int      face_no,
-                                    std::vector<Point<3> > &q_points);
+void
+QProjector<3>::project_to_face (const Quadrature<2>    &quadrature,
+                                const unsigned int      face_no,
+                                std::vector<Point<3> > &q_points);
 
 template <>
-Quadrature<1> QProjector<1>::project_to_all_faces (const Quadrature<0> &quadrature);
+Quadrature<1>
+QProjector<1>::project_to_all_faces (const Quadrature<0> &quadrature);
 
 
 template <>
-void QProjector<1>::project_to_subface (const Quadrature<0> &,
-                                       const unsigned int,
-                                       const unsigned int,
-                                       std::vector<Point<1> > &);
+void
+QProjector<1>::project_to_subface (const Quadrature<0> &,
+                                   const unsigned int,
+                                   const unsigned int,
+                                   std::vector<Point<1> > &);
 template <>
-void QProjector<2>::project_to_subface (const Quadrature<1>    &quadrature,
-                                       const unsigned int      face_no,
-                                       const unsigned int      subface_no,
-                                       std::vector<Point<2> > &q_points);
+void
+QProjector<2>::project_to_subface (const Quadrature<1>    &quadrature,
+                                   const unsigned int      face_no,
+                                   const unsigned int      subface_no,
+                                   std::vector<Point<2> > &q_points);
 template <>
-void QProjector<3>::project_to_subface (const Quadrature<2>    &quadrature,
-                                       const unsigned int      face_no,
-                                       const unsigned int      subface_no,
-                                       std::vector<Point<3> > &q_points);
+void
+QProjector<3>::project_to_subface (const Quadrature<2>    &quadrature,
+                                   const unsigned int      face_no,
+                                   const unsigned int      subface_no,
+                                   std::vector<Point<3> > &q_points);
 
 template <>
-Quadrature<1> QProjector<1>::project_to_all_subfaces (const Quadrature<0> &quadrature);
+Quadrature<1>
+QProjector<1>::project_to_all_subfaces (const Quadrature<0> &quadrature);
 
 
 template <>
 bool
 QIterated<1>::uses_both_endpoints (const Quadrature<1> &base_quadrature);
+
 template <>
 QIterated<1>::QIterated (const Quadrature<1> &base_quadrature,
                         const unsigned int   n_copies);
index ef32d229f39edf00a7fd4fc7c35c759797ecd6df..3434c8d5ca8aebef0bcfa3ab294aabb14843bb85 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002 by the deal authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -57,29 +57,33 @@ class SmartPointer
                                     /**
                                      * Standard constructor for null pointer.
                                      */
-    SmartPointer();
+    SmartPointer ();
 
                                     /*
                                      * Copy constructor for
-                                     * @p{SmartPointer}. We copy the object
-                                     * subscribed to from @p{tt}, but subscribe
+                                     * @p{SmartPointer}. We do now
+                                     * copy the object subscribed to
+                                     * from @p{tt}, but subscribe
                                      * ourselves to it again.
                                      */
     SmartPointer (const SmartPointer<T> &tt);
 
                                     /**
-                                     * Constructor taking a normal pointer.
-                                     * If possible, i.e. if the pointer
-                                     * is not a null pointer, the constructor
-                                     * subscribes to the given object to
-                                     * lock it, i.e. to prevent its
-                                     * destruction before the end of its use.
+                                     * Constructor taking a normal
+                                     * pointer.  If possible, i.e. if
+                                     * the pointer is not a null
+                                     * pointer, the constructor
+                                     * subscribes to the given object
+                                     * to lock it, i.e. to prevent
+                                     * its destruction before the end
+                                     * of its use.
                                      */
     SmartPointer (T *t);
 
 
                                     /**
-                                     * Destructor, removing the subscription.
+                                     * Destructor, removing the
+                                     * subscription.
                                      */
     ~SmartPointer();
     
@@ -183,14 +187,16 @@ class SmartPointer
 
 
 template <typename T>
-SmartPointer<T>::SmartPointer () :
+SmartPointer<T>::SmartPointer ()
+                :
                t (0)
 {}
 
 
 
 template <typename T>
-SmartPointer<T>::SmartPointer (T *t) :
+SmartPointer<T>::SmartPointer (T *t)
+                :
                t (t)
 {
   if (t != 0)
@@ -200,7 +206,8 @@ SmartPointer<T>::SmartPointer (T *t) :
 
 
 template <typename T>
-SmartPointer<T>::SmartPointer (const SmartPointer<T> &tt) :
+SmartPointer<T>::SmartPointer (const SmartPointer<T> &tt)
+                :
                t (tt.t)
 {
   if (t != 0)
@@ -237,7 +244,8 @@ SmartPointer<T> & SmartPointer<T>::operator = (T *tt)
 
 
 template <typename T>
-SmartPointer<T> & SmartPointer<T>::operator = (const SmartPointer<T>& tt)
+SmartPointer<T> &
+SmartPointer<T>::operator = (const SmartPointer<T>& tt)
 {
                                   // if objects on the left and right
                                   // hand side of the operator= are
index 38d6c8a4c873217f3b6ade893b7e91009e142fc0..ed484154066aaf722dda3f1e5ae7e75a1cbf41f6 100644 (file)
@@ -409,6 +409,7 @@ namespace internal
                                           * private.
                                           */
         Accessor ();
+        
                                          /**
                                           * Copy constructor. Not
                                           * needed, and invisible, so
@@ -897,7 +898,7 @@ class TableBase : public Subscriptor
                                       * don't know here whether
                                       * copying is expensive or not.
                                       */
-    const T & el (const TableIndices<N> &indices) const;    
+    const T & el (const TableIndices<N> &indices) const;
   
                                      /**
                                       * Direct read-only access to
@@ -953,7 +954,8 @@ class TableBase : public Subscriptor
  */
 template <int N,typename T>
 class Table : public TableBase<N,T>
-{};
+{
+};
 
 
 /**
@@ -1575,7 +1577,7 @@ TableIndices<1>::TableIndices ()
 
 
 inline
-TableIndices<1>::TableIndices (const unsigned int index1) 
+TableIndices<1>::TableIndices (const unsigned int index1)
 {
   this->indices[0] = index1;
 }
@@ -1620,7 +1622,7 @@ TableIndices<3>::TableIndices (const unsigned int index1,
 
 
 inline
-TableIndices<4>::TableIndices () 
+TableIndices<4>::TableIndices ()
 {
   this->indices[0] = this->indices[1] = this->indices[2] = this->indices[3] = 0;
 }
@@ -1642,9 +1644,12 @@ TableIndices<4>::TableIndices (const unsigned int index1,
 
 
 inline
-TableIndices<5>::TableIndices () 
+TableIndices<5>::TableIndices ()
 {
-  this->indices[0] = this->indices[1] = this->indices[2] = this->indices[3] = this->indices[4] = 0;
+  this->indices[0] = this->indices[1]
+                   = this->indices[2]
+                   = this->indices[3]
+                   = this->indices[4] = 0;
 }
 
 
@@ -1666,7 +1671,7 @@ TableIndices<5>::TableIndices (const unsigned int index1,
 
 
 inline
-TableIndices<6>::TableIndices () 
+TableIndices<6>::TableIndices ()
 {
   this->indices[0] = this->indices[1] = this->indices[2]
                   = this->indices[3] = this->indices[4] = 0;
@@ -1788,7 +1793,7 @@ namespace internal
     template <int N, typename T, bool C, unsigned int P>
     inline
     Accessor<N,T,C,P-1>
-    Accessor<N,T,C,P>::operator [] (const unsigned int i) const 
+    Accessor<N,T,C,P>::operator [] (const unsigned int i) const
     {
       Assert (i < table.size()[N-P],
               ExcIndexRange (i, 0, table.size()[N-P]));
@@ -1857,7 +1862,7 @@ namespace internal
     template <int N, typename T, bool C>
     inline
     typename Accessor<N,T,C,1>::reference
-    Accessor<N,T,C,1>::operator [] (const unsigned int i) const 
+    Accessor<N,T,C,1>::operator [] (const unsigned int i) const
     {
       Assert (i < table.size()[N-1],
               ExcIndexRange (i, 0, table.size()[N-1]));
@@ -1910,7 +1915,7 @@ TableBase<N,T>::~TableBase ()
 
 template <int N, typename T>
 TableBase<N,T>&
-TableBase<N,T>::operator = (const TableBase<N,T>& m) 
+TableBase<N,T>::operator = (const TableBase<N,T>& m)
 {
   reinit (m.size());
   if (!empty())
@@ -1924,7 +1929,7 @@ TableBase<N,T>::operator = (const TableBase<N,T>& m)
 template <int N, typename T>
 template <typename T2>
 TableBase<N,T>&
-TableBase<N,T>::operator = (const TableBase<N,T2>& m) 
+TableBase<N,T>::operator = (const TableBase<N,T2>& m)
 {
   reinit (m.size());
   if (!empty())
@@ -2047,7 +2052,7 @@ TableBase<N,T>::memory_consumption () const
 template <int N, typename T>
 inline
 unsigned int
-TableBase<N,T>::position (const TableIndices<N> &indices) const 
+TableBase<N,T>::position (const TableIndices<N> &indices) const
 {
                                    // specialize this for the
                                    // different numbers of dimensions,
@@ -2130,7 +2135,7 @@ TableBase<N,T>::data () const
 
 
 template <typename T>
-Table<1,T>::Table () 
+Table<1,T>::Table ()
 {}
 
 
@@ -2316,7 +2321,7 @@ Table<2,T>::n_cols () const
 
 
 template <typename T>
-Table<3,T>::Table () 
+Table<3,T>::Table ()
 {}
 
 
@@ -2402,7 +2407,7 @@ Table<3,T>::operator () (const unsigned int i,
 
 
 template <typename T>
-Table<4,T>::Table () 
+Table<4,T>::Table ()
 {}
 
 
@@ -2500,7 +2505,7 @@ Table<4,T>::operator () (const unsigned int i,
 
 
 template <typename T>
-Table<5,T>::Table () 
+Table<5,T>::Table ()
 {}
 
 
@@ -2608,7 +2613,7 @@ Table<5,T>::operator () (const unsigned int i,
 
 
 template <typename T>
-Table<6,T>::Table () 
+Table<6,T>::Table ()
 {}
 
 
index 571eef05e7fce184a13997f4a0505847f0558c21..0cebb6f1b5df7d2a3bdf33d8e4bd723c0654d8ea 100644 (file)
@@ -54,12 +54,12 @@ class TableEntryBase
                                      /**
                                       * Write the table entry as text.
                                       */
-    virtual void write_text (std::ostream &) const =0;
+    virtual void write_text (std::ostream &) const = 0;
     
                                      /**
                                       * Write the table entry in tex format.
                                       */
-    virtual void write_tex (std::ostream &) const =0;
+    virtual void write_tex (std::ostream &) const = 0;
 };
 
 
@@ -81,7 +81,7 @@ class TableEntry : public TableEntryBase
                                      /**
                                       * Constructor.
                                       */
-    TableEntry(const value_type value);
+    TableEntry (const value_type value);
 
                                      /**
                                       * Destructor.
@@ -92,7 +92,7 @@ class TableEntry : public TableEntryBase
                                       * Returns the value of this
                                       * table entry.
                                       */
-    value_type value() const;
+    value_type value () const;
 
                                      /**
                                       * Write the table entry as text.
@@ -180,7 +180,7 @@ class TableHandler
                                      /**
                                       * Constructor.
                                       */
-    TableHandler();
+    TableHandler ();
     
                                      /**
                                       * Adds a column (if not yet
index 8bafe9c136ea66222ec53c6ececc1f7faaa0bedd..b39b5ee528d8751402730721d5cbaa5a92791be9 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002 by the deal authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -209,7 +209,7 @@ class Tensor
                                     // also, it would be sufficient to make
                                     // the function unroll_loops a friend,
                                     // but that seems to be impossible as well.
-    template<int otherrank, int otherdim> friend class Tensor;
+    template <int otherrank, int otherdim> friend class Tensor;
 };
 
 
@@ -706,7 +706,7 @@ void contract (Tensor<4,dim>       &dest,
 template <int dim>
 void outer_product (Tensor<2,dim>       &dst,
                    const Tensor<1,dim> &src1,
-                   const Tensor<1,dim> &src2) 
+                   const Tensor<1,dim> &src2)
 {
   for (unsigned int i=0; i<dim; ++i)
     for (unsigned int j=0; j<dim; ++j)
@@ -724,7 +724,7 @@ void outer_product (Tensor<2,dim>       &dst,
 template <int dim>
 void outer_product (Tensor<3,dim>       &dst,
                    const Tensor<1,dim> &src1,
-                   const Tensor<2,dim> &src2) 
+                   const Tensor<2,dim> &src2)
 {
   for (unsigned int i=0; i<dim; ++i)
     for (unsigned int j=0; j<dim; ++j)
@@ -743,7 +743,7 @@ void outer_product (Tensor<3,dim>       &dst,
 template <int dim>
 void outer_product (Tensor<3,dim>       &dst,
                    const Tensor<2,dim> &src1,
-                   const Tensor<1,dim> &src2) 
+                   const Tensor<1,dim> &src2)
 {
   for (unsigned int i=0; i<dim; ++i)
     for (unsigned int j=0; j<dim; ++j)
@@ -766,7 +766,7 @@ void outer_product (Tensor<3,dim>       &dst,
 template <int dim>
 void outer_product (Tensor<1,dim>       &dst,
                    const double         src1,
-                   const Tensor<1,dim> &src2) 
+                   const Tensor<1,dim> &src2)
 {
   for (unsigned int i=0; i<dim; ++i)
     dst[i] = src1 * src2[i];
@@ -787,7 +787,7 @@ void outer_product (Tensor<1,dim>       &dst,
 template <int dim>
 void outer_product (Tensor<1,dim>       &dst,
                    const Tensor<1,dim>  src1,
-                   const double         src2) 
+                   const double         src2)
 {
   for (unsigned int i=0; i<dim; ++i)
     dst[i] = src1[i] * src2;
@@ -1018,7 +1018,7 @@ invert (const Tensor<2,dim> &t)
 template <int dim>
 inline
 Tensor<2,dim>
-transpose (const Tensor<2,dim> &t) 
+transpose (const Tensor<2,dim> &t)
 {
   Tensor<2,dim> tt = t;
   for (unsigned int i=0; i<dim; ++i)
@@ -1041,7 +1041,7 @@ transpose (const Tensor<2,dim> &t)
  */
 inline
 Tensor<2,1>
-transpose (const Tensor<2,1> &t) 
+transpose (const Tensor<2,1> &t)
 {
   return t;
 }
@@ -1057,7 +1057,7 @@ transpose (const Tensor<2,1> &t)
  */
 inline
 Tensor<2,2>
-transpose (const Tensor<2,2> &t) 
+transpose (const Tensor<2,2> &t)
 {
   const double x[2][2] = {{t[0][0], t[1][0]}, {t[0][1], t[1][1]}};
   return Tensor<2,2>(x);
@@ -1074,7 +1074,7 @@ transpose (const Tensor<2,2> &t)
  */
 inline
 Tensor<2,3>
-transpose (const Tensor<2,3> &t) 
+transpose (const Tensor<2,3> &t)
 {
   const double x[3][3] = {{t[0][0], t[1][0], t[2][0]},
                           {t[0][1], t[1][1], t[2][1]},
index 7030178acc58682bbda2dcd9945fb9fdd0f8e883..1e0429ff47b21380e2e9eab2328eb0b4fd43a089 100644 (file)
@@ -147,12 +147,14 @@ class Tensor<1,dim>
     Tensor<1,dim> & operator = (const Tensor<1,dim> &);
 
                                     /**
-                                     * Test for equality of two points.
+                                     * Test for equality of two
+                                     * tensors.
                                      */
     bool operator == (const Tensor<1,dim> &) const;
 
                                     /**
-                                     * Test for inequality of two points.
+                                     * Test for inequality of two
+                                     * tensors.
                                      */
     bool operator != (const Tensor<1,dim> &) const;
 
@@ -162,6 +164,7 @@ class Tensor<1,dim>
                                      * offset.
                                      */
     Tensor<1,dim> & operator += (const Tensor<1,dim> &);
+    
                                     /**
                                      * Subtract another vector.
                                      */
@@ -281,7 +284,7 @@ class Tensor<1,dim>
                                      * that seems to be impossible as
                                      * well.
                                      */
-    template<int otherrank, int otherdim>  friend class Tensor;
+    template <int otherrank, int otherdim>  friend class Tensor;
 
                                     /**
                                      * Point is allowed access to
index e61d0ad6b32d8e9f5b36b68a4ab76a07acb0a624..d1fc171c5b6692514fe9b4d1ce2d15f5240f1f34 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002 by the deal authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
index 03383d5478fbc54f951c215427f8f2be7e90f735..636e91961e64379d3f8c17b32288060a31a7697c 100644 (file)
@@ -704,7 +704,7 @@ namespace Threads
                                       * for use in the template
                                       * functions below.
                                       */
-    void deregister_thread ();  
+    void deregister_thread ();
   }
   
 }   // end declarations of namespace Threads
index a4551d4fb5a0e3ecbd9cd5afe1b5aa2a09a7bf00..427b025ff9eed862b321819fb7efdd7e1f194422 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002 by the deal authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
index af219f9f4009f8259851ae5179b322f3da3001e5..5527ffea404442cdfa9b01b0da695e7c18fbdc2a 100644 (file)
 
 
 template <int dim>
-AutoDerivativeFunction<dim>::AutoDerivativeFunction (const double hh,
-                                                    const unsigned int n_components,
-                                                    const double       initial_time):
+AutoDerivativeFunction<dim>::
+AutoDerivativeFunction (const double hh,
+                        const unsigned int n_components,
+                        const double       initial_time)
+                :
                Function<dim>(n_components, initial_time),
                h(1),
                ht(dim),
@@ -111,8 +113,10 @@ AutoDerivativeFunction<dim>::gradient (const Point<dim>   &p,
 
 
 template <int dim>
-void AutoDerivativeFunction<dim>::vector_gradient (const Point<dim>       &p,
-                                                  std::vector<Tensor<1,dim> > &gradients) const
+void
+AutoDerivativeFunction<dim>::
+vector_gradient (const Point<dim>            &p,
+                 std::vector<Tensor<1,dim> > &gradients) const
 {
   Assert (gradients.size() == this->n_components,
          ExcDimensionMismatch(gradients.size(), this->n_components));
@@ -135,6 +139,7 @@ void AutoDerivativeFunction<dim>::vector_gradient (const Point<dim>       &p,
          }
        break;
       }
+
       case Euler:
       {
        Point<dim> q1, q2;
@@ -152,6 +157,7 @@ void AutoDerivativeFunction<dim>::vector_gradient (const Point<dim>       &p,
          }
        break;
       }
+
       case FourthOrder:
       {
        Point<dim> q1, q2, q3, q4;
@@ -175,6 +181,7 @@ void AutoDerivativeFunction<dim>::vector_gradient (const Point<dim>       &p,
          }
        break;
       }
+
       default:
            Assert(false, ExcInvalidFormula());
     }
@@ -182,9 +189,11 @@ void AutoDerivativeFunction<dim>::vector_gradient (const Point<dim>       &p,
 
 
 template <int dim>
-void AutoDerivativeFunction<dim>::gradient_list (const std::vector<Point<dim> > &points,
-                                                std::vector<Tensor<1,dim> >    &gradients,
-                                                const unsigned int              comp) const
+void
+AutoDerivativeFunction<dim>::
+gradient_list (const std::vector<Point<dim> > &points,
+               std::vector<Tensor<1,dim> >    &gradients,
+               const unsigned int              comp) const
 {
   Assert (gradients.size() == points.size(),
          ExcDimensionMismatch(gradients.size(), points.size()));
@@ -202,6 +211,7 @@ void AutoDerivativeFunction<dim>::gradient_list (const std::vector<Point<dim> >
            }
        break;
       }
+
       case Euler:
       {
        Point<dim> q1, q2;
@@ -214,6 +224,7 @@ void AutoDerivativeFunction<dim>::gradient_list (const std::vector<Point<dim> >
            }
        break;
       }
+
       case FourthOrder:
       {
        Point<dim> q1, q2, q3, q4;
@@ -231,6 +242,7 @@ void AutoDerivativeFunction<dim>::gradient_list (const std::vector<Point<dim> >
            }
        break;
       }
+
       default:
            Assert(false, ExcInvalidFormula());
     }
@@ -264,6 +276,7 @@ vector_gradient_list (const std::vector<Point<dim> >            &points,
          }
        break;
       }
+
       case Euler:
       {
        Point<dim> q1, q2;
@@ -273,10 +286,12 @@ vector_gradient_list (const std::vector<Point<dim> >            &points,
              q1=points[p]+ht[i];
              q2=points[p]-ht[i];
              for (unsigned int comp=0; comp<this->n_components; ++comp)
-               gradients[p][comp][i]=(this->value(q1, comp)-this->value(q2, comp))/(2*h);
+               gradients[p][comp][i]=(this->value(q1, comp) -
+                                       this->value(q2, comp))/(2*h);
            }
        break;
       }
+
       case FourthOrder:
       {
        Point<dim> q1, q2, q3, q4;
@@ -295,6 +310,7 @@ vector_gradient_list (const std::vector<Point<dim> >            &points,
            }
        break;
       }
+      
       default:
            Assert(false, ExcInvalidFormula());
     }
index f9d17e4036cda6fa6bbfd01f87e7beb9119c1283..6b7ac3291758a096f283d3016d7198c065da8f01 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002 by the deal authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -16,7 +16,7 @@
 #include <cmath>
 
 
-ConvergenceTable::ConvergenceTable()  
+ConvergenceTable::ConvergenceTable()
 {}
 
 
@@ -77,8 +77,10 @@ void ConvergenceTable::evaluate_convergence_rates(const std::string &data_column
 }
 
 
-void ConvergenceTable::evaluate_convergence_rates(const std::string &data_column_key,
-                                                 const RateMode     rate_mode)
+
+void
+ConvergenceTable::evaluate_convergence_rates(const std::string &data_column_key,
+                                             const RateMode     rate_mode)
 {
   Assert(columns.count(data_column_key), ExcColumnNotExistent(data_column_key));
   
@@ -106,6 +108,7 @@ void ConvergenceTable::evaluate_convergence_rates(const std::string &data_column
     {
       case none:
            break;
+
       case reduction_rate:
            rate_key+="red.rate";
            Assert(columns.count(rate_key)==0, ExcRateColumnAlreadyExists(rate_key));
@@ -115,6 +118,7 @@ void ConvergenceTable::evaluate_convergence_rates(const std::string &data_column
            for (unsigned int i=1; i<n; ++i)
              add_value(rate_key, values[i-1]/values[i]);
            break;
+
       case reduction_rate_log2:
            rate_key+="red.rate.log2";
            Assert(columns.count(rate_key)==0, ExcRateColumnAlreadyExists(rate_key));
@@ -124,6 +128,7 @@ void ConvergenceTable::evaluate_convergence_rates(const std::string &data_column
            for (unsigned int i=1; i<n; ++i)
              add_value(rate_key, std::log(values[i-1]/values[i])/std::log(2.0));
            break;
+
       default:
            ExcNotImplemented();  
     }
@@ -145,7 +150,10 @@ void ConvergenceTable::evaluate_convergence_rates(const std::string &data_column
   add_column_to_supercolumn(rate_key, superkey);
 }
 
-void ConvergenceTable::omit_column_from_convergence_rate_evaluation(const std::string &key)
+
+
+void
+ConvergenceTable::omit_column_from_convergence_rate_evaluation(const std::string &key)
 {
   Assert(columns.count(key), ExcColumnNotExistent(key));
   
@@ -154,8 +162,10 @@ void ConvergenceTable::omit_column_from_convergence_rate_evaluation(const std::s
 }
 
 
-void ConvergenceTable::evaluate_all_convergence_rates(const std::string &reference_column_key,
-                                                     const RateMode rate_mode)
+
+void
+ConvergenceTable::evaluate_all_convergence_rates(const std::string &reference_column_key,
+                                                 const RateMode rate_mode)
 {
   for (std::map<std::string, Column>::const_iterator col_iter=columns.begin();
        col_iter!=columns.end(); ++col_iter)
@@ -164,7 +174,9 @@ void ConvergenceTable::evaluate_all_convergence_rates(const std::string &referen
 }
 
 
-void ConvergenceTable::evaluate_all_convergence_rates(const RateMode rate_mode)
+
+void
+ConvergenceTable::evaluate_all_convergence_rates(const RateMode rate_mode)
 {
   for (std::map<std::string, Column>::const_iterator col_iter=columns.begin();
        col_iter!=columns.end(); ++col_iter)
index 37feed730bc1159584ebd6c41716543cee1c694a..677801b0b75946e66191e274c061f9ab1aaab5b3 100644 (file)
@@ -28,10 +28,10 @@ const unsigned int DataOutBase::Patch<dim,spacedim>::no_neighbor;
 
 
 template <int dim, int spacedim>
-DataOutBase::Patch<dim,spacedim>::Patch () :
+DataOutBase::Patch<dim,spacedim>::Patch ()
+                :
                patch_index(no_neighbor),
                n_subdivisions (1)
-
                                   // all the other data has a
                                   // constructor of its own, except
                                   // for the "neighbors" field, which
@@ -60,7 +60,8 @@ DataOutBase::Patch<dim,spacedim>::memory_consumption () const
 
 
 
-DataOutBase::UcdFlags::UcdFlags (const bool write_preamble) :
+DataOutBase::UcdFlags::UcdFlags (const bool write_preamble)
+                :
                write_preamble (write_preamble)
 {}
 
@@ -68,7 +69,8 @@ DataOutBase::UcdFlags::UcdFlags (const bool write_preamble) :
 
 DataOutBase::PovrayFlags::PovrayFlags (const bool smooth,
                                       const bool bicubic_patch,
-                                      const bool external_data) :
+                                      const bool external_data)
+                :
                smooth (smooth),
                bicubic_patch(bicubic_patch),
                external_data(external_data)
@@ -76,7 +78,8 @@ DataOutBase::PovrayFlags::PovrayFlags (const bool smooth,
 
 
 DataOutBase::DXFlags::DXFlags (const bool write_multigrid,
-                              const bool write_neighbors) :
+                              const bool write_neighbors)
+                :
                write_multigrid(write_multigrid),
                write_neighbors(write_neighbors)
 {}
@@ -195,7 +198,8 @@ DataOutBase::EpsFlags::EpsFlags (const unsigned int  height_vector,
                                 const bool          draw_mesh,
                                 const bool          draw_cells,
                                 const bool          shade_cells,
-                                const ColorFunction color_function) :
+                                const ColorFunction color_function)
+                :
                height_vector(height_vector),
                color_vector(color_vector),
                size_type(size_type),
@@ -425,10 +429,11 @@ DataOutBase::GmvFlags::memory_consumption () const
 
 
 
-DataOutBase::TecplotFlags::TecplotFlags (const char* tecplot_binary_file_name) :
-  tecplot_binary_file_name(tecplot_binary_file_name)
-{
-}
+DataOutBase::TecplotFlags::
+TecplotFlags (const char* tecplot_binary_file_name)
+                :
+                tecplot_binary_file_name(tecplot_binary_file_name)
+{}
 
 
 
@@ -3082,21 +3087,21 @@ namespace
 {
   class TecplotMacros
   {
-  public:
-    TecplotMacros(const unsigned int n_nodes = 0,
-                 const unsigned int n_vars = 0,
-                 const unsigned int n_cells = 0,
-                 const unsigned int n_vert = 0);
-    ~TecplotMacros();
-    float & nd(const unsigned int i, const unsigned int j);
-    int   & cd(const unsigned int i, const unsigned int j);
-    std::vector<float> nodalData;
-    std::vector<int>   connData;
-  private:
-    unsigned int n_nodes;
-    unsigned int n_vars;
-    unsigned int n_cells;
-    unsigned int n_vert;
+    public:
+      TecplotMacros(const unsigned int n_nodes = 0,
+                    const unsigned int n_vars = 0,
+                    const unsigned int n_cells = 0,
+                    const unsigned int n_vert = 0);
+      ~TecplotMacros();
+      float & nd(const unsigned int i, const unsigned int j);
+      int   & cd(const unsigned int i, const unsigned int j);
+      std::vector<float> nodalData;
+      std::vector<int>   connData;
+    private:
+      unsigned int n_nodes;
+      unsigned int n_vars;
+      unsigned int n_cells;
+      unsigned int n_vert;
   };
 
 
@@ -3104,7 +3109,8 @@ namespace
   TecplotMacros::TecplotMacros(const unsigned int n_nodes,
                               const unsigned int n_vars,
                               const unsigned int n_cells,
-                              const unsigned int n_vert) :
+                              const unsigned int n_vert)
+                  :
                  n_nodes(n_nodes),
                  n_vars(n_vars),
                  n_cells(n_cells),
@@ -3118,13 +3124,13 @@ namespace
 
   inline
   TecplotMacros::~TecplotMacros()
-  {
-  }
+  {}
 
 
 
   inline
-  float & TecplotMacros::nd(const unsigned int i, const unsigned int j)
+  float & TecplotMacros::nd (const unsigned int i,
+                             const unsigned int j)
   {
     return nodalData[(i)*(n_nodes) + (j)]; 
   }
@@ -3132,7 +3138,8 @@ namespace
 
 
   inline
-  int & TecplotMacros::cd(const unsigned int i, const unsigned int j)
+  int & TecplotMacros::cd (const unsigned int i,
+                           const unsigned int j)
   {
     return connData[(i) + (j)*(n_vert)]; 
   }
@@ -3154,22 +3161,27 @@ void DataOutBase::write_tecplot_binary (const std::vector<Patch<dim,spacedim> >
   
 #ifndef DEAL_II_HAVE_TECPLOT
   
-  // simply call the ASCII output function if the Tecplot API isn't present
+                                   // simply call the ASCII output
+                                   // function if the Tecplot API
+                                   // isn't present
   write_tecplot (patches, data_names, flags, out);
   return;
   
 #else
   
-  // Tecplot binary output only good for 2D & 3D
+                                   // Tecplot binary output only good
+                                   // for 2D & 3D
   if (dim == 1)
     {
       write_tecplot (patches, data_names, flags, out);
       return;
     };
 
-  // if the user hasn't specified a file name we should
-  // call the ASCII function and use the ostream @p{out}
-  // instead of doing something silly later
+                                   // if the user hasn't specified a
+                                   // file name we should call the
+                                   // ASCII function and use the
+                                   // ostream @p{out} instead of doing
+                                   // something silly later
   char* file_name = (char*) flags.tecplot_binary_file_name;
 
   if (file_name == NULL)
@@ -3197,7 +3209,8 @@ void DataOutBase::write_tecplot_binary (const std::vector<Patch<dim,spacedim> >
                                   // and cells for later use
   unsigned int n_cells = 0,
                n_nodes = 0;
-  for (typename std::vector<Patch<dim,spacedim> >::const_iterator patch=patches.begin();
+  for (typename std::vector<Patch<dim,spacedim> >::const_iterator
+         patch=patches.begin();
        patch!=patches.end(); ++patch)
     switch (dim)
       {
@@ -4078,16 +4091,18 @@ DataOutInterface<dim,spacedim>::write (std::ostream &out,
 
 
 template <int dim, int spacedim>
-void DataOutInterface<dim,spacedim>::set_default_format(const OutputFormat fmt)
+void
+DataOutInterface<dim,spacedim>::set_default_format(const OutputFormat fmt)
 {
-  Assert(fmt != default_format, ExcNotImplemented());
+  Assert (fmt != default_format, ExcNotImplemented());
   default_fmt = fmt;
 }
 
 
 
 template <int dim, int spacedim>
-void DataOutInterface<dim,spacedim>::set_flags (const DXFlags &flags) 
+void
+DataOutInterface<dim,spacedim>::set_flags (const DXFlags &flags) 
 {
   dx_flags = flags;
 }
@@ -4095,7 +4110,8 @@ void DataOutInterface<dim,spacedim>::set_flags (const DXFlags &flags)
 
 
 template <int dim, int spacedim>
-void DataOutInterface<dim,spacedim>::set_flags (const UcdFlags &flags) 
+void
+DataOutInterface<dim,spacedim>::set_flags (const UcdFlags &flags) 
 {
   ucd_flags = flags;
 }
@@ -4103,7 +4119,8 @@ void DataOutInterface<dim,spacedim>::set_flags (const UcdFlags &flags)
 
 
 template <int dim, int spacedim>
-void DataOutInterface<dim,spacedim>::set_flags (const GnuplotFlags &flags) 
+void
+DataOutInterface<dim,spacedim>::set_flags (const GnuplotFlags &flags) 
 {
   gnuplot_flags = flags;
 }
@@ -4111,7 +4128,8 @@ void DataOutInterface<dim,spacedim>::set_flags (const GnuplotFlags &flags)
 
 
 template <int dim, int spacedim>
-void DataOutInterface<dim,spacedim>::set_flags (const PovrayFlags &flags) 
+void
+DataOutInterface<dim,spacedim>::set_flags (const PovrayFlags &flags) 
 {
   povray_flags = flags;
 }
@@ -4119,7 +4137,8 @@ void DataOutInterface<dim,spacedim>::set_flags (const PovrayFlags &flags)
 
 
 template <int dim, int spacedim>
-void DataOutInterface<dim,spacedim>::set_flags (const EpsFlags &flags) 
+void
+DataOutInterface<dim,spacedim>::set_flags (const EpsFlags &flags) 
 {
   eps_flags = flags;
 }
@@ -4127,7 +4146,8 @@ void DataOutInterface<dim,spacedim>::set_flags (const EpsFlags &flags)
 
 
 template <int dim, int spacedim>
-void DataOutInterface<dim,spacedim>::set_flags (const GmvFlags &flags) 
+void
+DataOutInterface<dim,spacedim>::set_flags (const GmvFlags &flags) 
 {
   gmv_flags = flags;
 }
@@ -4135,7 +4155,8 @@ void DataOutInterface<dim,spacedim>::set_flags (const GmvFlags &flags)
 
 
 template <int dim, int spacedim>
-void DataOutInterface<dim,spacedim>::set_flags (const TecplotFlags &flags) 
+void
+DataOutInterface<dim,spacedim>::set_flags (const TecplotFlags &flags) 
 {
   tecplot_flags = flags;
 }
@@ -4143,7 +4164,8 @@ void DataOutInterface<dim,spacedim>::set_flags (const TecplotFlags &flags)
 
 
 template <int dim, int spacedim>
-void DataOutInterface<dim,spacedim>::set_flags (const VtkFlags &flags) 
+void
+DataOutInterface<dim,spacedim>::set_flags (const VtkFlags &flags) 
 {
   vtk_flags = flags;
 }
@@ -4152,7 +4174,8 @@ void DataOutInterface<dim,spacedim>::set_flags (const VtkFlags &flags)
 
 template <int dim, int spacedim>
 std::string
-DataOutInterface<dim,spacedim>::default_suffix (const OutputFormat output_format_) const
+DataOutInterface<dim,spacedim>::
+default_suffix (const OutputFormat output_format_) const
 {
   OutputFormat output_format = output_format_;
   if (output_format == default_format)
@@ -4197,7 +4220,8 @@ DataOutInterface<dim,spacedim>::default_suffix (const OutputFormat output_format
 
 template <int dim, int spacedim>
 typename DataOutInterface<dim,spacedim>::OutputFormat
-DataOutInterface<dim,spacedim>::parse_output_format (const std::string &format_name)
+DataOutInterface<dim,spacedim>::
+parse_output_format (const std::string &format_name)
 {
   if (format_name == "dx")
     return dx;
@@ -4235,7 +4259,8 @@ DataOutInterface<dim,spacedim>::parse_output_format (const std::string &format_n
 
 
 template <int dim, int spacedim>
-std::string DataOutInterface<dim,spacedim>::get_output_format_names ()
+std::string
+DataOutInterface<dim,spacedim>::get_output_format_names ()
 {
   return "dx|ucd|gnuplot|povray|eps|gmv|tecplot|vtk";
 }
@@ -4243,7 +4268,8 @@ std::string DataOutInterface<dim,spacedim>::get_output_format_names ()
 
 
 template <int dim, int spacedim>
-void DataOutInterface<dim,spacedim>::declare_parameters (ParameterHandler &prm) 
+void
+DataOutInterface<dim,spacedim>::declare_parameters (ParameterHandler &prm) 
 {
   prm.declare_entry ("Output format", "gnuplot",
                     Patterns::Selection (get_output_format_names ()));
@@ -4284,7 +4310,8 @@ void DataOutInterface<dim,spacedim>::declare_parameters (ParameterHandler &prm)
 
 
 template <int dim, int spacedim>
-void DataOutInterface<dim,spacedim>::parse_parameters (ParameterHandler &prm) 
+void
+DataOutInterface<dim,spacedim>::parse_parameters (ParameterHandler &prm) 
 {
   const std::string& output_name = prm.get ("Output format");
   default_fmt = parse_output_format (output_name);
@@ -4323,6 +4350,7 @@ void DataOutInterface<dim,spacedim>::parse_parameters (ParameterHandler &prm)
 }
 
 
+
 template <int dim, int spacedim>
 unsigned int 
 DataOutInterface<dim,spacedim>::memory_consumption () const
index da7bcf9473530d986d25d002bdd35bb2fd64d7f4..cd706d3c06833a95d8f9679b12e6d35b095f37ac 100644 (file)
 
 
 
-ExceptionBase::ExceptionBase () :
+ExceptionBase::ExceptionBase ()
+                :
                file(""), line(0), function(""), cond(""), exc("")
 {}
 
 
 
 ExceptionBase::ExceptionBase (const char* f, const int l, const char *func,
-                             const char* c, const char *e) :
+                             const char* c, const char *e)
+                :
                file(f), line(l), function(func), cond(c), exc(e)
 {}
 
index fd46e617930e4886021727d3678df036f1dcf0b9..fdad456fd2f34ce9300fdd164299840a8a90d2c8 100644 (file)
@@ -20,7 +20,8 @@
 
 template <int dim>
 Function<dim>::Function (const unsigned int n_components,
-                        const double       initial_time) :
+                        const double       initial_time)
+                :
                FunctionTime(initial_time),
                n_components(n_components)
 {}
@@ -200,7 +201,8 @@ Function<dim>::memory_consumption () const
 //------------------------------------------------------------//
 
 template <int dim>
-ZeroFunction<dim>::ZeroFunction (const unsigned int n_components) :
+ZeroFunction<dim>::ZeroFunction (const unsigned int n_components)
+                :
                Function<dim> (n_components)
 {}
 
@@ -309,16 +311,19 @@ void ZeroFunction<dim>::vector_gradient_list (const std::vector<Point<dim> >
 
 template <int dim>
 ConstantFunction<dim>::ConstantFunction (const double value,
-                                        const unsigned int n_components) :
+                                        const unsigned int n_components)
+                :
                ZeroFunction<dim> (n_components),
   function_value    (value)
 {}
 
 
+
 template <int dim>
 ConstantFunction<dim>::~ConstantFunction () {}
 
 
+
 template <int dim>
 double ConstantFunction<dim>::value (const Point<dim> &,
                                     const unsigned int) const
@@ -327,6 +332,7 @@ double ConstantFunction<dim>::value (const Point<dim> &,
 }
 
 
+
 template <int dim>
 void ConstantFunction<dim>::vector_value (const Point<dim> &,
                                          Vector<double>   &return_value) const
@@ -338,6 +344,7 @@ void ConstantFunction<dim>::vector_value (const Point<dim> &,
 }
 
 
+
 template <int dim>
 void ConstantFunction<dim>::value_list (const std::vector<Point<dim> > &points,
                                        std::vector<double>            &values,
@@ -350,6 +357,7 @@ void ConstantFunction<dim>::value_list (const std::vector<Point<dim> > &points,
 }
 
 
+
 template <int dim>
 void ConstantFunction<dim>::vector_value_list (const std::vector<Point<dim> > &points,
                                               std::vector<Vector<double> >   &values) const
@@ -379,15 +387,17 @@ ConstantFunction<dim>::memory_consumption () const
 //------------------------------------------------------------//
 
 template <int dim>
-ComponentSelectFunction<dim>::ComponentSelectFunction (const unsigned int selected,
-                                                      const double value,
-                                                      const unsigned int n_components)
+ComponentSelectFunction<dim>::
+ComponentSelectFunction (const unsigned int selected,
+                         const double value,
+                         const unsigned int n_components)
                :
                ConstantFunction<dim> (value, n_components),
-  selected(selected)
+                selected(selected)
 {}
 
 
+
 template <int dim>
 void ComponentSelectFunction<dim>::vector_value (const Point<dim> &,
                                                 Vector<double>   &return_value) const
index e7e7555cbadaa1485813da092db781669feef38f..b1892e52d68cb73097488b4386b75832d4d9c4e7 100644 (file)
@@ -1321,8 +1321,8 @@ namespace Functions
                  const std::vector<double>      &weights)
                  :
                  Function<dim> (1),
-    fourier_coefficients (fourier_coefficients),
-    weights (weights)
+                  fourier_coefficients (fourier_coefficients),
+                  weights (weights)
   {
     Assert (fourier_coefficients.size() > 0, ExcInvalidArraySize());
     Assert (fourier_coefficients.size() == weights.size(),
index 9c9e2ae8fc9a2958986b95226f0077643b3f3834..8010871e32c26d2d2b5669b64c6e6f30d63eaea3 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002 by the deal authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
index 365acfae911710e60639d06ea134750ffeb8cc8d..00aadd4aec520925637fd7a7daf90e1015e965f0 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002 by the deal authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -20,15 +20,21 @@ FunctionTime::FunctionTime(const double initial_time)
                time(initial_time)
 {}
 
+
+
 FunctionTime::~FunctionTime()
 {}
 
+
+
 void
 FunctionTime::set_time (const double new_time)
 {
   time = new_time;
 }
 
+
+
 void
 FunctionTime::advance_time (const double delta_t)
 {
index 2e4a77dc7a4cc8f3261daa858d4a9980c26dc381..c3dbdf7d28f6e879bb096ae84456de045d9fa620 100644 (file)
@@ -112,28 +112,28 @@ void LogStream::pop ()
 
 
 void
-LogStream::depth_console(unsigned n)
+LogStream::depth_console (const unsigned n)
 {
   std_depth = n;
 }
 
 
 void
-LogStream::depth_file(unsigned n)
+LogStream::depth_file (const unsigned n)
 {
   file_depth = n;
 }
 
 
 void
-LogStream::log_execution_time (bool flag)
+LogStream::log_execution_time (const bool flag)
 {
   print_utime = flag;
 }
 
 
 void
-LogStream::log_time_differences (bool flag)
+LogStream::log_time_differences (const bool flag)
 {
   diff_utime = flag;
 }
index 37550c2ba3c60163d6bc2fb676bb1b2b931c50cf..5fc7a3e49a26de9ba4084409dc7bff1b0297629c 100644 (file)
@@ -101,7 +101,8 @@ unsigned int MultithreadInfo::get_n_cpus()
 #endif
 
 
-MultithreadInfo::MultithreadInfo () :
+MultithreadInfo::MultithreadInfo ()
+                :
                 n_cpus (get_n_cpus()),
                 n_default_threads (n_cpus)
 {}
index 22f3f6a1e3852e83f0ece7bc6a19e7f6094e0f5b..d2e7729c02f7fac5bfe11529e74f473bbc6ad84c 100644 (file)
@@ -31,6 +31,7 @@
 # include <limits>
 #endif
 
+//TODO[WB]: Remove the "status" flag -- it's only a remnant of the DiffPack past of this class and doesn't really serve any useful purpose any more
 
 namespace Patterns
 {  
@@ -72,7 +73,8 @@ namespace Patterns
   
 
   Integer::Integer (const int lower_bound,
-                   const int upper_bound) :
+                   const int upper_bound)
+                  :
                  lower_bound (lower_bound),
                  upper_bound (upper_bound)
   {}
@@ -157,7 +159,8 @@ namespace Patterns
 #endif
   
   Double::Double (const double lower_bound,
-                 const double upper_bound) :
+                 const double upper_bound)
+                  :
                  lower_bound (lower_bound),
                  upper_bound (upper_bound)
   {}
@@ -495,7 +498,8 @@ namespace Patterns
 
 
 
-  Bool::Bool () :
+  Bool::Bool ()
+                  :
                  Selection ("true|false")
   {}
 
@@ -538,7 +542,8 @@ namespace Patterns
 
 
 
-ParameterHandler::ParameterHandler () :
+ParameterHandler::ParameterHandler ()
+                :
                status(true) {}
 
 
@@ -668,10 +673,12 @@ bool ParameterHandler::declare_entry (const std::string           &entry,
   if (p->entries.find (entry) != p->entries.end())
     return false;
 
+                                   // entry doesn't yet exist, but
+                                   // map::operator[] will create it
   p->entries[entry] = std::make_pair(default_value, pattern.clone());
 
-                                  // check whether default answer matches
-                                  // the pattern
+                                  // check whether default answer
+                                  // matches the pattern
   if (!pattern.match(default_value))
     return false;
 
@@ -738,8 +745,8 @@ const std::string & ParameterHandler::get (const std::string &entry_string) cons
     };
 
 
-// entry exists; now find out whether
-                                  // it was changed:
+                                   // entry exists; now find out
+                                  // whether it was changed:
   Section::EntryType::const_iterator ptr;
   ptr = pc->entries.find (entry_string);
   if (ptr != pc->entries.end())
@@ -840,7 +847,8 @@ void
 ParameterHandler::log_parameters (LogStream &out)
 {
   out.push("parameters");
-                                  // dive recursively into the subsections
+                                  // dive recursively into the
+                                  // subsections
   log_parameters_section (out);
 
   out.pop();
@@ -1261,7 +1269,8 @@ ParameterHandler::Section::memory_consumption () const
 
 
 
-MultipleParameterLoop::MultipleParameterLoop() :
+MultipleParameterLoop::MultipleParameterLoop()
+                :
                n_branches(0)
 {}
 
@@ -1484,7 +1493,8 @@ MultipleParameterLoop::memory_consumption () const
 
 MultipleParameterLoop::Entry::Entry (const std::vector<std::string> &ssp,
                                     const std::string              &Name,
-                                    const std::string              &Value) :
+                                    const std::string              &Value)
+                :
                subsection_path (ssp), entry_name(Name), entry_value(Value)
 {}
 
index 1792e49d270d839d9b443a468076db7f1c27a049..9266f9f3bf4c7dd3d57333f9f71a067156143e5f 100644 (file)
@@ -121,8 +121,8 @@ namespace Polynomials
 
   template <typename number>
   void
-  Polynomial<number>::scale(std::vector<number> &coefficients,
-                            const number         factor)
+  Polynomial<number>::scale (std::vector<number> &coefficients,
+                             const number         factor)
   {
     double f = 1.;
     for (typename std::vector<number>::iterator c = coefficients.begin();
@@ -137,7 +137,7 @@ namespace Polynomials
 
   template <typename number>
   void
-  Polynomial<number>::scale(const number factor)
+  Polynomial<number>::scale (const number factor)
   {
     scale (coefficients, factor);
   }
@@ -146,8 +146,8 @@ namespace Polynomials
 
   template <typename number>
   void
-  Polynomial<number>::multiply(std::vector<number> &coefficients,
-                               const number         factor)
+  Polynomial<number>::multiply (std::vector<number> &coefficients,
+                                const number         factor)
   {
     for (typename std::vector<number>::iterator c = coefficients.begin();
          c != coefficients.end(); ++c)
@@ -230,15 +230,15 @@ namespace Polynomials
   template <typename number>
   template <typename number2>
   void
-  Polynomial<number>::shift(const number2 offset)
+  Polynomial<number>::shift (const number2 offset)
   {
-    shift(coefficients, offset);
+    shift (coefficients, offset);
   }
 
 
   template <typename number>
   void
-  Polynomial<number>::print(std::ostream& out) const
+  Polynomial<number>::print (std::ostream& out) const
   {
     for (int i=degree();i>=0;--i)
       {
index 1debf26e89aafbceeb58242f3f28bdad35d3e52e..0fdfb58753c56882b503a9ba5ea6ab6921f88998 100644 (file)
@@ -34,8 +34,8 @@ PolynomialSpace<dim>::compute_n_pols (const unsigned int n)
 template <>
 void
 PolynomialSpace<1>::
-compute_index(const unsigned int n,
-              unsigned int      (&index)[1]) const
+compute_index (const unsigned int n,
+               unsigned int      (&index)[1]) const
 {
   index[0] = n;
 }
index 49401242e893f1f632ecfd241be50e89df127896..8f303403dd00be2c13c1ce515462ae85c83ab6e0 100644 (file)
@@ -42,14 +42,14 @@ Quadrature<dim>::Quadrature (const unsigned int n_q) :
 
 template <int dim>
 Quadrature<dim>::Quadrature (const std::vector<Point<dim> > &points,
-                            const std::vector<double>               &weights)
+                            const std::vector<double>      &weights)
                :
                n_quadrature_points(points.size()),
                quadrature_points(points),
                weights(weights)
 {
-  Assert(weights.size() == points.size(),
-        ExcDimensionMismatch(weights.size(), points.size()));
+  Assert (weights.size() == points.size(),
+          ExcDimensionMismatch(weights.size(), points.size()));
 }
 
 
@@ -68,7 +68,8 @@ Quadrature<dim>::Quadrature (const std::vector<Point<dim> > &points)
 
 
 template <int dim>
-Quadrature<dim>::Quadrature (const Point<dim> &point):
+Quadrature<dim>::Quadrature (const Point<dim> &point)
+                :
                n_quadrature_points(1),
                quadrature_points(std::vector<Point<dim> > (1, point)),
                weights(std::vector<double> (1, 1.))
@@ -78,7 +79,8 @@ Quadrature<dim>::Quadrature (const Point<dim> &point):
 template <>
 Quadrature<0>::Quadrature (const Quadrature<-1> &,
                           const Quadrature<1> &)
-               : n_quadrature_points (0)
+               :
+                n_quadrature_points (0)
 {
   Assert (false, ExcInternalError());
 }
@@ -87,7 +89,8 @@ Quadrature<0>::Quadrature (const Quadrature<-1> &,
 
 template <>
 Quadrature<1>::Quadrature (const Quadrature<0> &,
-                          const Quadrature<1> &) :
+                          const Quadrature<1> &)
+                :
                 n_quadrature_points (0)
 {
   Assert (false, ExcInternalError());
@@ -97,7 +100,7 @@ Quadrature<1>::Quadrature (const Quadrature<0> &,
 
 template <int dim>
 Quadrature<dim>::Quadrature (const SubQuadrature &q1,
-                            const Quadrature<1>                     &q2)
+                            const Quadrature<1> &q2)
                :
                n_quadrature_points (q1.n_quadrature_points *
                                     q2.n_quadrature_points),
@@ -142,7 +145,8 @@ Quadrature<dim>::~Quadrature ()
 
 
 template <>
-const Point<0> & Quadrature<0>::point (const unsigned int) const
+const Point<0> &
+Quadrature<0>::point (const unsigned int) const
 {
   Assert (false, ExcInternalError());
   static const Point<0> dummy;
@@ -152,7 +156,8 @@ const Point<0> & Quadrature<0>::point (const unsigned int) const
 
 
 template <int dim>
-const Point<dim> & Quadrature<dim>::point (const unsigned int i) const
+const Point<dim> &
+Quadrature<dim>::point (const unsigned int i) const
 {
   Assert (i<n_quadrature_points, ExcIndexRange(i, 0, n_quadrature_points));
   return quadrature_points[i];
@@ -161,7 +166,8 @@ const Point<dim> & Quadrature<dim>::point (const unsigned int i) const
 
 
 template <>
-const std::vector<Point<0> > & Quadrature<0>::get_points () const
+const std::vector<Point<0> > &
+Quadrature<0>::get_points () const
 {
   Assert (false, ExcInternalError());
   return quadrature_points;
@@ -170,7 +176,8 @@ const std::vector<Point<0> > & Quadrature<0>::get_points () const
 
 
 template <int dim>
-const std::vector<Point<dim> > & Quadrature<dim>::get_points () const
+const std::vector<Point<dim> > &
+Quadrature<dim>::get_points () const
 {
   return quadrature_points;
 }
@@ -178,7 +185,8 @@ const std::vector<Point<dim> > & Quadrature<dim>::get_points () const
 
 
 template <>
-double Quadrature<0>::weight (const unsigned int) const
+double
+Quadrature<0>::weight (const unsigned int) const
 {
   Assert (false, ExcInternalError());
   return 0;
@@ -187,7 +195,8 @@ double Quadrature<0>::weight (const unsigned int) const
 
 
 template <int dim>
-double Quadrature<dim>::weight (const unsigned int i) const
+double
+Quadrature<dim>::weight (const unsigned int i) const
 {
   Assert (i<n_quadrature_points, ExcIndexRange(i, 0, n_quadrature_points));
   return weights[i];
@@ -196,7 +205,8 @@ double Quadrature<dim>::weight (const unsigned int i) const
 
 
 template <int dim>
-const std::vector<double> & Quadrature<dim>::get_weights () const
+const std::vector<double> &
+Quadrature<dim>::get_weights () const
 {
   return weights;
 }
@@ -204,7 +214,8 @@ const std::vector<double> & Quadrature<dim>::get_weights () const
 
 
 template <>
-const std::vector<double> & Quadrature<0>::get_weights () const
+const std::vector<double> &
+Quadrature<0>::get_weights () const
 {
   Assert (false, ExcInternalError());
   return weights;
@@ -225,19 +236,21 @@ Quadrature<dim>::memory_consumption () const
 
 
 template <>
-void QProjector<1>::project_to_face (const Quadrature<0> &,
-                                    const unsigned int,
-                                    std::vector<Point<1> > &)
+void
+QProjector<1>::project_to_face (const Quadrature<0> &,
+                                const unsigned int,
+                                std::vector<Point<1> > &)
 {
-  Assert(false, ExcNotImplemented());
+  Assert (false, ExcNotImplemented());
 }
 
 
 
 template <>
-void QProjector<2>::project_to_face (const Quadrature<1>      &quadrature,
-                                    const unsigned int        face_no,
-                                    std::vector<Point<2> >   &q_points)
+void
+QProjector<2>::project_to_face (const Quadrature<1>      &quadrature,
+                                const unsigned int        face_no,
+                                std::vector<Point<2> >   &q_points)
 {
   const unsigned int dim=2;
   Assert (face_no<2*dim, ExcIndexRange (face_no, 0, 2*dim));
@@ -267,9 +280,10 @@ void QProjector<2>::project_to_face (const Quadrature<1>      &quadrature,
 
 
 template <>
-void QProjector<3>::project_to_face (const Quadrature<2>    &quadrature,
-                                    const unsigned int      face_no,
-                                    std::vector<Point<3> > &q_points)
+void
+QProjector<3>::project_to_face (const Quadrature<2>    &quadrature,
+                                const unsigned int      face_no,
+                                std::vector<Point<3> > &q_points)
 {
   const unsigned int dim=3;
   Assert (face_no<2*dim, ExcIndexRange (face_no, 0, 2*dim));
@@ -318,10 +332,11 @@ void QProjector<3>::project_to_face (const Quadrature<2>    &quadrature,
 
 
 template <>
-void QProjector<1>::project_to_subface (const Quadrature<0> &,
-                                       const unsigned int,
-                                       const unsigned int,
-                                       std::vector<Point<1> > &)
+void
+QProjector<1>::project_to_subface (const Quadrature<0> &,
+                                   const unsigned int,
+                                   const unsigned int,
+                                   std::vector<Point<1> > &)
 {
   Assert(false, ExcNotImplemented());
 }
@@ -329,10 +344,11 @@ void QProjector<1>::project_to_subface (const Quadrature<0> &,
 
   
 template <>
-void QProjector<2>::project_to_subface (const Quadrature<1>    &quadrature,
-                                       const unsigned int      face_no,
-                                       const unsigned int      subface_no,
-                                       std::vector<Point<2> > &q_points)
+void
+QProjector<2>::project_to_subface (const Quadrature<1>    &quadrature,
+                                   const unsigned int      face_no,
+                                   const unsigned int      subface_no,
+                                   std::vector<Point<2> > &q_points)
 {
   const unsigned int dim=2;
   Assert (face_no<2*dim, ExcIndexRange (face_no, 0, 2*dim));
@@ -405,10 +421,11 @@ void QProjector<2>::project_to_subface (const Quadrature<1>    &quadrature,
 
 
 template <>
-void QProjector<3>::project_to_subface (const Quadrature<2>    &quadrature,
-                                       const unsigned int      face_no,
-                                       const unsigned int      subface_no,
-                                       std::vector<Point<3> > &q_points)
+void
+QProjector<3>::project_to_subface (const Quadrature<2>    &quadrature,
+                                   const unsigned int      face_no,
+                                   const unsigned int      subface_no,
+                                   std::vector<Point<3> > &q_points)
 {
   const unsigned int dim=3;
   Assert (face_no<2*dim, ExcIndexRange (face_no, 0, 2*dim));
@@ -728,7 +745,8 @@ QIterated<1>::uses_both_endpoints (const Quadrature<1> &base_quadrature)
 
 template <>
 QIterated<1>::QIterated (const Quadrature<1> &base_quadrature,
-                        const unsigned int   n_copies) :
+                        const unsigned int   n_copies)
+                :
                Quadrature<1> (uses_both_endpoints(base_quadrature) ?
                               (base_quadrature.n_quadrature_points-1) * n_copies + 1 :
                               base_quadrature.n_quadrature_points * n_copies) 
@@ -831,7 +849,8 @@ QIterated<1>::QIterated (const Quadrature<1> &base_quadrature,
 // of lower dimensional iterated quadrature formulae
 template <int dim>
 QIterated<dim>::QIterated (const Quadrature<1> &base_quadrature,
-                          const unsigned int   N) :
+                          const unsigned int   N)
+                :
                Quadrature<dim> (QIterated<dim-1>(base_quadrature, N),
                                 QIterated<1>(base_quadrature, N))
 {}
index edc0d64e646281caeb2cb2fa1bd774902c2317f5..86d92a06e2532dfad8cd32e320d1d907b637b151 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002 by the deal authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -30,7 +30,8 @@ number abs (const number a)
 
 template <>
 QGauss<1>::QGauss (const unsigned int n)
-  : Quadrature<1> (n)
+                :
+                Quadrature<1> (n)
 {
   const unsigned int m = (n+1)/2;
   long double z;
@@ -71,7 +72,8 @@ QGauss<1>::QGauss (const unsigned int n)
 
 
 template <>
-QGauss2<1>::QGauss2 () :
+QGauss2<1>::QGauss2 ()
+                :
                Quadrature<1> (2)
 {
                                   // points on [-1,1]
@@ -93,8 +95,10 @@ QGauss2<1>::QGauss2 () :
 }
 
 
+
 template <>
-QGauss3<1>::QGauss3 () :
+QGauss3<1>::QGauss3 ()
+                :
                Quadrature<1> (3)
 {
                                   // points on [-1,1]
@@ -122,8 +126,10 @@ QGauss3<1>::QGauss3 () :
 }
 
 
+
 template <>
-QGauss4<1>::QGauss4 () :
+QGauss4<1>::QGauss4 () 
+                :
                Quadrature<1> (4)
 {
                                   // points on [-1,1]
@@ -156,7 +162,8 @@ QGauss4<1>::QGauss4 () :
 
 
 template <>
-QGauss5<1>::QGauss5 () :
+QGauss5<1>::QGauss5 ()
+                :
                Quadrature<1> (5)
 {
                                   // points on [-1,1]
@@ -192,8 +199,10 @@ QGauss5<1>::QGauss5 () :
 }
 
 
+
 template <>
-QGauss6<1>::QGauss6 () :
+QGauss6<1>::QGauss6 ()
+                :
                Quadrature<1> (6)
 {
                                   // points on [-1,1]
@@ -233,8 +242,10 @@ QGauss6<1>::QGauss6 () :
 }
 
 
+
 template <>
-QGauss7<1>::QGauss7 () :
+QGauss7<1>::QGauss7 ()
+                :
                Quadrature<1> (7)
 {
                                   // points on [-1,1]
@@ -280,7 +291,8 @@ QGauss7<1>::QGauss7 () :
 
 
 template <>
-QMidpoint<1>::QMidpoint () :
+QMidpoint<1>::QMidpoint ()
+                :
                Quadrature<1>(1)
 {
   this->quadrature_points[0] = Point<1>(0.5);
@@ -288,8 +300,10 @@ QMidpoint<1>::QMidpoint () :
 }
 
 
+
 template <>
-QTrapez<1>::QTrapez () :
+QTrapez<1>::QTrapez ()
+                :
                Quadrature<1> (2)
 {
   static const double xpts[] = { 0.0, 1.0 };
@@ -303,8 +317,10 @@ QTrapez<1>::QTrapez () :
 }
 
 
+
 template <>
-QSimpson<1>::QSimpson () :
+QSimpson<1>::QSimpson ()
+                :
                Quadrature<1> (3)
 {
   static const double xpts[] = { 0.0, 0.5, 1.0 };
@@ -318,8 +334,10 @@ QSimpson<1>::QSimpson () :
 }
 
 
+
 template <>
-QMilne<1>::QMilne () :
+QMilne<1>::QMilne ()
+                :
                Quadrature<1> (5)
 {
   static const double xpts[] = { 0.0, .25, .5, .75, 1.0 };
@@ -333,8 +351,10 @@ QMilne<1>::QMilne () :
 }
 
 
+
 template <>
-QWeddle<1>::QWeddle () :
+QWeddle<1>::QWeddle ()
+                :
                Quadrature<1> (7)
 {
   static const double xpts[] = { 0.0, 1./6., 1./3., .5, 2./3., 5./6., 1.0 };
@@ -362,48 +382,89 @@ QGauss<dim>::QGauss (const unsigned int n)
 
 
 template <int dim>
-QGauss2<dim>::QGauss2 () :  Quadrature<dim> (QGauss2<dim-1>(), QGauss2<1>())  {}
+QGauss2<dim>::QGauss2 ()
+                :
+                Quadrature<dim> (QGauss2<dim-1>(), QGauss2<1>())
+{}
+
+
 
 template <int dim>
-QGauss3<dim>::QGauss3 () :
-               Quadrature<dim> (QGauss3<dim-1>(), QGauss3<1>()){}
+QGauss3<dim>::QGauss3 ()
+                :
+               Quadrature<dim> (QGauss3<dim-1>(), QGauss3<1>())
+{}
+
+
 
 template <int dim>
-QGauss4<dim>::QGauss4 () :
-               Quadrature<dim> (QGauss4<dim-1>(), QGauss4<1>()){}
+QGauss4<dim>::QGauss4 ()
+                :
+               Quadrature<dim> (QGauss4<dim-1>(), QGauss4<1>())
+{}
+
+
 
 template <int dim>
-QGauss5<dim>::QGauss5 () :
-               Quadrature<dim> (QGauss5<dim-1>(), QGauss5<1>()){}
+QGauss5<dim>::QGauss5 ()
+                :
+               Quadrature<dim> (QGauss5<dim-1>(), QGauss5<1>())
+{}
+
+
 
 template <int dim>
-QGauss6<dim>::QGauss6 () :
-               Quadrature<dim> (QGauss6<dim-1>(), QGauss6<1>()){}
+QGauss6<dim>::QGauss6 ()
+                :
+               Quadrature<dim> (QGauss6<dim-1>(), QGauss6<1>())
+{}
+
+
 
 template <int dim>
-QGauss7<dim>::QGauss7 () :
-               Quadrature<dim> (QGauss7<dim-1>(), QGauss7<1>()){}
+QGauss7<dim>::QGauss7 ()
+                :
+               Quadrature<dim> (QGauss7<dim-1>(), QGauss7<1>())
+{}
+
 
 
 template <int dim>
-QMidpoint<dim>::QMidpoint () :
-               Quadrature<dim> (QMidpoint<dim-1>(), QMidpoint<1>()){}
+QMidpoint<dim>::QMidpoint ()
+                :
+               Quadrature<dim> (QMidpoint<dim-1>(), QMidpoint<1>())
+{}
+
+
 
 template <int dim>
-QTrapez<dim>::QTrapez () :
-               Quadrature<dim> (QTrapez<dim-1>(), QTrapez<1>()){}
+QTrapez<dim>::QTrapez ()
+                :
+               Quadrature<dim> (QTrapez<dim-1>(), QTrapez<1>())
+{}
+
+
 
 template <int dim>
-QSimpson<dim>::QSimpson () :
-               Quadrature<dim> (QSimpson<dim-1>(), QSimpson<1>()){}
+QSimpson<dim>::QSimpson ()
+                :
+               Quadrature<dim> (QSimpson<dim-1>(), QSimpson<1>())
+{}
+
+
 
 template <int dim>
-QMilne<dim>::QMilne () :
-               Quadrature<dim> (QMilne<dim-1>(), QMilne<1>()){}
+QMilne<dim>::QMilne ()
+                :
+               Quadrature<dim> (QMilne<dim-1>(), QMilne<1>())
+{}
+
 
 template <int dim>
-QWeddle<dim>::QWeddle () :
-               Quadrature<dim> (QWeddle<dim-1>(), QWeddle<1>()){}
+QWeddle<dim>::QWeddle ()
+                :
+               Quadrature<dim> (QWeddle<dim-1>(), QWeddle<1>())
+{}
 
 
 // explicit specialization
index 52162862515f6e165acd901fd6bc8a80ef243cf7..23c64e124044f86f34cb5b3ba992c9a6015691d4 100644 (file)
@@ -84,13 +84,15 @@ ActiveObjectMonitor::deregister_object (const Subscriptor *p)
 
 
 
-Subscriptor::Subscriptor () :
+Subscriptor::Subscriptor ()
+                :
                counter (0),
                object_info (0)
 {}
 
 
-Subscriptor::Subscriptor (const Subscriptor &) :
+Subscriptor::Subscriptor (const Subscriptor &)
+                :
                counter (0),
                object_info (0)
 {}
@@ -134,14 +136,15 @@ void Subscriptor::subscribe () const
 }
 
 
-void Subscriptor::unsubscribe () const {
+void Subscriptor::unsubscribe () const
+{
   Assert (counter>0, ExcNotUsed());
   Threads::ThreadMutex::ScopedLock lock (subscription_lock);
   --counter;
 }
 
 
-unsigned int Subscriptor::n_subscriptions () const 
+unsigned int Subscriptor::n_subscriptions () const
 {
   return counter;
 }
index fd13f4d6cb9a226eb3ec3efb6b41ace1295c5aa8..fe02a13d88a2ddff01985acd17c081477a2d1b2f 100644 (file)
@@ -37,13 +37,19 @@ TableEntryBase::~TableEntryBase ()
 /*---------------------------------------------------------------------*/
 
 template <typename number>
-TableEntry<number>::TableEntry(const number value):
-               val(value)  {}
+TableEntry<number>::TableEntry(const number value)
+                :
+               val(value)
+{}
+
+
 
 template <typename number>
 TableEntry<number>::~TableEntry()
 {}
 
+
+
 template <typename number>
 number TableEntry<number>::value() const
 {
@@ -51,12 +57,15 @@ number TableEntry<number>::value() const
 }
 
 
+
 template <typename number>
 void TableEntry<number>::write_tex (std::ostream &out) const
 {
   out << val;
 }
 
+
+
 template <typename number>
 void TableEntry<number>::write_text (std::ostream &out) const
 {
@@ -66,7 +75,8 @@ void TableEntry<number>::write_text (std::ostream &out) const
 
 /*---------------------------------------------------------------------*/
 
-TableHandler::Column::Column(const std::string &tex_caption):
+TableHandler::Column::Column(const std::string &tex_caption)
+                :
                tex_caption(tex_caption),
                tex_format("c"),
                precision(4),
@@ -75,7 +85,9 @@ TableHandler::Column::Column(const std::string &tex_caption):
 {}
 
 
-TableHandler::Column::Column():
+
+TableHandler::Column::Column()
+                :
                tex_caption(),
                tex_format("c"),
                precision(4),
@@ -84,6 +96,7 @@ TableHandler::Column::Column():
 {}
 
 
+
 TableHandler::Column::~Column()
 {
   for (unsigned int i=0; i<entries.size(); ++i)
@@ -331,7 +344,7 @@ void TableHandler::write_text(std::ostream &out) const
 }
 
 
-void TableHandler::write_tex(std::ofstream &out) const
+void TableHandler::write_tex (std::ofstream &out) const
 {
   AssertThrow (out, ExcIO());
   
@@ -451,7 +464,7 @@ void TableHandler::write_tex(std::ofstream &out) const
 unsigned int TableHandler::n_rows() const
 {
   std::map<std::string, Column>::const_iterator col_iter=columns.begin();
-  unsigned int n=col_iter->second.entries.size();
+  unsigned int n = col_iter->second.entries.size();
   std::string first_name=col_iter->first;
 
   for (++col_iter; col_iter!=columns.end(); ++col_iter)
index fb74015428ec8db905f08dec4313fdf5bdbd9322..b27907abdd446e123cdda66e3557df5a49f40176 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002 by the deal authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
index e12ae0f97c95ca13a71b5e1198ce63a67611cdc6..5c877df8810e512e099eec5f0d92d24832a960b7 100644 (file)
@@ -30,6 +30,7 @@ TensorFunction<rank, dim>::TensorFunction (const double initial_time)
 {}
 
 
+
 template <int rank, int dim>
 TensorFunction<rank, dim>::~TensorFunction ()
 {}
index 3702ea2e67e50836e8896e780f15b33ba9f4bd5b..7708913b298f8fecd62ca2fe30866a5b46a6ccb1 100644 (file)
@@ -37,7 +37,7 @@ namespace Threads
     ThreadMutex  n_existing_threads_mutex;
 
   
-    void register_thread () 
+    void register_thread ()
     {
       ThreadMutex::ScopedLock lock (n_existing_threads_mutex);
       ++n_existing_threads_counter;
@@ -102,8 +102,7 @@ namespace Threads
   unsigned int n_existing_threads () 
   {
     ThreadMutex::ScopedLock lock (internal::n_existing_threads_mutex);
-    const unsigned int n = internal::n_existing_threads_counter;
-    return n;
+    return internal::n_existing_threads_counter;
   }
   
   
index cfb226921e6450ac16ecb6459ccc577946ccdd98..f30e469185b67bbad4c7b2aeb5f2bc4799f6d1d5 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002 by the deal authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -33,7 +33,8 @@ extern "C" {
 
 
 Timer::Timer()
-  : cumulative_time (0.)
+                :
+                cumulative_time (0.)
 {
   start();
 }
index aed95465c9c4b52cfd1e46e61f709b96195b6116..533de086c57edda2c172cd901075f909e9232ce1 100644 (file)
@@ -50,7 +50,7 @@
 #endif
 //
 // last known and checked version is 3.4:
-#if (__GNUC__ > 3) || ((__GNUC__ == 3) && (__GNUC_MINOR__ > 4))
+#if (__GNUC__ > 3) || ((__GNUC__ == 3) && (__GNUC_MINOR__ > 5))
 #  if defined(BOOST_ASSERT_CONFIG)
 #     error "Unknown compiler version - please run the configure tests and report the results"
 #  else
index ec11df6aa12556dacaaddce2cd5b4b0c49d89398..028897130a01800b2da530625e06a402d6444afd 100644 (file)
@@ -84,7 +84,8 @@ class TriaAccessor
     TriaAccessor (const Triangulation<dim> *parent =  0,
                  const int                 level  = -1,
                  const int                 index  = -1,
-                 const AccessorData       *       =  0) :
+                 const AccessorData       *       =  0)
+                    :
                    present_level (level),
                    present_index (index),
                    tria (parent) {};
@@ -124,7 +125,7 @@ class TriaAccessor
                                     /**
                                      *  Same as above.
                                      */
-    TriaAccessor &operator = (const TriaAccessor &);
+    TriaAccessor & operator = (const TriaAccessor &);
 
   protected:
     
@@ -296,8 +297,10 @@ class TriaObjectAccessor :  public TriaAccessor<dim>
     TriaObjectAccessor (const Triangulation<dim> *parent     =  0,
                        const int                 level      = -1,
                        const int                 index      = -1,
-                       const AccessorData       *local_data =  0) :
-                   TriaAccessor<dim> (parent, level, index, local_data) {};
+                       const AccessorData       *local_data =  0)
+                    :
+                   TriaAccessor<dim> (parent, level, index, local_data)
+      {};
 
                                     /**
                                      * Copy the data of a line. Only
@@ -343,7 +346,8 @@ class TriaObjectAccessor :  public TriaAccessor<dim>
                                      *
                                      * Implemented only for @p{celldim>1}.
                                      */
-    TriaIterator<dim,TriaObjectAccessor<1, dim> > line (const unsigned int i) const;
+    TriaIterator<dim,TriaObjectAccessor<1, dim> >
+    line (const unsigned int i) const;
 
                                     /**
                                      * Line index of the @p{i}th
@@ -361,7 +365,8 @@ class TriaObjectAccessor :  public TriaAccessor<dim>
                                      *
                                      * Implemented only for @p{celldim>2}.
                                      */
-    TriaIterator<dim,TriaObjectAccessor<2, dim> > quad (const unsigned int i) const;
+    TriaIterator<dim,TriaObjectAccessor<2, dim> >
+    quad (const unsigned int i) const;
 
                                     /**
                                      * Quad index of the @p{i}th
@@ -583,12 +588,12 @@ class TriaObjectAccessor :  public TriaAccessor<dim>
                                      * you are doing!  This value is
                                      * reserved for another purpose
                                      * and algorithms may not work if
-                                     * boundary cells have have this
+                                     * boundary cells have this
                                      * boundary indicator or if
                                      * interior cells have boundary
                                      * indicators other than 255.
                                      */
-    void set_boundary_indicator (unsigned char) const;
+    void set_boundary_indicator (const unsigned char) const;
 
                                     /**
                                      * Return whether this object is
@@ -751,7 +756,8 @@ class TriaObjectAccessor<0, dim> : public TriaAccessor<dim>
     TriaObjectAccessor (const Triangulation<dim> *parent     =  0,
                        const int                 level      = -1,
                        const int                 index      = -1,
-                       const AccessorData       *local_data =  0) :
+                       const AccessorData       *local_data =  0)
+                    :
                    TriaAccessor<dim> (parent, level, index, local_data)
       {
        Assert (false, ExcInternalError());
@@ -785,8 +791,10 @@ class TriaObjectAccessor<1, dim> :  public TriaAccessor<dim>
     TriaObjectAccessor (const Triangulation<dim> *parent     =  0,
                        const int                 level      = -1,
                        const int                 index      = -1,
-                       const AccessorData       *local_data =  0) :
-                   TriaAccessor<dim> (parent, level, index, local_data) {};
+                       const AccessorData       *local_data =  0)
+                    :
+                   TriaAccessor<dim> (parent, level, index, local_data)
+      {};
 
                                     /**
                                      *  Copy the data of the given
@@ -1024,7 +1032,7 @@ class TriaObjectAccessor<1, dim> :  public TriaAccessor<dim>
                                      * interior cells have boundary
                                      * indicators other than 255.
                                      */
-    void set_boundary_indicator (unsigned char) const;
+    void set_boundary_indicator (const unsigned char) const;
 
                                     /**
                                      * Return whether this line is at
@@ -1196,8 +1204,10 @@ class TriaObjectAccessor<2, dim> :  public TriaAccessor<dim>
     TriaObjectAccessor (const Triangulation<dim> *parent     =  0,
                        const int                 level      = -1,
                        const int                 index      = -1,
-                       const AccessorData       *local_data =  0) :
-                   TriaAccessor<dim> (parent, level, index, local_data) {};
+                       const AccessorData       *local_data =  0)
+                    :
+                   TriaAccessor<dim> (parent, level, index, local_data)
+      {};
 
                                     /**
                                      *  Copy the data of the given quad.
@@ -1456,12 +1466,12 @@ class TriaObjectAccessor<2, dim> :  public TriaAccessor<dim>
                                      * doing, since this value is
                                      * reserved for another purpose
                                      * and algorithms may not work if
-                                     * boundary cells have have this
+                                     * boundary cells have this
                                      * boundary indicator or if
                                      * interior cells have boundary
                                      * indicators other than 255.
                                      */
-    void set_boundary_indicator (unsigned char) const;
+    void set_boundary_indicator (const unsigned char) const;
 
                                     /**
                                      * Return whether this quad is at
@@ -1660,8 +1670,10 @@ class TriaObjectAccessor<3, dim> :  public TriaAccessor<dim>
     TriaObjectAccessor (const Triangulation<dim> *parent     =  0,
                        const int                 level      = -1,
                        const int                 index      = -1,
-                       const AccessorData       *local_data =  0) :
-                   TriaAccessor<dim> (parent, level, index, local_data) {};
+                       const AccessorData       *local_data =  0)
+                    :
+                   TriaAccessor<dim> (parent, level, index, local_data)
+      {};
 
                                     /**
                                      *  Copy the data of the given
@@ -1938,12 +1950,12 @@ class TriaObjectAccessor<3, dim> :  public TriaAccessor<dim>
                                      * doing, since this value is
                                      * reserved for another purpose
                                      * and algorithms may not work if
-                                     * boundary cells have have this
+                                     * boundary cells have this
                                      * boundary indicator or if
                                      * interior cells have boundary
                                      * indicators other than 255.
                                      */
-    void set_boundary_indicator (unsigned char) const;
+    void set_boundary_indicator (const unsigned char) const;
 
                                     /**
                                      * Return whether this hex is at
@@ -2128,8 +2140,10 @@ class CellAccessor :  public TriaObjectAccessor<dim,dim>
     CellAccessor (const Triangulation<dim> *parent     =  0,
                  const int                 level      = -1,
                  const int                 index      = -1,
-                 const AccessorData       *local_data =  0) :
-                   TriaObjectAccessor<dim,dim> (parent, level, index, local_data) {};
+                 const AccessorData       *local_data =  0)
+                    :
+                   TriaObjectAccessor<dim,dim> (parent, level, index, local_data)
+      {};
 
                                     /**
                                      *  Return a pointer to the
index f96000286855110fdd0c5fa85249f502bd9920db..c8762ebdc7b28b2d3e3dc271512651e74332a4e0 100644 (file)
@@ -181,7 +181,7 @@ unsigned char TriaObjectAccessor<1, dim>::boundary_indicator () const
 
 
 template <int dim>
-void TriaObjectAccessor<1, dim>::set_boundary_indicator (unsigned char boundary_ind) const
+void TriaObjectAccessor<1, dim>::set_boundary_indicator (const unsigned char boundary_ind) const
 {
   Assert (dim>=2, typename TriaAccessor<dim>::ExcNotUsefulForThisDimension());
   Assert (used(), typename TriaAccessor<dim>::ExcCellNotUsed());
@@ -410,7 +410,7 @@ unsigned char TriaObjectAccessor<2, dim>::boundary_indicator () const
 
 
 template <int dim>
-void TriaObjectAccessor<2, dim>::set_boundary_indicator (unsigned char boundary_ind) const
+void TriaObjectAccessor<2, dim>::set_boundary_indicator (const unsigned char boundary_ind) const
 {
   Assert (dim>=3, typename TriaAccessor<dim>::ExcNotUsefulForThisDimension());
   Assert (used(), typename TriaAccessor<dim>::ExcCellNotUsed());
@@ -830,7 +830,7 @@ unsigned char TriaObjectAccessor<3, dim>::boundary_indicator () const
 
 
 template <int dim>
-void TriaObjectAccessor<3, dim>::set_boundary_indicator (unsigned char boundary_ind) const
+void TriaObjectAccessor<3, dim>::set_boundary_indicator (const unsigned char boundary_ind) const
 {
   Assert (dim<4, typename TriaAccessor<dim>::ExcNotUsefulForThisDimension());
   Assert (used(), typename TriaAccessor<dim>::ExcCellNotUsed());
index 186ce42a925aeea5a94358ce0953748bf7e90679..1e4f3f89a8861e7585ddaf12960c8e806a94b62d 100644 (file)
@@ -14,7 +14,7 @@ target = $(basename $(shell echo step-*.cc))
 # run-time checking of parameters and internal states is performed, so
 # you should set this value to `on' while you develop your program,
 # and to `off' when running production computations.
-debug-mode = on
+debug-mode = off
 
 
 # As third field, we need to give the path to the top-level deal.II
index c7018adb45c6f3af6f94c2a3932b29d77135544a..1f4c9743522cfcee6368e0d895705cde34baa8dc 100644 (file)
@@ -206,7 +206,7 @@ BlockIndices::reinit (const std::vector<unsigned int> &n)
 
 inline
 std::pair<unsigned int,unsigned int>
-BlockIndices::global_to_local (const unsigned int i) const
+BlockIndices::global_to_local (const unsigned int i) const throw()
 {
   Assert (i<total_size(), ExcIndexRange(i, 0, total_size()));
 
@@ -221,7 +221,7 @@ BlockIndices::global_to_local (const unsigned int i) const
 inline
 unsigned int
 BlockIndices::local_to_global (const unsigned int block,
-                              const unsigned int index) const
+                              const unsigned int index) const throw()
 {
   Assert (block < n_blocks, ExcIndexRange(block, 0, n_blocks));
   Assert (index < start_indices[block+1]-start_indices[block],
@@ -233,7 +233,7 @@ BlockIndices::local_to_global (const unsigned int block,
 
 inline
 unsigned int
-BlockIndices::size () const
+BlockIndices::size () const throw()
 {
   return n_blocks;
 }
@@ -242,7 +242,7 @@ BlockIndices::size () const
 
 inline
 unsigned int
-BlockIndices::total_size () const
+BlockIndices::total_size () const throw()
 {
   return start_indices[n_blocks];
 }
@@ -251,7 +251,7 @@ BlockIndices::total_size () const
 
 inline
 unsigned int
-BlockIndices::block_size (const unsigned int block) const
+BlockIndices::block_size (const unsigned int block) const throw()
 {
   Assert (block < n_blocks, ExcIndexRange(block, 0, n_blocks));
   return start_indices[block+1]-start_indices[block];
@@ -272,7 +272,7 @@ BlockIndices::operator = (const BlockIndices &b)
 
 inline
 bool
-BlockIndices::operator == (const BlockIndices &b) const
+BlockIndices::operator == (const BlockIndices &b) const throw()
 {
   if (n_blocks != b.n_blocks)
     return false;
@@ -288,7 +288,7 @@ BlockIndices::operator == (const BlockIndices &b) const
 
 inline
 void
-BlockIndices::swap (BlockIndices &b)
+BlockIndices::swap (BlockIndices &b) throw()
 {
   Assert (n_blocks == b.n_blocks,
          ExcDimensionMismatch(n_blocks, b.n_blocks));  
@@ -301,7 +301,7 @@ BlockIndices::swap (BlockIndices &b)
 
 inline
 unsigned int
-BlockIndices::memory_consumption () const
+BlockIndices::memory_consumption () const throw()
 {
   return (sizeof(*this) + 
          start_indices.size() * sizeof(start_indices[0]));
@@ -320,7 +320,7 @@ BlockIndices::memory_consumption () const
  * @author Wolfgang Bangerth, 2000
  */
 inline
-void swap (BlockIndices &u, BlockIndices &v)
+void swap (BlockIndices &u, BlockIndices &v) throw()
 {
   u.swap (v);
 }
index ca309468293599eadbab762ef9c100ec3e989b1b..9bd5bcecc6c7fc35ed6768e89a2a079d448b7ffc 100644 (file)
 template <class MATRIX>
 class BlockDiagonalMatrix : public Subscriptor
 {
-public:
-                                  /**
-                                   * Constructor for an @p{n_blocks}
-                                   * by @p{n_blocks} matrix with
-                                   * diagonal blocks @p{M}.
-                                   */
-  BlockDiagonalMatrix (const MATRIX& M,
-                      unsigned int n_blocks);
-
-                                  /**
-                                   * Matrix-vector-multiplication.
-                                   */
-  template <typename number1, typename number2>
-  void vmult (BlockVector<number1>& dst,
-             const BlockVector<number2>& src) const;
+  public:
+                                     /**
+                                      * Constructor for an @p{n_blocks}
+                                      * by @p{n_blocks} matrix with
+                                      * diagonal blocks @p{M}.
+                                      */
+    BlockDiagonalMatrix (const MATRIX       &M,
+                         const unsigned int  n_blocks);
+
+                                     /**
+                                      * Matrix-vector-multiplication.
+                                      */
+    template <typename number1, typename number2>
+    void vmult (BlockVector<number1>& dst,
+                const BlockVector<number2>& src) const;
   
-                                  /**
-                                   * Transposed matrix-vector-multiplication.
-                                   */
-  template <typename number1, typename number2>
-  void Tvmult (BlockVector<number1>& dst,
-              const BlockVector<number2>& src) const;
-private:
-                                  /**
-                                   * Number of blocks.
-                                   */
-  unsigned int num_blocks;
-                                  /**
-                                   * Diagonal entry.
-                                   */
-  SmartPointer<const MATRIX> matrix;
+                                     /**
+                                      * Transposed matrix-vector-multiplication.
+                                      */
+    template <typename number1, typename number2>
+    void Tvmult (BlockVector<number1>& dst,
+                 const BlockVector<number2>& src) const;
+  private:
+                                     /**
+                                      * Number of blocks.
+                                      */
+    unsigned int num_blocks;
+
+                                     /**
+                                      * Diagonal entry.
+                                      */
+    SmartPointer<const MATRIX> matrix;
 };
 
 
@@ -73,9 +74,10 @@ private:
 
 template <class MATRIX>
 BlockDiagonalMatrix<MATRIX>::BlockDiagonalMatrix (const MATRIX& M,
-                                                 unsigned int num_blocks)
-  : num_blocks (num_blocks),
-    matrix(&M)
+                                                 const unsigned int num_blocks)
+                :
+                num_blocks (num_blocks),
+                matrix(&M)
 {}
 
 
@@ -83,7 +85,7 @@ template <class MATRIX>
 template <typename number1, typename number2>
 void
 BlockDiagonalMatrix<MATRIX>::vmult (BlockVector<number1>& dst,
-                                    const BlockVector<number2>& src) const
+                                    const BlockVector<number2>& src) const
 {
   Assert (dst.n_blocks()==num_blocks,
          ExcDimensionMismatch(dst.n_blocks(),num_blocks));
index c3b329c849b8c3ed84dde3fa4b9c81efe6c700c7..d68db16ec2c4330aee6a9f1c5379d8ff17f8de1b 100644 (file)
@@ -79,9 +79,9 @@ class BlockSparseMatrix : public Subscriptor
                                          * access, a const matrix
                                          * pointer is sufficient.
                                          */
-       Accessor (const BlockSparseMatrix<number>*,
-                 unsigned int row,
-                 unsigned short index);
+       Accessor (const BlockSparseMatrix<number> *m,
+                 const unsigned int               row,
+                 const unsigned short             index);
        
                                         /**
                                          * Row number of the element
@@ -744,17 +744,18 @@ class BlockSparseMatrix : public Subscriptor
 
 template <typename number>
 inline
-BlockSparseMatrix<number>::Accessor::Accessor (
-  const BlockSparseMatrix<number>* matrix,
-  unsigned int r,
-  unsigned short i)
-               : matrix(matrix),
-                 base_iterator(matrix->block(0,0).begin()),
-                 block_row(0),
-                 row_start(0),
-                 block_col(0),
-                 col_start(0),
-                 a_index(0)
+BlockSparseMatrix<number>::Accessor::
+Accessor (const BlockSparseMatrix<number> *matrix,
+          const unsigned int               r,
+          const unsigned short             i)
+               :
+                matrix(matrix),
+                base_iterator(matrix->block(0,0).begin()),
+               block_row(0),
+               row_start(0),
+               block_col(0),
+               col_start(0),
+               a_index(0)
 {
   Assert (i==0, ExcNotImplemented());
 
@@ -816,11 +817,12 @@ BlockSparseMatrix<number>::Accessor::value () const
 
 template <typename number>
 inline
-BlockSparseMatrix<number>::const_iterator::const_iterator(
-  const BlockSparseMatrix<number>* m,
-  unsigned int r,
-  unsigned short i)
-               : BlockSparseMatrix<number>::Accessor(m, r, i)
+BlockSparseMatrix<number>::const_iterator::
+const_iterator(const BlockSparseMatrix<number>* m,
+               unsigned int r,
+               unsigned short i)
+               :
+                BlockSparseMatrix<number>::Accessor(m, r, i)
 {}
 
 
@@ -919,7 +921,8 @@ BlockSparseMatrix<number>::const_iterator::operator-> () const
 template <typename number>
 inline
 bool
-BlockSparseMatrix<number>::const_iterator::operator == (const const_iterator& i) const
+BlockSparseMatrix<number>::const_iterator::
+operator == (const const_iterator& i) const
 {
   if (this->matrix != i->matrix)
     return false;
@@ -936,7 +939,8 @@ BlockSparseMatrix<number>::const_iterator::operator == (const const_iterator& i)
 template <typename number>
 inline
 bool
-BlockSparseMatrix<number>::const_iterator::operator != (const const_iterator& i) const
+BlockSparseMatrix<number>::const_iterator::
+operator != (const const_iterator& i) const
 {
   return !(*this == i);
 }
@@ -946,7 +950,8 @@ BlockSparseMatrix<number>::const_iterator::operator != (const const_iterator& i)
 template <typename number>
 inline
 bool
-BlockSparseMatrix<number>::const_iterator::operator < (const const_iterator& i) const
+BlockSparseMatrix<number>::const_iterator::
+operator < (const const_iterator& i) const
 {
   if (this->block_row<i->block_row)
     return true;
@@ -1357,9 +1362,9 @@ template <typename number>
 template <typename somenumber>
 void
 BlockSparseMatrix<number>::
-precondition_Jacobi (BlockVector<somenumber>          &dst,
+precondition_Jacobi (BlockVector<somenumber>       &dst,
                     const BlockVector<somenumber> &src,
-                    const number                           omega) const
+                    const number                  omega) const
 {
   Assert (rows == columns, ExcMatrixNotBlockSquare());
   Assert (dst.n_blocks() == rows,
@@ -1374,6 +1379,7 @@ precondition_Jacobi (BlockVector<somenumber>          &dst,
 }
 
 
+
 template <typename number>
 inline
 typename BlockSparseMatrix<number>::const_iterator
@@ -1382,6 +1388,8 @@ BlockSparseMatrix<number>::begin () const
   return const_iterator(this, 0, 0);
 }
 
+
+
 template <typename number>
 inline
 typename BlockSparseMatrix<number>::const_iterator
@@ -1390,6 +1398,8 @@ BlockSparseMatrix<number>::end () const
   return const_iterator(this, m(), 0);
 }
 
+
+
 template <typename number>
 inline
 typename BlockSparseMatrix<number>::const_iterator
@@ -1399,6 +1409,8 @@ BlockSparseMatrix<number>::begin (unsigned int r) const
   return const_iterator(this, r, 0);
 }
 
+
+
 template <typename number>
 inline
 typename BlockSparseMatrix<number>::const_iterator
index a3241e71b69cb8c85d164a83ed8169f9bc020bd8..2cff410b1787ded50f8194890910f76a2b13c351 100644 (file)
@@ -127,7 +127,7 @@ class BlockSparseMatrixEZ : public Subscriptor
                                      */
     const SparseMatrixEZ<Number>&
     block (const unsigned int row,
-          const unsigned int column) const;    
+          const unsigned int column) const;
 
                                     /**
                                      * Return the number of blocks in a
@@ -431,11 +431,9 @@ BlockSparseMatrixEZ<Number>::vmult (BlockVector<somenumber>       &dst,
   dst = 0.;
   
   for (unsigned int row=0; row<n_block_rows(); ++row)
-    {
-      for (unsigned int col=0; col<n_block_cols(); ++col)
-       block(row,col).vmult_add (dst.block(row),
-                                 src.block(col));
-    };
+    for (unsigned int col=0; col<n_block_cols(); ++col)
+      block(row,col).vmult_add (dst.block(row),
+                                src.block(col));
 }
 
 
@@ -443,9 +441,9 @@ BlockSparseMatrixEZ<Number>::vmult (BlockVector<somenumber>       &dst,
 template <typename Number>
 template <typename somenumber>
 void
-BlockSparseMatrixEZ<Number>::vmult_add (
-  BlockVector<somenumber>       &dst,
-  const BlockVector<somenumber> &src) const
+BlockSparseMatrixEZ<Number>::
+vmult_add (BlockVector<somenumber>       &dst,
+           const BlockVector<somenumber> &src) const
 {
   Assert (dst.n_blocks() == n_block_rows(),
          ExcDimensionMismatch(dst.n_blocks(), n_block_rows()));
@@ -453,11 +451,9 @@ BlockSparseMatrixEZ<Number>::vmult_add (
          ExcDimensionMismatch(src.n_blocks(), n_block_cols()));
 
   for (unsigned int row=0; row<n_block_rows(); ++row)
-    {
-      for (unsigned int col=0; col<n_block_cols(); ++col)
-       block(row,col).vmult_add (dst.block(row),
-                                 src.block(col));
-    };
+    for (unsigned int col=0; col<n_block_cols(); ++col)
+      block(row,col).vmult_add (dst.block(row),
+                                src.block(col));
 }
 
 
@@ -466,9 +462,9 @@ BlockSparseMatrixEZ<Number>::vmult_add (
 template <typename Number>
 template <typename somenumber>
 void
-BlockSparseMatrixEZ<Number>::Tvmult (
-  BlockVector<somenumber>   &dst,
-  const BlockVector<somenumber>& src) const
+BlockSparseMatrixEZ<Number>::
+Tvmult (BlockVector<somenumber>       &dst,
+        const BlockVector<somenumber> &src) const
 {
   Assert (dst.n_blocks() == n_block_cols(),
          ExcDimensionMismatch(dst.n_blocks(), n_block_cols()));
@@ -478,11 +474,9 @@ BlockSparseMatrixEZ<Number>::Tvmult (
   dst = 0.;
   
   for (unsigned int row=0; row<n_block_rows(); ++row)
-    {
-      for (unsigned int col=0; col<n_block_cols(); ++col)
-       block(row,col).Tvmult_add (dst.block(col),
-                                  src.block(row));
-    };
+    for (unsigned int col=0; col<n_block_cols(); ++col)
+      block(row,col).Tvmult_add (dst.block(col),
+                                 src.block(row));
 }
 
 
@@ -490,9 +484,9 @@ BlockSparseMatrixEZ<Number>::Tvmult (
 template <typename Number>
 template <typename somenumber>
 void
-BlockSparseMatrixEZ<Number>::Tvmult_add (
-  BlockVector<somenumber>    &dst,
-  const BlockVector<somenumber> &src) const
+BlockSparseMatrixEZ<Number>::
+Tvmult_add (BlockVector<somenumber>       &dst,
+            const BlockVector<somenumber> &src) const
 {
   Assert (dst.n_blocks() == n_block_cols(),
          ExcDimensionMismatch(dst.n_blocks(), n_block_cols()));
@@ -512,7 +506,7 @@ template <typename number>
 template <class STREAM>
 inline
 void
-BlockSparseMatrixEZ<number>::print_statistics(STREAM& out, bool full)
+BlockSparseMatrixEZ<number>::print_statistics (STREAM& out, bool full)
 {
   unsigned int used_total = 0;
   unsigned int allocated_total = 0;
index d741be8839632df8849eab7d6a4d90dd10d418de..527de59e91a6c0f03d3bb7f9c85ed78baad89eae 100644 (file)
@@ -505,7 +505,7 @@ template <class SparsityPatternBase>
 inline
 SparsityPatternBase &
 BlockSparsityPatternBase<SparsityPatternBase>::block (const unsigned int row,
-                                                 const unsigned int column)
+                                                      const unsigned int column)
 {
   Assert (row<rows, ExcIndexRange(row,0,rows));
   Assert (column<columns, ExcIndexRange(column,0,columns));
@@ -518,7 +518,7 @@ template <class SparsityPatternBase>
 inline
 const SparsityPatternBase &
 BlockSparsityPatternBase<SparsityPatternBase>::block (const unsigned int row,
-                            const unsigned int column) const
+                                                      const unsigned int column) const
 {
   Assert (row<rows, ExcIndexRange(row,0,rows));
   Assert (column<columns, ExcIndexRange(column,0,columns));
index 2b5cdf9fadaceb4dbab5d6f7ab2a5b2de24c798f..2da069eae754c9103431c77724b59beafacce791 100644 (file)
@@ -198,7 +198,7 @@ namespace internal
         public std::iterator<std::random_access_iterator_tag,
                              typename Types<number,constness>::NumberType>
 #else
-    random_access_iterator<typename Types<number,constness>::NumberType,int>
+        random_access_iterator<typename Types<number,constness>::NumberType,int>
 #endif      
     {
       private:
@@ -592,6 +592,8 @@ namespace internal
   }  // namespace BlockVectorIterators
 }  // namespace internal
 
+
+
 /**
  * A vector composed of several blocks each representing a vector of
  * its own.
@@ -1306,7 +1308,7 @@ Number& BlockVector<Number>::operator() (const unsigned int i)
 template <typename Number>
 inline
 BlockVector<Number> &
-BlockVector<Number>::operator *= (const Number factor) 
+BlockVector<Number>::operator *= (const Number factor)
 {
   scale (factor);
   return *this;
@@ -1317,7 +1319,7 @@ BlockVector<Number>::operator *= (const Number factor)
 template <typename Number>
 inline
 BlockVector<Number> &
-BlockVector<Number>::operator /= (const Number factor) 
+BlockVector<Number>::operator /= (const Number factor)
 {
   scale (1./factor);
   return *this;
@@ -1551,7 +1553,8 @@ namespace internal
     template <typename number, bool constness>
     inline
     bool
-    Iterator<number,constness>::operator == (const Iterator &i) const
+    Iterator<number,constness>::
+    operator == (const Iterator &i) const
     {
       Assert (parent == i.parent, ExcPointerToDifferentVectors());
 
@@ -1563,7 +1566,8 @@ namespace internal
     template <typename number, bool constness>
     inline
     bool
-    Iterator<number,constness>::operator == (const InverseConstnessIterator &i) const
+    Iterator<number,constness>::
+    operator == (const InverseConstnessIterator &i) const
     {
       Assert (parent == i.parent, ExcPointerToDifferentVectors());
 
@@ -1575,7 +1579,8 @@ namespace internal
     template <typename number, bool constness>
     inline
     bool
-    Iterator<number,constness>::operator != (const Iterator &i) const
+    Iterator<number,constness>::
+    operator != (const Iterator &i) const
     {
       Assert (parent == i.parent, ExcPointerToDifferentVectors());
 
@@ -1587,7 +1592,8 @@ namespace internal
     template <typename number, bool constness>
     inline
     bool
-    Iterator<number,constness>::operator != (const InverseConstnessIterator &i) const
+    Iterator<number,constness>::
+    operator != (const InverseConstnessIterator &i) const
     {
       Assert (parent == i.parent, ExcPointerToDifferentVectors());
 
@@ -1599,7 +1605,8 @@ namespace internal
     template <typename number, bool constness>
     inline
     bool
-    Iterator<number,constness>::operator < (const Iterator &i) const
+    Iterator<number,constness>::
+    operator < (const Iterator &i) const
     {
       Assert (parent == i.parent, ExcPointerToDifferentVectors());
 
@@ -1611,7 +1618,8 @@ namespace internal
     template <typename number, bool constness>
     inline
     bool
-    Iterator<number,constness>::operator < (const InverseConstnessIterator &i) const
+    Iterator<number,constness>::
+    operator < (const InverseConstnessIterator &i) const
     {
       Assert (parent == i.parent, ExcPointerToDifferentVectors());
 
@@ -1623,7 +1631,8 @@ namespace internal
     template <typename number, bool constness>
     inline
     bool
-    Iterator<number,constness>::operator <= (const Iterator &i) const
+    Iterator<number,constness>::
+    operator <= (const Iterator &i) const
     {
       Assert (parent == i.parent, ExcPointerToDifferentVectors());
 
@@ -1635,7 +1644,8 @@ namespace internal
     template <typename number, bool constness>
     inline
     bool
-    Iterator<number,constness>::operator <= (const InverseConstnessIterator &i) const
+    Iterator<number,constness>::
+    operator <= (const InverseConstnessIterator &i) const
     {
       Assert (parent == i.parent, ExcPointerToDifferentVectors());
 
@@ -1647,7 +1657,8 @@ namespace internal
     template <typename number, bool constness>
     inline
     bool
-    Iterator<number,constness>::operator > (const Iterator &i) const
+    Iterator<number,constness>::
+    operator > (const Iterator &i) const
     {
       Assert (parent == i.parent, ExcPointerToDifferentVectors());
 
@@ -1659,7 +1670,8 @@ namespace internal
     template <typename number, bool constness>
     inline
     bool
-    Iterator<number,constness>::operator > (const InverseConstnessIterator &i) const
+    Iterator<number,constness>::
+    operator > (const InverseConstnessIterator &i) const
     {
       Assert (parent == i.parent, ExcPointerToDifferentVectors());
 
@@ -1671,7 +1683,8 @@ namespace internal
     template <typename number, bool constness>
     inline
     bool
-    Iterator<number,constness>::operator >= (const Iterator &i) const
+    Iterator<number,constness>::
+    operator >= (const Iterator &i) const
     {
       Assert (parent == i.parent, ExcPointerToDifferentVectors());
 
@@ -1683,7 +1696,8 @@ namespace internal
     template <typename number, bool constness>
     inline
     bool
-    Iterator<number,constness>::operator >= (const InverseConstnessIterator &i) const
+    Iterator<number,constness>::
+    operator >= (const InverseConstnessIterator &i) const
     {
       Assert (parent == i.parent, ExcPointerToDifferentVectors());
 
@@ -1695,7 +1709,8 @@ namespace internal
     template <typename number, bool constness>
     inline
     typename Iterator<number,constness>::difference_type
-    Iterator<number,constness>::operator - (const Iterator &i) const
+    Iterator<number,constness>::
+    operator - (const Iterator &i) const
     {
       Assert (parent == i.parent, ExcPointerToDifferentVectors());
 
@@ -1708,7 +1723,8 @@ namespace internal
     template <typename number, bool constness>
     inline
     typename Iterator<number,constness>::difference_type
-    Iterator<number,constness>::operator - (const InverseConstnessIterator &i) const
+    Iterator<number,constness>::
+    operator - (const InverseConstnessIterator &i) const
     {
       Assert (parent == i.parent, ExcPointerToDifferentVectors());
 
@@ -1721,7 +1737,8 @@ namespace internal
     template <typename number, bool constness>
     inline
     Iterator<number,constness>
-    Iterator<number,constness>::operator + (const difference_type &d) const
+    Iterator<number,constness>::
+    operator + (const difference_type &d) const
     {
                                        // if the index pointed to is
                                        // still within the block we
@@ -1745,7 +1762,8 @@ namespace internal
     template <typename number, bool constness>
     inline
     Iterator<number,constness>
-    Iterator<number,constness>::operator - (const difference_type &d) const
+    Iterator<number,constness>::
+    operator - (const difference_type &d) const
     {
                                        // if the index pointed to is
                                        // still within the block we
@@ -1769,7 +1787,8 @@ namespace internal
     template <typename number, bool constness>
     inline
     Iterator<number,constness> &
-    Iterator<number,constness>::operator += (const difference_type &d)
+    Iterator<number,constness>::
+    operator += (const difference_type &d)
     {
                                        // if the index pointed to is
                                        // still within the block we
@@ -1796,7 +1815,8 @@ namespace internal
     template <typename number, bool constness>
     inline
     Iterator<number,constness> &
-    Iterator<number,constness>::operator -= (const difference_type &d)
+    Iterator<number,constness>::
+    operator -= (const difference_type &d)
     {
                                        // if the index pointed to is
                                        // still within the block we
index 433d5e66891427adb35559208ec66ea9d8639806..e83e986bbb1b3d1102ab4056cdca54b4871b5fe3 100644 (file)
@@ -220,9 +220,10 @@ class CompressedSparsityPattern : public Subscriptor
                                      * Check if a value at a certain
                                      * position may be non-zero.
                                      */
-    bool exists (const unsigned int i, const unsigned int j) const;
+    bool exists (const unsigned int i,
+                 const unsigned int j) const;
     
-                            /**
+                                     /**
                                      * Make the sparsity pattern
                                      * symmetric by adding the
                                      * sparsity pattern of the
index 2fb54ab4ef1cb84e84510ed0185f550862f3c8ce..cb979f41ecb9cb7f2a4a4cc0be40a14bc4eb6afa 100644 (file)
@@ -62,7 +62,8 @@ class EigenPower : private Solver<VECTOR>
                                         /**
                                          * Constructor. Set the shift parameter.
                                          */
-       AdditionalData (const double shift = 0.):
+       AdditionalData (const double shift = 0.)
+                        :
                        shift(shift)
          {}
        
@@ -160,11 +161,12 @@ class EigenInverse : private Solver<VECTOR>
                                          */
        AdditionalData (double relaxation = 1.,
                        unsigned int start_adaption = 6,
-                       bool use_residual = true):
-         relaxation(relaxation),
-         start_adaption(start_adaption),
-         use_residual(use_residual)
-      {}
+                       bool use_residual = true)
+                        :
+                        relaxation(relaxation),
+                        start_adaption(start_adaption),
+                        use_residual(use_residual)
+          {}
        
     };
     
@@ -213,17 +215,20 @@ class EigenInverse : private Solver<VECTOR>
 template <class VECTOR>
 EigenPower<VECTOR>::EigenPower (SolverControl &cn,
                                VectorMemory<VECTOR> &mem,
-                               const AdditionalData &data):
+                               const AdditionalData &data)
+                :
                Solver<VECTOR>(cn, mem),
                additional_data(data)
 {}
 
 
+
 template <class VECTOR>
 EigenPower<VECTOR>::~EigenPower ()
 {}
 
 
+
 template <class VECTOR>
 template <class MATRIX>
 void
@@ -300,17 +305,20 @@ EigenPower<VECTOR>::solve (double       &value,
 template <class VECTOR>
 EigenInverse<VECTOR>::EigenInverse (SolverControl &cn,
                                    VectorMemory<VECTOR> &mem,
-                                   const AdditionalData &data):
+                                   const AdditionalData &data)
+                :
                Solver<VECTOR>(cn, mem),
                additional_data(data)
 {}
 
 
+
 template <class VECTOR>
 EigenInverse<VECTOR>::~EigenInverse ()
 {}
 
 
+
 template <class VECTOR>
 template <class MATRIX>
 void
index fc85475dc3b9d141980d6ce5d64c79a608d19eb3..54318c691abe10e5d1136ffac1adba963835a8ac 100644 (file)
@@ -24,8 +24,7 @@
 
 
 template <class MATRIX, class VECTOR>
-FilteredMatrix<MATRIX,VECTOR>::
-FilteredMatrix () 
+FilteredMatrix<MATRIX,VECTOR>::FilteredMatrix ()
 {}
 
 
@@ -75,7 +74,7 @@ set_referenced_matrix (const MATRIX &m)
 
 template <class MATRIX, class VECTOR>
 void
-FilteredMatrix<MATRIX,VECTOR>::clear_constraints () 
+FilteredMatrix<MATRIX,VECTOR>::clear_constraints ()
 {
                                   // swap vectors to release memory
   std::vector<IndexValuePair> empty;
index 9efe89c77345f443795184eac455fbb040ea227b..1c854d37e1ae7af5fe41be04fbff8f566f363e68 100644 (file)
@@ -744,8 +744,6 @@ class FullMatrix : public Table<2,number>
                    int,
                    << "This function is not implemented for the given"
                    << " matrix dimension " << arg1);
-
-  private:
 };
 
 
@@ -787,9 +785,10 @@ FullMatrix<number>::copy_from (const MATRIX& M)
   reinit (M.m(), M.n());
   typename MATRIX::const_iterator entry;
   const typename MATRIX::const_iterator end = M.end();
-  for (entry = M.begin();entry != end;++entry)
+  for (entry = M.begin(); entry != end; ++entry)
     el(entry->row(), entry->column()) = entry->value();
- }
+}
+
 
 /*----------------------------   fullmatrix.h     ---------------------------*/
 
index 7d15b5c348479d039b694827a6eeddad2bece4d1..1df2423dc2ab43285ce9e20e12859fc244723f82 100644 (file)
@@ -51,7 +51,8 @@ FullMatrix<number>::FullMatrix (const unsigned int m,
 
 
 template <typename number>
-FullMatrix<number>::FullMatrix (const FullMatrix &m) :
+FullMatrix<number>::FullMatrix (const FullMatrix &m)
+                :
                Table<2,number> (m)
 {}
 
@@ -480,7 +481,7 @@ 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)
+                                  const unsigned int j)
 {
   Assert (this->data() != 0, ExcEmptyMatrix());
   
@@ -491,8 +492,8 @@ 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,
-                       const unsigned int k)
+                                  const unsigned int j, const number t,
+                                  const unsigned int k)
 {
   Assert (this->data() != 0, ExcEmptyMatrix());
   
@@ -501,6 +502,7 @@ 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)
 {
@@ -1143,8 +1145,6 @@ FullMatrix<number>::determinant () const
   
   Assert (this->n_cols() == this->n_rows(),
          ExcDimensionMismatch(this->n_cols(), this->n_rows()));
-  Assert ((this->n_cols()>=1) && (this->n_cols()<=3),
-         ExcNotImplemented(this->n_cols()));
   
   switch (this->n_cols()) 
     {
@@ -1160,6 +1160,8 @@ FullMatrix<number>::determinant () const
                      +this->el(2,0)*this->el(0,1)*this->el(1,2)
                      -this->el(2,0)*this->el(0,2)*this->el(1,1));
       default:
+            Assert (false,
+                    ExcNotImplemented(this->n_cols()));
             return 0;
     };
 }
@@ -1178,6 +1180,7 @@ FullMatrix<number>::norm2 () const
 }
 
 
+
 template <typename number>
 number
 FullMatrix<number>::relative_symmetry_norm2 () const
@@ -1198,6 +1201,7 @@ FullMatrix<number>::relative_symmetry_norm2 () const
 }
 
 
+
 template <typename number>
 template <typename number2>
 void
@@ -1435,7 +1439,7 @@ FullMatrix<number>::print_formatted (std::ostream       &out,
 
 template <typename number>
 void
-FullMatrix<number>::gauss_jordan()
+FullMatrix<number>::gauss_jordan ()
 {
   Assert (this->data() != 0, ExcEmptyMatrix());  
   Assert (this->n_cols() == this->n_rows(), ExcNotQuadratic());
@@ -1570,7 +1574,8 @@ FullMatrix<number>::householder(Vector<number2>& src)
 template <typename number>
 template <typename number2>
 double
-FullMatrix<number>::least_squares(Vector<number2>& dst, Vector<number2>& src)
+FullMatrix<number>::least_squares (Vector<number2>& dst,
+                                   Vector<number2>& src)
 {
   Assert (this->data() != 0, ExcEmptyMatrix());
   
index 13a4e4c9c8a295de4734daad9c4b6c5c3482c95d..c6e686b29cbb864c1af555da51645fc6f5a5e9ca 100644 (file)
@@ -94,9 +94,8 @@ class MatrixOut : public DataOutInterface<2,2>
                                      * modify the output of the
                                      * @ref{MatrixOut} class.
                                      */
-    class Options 
+    struct Options 
     {
-      public:
                                         /**
                                          * If @p{true}, only show the
                                          * absolute values of the
@@ -298,7 +297,7 @@ class MatrixOut : public DataOutInterface<2,2>
     static double get_gridpoint_value (const Matrix       &matrix,
                                       const unsigned int  i,
                                       const unsigned int  j,
-                                      const Options      &options);    
+                                      const Options      &options);
 };
 
 
index eae0912c0611296ccbbba27bca0bf126a952eb66..24f4408e7d8fbc8b0b7849e35aa835ab1cb0c462 100644 (file)
@@ -1591,6 +1591,8 @@ SparseMatrix<number>::begin (const unsigned int r) const
   return const_iterator(this, r, 0);
 }
 
+
+
 template <typename number>
 inline
 typename SparseMatrix<number>::const_iterator
index 67c396cb431a2322f749c4b1c227d6706eddf951..ac672ea4bf5da4fd8f30fa5d7724b6db6809f18c 100644 (file)
@@ -41,7 +41,8 @@
 
 
 template <typename number>
-SparseMatrix<number>::SparseMatrix () :
+SparseMatrix<number>::SparseMatrix ()
+                :
                cols(0),
                val(0),
                max_len(0)
@@ -50,7 +51,8 @@ SparseMatrix<number>::SparseMatrix () :
 
 
 template <typename number>
-SparseMatrix<number>::SparseMatrix (const SparseMatrix &m) :
+SparseMatrix<number>::SparseMatrix (const SparseMatrix &m)
+                :
                Subscriptor (m),
                cols(0),
                val(0),
@@ -372,7 +374,8 @@ SparseMatrix<number>::threaded_vmult (Vector<somenumber>       &dst,
 template <typename number>
 template <typename somenumber>
 void
-SparseMatrix<number>::Tvmult (Vector<somenumber>& dst, const Vector<somenumber>& src) const
+SparseMatrix<number>::Tvmult (Vector<somenumber>& dst,
+                              const Vector<somenumber>& src) const
 {
   Assert (val != 0, ExcMatrixNotInitialized());
   Assert (cols != 0, ExcMatrixNotInitialized());
@@ -395,7 +398,8 @@ SparseMatrix<number>::Tvmult (Vector<somenumber>& dst, const Vector<somenumber>&
 template <typename number>
 template <typename somenumber>
 void
-SparseMatrix<number>::vmult_add (Vector<somenumber>& dst, const Vector<somenumber>& src) const
+SparseMatrix<number>::vmult_add (Vector<somenumber>& dst,
+                                 const Vector<somenumber>& src) const
 {
   Assert (cols != 0, ExcMatrixNotInitialized());
   Assert (val != 0, ExcMatrixNotInitialized());
@@ -420,7 +424,8 @@ SparseMatrix<number>::vmult_add (Vector<somenumber>& dst, const Vector<somenumbe
 template <typename number>
 template <typename somenumber>
 void
-SparseMatrix<number>::Tvmult_add (Vector<somenumber>& dst, const Vector<somenumber>& src) const
+SparseMatrix<number>::Tvmult_add (Vector<somenumber>& dst,
+                                  const Vector<somenumber>& src) const
 {
   Assert (val != 0, ExcMatrixNotInitialized());
   Assert (cols != 0, ExcMatrixNotInitialized());
@@ -522,10 +527,11 @@ SparseMatrix<number>::matrix_norm_square (const Vector<somenumber>& v) const
 template <typename number>
 template <typename somenumber>
 void
-SparseMatrix<number>::threaded_matrix_norm_square (const Vector<somenumber> &v,
-                                                  const unsigned int        begin_row,
-                                                  const unsigned int        end_row,
-                                                  somenumber               *partial_sum) const
+SparseMatrix<number>::
+threaded_matrix_norm_square (const Vector<somenumber> &v,
+                             const unsigned int        begin_row,
+                             const unsigned int        end_row,
+                             somenumber               *partial_sum) const
 {
                                   // this function should not be called
                                   // when not in parallel mode.
@@ -634,11 +640,12 @@ SparseMatrix<number>::matrix_scalar_product (const Vector<somenumber>& u,
 template <typename number>
 template <typename somenumber>
 void
-SparseMatrix<number>::threaded_matrix_scalar_product (const Vector<somenumber> &u,
-                                                     const Vector<somenumber> &v,
-                                                     const unsigned int        begin_row,
-                                                     const unsigned int        end_row,
-                                                     somenumber               *partial_sum) const
+SparseMatrix<number>::
+threaded_matrix_scalar_product (const Vector<somenumber> &u,
+                                const Vector<somenumber> &v,
+                                const unsigned int        begin_row,
+                                const unsigned int        end_row,
+                                somenumber               *partial_sum) const
 {
                                   // this function should not be called
                                   // when not in parallel mode.
@@ -921,8 +928,9 @@ SparseMatrix<number>::precondition_SSOR (Vector<somenumber>       &dst,
 template <typename number>
 template <typename somenumber>
 void
-SparseMatrix<number>::precondition_SOR (Vector<somenumber>& dst, const Vector<somenumber>& src,
-                           const number om) const
+SparseMatrix<number>::precondition_SOR (Vector<somenumber>& dst,
+                                        const Vector<somenumber>& src,
+                                        const number om) const
 {
   Assert (cols != 0, ExcMatrixNotInitialized());
   Assert (val != 0, ExcMatrixNotInitialized());
@@ -936,8 +944,9 @@ SparseMatrix<number>::precondition_SOR (Vector<somenumber>& dst, const Vector<so
 template <typename number>
 template <typename somenumber>
 void
-SparseMatrix<number>::precondition_TSOR (Vector<somenumber>& dst, const Vector<somenumber>& src,
-                           const number om) const
+SparseMatrix<number>::precondition_TSOR (Vector<somenumber>& dst,
+                                         const Vector<somenumber>& src,
+                                         const number om) const
 {
   Assert (cols != 0, ExcMatrixNotInitialized());
   Assert (val != 0, ExcMatrixNotInitialized());
@@ -951,7 +960,8 @@ SparseMatrix<number>::precondition_TSOR (Vector<somenumber>& dst, const Vector<s
 template <typename number>
 template <typename somenumber>
 void
-SparseMatrix<number>::SOR (Vector<somenumber>& dst, const number om) const
+SparseMatrix<number>::SOR (Vector<somenumber>& dst,
+                           const number om) const
 {
   Assert (cols != 0, ExcMatrixNotInitialized());
   Assert (val != 0, ExcMatrixNotInitialized());
@@ -976,7 +986,8 @@ SparseMatrix<number>::SOR (Vector<somenumber>& dst, const number om) const
 template <typename number>
 template <typename somenumber>
 void
-SparseMatrix<number>::TSOR (Vector<somenumber>& dst, const number om) const
+SparseMatrix<number>::TSOR (Vector<somenumber>& dst,
+                            const number om) const
 {
   Assert (cols != 0, ExcMatrixNotInitialized());
   Assert (val != 0, ExcMatrixNotInitialized());
@@ -999,11 +1010,10 @@ SparseMatrix<number>::TSOR (Vector<somenumber>& dst, const number om) const
 template <typename number>
 template <typename somenumber>
 void
-SparseMatrix<number>::PSOR (
-  Vector<somenumber>& dst,
-  const std::vector<unsigned int>& permutation,
-  const std::vector<unsigned int>& inverse_permutation,
-  const number om) const
+SparseMatrix<number>::PSOR (Vector<somenumber>& dst,
+                            const std::vector<unsigned int>& permutation,
+                            const std::vector<unsigned int>& inverse_permutation,
+                            const number om) const
 {
   Assert (cols != 0, ExcMatrixNotInitialized());
   Assert (val != 0, ExcMatrixNotInitialized());
@@ -1039,11 +1049,10 @@ SparseMatrix<number>::PSOR (
 template <typename number>
 template <typename somenumber>
 void
-SparseMatrix<number>::TPSOR (
-  Vector<somenumber>& dst,
-  const std::vector<unsigned int>& permutation,
-  const std::vector<unsigned int>& inverse_permutation,
-  const number om) const
+SparseMatrix<number>::TPSOR (Vector<somenumber>& dst,
+                             const std::vector<unsigned int>& permutation,
+                             const std::vector<unsigned int>& inverse_permutation,
+                             const number om) const
 {
   Assert (cols != 0, ExcMatrixNotInitialized());
   Assert (val != 0, ExcMatrixNotInitialized());
@@ -1071,12 +1080,13 @@ SparseMatrix<number>::TPSOR (
 }
 
 
+
 template <typename number>
 template <typename somenumber>
 void
 SparseMatrix<number>::SOR_step (Vector<somenumber> &v,
-              const Vector<somenumber> &b,
-              const number        om) const
+                                const Vector<somenumber> &b,
+                                const number        om) const
 {
   Assert (cols != 0, ExcMatrixNotInitialized());
   Assert (val != 0, ExcMatrixNotInitialized());
@@ -1095,12 +1105,14 @@ SparseMatrix<number>::SOR_step (Vector<somenumber> &v,
     }
 }
 
+
+
 template <typename number>
 template <typename somenumber>
 void
 SparseMatrix<number>::TSOR_step (Vector<somenumber> &v,
-              const Vector<somenumber> &b,
-              const number        om) const
+                                 const Vector<somenumber> &b,
+                                 const number        om) const
 {
   Assert (cols != 0, ExcMatrixNotInitialized());
   Assert (val != 0, ExcMatrixNotInitialized());
@@ -1119,22 +1131,26 @@ SparseMatrix<number>::TSOR_step (Vector<somenumber> &v,
     }
 }
 
+
+
 template <typename number>
 template <typename somenumber>
 void
 SparseMatrix<number>::SSOR_step (Vector<somenumber> &v,
-              const Vector<somenumber> &b,
-              const number        om) const
+                                 const Vector<somenumber> &b,
+                                 const number        om) const
 {
   SOR_step(v,b,om);
   TSOR_step(v,b,om);
 }
 
 
+
 template <typename number>
 template <typename somenumber>
 void
-SparseMatrix<number>::SSOR (Vector<somenumber>& dst, const number om) const
+SparseMatrix<number>::SSOR (Vector<somenumber>& dst,
+                            const number om) const
 {
   Assert (cols != 0, ExcMatrixNotInitialized());
   Assert (val != 0, ExcMatrixNotInitialized());
index 057e390dd7b738b5b72aaf5f946f8ee7fe52531e..ba8b2e142059d3d2699cdb1d6b96c98f10ed3f3d 100644 (file)
@@ -103,8 +103,9 @@ class Vector
 //     Vector (const Vector<OtherNumber> &v);
     
                                     /**
-                                     * Constructor. Set dimension to @p{n} and
-                                     * initialize all elements with zero.
+                                     * Constructor. Set dimension to
+                                     * @p{n} and initialize all
+                                     * elements with zero.
                                      */
     explicit Vector (const unsigned int n);
 
@@ -456,12 +457,13 @@ class Vector
                                      * the cellwise ratio of true to estimated
                                      * error.
                                      *
-                                     * This vector is appropriately scaled to
-                                     * hold the result.
+                                     * This vector is appropriately
+                                     * scaled to hold the result.
                                      *
-                                     * If any of the @p{b[i]} is zero, the result
-                                     * is undefined. No attempt is made to
-                                     * catch such situations.
+                                     * If any of the @p{b[i]} is
+                                     * zero, the result is
+                                     * undefined. No attempt is made
+                                     * to catch such situations.
                                      */
     void ratio (const Vector<Number> &a,
                const Vector<Number> &b);
@@ -585,7 +587,8 @@ class Vector
 
 template <typename Number>
 inline
-Vector<Number>::Vector () :
+Vector<Number>::Vector ()
+                :
                dim(0),
                maxdim(0),
                val(0)
@@ -612,7 +615,8 @@ Vector<Number>::Vector (const InputIterator first, const InputIterator last)
 
 template <typename Number>
 inline
-Vector<Number>::Vector (const unsigned int n) :
+Vector<Number>::Vector (const unsigned int n)
+                :
                dim(0),
                maxdim(0),
                val(0)
@@ -713,7 +717,7 @@ Vector<Number>::end ()
 template <typename Number>
 inline
 typename Vector<Number>::const_iterator
-Vector<Number>::end () const 
+Vector<Number>::end () const
 {
   return &val[dim];
 }
@@ -742,7 +746,7 @@ Number& Vector<Number>::operator() (const unsigned int i)
 
 template <typename Number>
 inline
-Vector<Number> & Vector<Number>::operator *= (const Number factor) 
+Vector<Number> & Vector<Number>::operator *= (const Number factor)
 {
   scale (factor);
   return *this;
@@ -752,7 +756,7 @@ Vector<Number> & Vector<Number>::operator *= (const Number factor)
 
 template <typename Number>
 inline
-Vector<Number> & Vector<Number>::operator /= (const Number factor) 
+Vector<Number> & Vector<Number>::operator /= (const Number factor)
 {
   scale (1./factor);
   return *this;
index e2af40dc79e73b1a7fc1da4c79c3e75e37634096..3ee37bacb582b4f097e1323d815430d4b01784b3 100644 (file)
@@ -39,7 +39,8 @@ static inline Number sqr (const Number x)
 
 
 template <typename Number>
-Vector<Number>::Vector (const Vector<Number>& v) :
+Vector<Number>::Vector (const Vector<Number>& v)
+                :
                dim(v.size()),
                maxdim(v.size()),
                val(0)
@@ -271,7 +272,8 @@ Number Vector<Number>::lp_norm (const Number p) const
 
 
 template <typename Number>
-Number Vector<Number>::linfty_norm () const {
+Number Vector<Number>::linfty_norm () const
+{
   Assert (dim!=0, ExcEmptyVector());
 
   Number max0=0.,
@@ -391,7 +393,8 @@ 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)
+void Vector<Number>::sadd (const Number x, const Number a,
+                           const Vector<Number>& v)
 {
   Assert (dim!=0, ExcEmptyVector());
   Assert (dim == v.dim, ExcDimensionMismatch(dim, v.dim));
@@ -405,7 +408,8 @@ 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)
+                          const Vector<Number>& v, const Number b,
+                           const Vector<Number>& w)
 {
   Assert (dim!=0, ExcEmptyVector());
   Assert (dim == v.dim, ExcDimensionMismatch(dim, v.dim));
@@ -422,7 +426,8 @@ 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,
-                          const Vector<Number>& w, const Number c, const Vector<Number>& y)
+                          const Vector<Number>& w, const Number c,
+                           const Vector<Number>& y)
 {
   Assert (dim!=0, ExcEmptyVector());
   Assert (dim == v.dim, ExcDimensionMismatch(dim, v.dim));
index 0497fc8f1915123b39cd6dda7f49f696c0eda491..8536379115c808f2bf4af0e830574e32faf33978 100644 (file)
 #include <numeric>
 #include <functional>
 
-CompressedSparsityPattern::CompressedSparsityPattern () :
+CompressedSparsityPattern::CompressedSparsityPattern ()
+                :
                rows(0),
                cols(0)
 {}
 
 
 
-CompressedSparsityPattern::CompressedSparsityPattern (const CompressedSparsityPattern &s) :
+CompressedSparsityPattern::
+CompressedSparsityPattern (const CompressedSparsityPattern &s)
+                :
                Subscriptor(),
                rows(0),
                cols(0)
@@ -41,8 +44,9 @@ CompressedSparsityPattern::CompressedSparsityPattern (const CompressedSparsityPa
 
 CompressedSparsityPattern::CompressedSparsityPattern (const unsigned int m,
                                                      const unsigned int n) 
-               : rows(0),
-                 cols(0)
+               :
+                rows(0),
+                cols(0)
 {
   reinit (m,n);
 }
@@ -50,8 +54,9 @@ CompressedSparsityPattern::CompressedSparsityPattern (const unsigned int m,
 
 
 CompressedSparsityPattern::CompressedSparsityPattern (const unsigned int n)
-               : rows(0),
-                 cols(0)
+               :
+                rows(0),
+                cols(0)
 {
   reinit (n,n);
 }
@@ -100,7 +105,7 @@ CompressedSparsityPattern::empty () const
 
 
 unsigned int
-CompressedSparsityPattern::max_entries_per_row () const 
+CompressedSparsityPattern::max_entries_per_row () const
 {
   unsigned int m = 0;
   for (unsigned int i=1; i<rows; ++i)
@@ -137,7 +142,7 @@ CompressedSparsityPattern::exists (const unsigned int i,
 
 
 void
-CompressedSparsityPattern::symmetrize () 
+CompressedSparsityPattern::symmetrize ()
 {
   Assert (rows==cols, ExcNotSquare());
 
index f2d340de6dc5d12ee85c186ad43e15c6fe402a0f..c30794a372789b9476b6574b797dc5168111b5ca 100644 (file)
 
 // TYPEMAT and TYPEVEC are defined in sparsematrix?.cc
 
-template void SparseMatrix<TYPEMAT>::vmult<TYPEVEC> (
-  Vector<TYPEVEC> &,
-  const Vector<TYPEVEC> &) const;
-template void SparseMatrix<TYPEMAT>::Tvmult<TYPEVEC> (
-  Vector<TYPEVEC> &,
-  const Vector<TYPEVEC> &) const;
-template void SparseMatrix<TYPEMAT>::vmult_add<TYPEVEC> (
-  Vector<TYPEVEC> &,
-  const Vector<TYPEVEC> &) const;
-template void SparseMatrix<TYPEMAT>::Tvmult_add<TYPEVEC> (
-  Vector<TYPEVEC> &,
-  const Vector<TYPEVEC> &) const;
+template void SparseMatrix<TYPEMAT>::
+vmult<TYPEVEC> (Vector<TYPEVEC> &,
+                const Vector<TYPEVEC> &) const;
+template void SparseMatrix<TYPEMAT>::
+Tvmult<TYPEVEC> (Vector<TYPEVEC> &,
+                 const Vector<TYPEVEC> &) const;
+template void SparseMatrix<TYPEMAT>::
+vmult_add<TYPEVEC> (Vector<TYPEVEC> &,
+                    const Vector<TYPEVEC> &) const;
+template void SparseMatrix<TYPEMAT>::
+Tvmult_add<TYPEVEC> (Vector<TYPEVEC> &,
+                     const Vector<TYPEVEC> &) const;
 
 template TYPEVEC
-SparseMatrix<TYPEMAT>::matrix_norm_square<TYPEVEC> (
-  const Vector<TYPEVEC> &) const;
+SparseMatrix<TYPEMAT>::
+matrix_norm_square<TYPEVEC> (const Vector<TYPEVEC> &) const;
 
 template TYPEVEC
-SparseMatrix<TYPEMAT>::matrix_scalar_product<TYPEVEC> (
-  const Vector<TYPEVEC> &,
-  const Vector<TYPEVEC> &) const;
+SparseMatrix<TYPEMAT>::
+matrix_scalar_product<TYPEVEC> (const Vector<TYPEVEC> &,
+                                const Vector<TYPEVEC> &) const;
 
-template TYPEVEC SparseMatrix<TYPEMAT>::residual<TYPEVEC> (
-  Vector<TYPEVEC> &,
-  const Vector<TYPEVEC> &,
-  const Vector<TYPEVEC> &) const;
+template TYPEVEC SparseMatrix<TYPEMAT>::
+residual<TYPEVEC> (Vector<TYPEVEC> &,
+                   const Vector<TYPEVEC> &,
+                   const Vector<TYPEVEC> &) const;
 
-template void SparseMatrix<TYPEMAT>::precondition_SSOR<TYPEVEC> (
-  Vector<TYPEVEC> &,
-  const Vector<TYPEVEC> &,
-  const TYPEMAT) const;
+template void SparseMatrix<TYPEMAT>::
+precondition_SSOR<TYPEVEC> (Vector<TYPEVEC> &,
+                            const Vector<TYPEVEC> &,
+                            const TYPEMAT) const;
 
-template void SparseMatrix<TYPEMAT>::precondition_SOR<TYPEVEC> (
-  Vector<TYPEVEC> &,
-  const Vector<TYPEVEC> &,
-  const TYPEMAT) const;
+template void SparseMatrix<TYPEMAT>::
+precondition_SOR<TYPEVEC> (Vector<TYPEVEC> &,
+                           const Vector<TYPEVEC> &,
+                           const TYPEMAT) const;
 
-template void SparseMatrix<TYPEMAT>::precondition_TSOR<TYPEVEC> (
-  Vector<TYPEVEC> &,
-  const Vector<TYPEVEC> &,
-  const TYPEMAT) const;
+template void SparseMatrix<TYPEMAT>::
+precondition_TSOR<TYPEVEC> (Vector<TYPEVEC> &,
+                            const Vector<TYPEVEC> &,
+                            const TYPEMAT) const;
 
-template void SparseMatrix<TYPEMAT>::precondition_Jacobi<TYPEVEC> (
-  Vector<TYPEVEC> &,
-  const Vector<TYPEVEC> &,
-  const TYPEMAT) const;
+template void SparseMatrix<TYPEMAT>::
+precondition_Jacobi<TYPEVEC> (Vector<TYPEVEC> &,
+                              const Vector<TYPEVEC> &,
+                              const TYPEMAT) const;
 
-template void SparseMatrix<TYPEMAT>::SOR<TYPEVEC> (
-  Vector<TYPEVEC> &,
-  const TYPEMAT) const;
-template void SparseMatrix<TYPEMAT>::TSOR<TYPEVEC> (
-  Vector<TYPEVEC> &,
-  const TYPEMAT) const;
-template void SparseMatrix<TYPEMAT>::SSOR<TYPEVEC> (
-  Vector<TYPEVEC> &,
-  const TYPEMAT) const;
-template void SparseMatrix<TYPEMAT>::PSOR<TYPEVEC> (
-  Vector<TYPEVEC> &,
-  const std::vector<unsigned int>&,
-  const std::vector<unsigned int>&,
-  const TYPEMAT) const;
-template void SparseMatrix<TYPEMAT>::TPSOR<TYPEVEC> (
-  Vector<TYPEVEC> &,
-  const std::vector<unsigned int>&,
-  const std::vector<unsigned int>&,
-  const TYPEMAT) const;
-template void SparseMatrix<TYPEMAT>::SOR_step<TYPEVEC> (
-  Vector<TYPEVEC> &,
-  const Vector<TYPEVEC> &,
-  const TYPEMAT) const;
-template void SparseMatrix<TYPEMAT>::TSOR_step<TYPEVEC> (
-  Vector<TYPEVEC> &,
-  const Vector<TYPEVEC> &,
-  const TYPEMAT) const;
-template void SparseMatrix<TYPEMAT>::SSOR_step<TYPEVEC> (
-  Vector<TYPEVEC> &,
-  const Vector<TYPEVEC> &, 
-  const TYPEMAT) const;
+template void SparseMatrix<TYPEMAT>::
+SOR<TYPEVEC> (Vector<TYPEVEC> &,
+              const TYPEMAT) const;
+template void SparseMatrix<TYPEMAT>::
+TSOR<TYPEVEC> (Vector<TYPEVEC> &,
+               const TYPEMAT) const;
+template void SparseMatrix<TYPEMAT>::
+SSOR<TYPEVEC> (Vector<TYPEVEC> &,
+               const TYPEMAT) const;
+template void SparseMatrix<TYPEMAT>::
+PSOR<TYPEVEC> (Vector<TYPEVEC> &,
+               const std::vector<unsigned int>&,
+               const std::vector<unsigned int>&,
+               const TYPEMAT) const;
+template void SparseMatrix<TYPEMAT>::
+TPSOR<TYPEVEC> (Vector<TYPEVEC> &,
+                const std::vector<unsigned int>&,
+                const std::vector<unsigned int>&,
+                const TYPEMAT) const;
+template void SparseMatrix<TYPEMAT>::
+SOR_step<TYPEVEC> (Vector<TYPEVEC> &,
+                   const Vector<TYPEVEC> &,
+                   const TYPEMAT) const;
+template void SparseMatrix<TYPEMAT>::
+TSOR_step<TYPEVEC> (Vector<TYPEVEC> &,
+                    const Vector<TYPEVEC> &,
+                    const TYPEMAT) const;
+template void SparseMatrix<TYPEMAT>::
+SSOR_step<TYPEVEC> (Vector<TYPEVEC> &,
+                    const Vector<TYPEVEC> &, 
+                    const TYPEMAT) const;

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.