Simplified and optimized gtVectorLab2Cryst and gtVectorCryst2Lab.

Signed-off-by: Yoann Guilhem <yoann.guilhem@esrf.fr>

zUtil_Cryst/gtVectorCryst2Lab.m

@@ -88,14 +88,17 @@ if M == 1 && N > 1
     Vs = squeeze(sum(Vc(:, [1 1 1], ones(N, 1)).*g, 1)).';
 
 elseif M >= 1 && N == 1
-    % Transpose g to get gT, reshape it and expand it along 1st axis
-    gT = reshape(g.', [1 3 3]);
-
-    % Reshape Vc
-    Vc = reshape(Vc, [M 1 3]);
+%    % Transpose g to get gT, reshape it and expand it along 1st axis
+%    gT = reshape(g.', [1 3 3]);
+%
+%    % Reshape Vc
+%    Vc = reshape(Vc, [M 1 3]);
+%
+%    % Multiply and sum along 3rd axis to reproduce Vs = gT . Vc
+%    Vs = sum(bsxfun(@times, gT, Vc), 3);
 
-    % Multiply and sum along 3rd axis to reproduce Vs = gT . Vc
-    Vs = sum(bsxfun(@times, gT, Vc), 3);
+    % Can be simply done with this formulae 
+    Vs = Vc * g;
 
 elseif M == N && diag
     % Transpose and reshape Vc

zUtil_Cryst/gtVectorLab2Cryst.m

@@ -88,14 +88,17 @@ if M == 1 && N > 1
     Vc = squeeze(sum(Vs([1 1 1], :, ones(N, 1)).*g, 2)).';
 
 elseif M >= 1 && N == 1
-    % Transpose g to get gT, reshape it and expand it along 1st axis
-    g = reshape(g, [1 3 3]);
-
-    % Reshape Vs
-    Vs = reshape(Vs, [M 1 3]);
+%    % Transpose g to get gT, reshape it and expand it along 1st axis
+%    g = reshape(g, [1 3 3]);
+%
+%    % Reshape Vs
+%    Vs = reshape(Vs, [M 1 3]);
+%
+%    % Multiply and sum along 3rd axis to reproduce Vc = g . Vs
+%    Vc = sum(bsxfun(@times, g, Vs), 3);
 
-    % Multiply and sum along 3rd axis to reproduce Vc = g . Vs
-    Vc = sum(bsxfun(@times, g, Vs), 3);
+    % Can be simply done with this formulae 
+    Vc = (g * Vs.').';
 
 elseif M == N && diag
     % Transpose and reshape Vc