perftest.cpp

Go to the documentation of this file.

#include <vector>
#include <algorithm>
#include <numeric>
#include <boost/bind.hpp>
#include <cmath>
#include <functional>
#include <boost/numeric/ublas/vector.hpp>
#include <boost/numeric/ublas/vector_expression.hpp>

namespace ublas = boost::numeric::ublas;
using namespace std;
using namespace gplib;

const int size = 1000000;
const int iterations = 1;

void f(vector<double> &v1, vector<double> &v2, vector<double> &v3, vector<
    double> &v4)
  {
    transform(v1.begin(), v1.end(), v2.begin(), v3.begin(), boost::bind(
        std::multiplies<double>(), boost::bind<double>(std::minus<double>(),
            _1, _2), boost::bind<double>(std::minus<double>(), _1, _2)));
    copy(v2.begin(), v2.end(), v4.begin());
    double value = accumulate(v3.begin(), v3.end(), 0.0);
  }

void g(vector<double> &v1, vector<double> &v2, vector<double> &v3, vector<
    double> &v4)
  {
    double value = 0;
    for (size_t i = 0; i < v1.size(); ++i)
      {
        v3.at(i) = pow(v1.at(i) - v2.at(i), 2);
        v4.at(i) = v2.at(i);
        value += v3.at(i);
        //v3.at(i) = v1.at(i) * v2.at(i);
      }
  }

void h(vector<double> &v1, vector<double> &v2, vector<double> &v3, vector<
    double> &v4)
  {
    typedef vector<double>::iterator tit;
    double value = 0;
    for (tit it1 = v1.begin(), it2 = v2.begin(), it3 = v3.begin(), it4 =
        v4.begin(); it1 != v1.end(); ++it1, ++it2, ++it3, ++it4)
      {
        *it3 = (*it1 - *it2) * (*it1 - *it2);
        *it4 = *it2;
        value += *it3;
      }
  }

void ublasf(ublas::vector<double> &v1, ublas::vector<double> &v2,
    ublas::vector<double> &v3, ublas::vector<double> &v4)
  {
    transform(v1.begin(), v1.end(), v2.begin(), v3.begin(), boost::bind(
        std::multiplies<double>(), boost::bind<double>(std::minus<double>(),
            _1, _2), boost::bind<double>(std::minus<double>(), _1, _2)));
    v2 = v4;
    double value = ublas::sum(v3);
  }

void ublasg(ublas::vector<double> &v1, ublas::vector<double> &v2,
    ublas::vector<double> &v3, ublas::vector<double> &v4)
  {
    double value = 0;
    for (size_t i = 0; i < v1.size(); ++i)
      {
        v3(i) = pow(v1(i) - v2(i), 2);
        v4(i) = v2(i);
        value += v3(i);
        //v3.at(i) = v1.at(i) * v2.at(i);
      }
  }

void ublash(ublas::vector<double> &v1, ublas::vector<double> &v2,
    ublas::vector<double> &v3, ublas::vector<double> &v4)
  {
    typedef ublas::vector<double>::iterator tit;
    double value = 0;
    for (tit it1 = v1.begin(), it2 = v2.begin(), it3 = v3.begin(), it4 =
        v4.begin(); it1 != v1.end(); ++it1, ++it2, ++it3, ++it4)
      {
        *it3 = (*it1 - *it2) * (*it1 - *it2);
        *it4 = *it2;
        value += *it3;
      }
  }
int main()
  {
    vector<double> V1(size, 2.), V2(size, 1.5), V3(size, 1.1), V4(size, 1.3);
    f(V1, V2, V3, V4);
    g(V1, V2, V3, V4);
    h(V1, V2, V3, V4);

    ublas::vector<double> U1(size), U2(size), U3(size), U4(size);

    ublasf(U1, U2, U3, U4);
    ublasg(U1, U2, U3, U4);
    ublash(U1, U2, U3, U4);
  }