From: Toby D. Young Date: Thu, 22 Jul 2010 13:34:15 +0000 (+0000) Subject: Bug fix: Relaxed tolerance on is_hermitian() and is_symmetric() operations on PETScMa... X-Git-Tag: v8.0.0~5777 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=81302404a21291e40d0e68de54543f05d10ae873;p=dealii.git Bug fix: Relaxed tolerance on is_hermitian() and is_symmetric() operations on PETScMatrixBase such that they can return 1 (true) for the deal.II default tolerance git-svn-id: https://svn.dealii.org/trunk@21556 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/lac/include/lac/petsc_matrix_base.h b/deal.II/lac/include/lac/petsc_matrix_base.h index b4f4d80426..857398860b 100644 --- a/deal.II/lac/include/lac/petsc_matrix_base.h +++ b/deal.II/lac/include/lac/petsc_matrix_base.h @@ -1091,19 +1091,25 @@ namespace PETScWrappers void transpose (); /** - * Test whether a matrix is symmetric. - * Default tolerance is zero. + * Test whether a matrix is + * symmetric. Default + * tolerance is + * $1000\times32$-bit machine + * precision. */ - PetscTruth is_symmetric (const double tol = 0.0); + PetscTruth is_symmetric (const double tolerance = 1.e-12); #if DEAL_II_PETSC_VERSION_GTE(2,3,0) /** - * Test whether a matrix is Hermitian, - * i.e. it is the complex conjugate - * of its transpose. Default tolerance - * is zero. + * Test whether a matrix is + * Hermitian, i.e. it is the + * complex conjugate of its + * transpose. Default + * tolerance is + * $1000\times32$-bit machine + * precision. */ - PetscTruth is_hermitian (const double tol = 0.0); + PetscTruth is_hermitian (const double tolerance = 1.e-12); #endif /*