Skip to content
Snippets Groups Projects
Commit 4bec7dbf authored by Laura Nervo's avatar Laura Nervo
Browse files

Bug fixed with gtCrystSignedHKLs arguments : spacegroup -> symm


Signed-off-by: default avatarLaura Nervo <laura.nervo@esrf.fr>
parent 7550c42e
No related branches found
No related tags found
No related merge requests found
...@@ -110,8 +110,10 @@ d0 = parameters.cryst.dsp; ...@@ -110,8 +110,10 @@ d0 = parameters.cryst.dsp;
% hkl = []; % hkl = [];
% thtype = []; % thtype = [];
% %
% symm = gtCrystGetSymmetryOperators([], spacegroup);
%
% for i = 1:length(hkltypes_used) % for i = 1:length(hkltypes_used)
% allhkls_i = gtCrystSignedHKLs(hkltypes_used(i,:),spacegroup); % allhkls_i = gtCrystSignedHKLs(hkltypes_used(i,:),symm);
% nhkls = length(allhkls_i)/2; % nhkls = length(allhkls_i)/2;
% hklsp = [hklsp; allhkls_i(nhkls+1:end,:)]; % hklsp = [hklsp; allhkls_i(nhkls+1:end,:)];
% hkl = [hkl; repmat(hkltypes_used(i,:),nhkls,1)]; % hkl = [hkl; repmat(hkltypes_used(i,:),nhkls,1)];
......
...@@ -65,6 +65,7 @@ intsc.Rvec = []; ...@@ -65,6 +65,7 @@ intsc.Rvec = [];
fzone_ext = gtCrystRodriguesFundZone(spacegroup,tol.sm); % fund. zone with safety margin fzone_ext = gtCrystRodriguesFundZone(spacegroup,tol.sm); % fund. zone with safety margin
fzone_acc = gtCrystRodriguesFundZone(spacegroup,0); % exact fund. zone fzone_acc = gtCrystRodriguesFundZone(spacegroup,0); % exact fund. zone
symm = gtCrystGetSymmetryOperators([], spacegroup);
%% Determine line coordinates in Rodrigues space %% Determine line coordinates in Rodrigues space
...@@ -91,7 +92,7 @@ else ...@@ -91,7 +92,7 @@ else
%[Rvec,isec1,isec2] = gtRodriguesVectors2(plnorm(i,:),hkl(i,:),spacegroup,latticepar,tol.sm); %[Rvec,isec1,isec2] = gtRodriguesVectors2(plnorm(i,:),hkl(i,:),spacegroup,latticepar,tol.sm);
% All signed hkl-s for given hkl family % All signed hkl-s for given hkl family
shkls = gtCrystSignedHKLs(hkl(i,:),spacegroup); shkls = gtCrystSignedHKLs(hkl(i,:), symm);
[Rvec,isec1,isec2] = gtCrystRodriguesVector(plnorm(i,:),shkls,Bmat,fzone_ext); [Rvec,isec1,isec2] = gtCrystRodriguesVector(plnorm(i,:),shkls,Bmat,fzone_ext);
......
...@@ -258,11 +258,11 @@ if (any(gtIndexAllGrainValues(grain1,'R_onedge',[],1,[])) || ~isempty(drot)) ...@@ -258,11 +258,11 @@ if (any(gtIndexAllGrainValues(grain1,'R_onedge',[],1,[])) || ~isempty(drot))
else else
error('Field parameters.cryst.hkl is not set correctly.') error('Field parameters.cryst.hkl is not set correctly.')
end end
symm = gtCrystGetSymmetryOperators([], spacegroup);
% All signed hkl-s for the given family (both hkl and -h-k-l) % All signed hkl-s for the given family (both hkl and -h-k-l)
shkl1 = gtCrystSignedHKLs(hklsp(:,1)',spacegroup)'; shkl1 = gtCrystSignedHKLs(hklsp(:,1)',symm)';
shkl2 = gtCrystSignedHKLs(hklsp(:,2)',spacegroup)'; shkl2 = gtCrystSignedHKLs(hklsp(:,2)',symm)';
shkl3 = gtCrystSignedHKLs(hklsp(:,3)',spacegroup)'; shkl3 = gtCrystSignedHKLs(hklsp(:,3)',symm)';
end end
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment