Skip to content
Snippets Groups Projects
Commit 52364b63 authored by Yoann Guilhem's avatar Yoann Guilhem Committed by Nicola Vigano
Browse files

Simplify symmetry operator extraction.

parent a7bd853e
No related branches found
No related tags found
No related merge requests found
......@@ -21,11 +21,11 @@ function hkl_reflections = gtGetReflections(hkl, spacegroup)
% from Risoe group
% Get symmetry operators
tmp = gtCrystGetSymmetryOperators([], spacegroup);
symm = gtCrystGetSymmetryOperators([], spacegroup);
if size(hkl,2) == 4
symm = {tmp.g};
symm = {symm.g};
elseif size(hkl,2) == 3
symm = {tmp.g3};
symm = {symm.g3};
end
% Apply to input reflection
......
......@@ -41,11 +41,11 @@ load('parameters.mat');
spacegroup = parameters.cryst.spacegroup;
% Get symmetry operators
tmp = gtCrystGetSymmetryOperators([], spacegroup);
if size(hkl,2) == 4
symm = {tmp.g};
elseif size(hkl,2) == 3
symm = {tmp.g3};
symm = gtCrystGetSymmetryOperators([], spacegroup);
if size(hkl, 2) == 4
symm = {symm.g};
elseif size(hkl, 2) == 3
symm = {symm.g3};
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
......
......@@ -19,8 +19,8 @@ for i=1:length(boundaries_structure)
end
% Get symmetry operators
tmp = gtCrystGetSymmetryOperators([], spacegroup);
symm = {tmp.g3};
symm = gtCrystGetSymmetryOperators([], spacegroup);
symm = {symm.g3};
% loop through the boundaries
for ii=1:size(boundaries_structure, 2)
......
......@@ -23,11 +23,11 @@ function [allshkls, allhklinds, mult] = gtCrystSignedHKLs(hkl,spacegroup)
%
% Get symmetry operators
tmp = gtCrystGetSymmetryOperators([], spacegroup);
symm = gtCrystGetSymmetryOperators([], spacegroup);
if size(hkl,2) == 4
symm = {tmp.g};
symm = {symm.g};
elseif size(hkl,2) == 3
symm = {tmp.g3};
symm = {symm.g3};
end
allshkls = [];
......
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