From affd383ef8aa0338479b7a9247253c6c810a9453 Mon Sep 17 00:00:00 2001 From: Philip Eisenlohr Date: Thu, 14 Apr 2011 19:09:44 +0000 Subject: [PATCH] simplified disorientation calculation. performed FEM check that "grainrotation" reflects the *active* rotation from the starting orientation to the current one with respect to the lab coord system (as expected). --- code/crystallite.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/code/crystallite.f90 b/code/crystallite.f90 index 85fba5532..1356c2459 100644 --- a/code/crystallite.f90 +++ b/code/crystallite.f90 @@ -2936,9 +2936,9 @@ logical error else orientation = math_RtoQuaternion(transpose(R)) endif - crystallite_rotation(1:4,g,i,e) = math_QuaternionDisorientation(math_qConj(orientation), & - math_qConj(crystallite_orientation0(1:4,g,i,e)), & - 0_pInt ) ! we don't want symmetry here + crystallite_rotation(1:4,g,i,e) = math_QuaternionDisorientation(crystallite_orientation0(1:4,g,i,e), & ! active rotation from ori0 + orientation, & ! to current orientation + 0_pInt ) ! we don't want symmetry here crystallite_orientation(1:4,g,i,e) = orientation enddo enddo