diff --git a/grid/albertagrid/albertagrid.cc b/grid/albertagrid/albertagrid.cc index 3e35ee0476c8d19aa8933c992fb64e522b2ee344..deaa32bcd0fbd025b4aaeade370254d3b50bee19 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;