From: Denis Davydov <davydden@gmail.com>
Date: Tue, 3 Jan 2017 17:23:55 +0000 (+0100)
Subject: make an Assert in LA::distributed::Vector::compress_finish() less rigid
X-Git-Tag: v8.5.0-rc1~294^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F3727%2Fhead;p=dealii.git

make an Assert in LA::distributed::Vector::compress_finish() less rigid

Otherwise p::d::SolutionTransfer::interpolate() could trigger the
Assert due to round-off errors in certain scenarious.
---

diff --git a/include/deal.II/lac/la_parallel_vector.templates.h b/include/deal.II/lac/la_parallel_vector.templates.h
index 24a6ee0c16..37d8433162 100644
--- a/include/deal.II/lac/la_parallel_vector.templates.h
+++ b/include/deal.II/lac/la_parallel_vector.templates.h
@@ -680,7 +680,7 @@ namespace LinearAlgebra
                    j++, read_position++)
                 Assert(*read_position == Number() ||
                        std::abs(local_element(j) - *read_position) <=
-                       std::abs(local_element(j)) * 1000. *
+                       std::abs(local_element(j)) * 10000. *
                        std::numeric_limits<real_type>::epsilon(),
                        ExcNonMatchingElements(*read_position, local_element(j),
                                               part.this_mpi_process()));