MTL 带状矩阵乘法

原创 2005年03月03日 15:18:00

//MTL 带状矩阵乘法

//整理 by RobinKin
//

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

/*

  Sample Output

  Array A in packed form
  [
  [1,1,1,],
  [2,2,2,],
  [3,3,3,],
  [4,4,],
  [5,],
  ]
  x
  [1,2,3,4,5,]
  y
  [1,1,1,1,1,]
  Ax + y
  [15,39,75,75,53,]

  */

using namespace mtl;

typedef matrix< double, banded<>, banded<>, row_major>::type Matrix;
typedef dense1D<double> Vector;

int
main()
{
  const int M = 5;
  const int N = 5;
 
  Matrix A(M, N, 0, 2);

  Vector x(N);
  Vector y(M);

  //           1  1  1              1        1
  //
  //              2  2  2           2        1
  //
  //       A =       3  3  3    x = 3    y = 1
  //
  //                    4  4        4        1
  //
  //                       5        5        1

  int row = 1;
//赋值
  Matrix::iterator ri = A.begin();
  while (ri != A.end()) {
    Matrix::Row::iterator i = (*ri).begin();
    while (i != (*ri).end()) {
      *i = row;
      ++i;
    }
    ++row;
    ++ri;
  }

  for (int i = 0; i < N; ++i)
    x[i] = i + 1;

  mtl::set_value(y, 1);


  std::cout << "Array A in packed form" << std::endl;
  print_row(A);

  std::cout << "x" << std::endl;
  print_vector(x);
  std::cout << "y" << std::endl;
  print_vector(y);

//y=2*x*A+3*y
  mult(A, scaled(x, 2), scaled(y, 3), y);

  std::cout << "Ax + y" << std::endl;
  print_vector(y);

  return 0;
}

基于c语言实现的矩阵乘法

  • 2016年03月30日 10:21
  • 2.43MB
  • 下载

矩阵乘法.c

  • 2015年01月13日 15:05
  • 1KB
  • 下载

列向量和行向量看待矩阵乘法

声明: 仅个人小记 前言: 主要是引入一个新的看待矩阵乘法的角度觉得这个挺重要的,故做记录列向量角度,矩阵左乘AB = C 结合上图,我们可以知道,结果矩阵C中的第 j 列完全可以表示为矩阵A中列...

n维矩阵乘法

  • 2012年12月11日 19:02
  • 220KB
  • 下载

基于C语言的矩阵乘法

  • 2015年07月27日 14:08
  • 2KB
  • 下载

POJ 2778 AC自动机+矩阵乘法

这题要写的东西很多,以至于我不知道从哪里开始写了= =。 总之要正确理解AC自动机中的状态转移。时间比较晚了,早点回去...  我要注意的地方就是正确理解fail指针的作用。在失配时,fail指针...

linux多线程实现矩阵乘法

  • 2011年11月08日 17:08
  • 2KB
  • 下载

《多核程序设计》学习笔记:矩阵乘法并行化

矩阵乘法的并行化与实现

矩阵乘法运算

  • 2012年04月10日 10:10
  • 498KB
  • 下载
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:MTL 带状矩阵乘法
举报原因:
原因补充:

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