MTL 矩阵转置

// 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);
}