From: wolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Mon, 7 Feb 2000 12:29:11 +0000 (+0000)
Subject: Small clean-up.
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=7dde5f365fc44c4411ff9120d16ed91d15dbeb0b;p=dealii-svn.git

Small clean-up.


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

diff --git a/deal.II/lac/include/lac/sparse_vanka.h b/deal.II/lac/include/lac/sparse_vanka.h
index 71f0623d19..76b0d1ad81 100644
--- a/deal.II/lac/include/lac/sparse_vanka.h
+++ b/deal.II/lac/include/lac/sparse_vanka.h
@@ -499,14 +499,29 @@ class SparseBlockVanka : public SparseVanka<number>
 
 				     /**
 				      * In this field, we precompute
-				      * the first and the one after
-				      * the last index of each
-				      * block. This computation is
-				      * done in the constructor, to
-				      * avoid recomputing each time
-				      * the preconditioner is called.
+				      * for each block which degrees
+				      * of freedom belong to it. Thus,
+				      * if #dof_masks[i][j]==true#,
+				      * then DoF #j# belongs to block
+				      * #i#. Of course, no other
+				      * #dof_masks[l][j]# may be
+				      * #true# for #l!=i#. This
+				      * computation is done in the
+				      * constructor, to avoid
+				      * recomputing each time the
+				      * preconditioner is called.
 				      */
     vector<vector<bool> > dof_masks;
+
+				     /**
+				      * Compute the contents of the
+				      * field #dof_masks#. This
+				      * function is called from the
+				      * constructor.
+				      */
+    void compute_dof_masks (const SparseMatrix<number> &M,
+			    const vector<bool>         &selected,
+			    const BlockingStrategy      blocking_strategy);
 };
 
 
diff --git a/deal.II/lac/include/lac/sparse_vanka.templates.h b/deal.II/lac/include/lac/sparse_vanka.templates.h
index a2abf3c700..be154d391b 100644
--- a/deal.II/lac/include/lac/sparse_vanka.templates.h
+++ b/deal.II/lac/include/lac/sparse_vanka.templates.h
@@ -397,6 +397,17 @@ SparseBlockVanka<number>::SparseBlockVanka (const SparseMatrix<number> &M,
                 n_blocks (n_blocks),
                 dof_masks (n_blocks,
 		           vector<bool>(M.m(), false))
+{
+  compute_dof_masks (M, selected, blocking_strategy);
+};
+
+
+
+template <typename number>
+void
+SparseBlockVanka<number>::compute_dof_masks (const SparseMatrix<number> &M,
+					     const vector<bool>         &selected,
+					     const BlockingStrategy      blocking_strategy) 
 {
   Assert (n_blocks > 0, ExcInternalError());