Axis of rotation during twinning tutorial

[grandr...@gmail.com]

I am following along the twinning analysis tutorial in http://mtex-toolbox.github.io/files/doc/TwinningAnalysis.html  When I get to the step

isTwinning = angle(gB_MgMg.misorientation,twinning) < 5*degree

This seems to be checking if the misorientation is within 5 degrees of the twinning rotation.  
The problem is, this is insufficient for a twinning definition as the axis of rotation is also critical.
I understand that the relationship for twinning is defined earlier, my question is at what line of the code is this rotation axis checked?
The angle() command checks the angular relationship between two points, and does not appear to return an axis.

Regards, Jessica

[Ralf Hielscher]

Hi Jessica,

the command 

angle(mori1, mori2)

computes the angular deviation between two misorientations. This also takes a axis  into account, i.e., if both mori1 and mori2 have the misorientation angle 180 degree but their misorientation axis differ by 10 degree then the angle between mori1 and mori2 is 10 degree as well. 

You may investigate this by playing around with MTEX.

Ralf.