]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
fixed a small bug
authorkanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 30 Mar 2012 19:13:54 +0000 (19:13 +0000)
committerkanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 30 Mar 2012 19:13:54 +0000 (19:13 +0000)
git-svn-id: https://svn.dealii.org/trunk@25341 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/include/deal.II/meshworker/simple.h

index edab5351ff7ef866927aaaececcb872448f289dd..4268345b679579a3a0ad0fb734ca0b27d2043467 100644 (file)
@@ -1076,7 +1076,7 @@ namespace MeshWorker
        {
          for (unsigned int j=0; j<i1.size(); ++j)
            for (unsigned int k=0; k<i2.size(); ++k)
-             if (std::fabs(M(j,k)) >= 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)
        {

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.