gtHex2Cart : Transforms from direct hexagonal to cartesian coordinates and...

gtHex2Cart : Transforms from direct hexagonal to cartesian coordinates and normalize. Same as the old function gtCrystHKL2Cart

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

git-svn-id: https://svn.code.sf.net/p/dct/code/trunk@981 4c865b51-4357-4376-afb4-474e03ccb993

zUtil_Cryst/gtHex2Cart.m

-function cartesian = gtHex2Cart(hkil, latticepar)
+function normalised_hkls = gtHex2Cart(all_hkils, latticepar)
 % GTHEX2CART  
-%     cartesian = gtHex2Cart(hkil, latticepar)
-%     ----------------------------------------
+%     normalised_hkls = gtHex2Cart(all_hkils, latticepar)
+%     ---------------------------------------------------
 %     Transforms from direct hexagonal to cartesian coordinates
 %     and normalise
 %
 % Same as the old function gtCrystHKL2Cart
 
 
-for i=1:size(hkil,1)
-    cartesian(i,1) = hkil(i,1) + 0.5 * hkil(i,2);
-    cartesian(i,2) = 3^0.5/2 *hkil(i,2);
-    cartesian(i,3) = hkil(i,4);
-
-    % allow lattice parameters
-    cartesian(i,1) = cartesian(i,1)*2/(sqrt(3)*latticepar(1));
-    cartesian(i,2) = cartesian(i,2)*2/(sqrt(3)*latticepar(1));
-    cartesian(i,3) = cartesian(i,3)/latticepar(3);
-end
+% going to cartesian reciprocal space + normalisation of the cartesian
+% space
+all_hkls(:, 1) = all_hkils(:, 1) + 0.5 * all_hkils(:, 2);
+all_hkls(:, 2) = sqrt(3)/2 * all_hkils(:, 2);
+all_hkls(:, 3) = all_hkils(:, 4);
 
+all_hkls(:, 1) = all_hkls(:, 1)*2/(sqrt(3)*latticepar(1));
+all_hkls(:, 2) = all_hkls(:, 2)*2/(sqrt(3)*latticepar(1));
+all_hkls(:, 3) = all_hkls(:, 3)/latticepar(3);
 
-% normalise
-tmp = sqrt(sum((cartesian.*cartesian),2));
-cartesian = cartesian./(repmat(tmp,1,3));
+% normalise hkls vector
+tmp = sqrt(sum((all_hkls.*all_hkls),2));
+normalised_hkls = all_hkls./(repmat(tmp,1,3));
 
 end