正在生成编译时数组结构(c++17）

我有以下C++17代码，用于生成数组的编译时元组，其中零数组只是为了示例，在我的实现中，它们将是满的(使用-std=C++1z-fconcepts编译)。

#include <array>
#include <tuple>
#include <cmath>
template <std::size_t nrA, std::size_t ncA, std::size_t nrB, std::size_t ncB,
typename number=double>
constexpr auto operator *(const std::array<std::array<number,ncA>,nrA> & A,
const std::array<std::array<number,ncB>,nrB> & B)
{
std::array<std::array<number,ncB>,nrA> res{} ;
for (auto k=0u;k<ncB;++k)
for (auto i=0u;i<nrA;++i)
for (auto j=0u;j<nrB;++j)
res[i][k] += A[i][j]*B[j][k];
return res ;
}
constexpr auto logg2(const auto N)
{
auto res = 0;
auto n = N;
while (n != 0) 
{
n /= 2;
++res;
}
return res;
}
template <std::size_t N,typename number=double>
constexpr auto create_R()
{
return std::array<std::array<double,2*N>,N>{};
}
template <std::size_t N,typename number=double>
constexpr auto create_RT()
{
return std::array<std::array<double,N>,2*N>{};
}
template <std::size_t N,std::size_t ...Is>
constexpr auto make_impl(const std::index_sequence<Is...>)
{
return std::make_tuple(std::make_tuple(create_R<(N >> Is)>()...),
std::make_tuple(create_RT<(N >> Is)>()...));
}
template <std::size_t N,typename number=double>
constexpr auto make()
{
return make_impl<N/2>(std::make_index_sequence<logg2(N/2) - 1>());
}
int main(int argc, char *argv[])
{
const auto n = 4u;
const auto A = std::array<std::array<double,2*n>,2*n>{};
const auto [R,RT] = make<2*n>();
}

我想将make<>()修改为make<>(A)，并返回一个结构化绑定[R,RT,As]，其中As是一个元组，其中包含以下数组

A,
std::get<0>(R)*A*std::get<0>(RT),
std::get<1>(R)*std::get<0>(R)*A*std::get<0>(RT)*std::get<1>(RT)
...

我已经尝试了一段时间，但没有找到解决方案。

有什么想法吗？

编辑1

根据@MaxLanghof的要求，以下打印矩阵：

template <std::size_t nr, std::size_t nc, typename number=double>
constexpr auto print(const std::array<std::array<number,nc>,nr> & A)
{
for (auto i=0u;i<nr;++i)
{
for (auto j=0u;j<nc;++j)
std::cout << std::right << std::setw(12) << A[i][j];
std::cout << std::endl ;
}
std::cout << std::endl ;
}

并将以下行添加到main()

print(A);
print(std::get<0>(R)*A*std::get<0>(RT));
print(std::get<1>(R)*std::get<0>(R)*A*std::get<0>(RT)*std::get<1>(RT));

获得以下输出

0           0           0           0           0           0           0           0
0           0           0           0           0           0           0           0
0           0           0           0           0           0           0           0
0           0           0           0           0           0           0           0
0           0           0           0           0           0           0           0
0           0           0           0           0           0           0           0
0           0           0           0           0           0           0           0
0           0           0           0           0           0           0           0
0           0           0           0
0           0           0           0
0           0           0           0
0           0           0           0
0           0
0           0