gtPlaceSubVolume: restored previous performance of C++ functions

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

zUtil_Cxx/3D_ops/gtCxxPlaceSubVolumeAssign.cpp

@@ -15,16 +15,17 @@ void mexFunction( int nlhs, mxArray * plhs[], int nrhs, const mxArray * prhs[] )
     return;
   }
 
-  mxArray * const mat_output = (mxArray *)prhs[0];
   const mxArray * const mat_input = prhs[1];
   const mxArray * const mat_shifts_op = prhs[2];
   const mxArray * const mat_shifts_ip = prhs[3];
   const mxArray * const mat_dims = prhs[4];
 
-  /* Incoming volumes dimensions */
-
   /* Don't copy, just modify the incoming matrix */
-  mxUnshareArray(mat_output, true);
+  plhs[0] = mxCreateSharedDataCopy(prhs[0]);
+  mxArray * const mat_output = plhs[0];
+
+  //  mxArray * const mat_output = (mxArray *)prhs[0];
+  //  mxUnshareArray(mat_output, true);
 
   GT3D::GtMatlabSubVolumePlacer placer(mat_output, mat_input, mat_shifts_op,
       mat_shifts_ip, mat_dims);

zUtil_Cxx/3D_ops/gtCxxPlaceSubVolumeInterf.cpp

@@ -15,16 +15,17 @@ void mexFunction( int nlhs, mxArray * plhs[], int nrhs, const mxArray * prhs[] )
     return;
   }
 
-  mxArray * const mat_output = (mxArray *)prhs[0];
   const mxArray * const mat_input = prhs[1];
   const mxArray * const mat_shifts_op = prhs[2];
   const mxArray * const mat_shifts_ip = prhs[3];
   const mxArray * const mat_dims = prhs[4];
 
-  /* Incoming volumes dimensions */
-
   /* Don't copy, just modify the incoming matrix */
-  mxUnshareArray(mat_output, true);
+  plhs[0] = mxCreateSharedDataCopy(prhs[0]);
+  mxArray * const mat_output = plhs[0];
+
+  //  mxArray * const mat_output = (mxArray *)prhs[0];
+  //  mxUnshareArray(mat_output, true);
 
   GT3D::GtMatlabSubVolumePlacer placer(mat_output, mat_input, mat_shifts_op,
       mat_shifts_ip, mat_dims);

zUtil_Cxx/3D_ops/gtCxxPlaceSubVolumeSum.cpp

@@ -15,16 +15,17 @@ void mexFunction( int nlhs, mxArray * plhs[], int nrhs, const mxArray * prhs[] )
     return;
   }
 
-  mxArray * const mat_output = (mxArray *)prhs[0];
   const mxArray * const mat_input = prhs[1];
   const mxArray * const mat_shifts_op = prhs[2];
   const mxArray * const mat_shifts_ip = prhs[3];
   const mxArray * const mat_dims = prhs[4];
 
-  /* Incoming volumes dimensions */
-
   /* Don't copy, just modify the incoming matrix */
-  mxUnshareArray(mat_output, true);
+  plhs[0] = mxCreateSharedDataCopy(prhs[0]);
+  mxArray * const mat_output = plhs[0];
+
+  //  mxArray * const mat_output = (mxArray *)prhs[0];
+  //  mxUnshareArray(mat_output, true);
 
   GT3D::GtMatlabSubVolumePlacer placer(mat_output, mat_input, mat_shifts_op,
       mat_shifts_ip, mat_dims);

zUtil_Imaging/gtPlaceSubVolume.m

@@ -88,11 +88,11 @@ function output = gtPlaceSubVolume(output, input, shift, index, assign_op, use_c
 
         switch (assign_op)
             case {'zero', 'conflict', 'adaptive'} % Mark overlapping...
-                gtCxxPlaceSubVolumeInterf(output, input, shifts_op, shifts_ip, dims);
+                output = gtCxxPlaceSubVolumeInterf(output, input, shifts_op, shifts_ip, dims);
             case 'summed'
-                gtCxxPlaceSubVolumeSum(output, input, shifts_op, shifts_ip, dims);
+                output = gtCxxPlaceSubVolumeSum(output, input, shifts_op, shifts_ip, dims);
             case {'assign', 'parent'}
-                gtCxxPlaceSubVolumeAssign(output, input, shifts_op, shifts_ip, dims);
+                output = gtCxxPlaceSubVolumeAssign(output, input, shifts_op, shifts_ip, dims);
             otherwise
                 error('PLACE:wrong_argument', 'No option for "%s"', assign_op);
         end