Skip to content
Snippets Groups Projects
Commit 06aac2e9 authored by Nicola Vigano's avatar Nicola Vigano Committed by Nicola Vigano
Browse files

Omega prediction: do not compute stuff that is not emitted

parent 36deb818
No related branches found
No related tags found
No related merge requests found
......@@ -115,20 +115,22 @@ if (~isempty(omind))
% ROTATION MATRICES AND PLANE NORMALS IN LAB
% Get rotation matrix and multpily the input plane normals to get
% them in the diffracting position
rot = gtMathsRotationTensor(om, rotcomp);
if (nargout > 1)
% Get rotation matrix and multiply the input plane normals to get
% them in the diffracting position
rot = gtMathsRotationTensor(om, rotcomp);
% expand for multiplication
plt = reshape(pl, 1, 3, []);
plt = plt([1 1 1], :, :);
% expand for multiplication
plt = reshape(pl, 1, 3, []);
plt = plt([1 1 1], :, :);
plrot = rot .* plt;
plrot = sum(plrot, 2);
plrot = reshape(plrot, 3, []);
plrot = rot .* plt;
plrot = sum(plrot, 2);
plrot = reshape(plrot, 3, []);
plrot = ss3 .* plrot;
plsigned = ss3 .* pl;
plrot = ss3 .* plrot;
plsigned = ss3 .* pl;
end
else
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% All four omega indices
......
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