From: guido Date: Mon, 16 Dec 2002 20:55:44 +0000 (+0000) Subject: Empty blocks are ignored in vmult and Tvmult X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=f6fd3ff7ad08e27e825f11ff9fec954878032cb7;p=dealii-svn.git Empty blocks are ignored in vmult and Tvmult git-svn-id: https://svn.dealii.org/trunk@6832 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/lac/include/lac/block_matrix_array.h b/deal.II/lac/include/lac/block_matrix_array.h index 334d77af2e..ebf5b893ce 100644 --- a/deal.II/lac/include/lac/block_matrix_array.h +++ b/deal.II/lac/include/lac/block_matrix_array.h @@ -388,6 +388,8 @@ BlockMatrixArray::vmult_add (BlockVector& dst, for (; m != end ; ++m) { + if (m->matrix->empty()) + continue; if (m->prefix==1.) { if (m->transpose) @@ -445,6 +447,8 @@ BlockMatrixArray::Tvmult_add (BlockVector& dst, for (; m != end ; ++m) { + if (m->matrix->empty()) + continue; if (m->prefix==1.) { if (m->transpose)