From 53d38e232daf83e48e04aa47e98e787458a97b55 Mon Sep 17 00:00:00 2001 From: wolf Date: Thu, 18 May 2000 12:07:42 +0000 Subject: [PATCH] Add more assertions. git-svn-id: https://svn.dealii.org/trunk@2887 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/include/numerics/matrices.h | 4 ++++ deal.II/deal.II/source/numerics/matrices.cc | 7 +++++++ 2 files changed, 11 insertions(+) diff --git a/deal.II/deal.II/include/numerics/matrices.h b/deal.II/deal.II/include/numerics/matrices.h index 14aa6cdb3f..679948e1de 100644 --- a/deal.II/deal.II/include/numerics/matrices.h +++ b/deal.II/deal.II/include/numerics/matrices.h @@ -483,6 +483,10 @@ class MatrixTools : public MatrixCreator * Exception */ DeclException0 (ExcMatrixNotBlockSquare); + /** + * Exception + */ + DeclException0 (ExcBlocksDontMatch); }; diff --git a/deal.II/deal.II/source/numerics/matrices.cc b/deal.II/deal.II/source/numerics/matrices.cc index 9e64d9de02..fa62dc9371 100644 --- a/deal.II/deal.II/source/numerics/matrices.cc +++ b/deal.II/deal.II/source/numerics/matrices.cc @@ -698,6 +698,13 @@ MatrixTools::apply_boundary_values (const map &boundar Assert (matrix.get_sparsity_pattern().get_row_indices() == matrix.get_sparsity_pattern().get_column_indices(), ExcMatrixNotBlockSquare()); + Assert (matrix.get_sparsity_pattern().get_column_indices() == + solution.get_block_indices (), + ExcBlocksDontMatch ()); + Assert (matrix.get_sparsity_pattern().get_row_indices() == + right_hand_side.get_block_indices (), + ExcBlocksDontMatch ()); + // if no boundary values are to be applied // simply return -- 2.39.5