From d99f321adba7537c2b03f1c2ca41d2fb77508478 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 8 Mar 2012 14:37:28 +0000
Subject: [PATCH] Don't divide by the determinant in the gradient; add a
default constructor
[[Imported from SVN: r8562]]
---
dune/gfe/orthogonalmatrix.hh | 12 ++++++------
1 file changed, 6 insertions(+), 6 deletions(-)
diff --git a/dune/gfe/orthogonalmatrix.hh b/dune/gfe/orthogonalmatrix.hh
index b565c4e0..7e02ac7d 100644
--- a/dune/gfe/orthogonalmatrix.hh
+++ b/dune/gfe/orthogonalmatrix.hh
@@ -12,17 +12,17 @@ class OrthogonalMatrix
{
public:
+ /** \brief Default constructor -- leaves the matrix uninitialized */
+ OrthogonalMatrix()
+ {}
+
/** \brief Constructor from a general matrix
*
- * The input matrix gets normalized to det = 1. Since computing
- * the determinant is expensive you may not want to use this method
- * if you know your input matrix to be orthogonal anyways.
+ * It is not checked whether the matrix is really orthogonal
*/
explicit OrthogonalMatrix(const Dune::FieldMatrix<T,N,N>& matrix)
: data_(matrix)
- {
- data_ /= matrix.determinant();
- }
+ {}
/** \brief Const access to the matrix entries */
const Dune::FieldMatrix<T,N,N>& data() const
--
GitLab