From 76e7a0f48359d01bb9fb3544e088233b80086faa Mon Sep 17 00:00:00 2001
From: Gauthier Becker <gauthierbecker@gmail.com>
Date: Sat, 23 Jul 2011 09:28:35 +0000
Subject: [PATCH] Fix bug with new definition of fullMatrix::setAll() changed
in revision 9625. In fact fullMatrix::invert doesn't work anymore when
meshing with high order (degree 3). I resize the matrix to avoid fatal error
due to bad result matrix size.
---
Numeric/fullMatrix.cpp | 1 +
1 file changed, 1 insertion(+)
diff --git a/Numeric/fullMatrix.cpp b/Numeric/fullMatrix.cpp
index d35cf99734..b8ecb795b2 100644
--- a/Numeric/fullMatrix.cpp
+++ b/Numeric/fullMatrix.cpp
@@ -297,6 +297,7 @@ bool fullMatrix<double>::invert(fullMatrix<double> &result) const
{
int M = size1(), N = size2(), lda = size1(), info;
int *ipiv = new int[std::min(M, N)];
+ result.resize(M,N,false);
result.setAll(*this);
F77NAME(dgetrf)(&M, &N, result._data, &lda, ipiv, &info);
if(info == 0){
--
GitLab