matrix_functions.cpp

来自「矩阵运算源码最新版本」· C++ 代码 · 共 24 行

#include <iostream>
#include <boost/numeric/mtl/mtl.hpp>

int main(int argc, char* argv[])
{
    using namespace mtl;
    
    const unsigned n= 10;
    compressed2D<double>                           A(n, n);
    dense2D<float, matrix::parameters<col_major> > B(n, n);
    morton_dense<double, 0x55555555>               C(n, n);
    morton_dense<double, 0x555555f0>               D(n, n);

    matrix::hessian_setup(B, 1.0);
    matrix::hessian_setup(C, 2.0);
    matrix::hessian_setup(D, 3.0);
    
    std::cout << "one_norm(B) is " << one_norm(B)<< "\n";
    std::cout << "infinity_norm(B) is " << infinity_norm(B)<< "\n";
    std::cout << "frobenius_norm(B) is " << frobenius_norm(B)<< "\n";
    
    return 0;
}