From a039fa29fb240080180be69a45489f0595b03fa3 Mon Sep 17 00:00:00 2001
From: Denis Davydov <davydden@gmail.com>
Date: Tue, 16 Jun 2015 23:16:37 +0200
Subject: [PATCH] modify check of non-zero diagonal element in
 PETScWrappers::apply_boundary_values to complex case

---
 source/numerics/matrix_tools.cc | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

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;
                 }
 
-- 
2.39.5