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