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

MTL——C++实现的矩阵模板类

  • 2008年10月14日 09:45
  • 42KB
  • 下载

带状矩阵的压缩存储

对于n*n的方阵,若它的全部非零元素落在一个以主对角线为中心的带状区域中,这个带状区域包含主对角线 ,以及主对角线下面及上面各b条对角线上的元素,那么称该方阵为半带宽为b的带状矩阵。 带状矩阵的特点...
  • define_danmu_primer
  • define_danmu_primer
  • 2016-01-11 15:24:56
  • 1800

矩阵的压缩存储--特殊矩阵--带状矩阵

矩阵的压缩存储–特殊矩阵–带状矩阵带状矩阵指矩阵中所有的非零元素都集中在以对角线为中心的带状区域中本篇以最常见的三对角带状矩阵为例,示例代码:#include #include #define N...
  • y_universe
  • y_universe
  • 2017-12-11 18:12:58
  • 163

MTL 解下三角带状矩阵线形方程

//整理by RobinKin#include "mtl/matrix.h"#include "mtl/mtl.h"#include "mtl/utils.h"#include "mtl/linalg...
  • RobinKin
  • RobinKin
  • 2005-03-09 10:11:00
  • 1699

硬盘坏道镜像工具MTL

  • 2016年10月29日 11:02
  • 1.68MB
  • 下载

obj + mtl 格式说明

 OBJ(或 .OBJ)是一种开放的几何定义文件格式,最初由Wavefront Technologies公司开发,用以描述其Advanced Visualizer动画包。该格式已被其他3D图形应...
  • weiming8517
  • weiming8517
  • 2015-07-24 10:30:24
  • 17066

使用three.js加载obj+mtl文件

因在项目中运用到WEBGL的技术。所以研究了一下下。three.js库看来十分强大,此处为加载带材质的obj文件 1.下载相关js    https://threejs.org/     thee....
  • qq_37642205
  • qq_37642205
  • 2018-01-05 13:53:41
  • 1104

Obj与mtl 3D文件格式详解

OBJ文件       obj文件格式是Wavefront公司为它的一套基于工作站的3D建模和动画软件"AdvancedVisualizer"开发的一种文件格式。 OBJ文件是一种标准的3D模型文件...
  • z444_579
  • z444_579
  • 2016-08-17 15:24:47
  • 3745

.obj文件格式与.mtl文件格式

最近在学习obj文件格式,上网查了些资料,很难找到比较全面的文章,尤其是对.mtl文件的说明甚少。今天把最近搜索的资料整合了一下。这里的obj文件格式指的是Wavefront公司为它的一套基于工作站的...
  • u013467442
  • u013467442
  • 2015-07-07 17:55:23
  • 5349

.obj文件格式与.mtl文件格式(转)

.obj文件格式与.mtl文件格式
  • sinat_30751333
  • sinat_30751333
  • 2016-03-23 11:34:18
  • 1505
收藏助手
不良信息举报
您举报文章:MTL 带状矩阵乘法
举报原因:
原因补充:

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