tenporary fix for matrix.

git-svn-id: https://svn.dealii.org/branches/branch_petscscalar_complex@31035 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/source/numerics/matrix_tools.cc b/deal.II/source/numerics/matrix_tools.cc
index ad3c1936a24f8b6d3655a1836389da0ca6454f87..136fd5baeef2dbdc63cb67a91e1a09eb88ec9af8 100644 (file)
--- a/deal.II/source/numerics/matrix_tools.cc
+++ b/deal.II/source/numerics/matrix_tools.cc
@@ -2422,16 +2422,21 @@ namespace MatrixTools
         // find such an entry, take one
 #ifndef PETSC_USE_COMPLEX
         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)
+            {
+              average_nonzero_diagonal_entry = std::fabs(matrix.diag_element(i));
+              break;
+            }
 #else
         PetscScalar average_nonzero_diagonal_entry = (PetscScalar) std::complex<double> (1,1);
-#endif
-
         for (types::global_dof_index i=local_range.first; i<local_range.second; ++i)
-          if (matrix.diag_element(i) != std::complex<double> (0,0))
+          if (matrix.diag_element(i) != std::complex<double> (0,0)) // PETSC_NULL?
             {
               average_nonzero_diagonal_entry = std::fabs(matrix.diag_element(i));
               break;
             }
+#endif
 
         // figure out which rows of the matrix we
         // have to eliminate on this processor