Geo 4x4 tensors: fixed non-element-wise operations

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

zUtil_Geo/gtGeoLab2Sam.m

@@ -104,8 +104,8 @@ if (use_diffractometer)
         labXYZ_rot = labXYZ_rot(1:3, 1, :);
         labXYZ_rot = permute(labXYZ_rot, [3 1 2]);
     else
-        labXYZ_rot = permute(labXYZ_rot, [2 3 1]);
-        labXYZ_rot = labXYZ_rot(:, 1:3);
+        labXYZ_rot = labXYZ_rot(1:3, :, :);
+        labXYZ_rot = reshape(labXYZ_rot, 3, [])';
     end
 else
     ones_in = ones(size(labXYZ, 1), 1);

zUtil_Geo/gtGeoProjForReconstruction.m

@@ -105,7 +105,6 @@ function proj_geom = gtGeoProjForReconstruction(projvec_sam, omega, vol_center,.
     end
     projvec_rec = gtMathsNormalizeVectorsList(projvec_rec);
 
-
     if (strcmpi(rectype, 'ASTRA_grain') || strcmpi(rectype, 'ASTRA_absorption'))
         % Volume center in REC reference
         if (strcmpi(rectype, 'ASTRA_absorption'))

zUtil_Geo/gtGeoSam2Lab.m

@@ -105,8 +105,7 @@ if (use_diffractometer)
         labXYZ = permute(labXYZ, [3 1 2]);
     else
         labXYZ = labXYZ(1:3, :, :);
-        labXYZ = reshape(labXYZ, 3, []);
-        labXYZ = permute(labXYZ, [2 3 1]);
+        labXYZ = reshape(labXYZ, 3, [])';
     end
 else
     ones_in = ones(size(samXYZ, 1), 1);