MTL 矩阵转置

来源：互联网发布：数据库有些什么算法编辑：程序博客网时间：2024/05/18 03:23

// RobinKin 整理

#include <mtl/matrix.h>
#include <mtl/mtl.h>
#include <mtl/utils.h>
using namespace mtl;

//元素是double, 形状是矩形(还有对称,上三角,下三角,带状矩阵等形状), 元素紧密(不稀疏,不压缩),

//列主(还可以行主)

typedef matrix< double, rectangle<>, dense<>, column_major>::type MATRIX;

//元素压缩存储
typedef matrix< double, rectangle<>, compressed<>, column_major>::type MATRIX1;

//矩阵元素是int 类型
typedef matrix< int,rectangle<> , dense<>, column_major>::type MATRIX2;

int
main (int , char *[])
{
MATRIX A(3, 5), B(5, 3);
MATRIX1 A1(3, 5), B1(5, 3);
MATRIX2 A2(3, 5), B2(5, 3);
for (int j = 0; j < 5; ++j)
    for (int i = 0; i < 3; ++i)
    {
      A(i,j) = i * 5 + j;
      A1(i,j) = i * 5 + j;
      A2(i,j) = i * 5 + j;
    }

//转置后的结果是一样的

transpose(A, B);
transpose(A1, B1);
transpose(A2, B2);
print_all_matrix(A);
print_all_matrix(B);

std::cout<<std::endl;
print_all_matrix(A1);
print_all_matrix(B1);

std::cout<<std::endl;
print_all_matrix(A2);
print_all_matrix(B2);
exit (0);
}