From: Denis Davydov <davydden@gmail.com>
Date: Tue, 16 Jun 2015 21:16:37 +0000 (+0200)
Subject: modify check of non-zero diagonal element in PETScWrappers::apply_boundary_values... 
X-Git-Tag: v8.3.0-rc1~104^2~8
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=a039fa29fb240080180be69a45489f0595b03fa3;p=dealii.git

modify check of non-zero diagonal element in PETScWrappers::apply_boundary_values to complex case
---

diff --git a/source/numerics/matrix_tools.cc b/source/numerics/matrix_tools.cc
index 3f38c170fc..8a385ed620 100644
--- a/source/numerics/matrix_tools.cc
+++ b/source/numerics/matrix_tools.cc
@@ -2351,9 +2351,9 @@ namespace MatrixTools
             // find such an entry, take one
             PetscScalar average_nonzero_diagonal_entry = 1;
             for (types::global_dof_index i=local_range.first; i<local_range.second; ++i)
-              if (matrix.diag_element(i) != 0)
+              if (matrix.diag_element(i) != PetscScalar ())
                 {
-                  average_nonzero_diagonal_entry = std::fabs(matrix.diag_element(i));
+                  average_nonzero_diagonal_entry = std::abs(matrix.diag_element(i));
                   break;
                 }