From: David Wells <wellsd2@rpi.edu>
Date: Thu, 15 Mar 2018 20:03:12 +0000 (-0400)
Subject: Fix some (rectangular) FullMatrix methods.
X-Git-Tag: v9.0.0-rc1~322^2~5
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=724ac687914976005b5ec631d8c9fdcfdd29d063;p=dealii.git

Fix some (rectangular) FullMatrix methods.
---

diff --git a/include/deal.II/lac/full_matrix.templates.h b/include/deal.II/lac/full_matrix.templates.h
index b511693b1f..9859a8992e 100644
--- a/include/deal.II/lac/full_matrix.templates.h
+++ b/include/deal.II/lac/full_matrix.templates.h
@@ -255,10 +255,10 @@ number FullMatrix<number>::residual (Vector<number2> &dst,
   number res = 0.;
   const size_type size_m = m(),
                   size_n = n();
-  for (size_type i=0; i<size_n; ++i)
+  for (size_type i=0; i<size_m; ++i)
     {
       number s = number(right(i));
-      for (size_type j=0; j<size_m; ++j)
+      for (size_type j=0; j<size_n; ++j)
         s -= number(src(j)) * (*this)(i,j);
       dst(i) = s;
       res += s*s;
@@ -364,7 +364,7 @@ void FullMatrix<number>::add_row (const size_type i,
 {
   Assert (!this->empty(), ExcEmptyMatrix());
 
-  for (size_type k=0; k<m(); ++k)
+  for (size_type k=0; k<n(); ++k)
     (*this)(i,k) += s*(*this)(j,k);
 }
 
@@ -378,7 +378,7 @@ void FullMatrix<number>::add_row (const size_type i,
 {
   Assert (!this->empty(), ExcEmptyMatrix());
 
-  const size_type size_m = m();
+  const size_type size_m = n();
   for (size_type l=0; l<size_m; ++l)
     (*this)(i,l) += s*(*this)(j,l) + t*(*this)(k,l);
 }
@@ -390,7 +390,7 @@ void FullMatrix<number>::add_col (const size_type i, const number s,
 {
   Assert (!this->empty(), ExcEmptyMatrix());
 
-  for (size_type k=0; k<n(); ++k)
+  for (size_type k=0; k<m(); ++k)
     (*this)(k,i) += s*(*this)(k,j);
 }
 
@@ -402,7 +402,7 @@ void FullMatrix<number>::add_col (const size_type i, const number s,
 {
   Assert (!this->empty(), ExcEmptyMatrix());
 
-  for (size_t l=0; l<n(); ++l)
+  for (size_t l=0; l<m(); ++l)
     (*this)(l,i) += s*(*this)(l,j) + t*(*this)(l,k);
 }