From 2e28a5e39cf40ae2f69b40a124f11e15e4509a6d Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Sun, 3 Jun 2012 19:39:43 +0000
Subject: [PATCH] Do not check for zero entries when restricting occupation
patterns.
Because some people restrict patterns for uninitialized matrices.
In that case you get
a) valgrind warnings
b) crashes, if an entry happens to be zero, it is not restricted.
[[Imported from SVN: r6332]]
---
dune/solvers/transferoperators/genericmultigridtransfer.hh | 4 ----
1 file changed, 4 deletions(-)
diff --git a/dune/solvers/transferoperators/genericmultigridtransfer.hh b/dune/solvers/transferoperators/genericmultigridtransfer.hh
index 2388f0c..2ea769e 100644
--- a/dune/solvers/transferoperators/genericmultigridtransfer.hh
+++ b/dune/solvers/transferoperators/genericmultigridtransfer.hh
@@ -770,8 +770,6 @@ public:
for (; m!=mEnd; ++m)
{
- if (m->infinity_norm()==0)
- continue;
int w = m.index();
// Loop over all coarse grid vectors iv that have v in their support
@@ -831,8 +829,6 @@ public:
for (; m!=mEnd; ++m)
{
- if (m->infinity_norm()==0)
- continue;
int w = m.index();
int w_block = w/virtualBlockSize;
int w_inblock = w%virtualBlockSize;
--
GitLab