gtCrystSignedHKLs : added output hklsp and thtype

Signed-off-by: Laura Nervo <laura.nervo@esrf.fr>

zUtil_Cryst/gtCrystSignedHKLs.m

-function [allshkls, allhklinds, mult]  = gtCrystSignedHKLs(hkl, symm)
+function [allshkls, allhklinds, mult, hklsp, thtype]  = gtCrystSignedHKLs(hkl, symm)
 % GTCRYSTSIGNEDHKLS  Gives a complete list of signed hkl planes, by
 %     applying symmetry operators of the specified spacegroup to the input
 %     list of plane families.
@@ -6,21 +6,29 @@ function [allshkls, allhklinds, mult]  = gtCrystSignedHKLs(hkl, symm)
 %     'lattice_system' and 'crystal_system' instead of spacegroups should
 %     be considered.
 %
-%     [allshkls, allhklinds, mult] = gtCrystSignedHKLs(hkl,symm)
-%     ----------------------------------------------------------
+%     [allshkls, allhklinds, mult, hklsp, thtype] = gtCrystSignedHKLs(hkl, symm)
+%     --------------------------------------------------------------------------
 %
 %     INPUT:
-%       hkl        = <double>     plane families (one in a column)
+%       hkl        = <double>     plane families (one in a column) as in
+%                                 parameters.cryst
 %                                 (3,n) for cubic lattices 
 %                                 (4,n) for hexagonal lattices
 %       symm       = <struct>     structure containing .g3/.g with symmetry
 %                                 operators
 %
 %     OUTPUT:
-%       allshkls   =              complete list of signed hkl planes
+%       allshkls   = <double>     complete list of signed hkl planes
 %                                 (including Friedel pairs)
-%       allhklinds =              indices in allshkls according to input hkl 
-%       mult       = <int8>       multiplicity of hkl families (1,k)
+%       allhklinds = <double>     indices in allshkls according to input hkl 
+%       mult       = <double>     multiplicity of hkl families (1,k)
+%       hklsp      = <double>     list of all the hkls (friedel not considered)
+%       thtype     = <double>     theta types among the given reflections
+%                                 according to input hkl
+%
+%
+%     Version 003 21-08-2013 by LNervo
+%       Added output hklsp and thtype to simplify calculation
 %
 %     Version 002 10-07-2013 by LNervo
 %       Added check of class for hkl list
@@ -54,6 +62,9 @@ ncols = size(hkl,2);
 
 allshkls   = [];
 allhklinds = [];
+mult       = [];
+hklsp      = [];
+thtype     = [];
 
 % Loop through input hkl types
 for ii = 1:size(hkl,1)
@@ -70,12 +81,14 @@ for ii = 1:size(hkl,1)
     newshkls = unique(shkls, 'rows');
     % sort according to h > k > l
     newshkls = sortrows(newshkls, [-1 -2 -ncols]);
+    nhkls = size(newshkls,1)/2;
     
-	  allshkls = [allshkls; newshkls];
-    
-    mult(ii,1)  = size(newshkls,1);
-    
-    allhklinds = [allhklinds; repmat(ii,size(newshkls,1),1)];
+    mult       = [mult; nhkls*2];
+	  allshkls   = [allshkls; newshkls];
+    allhklinds = [allhklinds; repmat(ii,nhkls*2,1)];
+    %newshkls   = sortrows(newshkls, [ncols 1 2]); 
+    hklsp      = [hklsp; newshkls(nhkls+1:end,:)];
+    thtype     = [thtype; repmat(ii,nhkls,1)];
 
 end