From d14089c49b8ad75bc43c87893e41afac3d76ec55 Mon Sep 17 00:00:00 2001
From: Felix Gruber <gruber@igpm.rwth-aachen.de>
Date: Wed, 17 Aug 2016 16:08:47 +0200
Subject: [PATCH] [bugfix] fix Matrix::transpose()
Due to a confusion between the number of lines N and the number of
columns M, Matrix::transpose() gave wrong results for matrices with
N != M.
---
dune/istl/matrix.hh | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/dune/istl/matrix.hh b/dune/istl/matrix.hh
index 946d28ea4..e8f8e1563 100644
--- a/dune/istl/matrix.hh
+++ b/dune/istl/matrix.hh
@@ -738,9 +738,9 @@ namespace MatrixImp
/** \brief Return the transpose of the matrix */
Matrix transpose() const {
- Matrix out(N(), M());
- for (size_type i=0; i<M(); i++)
- for (size_type j=0; j<N(); j++)
+ Matrix out(M(), N());
+ for (size_type i=0; i<N(); i++)
+ for (size_type j=0; j<M(); j++)
out[j][i] = (*this)[i][j];
return out;
--
GitLab