vector_norm.cpp

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

// File: vector_norm.cpp

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

int main(int argc, char* argv[])
{
    using namespace mtl;

    typedef std::complex<double>  cdouble;
    dense_vector<cdouble>         v(10000);

    // Initialize vector
    for (int i= 0; i < size(v); i++)
	v[i]= cdouble(i+1, 10000-i);

    std::cout << "one_norm(v) is " << one_norm(v)<< "\n";
    
    std::cout << "two_norm(v) is " << two_norm(v)<< "\n";
    
    std::cout << "infinity_norm(v) is " << infinity_norm(v)<< "\n";
    
    // Unroll computation of two-norm to 6 independent statements
    std::cout << "two_norm<6>(v) is " << two_norm<6>(v)<< "\n";

    return 0;
}