Boost::ublasで行列計算する その３：ゼロ行列・単位行列

実際のところ

ゼロ行列

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

int main () {
    using namespace boost::numeric::ublas;
    zero_matrix<double> m (3, 3);
    std::cout << m << std::endl;
}

単位行列

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

int main () {
    using namespace boost::numeric::ublas;
    identity_matrix<double> m (3);
    std::cout << m << std::endl;
}

参考もと

線形代数 - boostjp
Matrix - 1.68.0
Matrix - 1.68.0