From aebb1f54ae76fb915a5b3806f2c65cca8f5713c8 Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Wed, 25 Aug 2010 19:56:55 +0000 Subject: [PATCH] Fix bug introduced in r21708 git-svn-id: https://svn.dealii.org/trunk@21715 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/lac/include/lac/constraint_matrix.templates.h | 9 +++++++++ 1 file changed, 9 insertions(+) diff --git a/deal.II/lac/include/lac/constraint_matrix.templates.h b/deal.II/lac/include/lac/constraint_matrix.templates.h index 042080ceb1..4717cfa796 100644 --- a/deal.II/lac/include/lac/constraint_matrix.templates.h +++ b/deal.II/lac/include/lac/constraint_matrix.templates.h @@ -1507,6 +1507,9 @@ namespace internals unsigned int * &col_ptr, number * &val_ptr) { + if (column_end == column_start) + return; + AssertIndexRange (column_end-1, global_cols.size()); const unsigned int loc_row = global_rows.local_row(i); @@ -1606,6 +1609,9 @@ namespace internals const FullMatrix &local_matrix, SparseMatrix *sparse_matrix) { + if (column_end == column_start) + return; + AssertIndexRange (column_end-1, global_rows.size()); const SparsityPattern & sparsity = sparse_matrix->get_sparsity_pattern(); #ifndef DEBUG @@ -1758,6 +1764,9 @@ namespace internals const Table<2,bool> &dof_mask, std::vector::iterator &col_ptr) { + if (column_end == column_start) + return; + const unsigned int loc_row = global_rows.local_row(i); // fast function if there are no indirect -- 2.39.5