Bye Bye Moore

PoCソルジャーな零細事業主が作業メモを残すブログ

Boost::ublasで行列計算する その5:行と列

実際のところ

matrix_row(行)とmatrix_column(列)

#include <iostream>
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/matrix_proxy.hpp>
#include <boost/numeric/ublas/io.hpp>

using namespace std;
using namespace boost::numeric::ublas;

int main()
{
    // double型を要素とする行列
    matrix<double> v(3,3);

    // 各要素の参照と代入
    v(0,0) = 3.0;  v(0,1) = 0.0; v(0,2) = 4.0; 
    v(1,0) = 3.0;  v(1,1) = 0.0; v(1,2) = 4.0; 
    v(2,0) = 3.0;  v(2,1) = 0.0; v(2,2) = 4.0; 

    matrix_row<matrix<double>> mr(v,0);
    for (unsigned j = 0; j < mr.size(); j++) {
      std::cout << "mr(" 
	        << j 
		<< "): " 
		<< mr(j) 
		<< " + "
		<< j
                << " = " 
		<< mr(j) + j
		<< std::endl;
    }
}