MTL 矩阵划分示例

原创 2005年03月03日 14:10:00

//整理 by RobinKin  

#include "mtl/matrix.h"
#include "mtl/linalg_vec.h"
#include "mtl/utils.h"
#include <list>

/*
  Sample output:

 Partitioning Vectors
 //向量划分
[0,1,2,3,4,5,6,7,8,9,]
split in half
[0,1,2,3,4,]
[5,6,7,8,9,]
split into thirds
[0,1,2,]
[3,4,5,6,]
[7,8,9,]
 
Partitioning Matrices
 //矩阵划分
4x4
[
[0,1,2,3],
[4,5,6,7],
[8,9,10,11],
[12,13,14,15]
]

划分成 2x2的 子矩阵(行2个 ,列两个)
split into four submatrices
Top mat: 2x2
2x2
[
[0,1],
[4,5]
]
2x2
[
[2,3],
[6,7]
]
 
2x2
[
[8,9],
[12,13]
]
2x2
[
[10,11],
[14,15]
]
 划分成 2x2的 子矩阵
split into six submatrics
Top mat: 3x2 (行3 ,列2)
1x2
[
[0,1]
]
1x2
[
[2,3]
]
 
2x2
[
[4,5],
[8,9]
]
2x2
[
[6,7],
[10,11]
]
 
1x2
[
[12,13]
]
1x2
[
[14,15]
]

 
  */

int
main()
{
  using namespace mtl;

  /* Partitioning Vectors */
  std::cout << "Partitioning Vectors" << std::endl << std::endl;
  const int N = 10;
  int dx[N];
  int i, j;

  for (i = 0; i < N; ++i)
    dx[i] = i;

  typedef external_vec<int> Vec;
  typedef dense1D<Vec::subrange_type> PVec;

  Vec x(dx, N);
  print_vector(x);

  std::cout << "split in half" << std::endl;
  PVec xp1(2); //分成两个
  subdivide(5, x, xp1);//每个子向量5个元素
  print_partitioned_vector(xp1);

  std::cout << "split into thirds" << std::endl;
  PVec xp2(3);//分3 份
  int splits[] = { 3,7  };
  partition(array_to_vec(splits), x, xp2);
  print_partitioned_vector(xp2);

  /* Partitioning Matrices */
  std::cout << std::endl;
  std::cout << "Partitioning Matrices" << std::endl << std::endl;

  typedef matrix<int, rectangle<>, dense<>, row_major>::type Matrix;

  typedef matrix<Matrix::submatrix_type,
                 rectangle<>, dense<>, row_major>::type PMatrix;

  const int M = 4;
  Matrix A(M, M);
  for (i = 0; i < M; ++i)
    for (j = 0; j < M; ++j)
      A(i,j) = i * M + j;

  print_all_matrix(A);

  std::cout << "split into four submatrices" << std::endl;
 
  PMatrix Ap1(2,2);//2x2个子矩阵
  subdivide(2, 2, A, Ap1);
  print_partitioned_matrix(Ap1);

  std::cout << "split into six submatrics" << std::endl;

 
  int row_splits[] = { 1, 3 };//按照1,3 划分行
  int col_splits[] = { 2 };
  PMatrix Ap2(3,2);//3x2的子矩阵
  partition(array_to_vec(row_splits),
     array_to_vec(col_splits),
     A, Ap2);
  print_partitioned_matrix(Ap2);

  return 0;
}

版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

MTL矩阵模板库

  • 2014-07-28 17:27
  • 1.36MB
  • 下载

mtl_矩阵模板库

  • 2016-10-09 19:38
  • 157KB
  • 下载

OpenCV—基本矩阵操作与示例

图像处理中对矩阵的操作非常重要,本文总结了使用OpenCV进行矩阵的创建、初始化以及基本矩阵操作,给出了示例代码。
  • iracer
  • iracer
  • 2016-05-02 16:30
  • 23179

设计模式-矩阵示例

  • 2014-08-07 23:14
  • 361KB
  • 下载

动态规划(dynamic programming)及示例(矩阵连乘、最长公共子序列、三角剖分)

动态规划是一个比较重要的算法思想,只要学了算法相关课程,一定会有这一章节。我会介绍动态规划的主要思想,也会在后面用Java实现几个经典算法(矩阵连乘、最长公共子序列、多边形三角剖分)作为示例加深理解。

c++计算矩阵乘积示例

  • 2014-01-09 18:08
  • 7.03MB
  • 下载

CUDA矩阵乘法——VS2010中使用CUDA示例

文章转自:http://blog.csdn.net/augusdi/article/details/12617177 各工具或库的版本: IDE:VS2008 VS2010 (使用MSVC编译器...

java矩阵运算包ujmp中的一些小示例和注意事项

本人最近在用ujmp包写一些程序,ujmp包是针对于超大数据量计算的矩阵的运算包,并且有图形显示的功能且支持多种文件格式的读取和输出,还支持连接数据库,matlab数据类型和weka数据类型,总体来说...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)