Skip to content
Snippets Groups Projects
Commit cb045e4e authored by Yoann Guilhem's avatar Yoann Guilhem
Browse files

Code formatting only.

parent c87e199d
No related branches found
No related tags found
No related merge requests found
function gtReadBoundaryProperties(phaseid, boundaries_structure, vol_out, r_vectors) function gtReadBoundaryProperties(phaseid, boundaries_structure, vol_out, r_vectors)
% gtReadBoundaryProperties(phaseid, boundaries_structure, vol_out, r_vectors) % gtReadBoundaryProperties(phaseid, boundaries_structure, vol_out, r_vectors)
disp('have changed this to axis convention matlab1 // instX , matlab2 // instY, matlab3 // instZ') disp('have changed this to axis convention matlab1 // instX , matlab2 // instY, matlab3 // instZ');
disp('AK, 11/8/2010') disp('AK, 11/8/2010');
pause(3) pause(3);
%continue to add to boundaries_structure.mat - orientations, %continue to add to boundaries_structure.mat - orientations,
%disorientation, r-vectors etc. %disorientation, r-vectors etc.
...@@ -23,9 +23,9 @@ symm = gtCrystGetSymmetryOperators([], spacegroup); ...@@ -23,9 +23,9 @@ symm = gtCrystGetSymmetryOperators([], spacegroup);
symm = {symm.g3}; symm = {symm.g3};
% loop through the boundaries % loop through the boundaries
for ii=1:size(boundaries_structure, 2) for ii = 1:size(boundaries_structure, 2)
if boundaries_structure(ii).count<4 %can still be useful to calc index planes where we have only one of the grains if (boundaries_structure(ii).count < 4) %can still be useful to calc index planes where we have only one of the grains
continue; continue;
else else
%Read in the r-vector data for the two grains %Read in the r-vector data for the two grains
...@@ -33,23 +33,23 @@ for ii=1:size(boundaries_structure, 2) ...@@ -33,23 +33,23 @@ for ii=1:size(boundaries_structure, 2)
grain2=boundaries_structure(ii).grain2; grain2=boundaries_structure(ii).grain2;
if grain1==0 if grain1==0
g1=[]; g1 = [];
g1equiv=[]; g1equiv = [];
elseif grain1==-1 elseif grain1==-1
continue; continue;
else else
R1= r_vectors(grain1,2:4); R1 = r_vectors(grain1, 2:4);
boundaries_structure(ii).grain1_R_vector=R1; boundaries_structure(ii).grain1_R_vector = R1;
g1 = gtMathsRod2OriMat(R1); g1 = gtMathsRod2OriMat(R1);
end end
if grain2==0 if grain2==0
g2=[]; g2 = [];
elseif grain2==-1 elseif grain2==-1
continue continue;
else else
R2= r_vectors(grain2,2:4); R2 = r_vectors(grain2,2:4);
boundaries_structure(ii).grain2_R_vector=R2; boundaries_structure(ii).grain2_R_vector = R2;
g2 = gtMathsRod2OriMat(R2); g2 = gtMathsRod2OriMat(R2);
end end
...@@ -58,7 +58,7 @@ for ii=1:size(boundaries_structure, 2) ...@@ -58,7 +58,7 @@ for ii=1:size(boundaries_structure, 2)
if ~isempty(g1) && ~isempty(g2) if ~isempty(g1) && ~isempty(g2)
% If we have both grains % If we have both grains
% Need to search symmetry equivalents for the minimum misorientation % Need to search symmetry equivalents for the minimum misorientation
rot_offset=[]; rot_offset = [];
% warning('g*symm maybe should be symm*g') % warning('g*symm maybe should be symm*g')
for jj = 1:length(symm) for jj = 1:length(symm)
g1equiv = g1*symm{jj}; g1equiv = g1*symm{jj};
...@@ -171,11 +171,11 @@ for ii=1:size(boundaries_structure, 2) ...@@ -171,11 +171,11 @@ for ii=1:size(boundaries_structure, 2)
boundaries_structure(ii).misorientation_axis_hex=hkil; boundaries_structure(ii).misorientation_axis_hex=hkil;
end end
else else
disp('crystallography not supported!') disp('crystallography not supported!');
end end
end % skip small boundaries end % skip small boundaries
end end
save('boundaries_structure.mat', 'boundaries_structure') save('boundaries_structure.mat', 'boundaries_structure');
end % end of function end % end of function
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