diff --git a/7_fed2/gtFedPredictOmegaMultiple.m b/7_fed2/gtFedPredictOmegaMultiple.m index 44ef0ea8156f9d5044f92ff4051b873e973e8fed..4f23664d0679a066c8bcff93cc55095191a43307 100644 --- a/7_fed2/gtFedPredictOmegaMultiple.m +++ b/7_fed2/gtFedPredictOmegaMultiple.m @@ -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