From f6fd3ff7ad08e27e825f11ff9fec954878032cb7 Mon Sep 17 00:00:00 2001 From: guido Date: Mon, 16 Dec 2002 20:55:44 +0000 Subject: [PATCH] Empty blocks are ignored in vmult and Tvmult git-svn-id: https://svn.dealii.org/trunk@6832 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/lac/include/lac/block_matrix_array.h | 4 ++++ 1 file changed, 4 insertions(+) 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) -- 2.39.5