如何有效地对非方阵进行转置？

问题描述

我做了一个矩阵类，我想实现一个换位方法：

template<typename T>
void Matrix<T>::Transpose()
{
    // for square matrices
    if(this->Width() == this->Height())
    {
        for(std::size_t y = 1; y < this->Height(); ++y)
        {
            for(std::size_t x = 0; x < y; ++x)
            {
                // the function operator is used to access the entries of the matrix
                std::swap((*this)(x, y), (*this)(y, x));
            }
        }
    }
    else
    {
        // TODO
    }
}

问题是如何实现非方阵的转置方法，而不分配一个全新的矩阵(该类用于大的稠密矩阵)，而不是原地分配。有什么办法吗？

推荐答案

转置矩阵的最有效方法是根本不转置矩阵。

通过另外存储行和列步长以及偏移量，以允许在同一数据缓冲区上定义子矩阵、切片或任何东西的方式设计矩阵类可能是最有效的。然后访问元素，您使用这些数据来计算指数。要转置，您只需操作这些步长值。

您可以查看OpenCV的Matrix实现(仅用于功能实现，不用于类设计！)

这篇关于如何有效地对非方阵进行转置？的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持！