vector_min_max_test.cpp

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

// Software License for MTL
// 
// Copyright (c) 2007 The Trustees of Indiana University. All rights reserved.
// Authors: Peter Gottschling and Andrew Lumsdaine
// 
// This file is part of the Matrix Template Library
// 
// See also license.mtl.txt in the distribution.

#include <iostream>
#include <cmath>
#include <boost/test/minimal.hpp>
#include <boost/numeric/mtl/concept/collection.hpp>
#include <boost/numeric/mtl/vector/dense_vector.hpp>
#include <boost/numeric/mtl/vector/parameter.hpp>
#include <boost/numeric/mtl/operation/operators.hpp>
#include <boost/numeric/mtl/operation/min.hpp>
#include <boost/numeric/mtl/operation/max.hpp>


using namespace mtl;
using namespace std;  
    

template <typename Vector>
void test(Vector& v, const char* name)
{
    typedef typename Collection<Vector>::value_type value_type;

    for (int i= 0; i < size(v); i++)
	v[i]= value_type(double(i+1) * pow(-1.0, i)); 

    std::cout << "\n" << name << "  --- v = " << v; std::cout.flush();

    std::cout << "min(v) = " << min(v) << "\n"; std::cout.flush();
    if (min(v) != -4.0) throw "min wrong";

    std::cout << "min<4>(v) = " << min<4>(v) << "\n"; std::cout.flush();
    if (min<4>(v) != -4.0) throw "min<4> wrong";

    std::cout << "max(v) = " << max(v) << "\n"; std::cout.flush();
    if (max(v) != 5.0) throw "max wrong";
}
 

int test_main(int argc, char* argv[])
{
    using mtl::vector::parameters;

    dense_vector<float>   u(5);
    dense_vector<double>  x(5);

    std::cout << "Testing vector operations\n";

    test(u, "test float");
    test(x, "test double");

    dense_vector<float, parameters<row_major> >   ur(5);
    test(ur, "test float in row vector");

    return 0;
}