diff --git a/Numeric/fullMatrix.cpp b/Numeric/fullMatrix.cpp
index f7d154e46f0af611d77f6785a0374be3ccceb46e..9875746c20fe061fd42178fc152cd6c0a7b6ffd4 100644
--- a/Numeric/fullMatrix.cpp
+++ b/Numeric/fullMatrix.cpp
@@ -311,6 +311,9 @@ void fullMatrix<double>::registerBindings(binding *b)
   cm = cb->addMethod("gemm_naive", &fullMatrix<double>::gemm_naive);
   cm->setArgNames("A","B","alpha","beta",NULL);
   cm->setDescription("this = beta*this + alpha * (A.B)");
+  cm = cb->addMethod("print", &fullMatrix<double>::print);
+  cm->setArgNames("name",NULL);
+  cm->setDescription("print the matrix");
   cm = cb->setConstructor<fullMatrix<double>,int,int>();
   cm->setDescription ("A new matrix of size 'nRows' x 'nColumns'");
   cm->setArgNames("nRows","nColumns",NULL);
diff --git a/Numeric/fullMatrix.h b/Numeric/fullMatrix.h
index 80addf83f6625d9330b4fc87aad8a83e869a90e6..f2a2cdefa928c05b9d0507adafc438809b336446 100644
--- a/Numeric/fullMatrix.h
+++ b/Numeric/fullMatrix.h
@@ -164,6 +164,17 @@ class fullMatrix
     _own_data = false;
     _data = original._data + c_start * _r;
   }
+  void setAsShapeProxy(fullMatrix<scalar> &original, int nbRow, int nbCol)
+  {
+    if(_data && _own_data)
+      delete [] _data;
+    _c = nbCol;
+    _r = nbRow;
+    if(_c*_r != original._c*original._r)
+      Msg::Error("trying to reshape a fullMatrix without conserving the total number of entries");
+    _own_data = false;
+    _data = original._data;
+  }
   inline int size1() const { return _r; }
   inline int size2() const { return _c; }
   fullMatrix<scalar> & operator = (const fullMatrix<scalar> &other)
@@ -356,9 +367,9 @@ class fullMatrix
   }
 #endif
   ;
-  void print(const char *name="") const 
+  void print(const std::string name = "") const 
   {
-    printf("Printing matrix %s:\n", name);
+    printf("Printing matrix %s:\n", name.c_str());
     int ni = size1();
     int nj = size2();
     for(int I = 0; I < ni; I++){