void mitk::AnisotropicRegistrationCommon::PropagateMatrices( const MatrixList &src, MatrixList &dst, const Rotation &rotation ) { const vnl_matrix_fixed < double, 3, 3 > rotationT = rotation.GetTranspose(); #pragma omp parallel for for ( size_t i = 0; i < src.size(); ++i ) { dst[i] = rotation * src[i] * rotationT; } }