Cxx/gtPlaceSubVol: safer dealing with matlab volumes

Signed-off-by: Nicola Vigano <nicola.vigano@esrf.fr>

zUtil_Cxx/include/GtMatlabVolumeLooper.h

@@ -10,6 +10,14 @@
 
 #include "mex.h"
 
+extern "C"
+{
+  bool mxUnshareArray(mxArray *array_ptr, bool noDeepCopy);
+  mxArray *mxCreateSharedCopy(const mxArray *pr);
+  mxArray *mxCreateSharedDataCopy(const mxArray *pr);
+  mxArray *mxUnreference(const mxArray *pr);
+}
+
 template<typename Type>
 class GtMatlabPointer {
   const mwSize multiplier;

zUtil_Cxx/internal_gtAssignGrainToVol.cpp

@@ -74,7 +74,7 @@ void mexFunction( int nlhs, mxArray * plhs[], int nrhs, const mxArray * prhs[] )
   /* -- No assert on correctness of parameters -- */
 
   /* Don't copy, just modify the incoming matrix */
-  plhs[0] = (mxArray *) prhs[0];
+  plhs[0] = mxCreateSharedDataCopy(prhs[0]);
 
   switch(mxGetClassID(prhs[1])) {
     case mxDOUBLE_CLASS: {

zUtil_Cxx/internal_gtAssignGrainToVol_interf.cpp

@@ -75,7 +75,7 @@ void mexFunction( int nlhs, mxArray * plhs[], int nrhs, const mxArray * prhs[] )
   /* -- No assert on correctness of parameters -- */
 
   /* Don't copy, just modify the incoming matrix */
-  plhs[0] = (mxArray *) prhs[0];
+  plhs[0] = mxCreateSharedDataCopy(prhs[0]);
 
   switch(mxGetClassID(prhs[1])) {
     case mxDOUBLE_CLASS: {

zUtil_Cxx/internal_gtAssignGrainToVol_sum.cpp

@@ -73,7 +73,7 @@ void mexFunction( int nlhs, mxArray * plhs[], int nrhs, const mxArray * prhs[] )
   /* -- No assert on correctness of parameters -- */
 
   /* Don't copy, just modify the incoming matrix */
-  plhs[0] = (mxArray *) prhs[0];
+  plhs[0] = mxCreateSharedDataCopy(prhs[0]);
 
   switch(mxGetClassID(prhs[1])) {
     case mxDOUBLE_CLASS: {

zUtil_Cxx/internal_gtDilateGrains_voxel.cpp

@@ -83,7 +83,7 @@ void mexFunction( int nlhs, mxArray * plhs[], int nrhs, const mxArray * prhs[] )
   /* -- No assert on correctness of parameters -- */
 
   /* Don't copy, just modify the incoming matrix */
-  plhs[0] = (mxArray *) prhs[0];
+  plhs[0] = mxCreateSharedDataCopy(prhs[0]);
 
   switch(mxGetClassID(prhs[1])) {
     case mxDOUBLE_CLASS: {