variant analysis and doubt

[wanghaijian217]

Discussed in #17

Hi, frank,

Thank you for sharing the code of variant analysis in lath martensite. One of the problems that has puzzled me is that if I draw the ratio of boundary length of V1-V2... V1-V24 variant pairs and the score of variant pairs that follow the K-S relationship.
Only V1-V2 is given in plotMap_variantPairs. m in ORtools; V1-V3 (5); V1-V4; V1-V6 information is not given for other variant pairs.

I wrote a part of it myself, but it didn't work out the same way as other organizations, so I'd like you to look at it or better to provide the relevant code. I would appreciate it if you could provide the variantpairs boundary length density and the variantpairs fraction code.

%% Determine the variant pairs
% Derive variant types 1 to 6 from variants 1 to 24
% variantGrains.prop.variantType = variantGrains.variantId - (variantGrains.packetId-1) * 24/4;
variantGrains.prop.variantType = variantGrains.variantId;
% Get the boundary segments and neighbouring variants
variantBoundary = variantGrains.boundary(variantGrains.CS,variantGrains.CS);
variantBoundaryIds = variantBoundary.grainId;

[c,~] = ismember(variantBoundaryIds,variantGrains.id);
variantBoundaryIds(any(~c,2),:) = [];
variantBoundary(any(~c,2)) = [];

%Identify any special boundaries, here V1-V2 : V1:V6
varTypes = variantGrains(variantGrains.id2ind(variantBoundaryIds)).variantType;
% cond(1,:) = any(ismember(varTypes,1),2) & any(ismember(varTypes,2),2) | ...
%     any(ismember(varTypes,3),2) & any(ismember(varTypes,4),2) | ...
%     any(ismember(varTypes,5),2) & any(ismember(varTypes,6),2);     %V1-V2
% cond(2,:) = any(ismember(varTypes,1),2) & any(ismember(varTypes,3),2) | ...
%     any(ismember(varTypes,1),2) & any(ismember(varTypes,5),2) | ...
%     any(ismember(varTypes,2),2) & any(ismember(varTypes,4),2) | ...
%     any(ismember(varTypes,2),2) & any(ismember(varTypes,6),2) | ...
%     any(ismember(varTypes,3),2) & any(ismember(varTypes,5),2) | ...
%     any(ismember(varTypes,4),2) & any(ismember(varTypes,6),2);     %V1-V3(V5)
% cond(3,:) = any(ismember(varTypes,1),2) & any(ismember(varTypes,6),2) | ...
%     any(ismember(varTypes,2),2) & any(ismember(varTypes,3),2) | ...
%     any(ismember(varTypes,4),2) & any(ismember(varTypes,5),2);     %V1-V6
% cond(4,:) = any(ismember(varTypes,1),2) & any(ismember(varTypes,4),2) | ...
%     any(ismember(varTypes,2),2) & any(ismember(varTypes,5),2) | ...
%     any(ismember(varTypes,3),2) & any(ismember(varTypes,6),2);     %V1-V4

%-------------------------------------------------------------------
vec1=any(ismember(varTypes,1),2);
vec2=any(ismember(varTypes,2),2);
vec3=any(ismember(varTypes,3),2);
vec4=any(ismember(varTypes,4),2);
vec5=any(ismember(varTypes,5),2);
vec6=any(ismember(varTypes,6),2);
vec7=any(ismember(varTypes,7),2);
vec8=any(ismember(varTypes,8),2);
vec9=any(ismember(varTypes,9),2);
vec10=any(ismember(varTypes,10),2);
vec11=any(ismember(varTypes,11),2);
vec12=any(ismember(varTypes,12),2);
vec13=any(ismember(varTypes,13),2);
vec14=any(ismember(varTypes,14),2);
vec15=any(ismember(varTypes,15),2);
vec16=any(ismember(varTypes,16),2);
vec17=any(ismember(varTypes,17),2);
vec18=any(ismember(varTypes,18),2);
vec19=any(ismember(varTypes,19),2);
vec20=any(ismember(varTypes,20),2);
vec21=any(ismember(varTypes,21),2);
vec22=any(ismember(varTypes,22),2);
vec23=any(ismember(varTypes,23),2);
vec24=any(ismember(varTypes,24),2);
%-------------------------------------------------------------
%======================================
cond(1,:) = vec1&vec2|vec3&vec4|vec5&vec6|vec7&vec8|vec9&vec10|vec11&vec12|vec13&vec14|vec15&vec16|vec17&vec18|vec19&vec20|vec21&vec22|vec23&vec24;

cond(2,:) = vec1&vec3|vec2&vec4|vec3&vec5|vec4&vec6|vec5&vec7|vec6&vec8|vec7&vec9|vec8&vec10|vec9&vec11|vec10&vec12|vec11&vec13|vec12&vec14|vec13&vec15|vec14&vec16|vec15&vec17|vec16&vec18|vec17&vec19|vec18&vec20|vec19&vec21|vec20&vec22|vec21&vec23|vec22&vec24|vec1&vec5|vec2&vec6|vec3&vec7|vec4&vec8|vec5&vec9|vec6&vec10|vec7&vec11|vec8&vec12|vec9&vec13|vec10&vec14|vec11&vec15|vec12&vec16|vec13&vec17|vec14&vec18|vec15&vec19|vec16&vec20|vec17&vec21|vec18&vec22|vec19&vec23|vec20&vec24;

cond(3,:) = vec1&vec4|vec2&vec5|vec3&vec6|vec4&vec7|vec5&vec8|vec6&vec9|vec7&vec10|vec8&vec11|vec9&vec12|vec10&vec13|vec11&vec14|vec12&vec15|vec13&vec16|vec14&vec17|vec15&vec18|vec16&vec19|vec17&vec20|vec18&vec21|vec19&vec22|vec20&vec23|vec21&vec24;

cond(4,:) = vec1&vec6|vec2&vec7|vec3&vec8|vec4&vec9|vec5&vec10|vec6&vec11|vec7&vec12|vec8&vec13|vec9&vec14|vec10&vec15|vec11&vec16|vec12&vec17|vec13&vec18|vec14&vec19|vec15&vec20|vec16&vec21|vec17&vec22|vec18&vec23|vec19&vec24;

%======================================

cond(5,:) = vec1&vec7|vec2&vec8|vec3&vec9|vec4&vec10|vec5&vec11|vec6&vec12|vec7&vec13|vec8&vec14|vec9&vec15|vec10&vec16|vec11&vec17|vec12&vec18|vec13&vec19|vec14&vec20|vec15&vec21|vec16&vec22|vec17&vec23|vec18&vec24;

cond(6,:) = vec1&vec8|vec2&vec9|vec3&vec10|vec4&vec11|vec5&vec12|vec6&vec13|vec7&vec14|vec8&vec15|vec9&vec16|vec10&vec17|vec11&vec18|vec12&vec19|vec13&vec20|vec14&vec21|vec15&vec22|vec16&vec23|vec17&vec24;

cond(7,:) = vec1&vec9|vec2&vec10|vec3&vec11|vec4&vec12|vec5&vec13|vec6&vec14|vec7&vec15|vec8&vec16|vec9&vec17|vec10&vec18|vec11&vec19|vec12&vec20|vec13&vec21|vec14&vec22|vec15&vec23|vec16&vec24|vec1&vec19|vec2&vec20|vec3&vec21|vec4&vec22|vec5&vec23|vec6&vec24;

cond(8,:) = vec1&vec10|vec2&vec11|vec3&vec12|vec4&vec13|vec5&vec14|vec6&vec15|vec7&vec16|vec8&vec17|vec9&vec18|vec10&vec19|vec11&vec20|vec12&vec21|vec13&vec22|vec14&vec23|vec15&vec24|vec1&vec14|vec2&vec15|vec3&vec16|vec4&vec17|vec5&vec18|vec6&vec19|vec7&vec20|vec8&vec21|vec9&vec22|vec10&vec23|vec11&vec24;

%========================================

cond(9,:) = vec1&vec11|vec2&vec12|vec3&vec13|vec4&vec14|vec5&vec15|vec6&vec16|vec7&vec17|vec8&vec18|vec9&vec19|vec10&vec20|vec11&vec21|vec12&vec22|vec13&vec23|vec14&vec24|vec1&vec13|vec2&vec14|vec3&vec15|vec4&vec16|vec5&vec17|vec6&vec18|vec7&vec19|vec8&vec20|vec9&vec21|vec10&vec22|vec11&vec23|vec12&vec24;

cond(10,:) = vec1&vec12|vec2&vec13|vec3&vec14|vec4&vec15|vec5&vec16|vec6&vec17|vec7&vec18|vec8&vec19|vec9&vec20|vec10&vec21|vec11&vec22|vec12&vec23|vec13&vec24|vec1&vec20|vec2&vec21|vec3&vec22|vec4&vec23|vec5&vec24;

cond(11,:) = vec1&vec15|vec2&vec16|vec3&vec17|vec4&vec18|vec5&vec19|vec6&vec20|vec7&vec21|vec8&vec22|vec9&vec23|vec10&vec24|vec1&vec23|vec2&vec24;

cond(12,:) = vec1&vec16|vec2&vec17|vec3&vec18|vec4&vec19|vec5&vec20|vec6&vec21|vec7&vec22|vec8&vec23|vec9&vec24;

%=========================================

cond(13,:) = vec1&vec17|vec2&vec18|vec3&vec19|vec4&vec20|vec5&vec21|vec6&vec22|vec7&vec23|vec8&vec24;

cond(14,:) = vec1&vec18|vec2&vec19|vec3&vec20|vec4&vec21|vec5&vec22|vec6&vec23|vec7&vec24|vec1&vec22|vec2&vec23|vec3&vec24;

cond(15,:) = vec1&vec21|vec2&vec22|vec3&vec23|vec4&vec24;

cond(16,:) = vec1&vec24|vec2&vec3|vec4&vec5|vec6&vec7|vec8&vec9|vec10&vec11|vec12&vec13|vec14&vec15|vec16&vec17|vec18&vec19|vec20&vec21|vec22&vec23;

for ii = 1:size(cond,1)
    variantPairs_boundary{ii} = variantBoundary(cond(ii,:));
end

[wanghaijian217]

If possible, the remaining V1-V7; can be given; V1-V8; V1-V9 (V19) and so on are 16 variant pairs, such as the variant pair boundary length fraction in the form of Figure 5 in this picture in DOI: 10.1016/j.actamat.2011. 12.018

[AzdiarGazder]

Please refer to ORTools_Example 10 for a detailed description and calculation of variant pairs.