vector_reduction.cpp

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

// File: vector_reduction.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(100);

    for (int i= 0; i < size(v); i++)
	v[i]= cdouble(i+1, 100-i);

    std::cout << "sum(v) is " << sum(v)<< "\n";
    
    std::cout << "product(v) is " <<  product(v)<< "\n";
    
    std::cout << "sum<6>(v) is " <<  sum<6>(v)<< "\n";

    return 0;
}