From: Wolfgang Bangerth Date: Wed, 13 Feb 2013 14:05:26 +0000 (+0000) Subject: Fix some confusion from yesterday. Not a bug here because matrices have to be X-Git-Tag: v8.0.0~1319 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=5836591b37a8bb51d0173dfd7b933da812c52b0b;p=dealii.git Fix some confusion from yesterday. Not a bug here because matrices have to be square, but still harder to read than necessary. git-svn-id: https://svn.dealii.org/trunk@28368 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/source/lac/sparse_direct.cc b/deal.II/source/lac/sparse_direct.cc index eb42a4d19e..8585f4cba4 100644 --- a/deal.II/source/lac/sparse_direct.cc +++ b/deal.II/source/lac/sparse_direct.cc @@ -2,7 +2,7 @@ // $Id$ // Version: $Name$ // -// Copyright (C) 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2009, 2010, 2011, 2012 by the deal.II authors +// Copyright (C) 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2009, 2010, 2011, 2012, 2013 by the deal.II authors // // This file is subject to QPL and may not be distributed // without copyright and license information. Please refer @@ -1728,7 +1728,7 @@ factorize (const Matrix &matrix) // loop over the elements of the matrix row by row, as suggested // in the documentation of the sparse matrix iterator class - for (unsigned int row = 0; row < matrix.n(); ++row) + for (unsigned int row = 0; row < matrix.m(); ++row) { for (typename Matrix::const_iterator p=matrix.begin(row); p!=matrix.end(row); ++p) @@ -1962,7 +1962,7 @@ void SparseDirectMUMPS::initialize_matrix (const Matrix &matrix) // loop over the elements of the matrix row by row, as suggested // in the documentation of the sparse matrix iterator class - for (unsigned int row = 0; row < matrix.n(); ++row) + for (unsigned int row = 0; row < matrix.m(); ++row) { for (typename Matrix::const_iterator ptr = matrix.begin (row); ptr != matrix.end (row); ++ptr)