C++ MTL 矩阵 scale 示例 （整理 by RobinKin from DevonIT.inc)

// -*- c++ -*-
//
// $COPYRIGHT$
//
//===========================================================================

#include <iostream>
using namespace std;

#include <mtl/matrix.h>
#include <mtl/mtl.h>
#include <mtl/utils.h>

/*

  Sample Output

  Before scaling row 2:
  3x3
  [
  [5,5.5,6],
  [2.5,3,3.5],
  [1,1.5,2]
  ]
  After scaling row 2:
  3x3
  [
  [5,5.5,6],
  [5,6,7],
  [1,1.5,2]
  ]
  
  */

using namespace mtl;

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

int
main()
{
  const Matrix::size_type N = 3;

  Matrix A(N, N);

  A(0,0) = 5;   A(0,1) = 4;   A(0,2) = 3;
  A(1,0) = 2.5; A(1,1) = 3;   A(1,2) = 3.5;
  A(2,0) = 1;   A(2,1) = 1.5; A(2,2) = 2;

  cout << endl;
  cout << "Before scaling row 2:" << endl;
  print_all_matrix(A);

  double scal = A(0,0) / A(1,0);

  copy(scaled(A[1], scal), A[1]);

  cout << "After scaling row 2:" << endl;
  print_all_matrix(A);

  return 0;
}