gtCalculateGrain: fixed to work with older and newer versions of grain structures

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

4_grains/gtCalculateGrain.m

@@ -159,8 +159,9 @@ St4(:, :, todel, :) = [];
 pl(:, todel)        = [];
 
 thetatypesp = cryst.thetatypesp(~todel);
+thetatypesp = reshape(thetatypesp, [], 1);
 hkl = cryst.hkl(:, thetatypesp)';
-hklsp  = cryst.hklsp(:, ~todel);
+hklsp  = cryst.hklsp(:, ~todel)';
 
 %%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 % % Following the convention of the matching output, the omega value
@@ -217,10 +218,10 @@ grain.allblobs.plorig = [pl'; pl'; pl'; pl'];
 grain.allblobs.pl = pls';
 grain.allblobs.pllab = pllab';
 grain.allblobs.hkl = [hkl; hkl; hkl; hkl];
-grain.allblobs.hklsp = [hklsp'; -hklsp'; hklsp'; -hklsp'];
+grain.allblobs.hklsp = [hklsp; -hklsp; hklsp; -hklsp];
 grain.allblobs.sintheta = [sinth; sinth; sinth; sinth];
 grain.allblobs.theta = [theta; theta; theta; theta];
-grain.allblobs.thetatype = [thetatypesp'; thetatypesp'; thetatypesp'; thetatypesp'];
+grain.allblobs.thetatype = [thetatypesp; thetatypesp; thetatypesp; thetatypesp];
 grain.allblobs.omega = om';
 grain.allblobs.omind = omind';
 grain.allblobs.dvec = dveclab';