matrix_functions2.cpp

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

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

int main(int argc, char* argv[])
{
    using namespace mtl;
    
    typedef std::complex<double>      cdouble;
    const unsigned                    xd= 2, yd= 5, n= xd * yd;
    compressed2D<cdouble>             A(n, n);
    matrix::laplacian_setup(A, xd, yd); 

    // Fill imaginary part of the matrix
    A*= cdouble(1, -1);
    std::cout << "A is\n" << with_format(A, 7, 1) << "\n";

    std::cout << "trace(A) is " << trace(A) << "\n\n";

    std::cout << "conj(A) is\n" << with_format(conj(A), 7, 1) << "\n";

    std::cout << "trans(A) is\n" << with_format(trans(A), 7, 1) << "\n";

    return 0;
}