C++ MTL 矩阵 行交换 示例 （整理 by RobinKin from DevonIT)

#include "mtl/matrix.h"
#include "mtl/mtl.h"
#include "mtl/utils.h"

/*
  Swaps the first row with the row that has the largest element in
  column 1. This is the kind of operation one would do inside a
  Gaussian elimintation algorithm.

  Sample Output

  3x3
  [
  [1,1.5,4.5],
  [3,2.5,9.5],
  [2,3.5,5.5]
  ]
  3x3
  [
  [3,2.5,9.5],
  [1,1.5,4.5],
  [2,3.5,5.5]
  ]

*/

using namespace mtl;
using namespace std;

int
main()
{
  //begin
  typedef matrix< double,
                  rectangle<>,
                  dense<external>,
                  column_major>::type Matrix;
  const Matrix::size_type N = 3;
  Matrix::size_type large;
  double dA[] = { 1, 3, 2, 1.5, 2.5