Plotting the twin relationships

[Elizaveta Kovaleva]

Good day all, 

I am trying to figure out the relationships between twins in monazite; 

I am following this script: https://mtex-toolbox.github.io/TwinningBoundaries.html

For some reason, MATLAB does not follow the command to restrict the grain boundary to a specific phase (Monazite-Monazite in my case):

gB_MonMon = gB('Monazite-(Ce)','Monazite-(Ce)')

Output: 

"Error using ischar

Not enough input arguments.

Error in grainBoundary/subsref (line 31)

    if ~ischar(gB.CSList{phId})"

See attached image. 

Any further operations do not go through either :(

I am trying to identify the angle and the rotation axis between the twins.

I would appreciate any help!

Thank you

Best wishes

Elizaveta

[Rüdiger Kilian]

Hi,
could you please just repeat your problem on the new forum/discussion list? https://github.com/mtex-toolbox/mtex/discussions
Cheers, Rüdiger

[Elizaveta Kovaleva]

Hi all, 

I have downloaded the newest version and all works there. Thanks! 

But I got quite lost here:

"We may determin the center of the cluster and check whether it is close to some special orientation relation ship

% determine the mean of the cluster 

mori_mean = mean(mori,'robust') 

  % determine the closest special orientation relation ship 

round2Miller(mori_mean)mori_mean = misorientation 

size: 1 x 1 

crystal symmetry : Magnesium (6/mmm, X||a*, Y||b, Z||c*) 

crystal symmetry : Magnesium (6/mmm, X||a*, Y||b, Z||c*) 

antipodal: true 

Bunge Euler angles in degree 

phi1 Phi phi2 Inv. 

90.5949 86.0962 269.995 0 

(10-1-2) || (-101-2) [1-210] || [1-210] error: 0.65°

Bases on the output above we may now define the special orientation relationship as

twinning = orientation.map(Miller(0,1,-1,-2,CS),Miller(0,-1,1,-2,CS),... Miller(2,-1,-1,0,CS),Miller(2,-1,-1,0,CS))  "

Where have these Miller indices in the last line came from? They are not the same as in the output? 

Thank you

Best wishes

Elizaveta

[Rüdiger]

Hi,

that would have been my suggestion. Could you please still ask in the
new forum? https://github.com/mtex-toolbox/mtex/discussions

Cheers,
Rüdiger