基于模板C++创建转置矩阵函数

Create transpose matrix function based on template C++

本文关键字:函数 转置 创建 于模板 C++      更新时间:2023-10-16

我正在尝试使用这个惊人的指南创建我的第一个模板,但我遇到了我不知道如何处理多个嵌套参数的麻烦:

我创建的原始函数实现如下:

QVector< QVector<double> > transpose(QVector< QVector<double> > &matrix)
{
    int pre_numcols = matrix.size();
    int pre_numrows = matrix[0].size();
    QVector< QVector<double> > transposed(pre_numrows);
    QVector<double> newcols(pre_numcols);
    qFill(newcols.begin(), newcols.end(), 0.0);
    qFill(transposed.begin(), transposed.end(), newcols);
    qDebug()<<transposed.size();
    qDebug()<<transposed[0].size();
    for (int i = 0; i < pre_numcols; ++i)
    {
        for (int j = 0; j < pre_numrows; ++j)
        {
            transposed[j][i] = matrix[i][j];
        }
    }
    return transposed;
}

考虑到我只是重新洗牌,许多类型都是可能的。

QVector可以被std::vector取代,double甚至可以被string取代。

我走到了这一步:

template<class T>
T transpose(T &matrix)
{
    int pre_numcols = matrix.size();
    int pre_numrows = matrix[0].size();
    T transposed(pre_numrows);
    QVector<double> newcols(pre_numcols);  // How to do this one?

如果我T读作"装有某种小类型的容器容器",鉴于它是一个子集,我该如何声明 newcols 变量?

注意:稍后我将编辑qFill()部分以适应其他情况,例如std::

您可以将模板用作模板参数:

template <
    template <typename> class Container,
    typename ValueType
>
void foo(Container<ValueType> const & c);

因此,您的模板将变为:

template <
    template <typename> class Container,
    typename ValueType
>
Container< Container<ValueType> > transpose(Container< Container<ValueType> > & matrix)
{
    int pre_numcols = matrix.size();
    int pre_numrows = matrix[0].size();
    Container< Container<ValueType> > transposed(pre_numrows);
    Container<ValueType> newcols(pre_numcols);
    // ...
    return transposed;
}

如果它像那样工作,那就太好了。但是,与现实中一样,新的问题出现了! std::vector不是template <typename T> class vector而是

template <
    typename T,
    typename Allocator = allocator<T>
>
class vector

所以我们必须将我们的函数更改为

template <
    template <typename, typename> class C,
    typename T,
    template <typename> class A = std::allocator,
    typename InnerType = C< T, A<T> >,
    typename OuterType = C< InnerType, A<InnerType> >
>
OuterType transpose(OuterType & matrix)
{
    int pre_numcols = matrix.size();
    int pre_numrows = matrix[0].size();
    OuterType transposed(pre_numrows);
    InnerType newcols(pre_numcols);
    // ...
    return transposed;
}

这要不那么整洁,我不确定Qt容器是否与此兼容。

既然您说 C++11 可以为您所用,那么您可以使用第一个函数(使用 template <typename> class Container )并使用标准容器的模板别名:

template <typename T> using vector = std::vector<T, std::allocator<T>>;
vector<vector<int>> v;
auto vprime = transpose<vector>(v);

这是使用辅助模板元函数的另一种可能的解决方案,我认为它比以前的解决方案更干净:

namespace
{
    template <typename T> struct get_inner_i {};
    template <template <typename> class T, typename Inner> 
    struct get_inner_i<T<Inner>> { typedef Inner type; };
    template <
        template <typename, typename> class T, 
        typename Inner,
        template <typename> class Allocator
    > struct get_inner_i<T<Inner, Allocator<Inner>>> { typedef Inner type; };
    template <typename T> using get_inner = typename get_inner_i<T>::type;
}
template <typename MatrixType>
MatrixType transpose(MatrixType const & matrix)
{
    auto const nrows = matrix.size();
    auto const ncols = nrows > 0 ? matrix[0].size() : 0;
    MatrixType transposed(ncols, get_inner<MatrixType>(nrows));
    for(auto k = 0; k < nrows; ++k)
        for(auto j = 0; j < ncols; ++j)
            transposed[j][k] = matrix[k][j];
    return transposed;
}

为什么这个更冗长的解决方案更优越?这不是关于我们必须编写多少代码,而是我们的函数对用户来说有多容易和直观。在此版本中,我们再次只有T作为模板参数,因此不需要为函数显式指定模板参数。

此外,我们的元函数可以推断它是否是具有自定义分配器参数的容器,用户不必执行上面提到的using模板别名技巧。


这是另一种解决方案,再次优于以前的解决方案,它使用 std::begin 来确定容器的值类型,而不考虑其模板参数,或者它是否是模板,只要它提供beginend迭代器,或者它是一个 C 样式数组:

namespace
{
    template <typename Container>
    struct value_type_i
    {
        typedef typename std::decay<
            decltype(*std::begin(std::declval<
                Container const &
            >()))
        >::type type;
    };
    template <typename Container>
    using value_type = typename value_type_i<Container>::type;
}
template <typename MatrixType>
MatrixType transpose(MatrixType const & matrix)
{
    auto const nrows = matrix.size();
    auto const ncols = nrows > 0 ? matrix[0].size() : 0;
    MatrixType transposed(ncols, value_type<MatrixType>(nrows));
    for(auto k = 0; k < nrows; ++k)
        for(auto j = 0; j < ncols; ++j)
            transposed[j][k] = matrix[k][j];
    return transposed;
}

您可以在此处看到它的工作原理:http://ideone.com/cCAyFD