Skip to content
Snippets Groups Projects
Commit 9a743134 authored by Nicola Vigano's avatar Nicola Vigano
Browse files

Cryst: fixed a few bugs for hkl lists with less than 2 families on the detector

parent 8987eaae
No related branches found
No related tags found
No related merge requests found
...@@ -83,9 +83,9 @@ formfactor(ind) = []; ...@@ -83,9 +83,9 @@ formfactor(ind) = [];
hkl = hkl(ind_tt, :); hkl = hkl(ind_tt, :);
d0 = d0(ind_tt); d0 = d0(ind_tt);
intensity = intensity(ind_tt); intensity = intensity(ind_tt);
thtype = (1:length(hkl))';
mult = mult(ind_tt); mult = mult(ind_tt);
formfactor = formfactor(ind_tt); formfactor = formfactor(ind_tt);
thtype = (1:numel(mult))';
disp(['Families on the detector: ' num2str(size(hkl,1))]) disp(['Families on the detector: ' num2str(size(hkl,1))])
disp('Now reflections are sorted by increasing theta values') disp('Now reflections are sorted by increasing theta values')
......
...@@ -58,7 +58,7 @@ end ...@@ -58,7 +58,7 @@ end
[~, ~, mult, hklsp, thtype] = gtCrystSignedHKLs(hkltypes_used, symm); [~, ~, mult, hklsp, thtype] = gtCrystSignedHKLs(hkltypes_used, symm);
intensity = []; intensity = [];
for ii = 1:length(hkltypes_used) for ii = 1:numel(mult)
nhkls = mult(ii)/2; nhkls = mult(ii)/2;
intensity = [intensity; repmat(cryst.int(ii),nhkls,1)]; intensity = [intensity; repmat(cryst.int(ii),nhkls,1)];
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