diff --git a/Numeric/fullMatrix.cpp b/Numeric/fullMatrix.cpp
index d35cf997345f7e0b82ef18d93edf9d6bf4afb0a7..b8ecb795b2fab9a381ef97906fbb10c9e661e806 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){