Twins: Fixed orientation prediction

Signed-off-by: Nicola Vigano <nicola.vigano@esrf.fr>

zUtil_Twins/gtTwinOrientations.m

-function [g_twins] = gtTwinOrientations(R_vector, twin_angle, twin_axis, spacegroup, lp, plot_variants)
+function g_twins = gtTwinOrientations(R_vector, twin_angle, twin_axis, spacegroup, lp, plot_variants)
 % GTTWINORIENTATIONS
-%     [g_twins] = gtTwinOrientations(R_vector, twin_angle, twin_axis, spacegroup, lp, [plot_variants])
+%     g_twins = gtTwinOrientations(R_vector, twin_angle, twin_axis, spacegroup, lp, [plot_variants])
 %     ------------------------------------------------------------------------------------------------
 %     By passing in the parent R_vector, return the g (sample to crystal transformation) plus
 %     all the twin variants
@@ -57,7 +57,7 @@ function [g_twins] = gtTwinOrientations(R_vector, twin_angle, twin_axis, spacegr
     for ii = num_axes:-1:1
         rotcomp = gtMathsRotationMatrixComp(twin_axes(ii, :)', 'col');
         cry2twins(:, :, ii) = gtMathsRotationTensor(twin_angle, rotcomp);
-        gs(:, :, ii) = sam2cry * cry2twins(:, :, ii);
+        gs(:, :, ii) = cry2twins(:, :, ii) * sam2cry;
     end
 
     orimats = gtMathsAngleAxis2OriMat(twin_angle(1, ones(num_axes, 1)), twin_axes');