From: kanschat Date: Fri, 30 Mar 2012 19:13:54 +0000 (+0000) Subject: fixed a small bug X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=da2ce5f43f38c966621122fd09ac7ee8e0cc7844;p=dealii-svn.git fixed a small bug git-svn-id: https://svn.dealii.org/trunk@25341 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/include/deal.II/meshworker/simple.h b/deal.II/include/deal.II/meshworker/simple.h index edab5351ff..4268345b67 100644 --- a/deal.II/include/deal.II/meshworker/simple.h +++ b/deal.II/include/deal.II/meshworker/simple.h @@ -1076,7 +1076,7 @@ namespace MeshWorker { for (unsigned int j=0; j= threshold) + if (std::fabs(M(k,j)) >= threshold) if(mg_constrained_dofs->at_refinement_edge(level, i1[j]) && !mg_constrained_dofs->at_refinement_edge(level, i2[k])) { @@ -1128,9 +1128,9 @@ namespace MeshWorker if(mg_constrained_dofs != 0) { assemble_in((*interface_in)[level], info.matrix(k,false).matrix, - info.indices, info.indices, level); + info.indices_by_block[row], info.indices_by_block[column], level); assemble_out((*interface_out)[level],info.matrix(k,false).matrix, - info.indices_by_block[row], info.indices_by_block[column], level); + info.indices_by_block[column], info.indices_by_block[row], level); } } } @@ -1144,7 +1144,6 @@ namespace MeshWorker { const unsigned int level1 = info1.cell->level(); const unsigned int level2 = info2.cell->level(); - if (local_indices.size() == 0) {