diff --git a/Framework/Geometry/inc/MantidGeometry/Crystal/SymmetryOperation.h b/Framework/Geometry/inc/MantidGeometry/Crystal/SymmetryOperation.h index 9842d8ab46f9eeb614b9e77275e9e84a798db986..2c3ec48e71d9e677059bbdc0831d47281b10e3da 100644 --- a/Framework/Geometry/inc/MantidGeometry/Crystal/SymmetryOperation.h +++ b/Framework/Geometry/inc/MantidGeometry/Crystal/SymmetryOperation.h @@ -161,7 +161,7 @@ protected: const V3R &vector) const; size_t m_order; - Kernel::IntMatrix m_inverseMatrix; + Kernel::IntMatrix m_transposedInverseMatrix; V3R m_reducedVector; std::string m_identifier; diff --git a/Framework/Geometry/src/Crystal/SymmetryOperation.cpp b/Framework/Geometry/src/Crystal/SymmetryOperation.cpp index 146eeec8666bf7560693fd8c956d81aad634c4c0..5dbf7e3c942b8a3a48afff78a73675bcd49d1a66 100644 --- a/Framework/Geometry/src/Crystal/SymmetryOperation.cpp +++ b/Framework/Geometry/src/Crystal/SymmetryOperation.cpp @@ -9,7 +9,7 @@ namespace Geometry { /// Default constructor, results in identity. SymmetryOperation::SymmetryOperation() - : m_order(1), m_inverseMatrix(Kernel::IntMatrix(3, 3, true)), + : m_order(1), m_transposedInverseMatrix(Kernel::IntMatrix(3, 3, true)), m_reducedVector(), m_identifier(), m_matrixVectorPair() { m_identifier = SymmetryOperationSymbolParser::getNormalizedIdentifier( m_matrixVectorPair.getMatrix(), m_matrixVectorPair.getVector()); @@ -52,9 +52,9 @@ void SymmetryOperation::init(const Kernel::IntMatrix &matrix, MatrixVectorPair<int, V3R>(matrix, getWrappedVector(vector)); // Inverse matrix for HKL operations. - m_inverseMatrix = Kernel::IntMatrix(matrix); - m_inverseMatrix.Invert(); - m_inverseMatrix = m_inverseMatrix.Transpose(); + m_transposedInverseMatrix = Kernel::IntMatrix(matrix); + m_transposedInverseMatrix.Invert(); + m_transposedInverseMatrix = m_transposedInverseMatrix.Transpose(); m_order = getOrderFromMatrix(m_matrixVectorPair.getMatrix()); m_identifier = SymmetryOperationSymbolParser::getNormalizedIdentifier( @@ -136,7 +136,7 @@ bool SymmetryOperation::hasTranslation() const { * @return :: Transformed index triplet. */ Kernel::V3D SymmetryOperation::transformHKL(const Kernel::V3D &hkl) const { - return m_inverseMatrix * hkl; + return m_transposedInverseMatrix * hkl; } /**