Skip to content
Snippets Groups Projects
gtAngularConsCos.m 296 B
function cos=gtAngularConsCos(spacegroup,f1,f2)

ref=parameters.cryst.hklsp;

cos=[];

v1=gtGetReflections(ref(f1,:)) ;
v2=gtGetReflections(ref(f2,:)) ;
for i=1:size(v1,1)
  for j=1:size(v2,1)
    cos=[cos; abs(v1(i,:)*v2(j,:)')/norm(v1)/norm(v2)] ;
  end
end
cos=unique(cos);
cos=[1; cos];

end