From e60f25d9c3854666178a5afcee74da115ace4be8 Mon Sep 17 00:00:00 2001
From: Guido Kanschat <dr.guido.kanschat@gmail.com>
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<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)
-- 
2.39.5