From: Martin Kronbichler Date: Wed, 25 Aug 2010 19:56:55 +0000 (+0000) Subject: Fix bug introduced in r21708 X-Git-Tag: v8.0.0~5683 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=aebb1f54ae76fb915a5b3806f2c65cca8f5713c8;p=dealii.git Fix bug introduced in r21708 git-svn-id: https://svn.dealii.org/trunk@21715 0785d39b-7218-0410-832d-ea1e28bc413d --- 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