From: Denis Davydov 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=b2a46fdc589a5233c16a1de809c90bb6b5728d6e;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::epsilon(), ExcNonMatchingElements(*read_position, local_element(j), part.this_mpi_process()));