From 4383058ea98ecb3bf60d549d1cbdb9a818018309 Mon Sep 17 00:00:00 2001 From: kanschat Date: Mon, 12 Mar 2012 22:56:34 +0000 Subject: [PATCH] add initialization of local blocks git-svn-id: https://svn.dealii.org/trunk@25265 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/include/deal.II/meshworker/simple.h | 23 +++++++++++++++++++-- 1 file changed, 21 insertions(+), 2 deletions(-) diff --git a/deal.II/include/deal.II/meshworker/simple.h b/deal.II/include/deal.II/meshworker/simple.h index 6e40780534..a7b3c6ef85 100644 --- a/deal.II/include/deal.II/meshworker/simple.h +++ b/deal.II/include/deal.II/meshworker/simple.h @@ -165,7 +165,7 @@ namespace MeshWorker * initialize_info() will * generate one local matrix * for each block row and - * column, whch will be + * column, which will be * numbered * lexicographically, row by * row. @@ -650,7 +650,26 @@ namespace MeshWorker inline void MatrixSimple::initialize_info(DOFINFO& info, bool face) const { - info.initialize_matrices(1, face); + const unsigned int n = local_indices.size(); + + if (n == 0) + info.initialize_matrices(1, face); + else + { + info.initialize_matrices(n*n, face); + unsigned int k=0; + for (unsigned int i=0;i