Fixed variables order for allblobs

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

4_grains/gtCalculateGrain_p.m

@@ -553,7 +553,8 @@ for det_ind = app.det_ind
         end
 
         for ii_gr = 1:num_grains
-            % Initialse output variables
+            allblobs = gtGrainAllblobsDefinition();
+
             allblobs.type = mode;
 
             allblobs.plcry = plcry(:, :, ii_gr);
@@ -598,6 +599,8 @@ for det_ind = app.det_ind
             grain{ii_gr}.allblobs(det_ind) = allblobs;
         end
     else
+        allblobs = gtGrainAllblobsDefinition();
+
         allblobs.type = mode;
         allblobs.plcry = plcry;
         allblobs.plorig = plorig;

7_fed2/gtFedPredictOmegaMultiple.m

@@ -48,10 +48,12 @@ function [om, plrot, plsigned, rot, omind] = gtFedPredictOmegaMultiple(...
 %    plrot    - the plane normals rotated in the diffraction positions;
 %               no. 3 and 4 are inverted after rotation (pointing oppposite
 %               to beam) to ease calculation at a later stage;
-%               size (3xn or 3xnx4)
+%               size (3xn or 3xnx4) (in LAB coordinates)
 %    plsigned - the original input plane normals with no. 3 and 4 inverted;
+%               in Sam coordinates
 %                 plrot(:,i,j) = rot(:,:,i,j)*plsigned(:,i,j) the  is where (3xn or 3xnx4)               
 %    rot      - rotation matrices (3x3xn or 3x3xnx4)
+%               Saml->Lab for column vectors, Lab->Sam for row vectors
 %    omind    - omega indices (1xn or 4xn) 
 %
 
@@ -165,7 +167,7 @@ else
     om([1 2], chom)    = om([2 1], chom);
     omind([1 2], chom) = omind([2 1], chom);
 
-    % ROTATION MATRICES AND PLANE NORMALS IN LAB
+    % ROTATION MATRICES AND PLANE NORMALS IN LAB (plsigned still in SAM)
 
     plsigned = zeros(3, nn, 4);
     plsigned(:, :, 1) = pl;