]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Empty blocks are ignored in vmult and Tvmult
authorguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 16 Dec 2002 20:55:44 +0000 (20:55 +0000)
committerguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 16 Dec 2002 20:55:44 +0000 (20:55 +0000)
git-svn-id: https://svn.dealii.org/trunk@6832 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/lac/include/lac/block_matrix_array.h

index 334d77af2ebb81abf7edbdaad0c5b4d0cb0cb398..ebf5b893ce137a59e273e8af9d48f8432b192dfe 100644 (file)
@@ -388,6 +388,8 @@ BlockMatrixArray<MATRIX>::vmult_add (BlockVector<number>& dst,
   
   for (; m != end ; ++m)
     {
+      if (m->matrix->empty())
+       continue;
       if (m->prefix==1.)
        {
          if (m->transpose)
@@ -445,6 +447,8 @@ BlockMatrixArray<MATRIX>::Tvmult_add (BlockVector<number>& dst,
   
   for (; m != end ; ++m)
     {
+      if (m->matrix->empty())
+       continue;
       if (m->prefix==1.)
        {
          if (m->transpose)

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.