From 85115d92c19c687a186282bfe85efdb39e2954c8 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Robert=20Kl=C3=B6fkorn?= <robertk@dune-project.org>
Date: Wed, 5 Oct 2005 19:01:11 +0000
Subject: [PATCH] jacobianInverse returns transposed inverse matrix.
[[Imported from SVN: r3041]]
---
grid/albertagrid/albertagrid.cc | 8 +-------
1 file changed, 1 insertion(+), 7 deletions(-)
diff --git a/grid/albertagrid/albertagrid.cc b/grid/albertagrid/albertagrid.cc
index 3e35ee047..deaa32bcd 100644
--- a/grid/albertagrid/albertagrid.cc
+++ b/grid/albertagrid/albertagrid.cc
@@ -470,13 +470,7 @@ namespace Dune
// Jinv = A^-1
assert( builtElMat_ == true );
- FieldMatrix< double , mydim , mydim > tmp;
- //elDet_ = std::abs( FMatrixHelp::invertMatrix(elMat_,Jinv_) );
- elDet_ = std::abs( FMatrixHelp::invertMatrix(elMat_,tmp) );
-
- for(int i=0; i<mydim; i++)
- for(int j=0 ; j<mydim ; j++)
- Jinv_[i][j] = tmp[j][i];
+ elDet_ = std::abs( FMatrixHelp::invertMatrix_retTransposed(elMat_,Jinv_) );
assert(elDet_ > 1.0E-25);
calcedDet_ = true;
--
GitLab