#include <vector>

#include <algorithm>

#include <numeric>

#include "cppunit/cppunit_proxy.h"

#if !defined (STLPORT) || defined(_STLP_USE_NAMESPACES)

using namespace std;

#endif

//

// TestCase class

//

class InnerprodTest : public CPPUNIT_NS::TestCase

{

  CPPUNIT_TEST_SUITE(InnerprodTest);

  CPPUNIT_TEST(inprod0);

  CPPUNIT_TEST(inprod1);

  CPPUNIT_TEST(inprod2);

  CPPUNIT_TEST_SUITE_END();

protected:

  void inprod0();

  void inprod1();

  void inprod2();

  static size_t add(size_t a_, size_t b_) {

    return a_ + b_;

  }

  static size_t mult(size_t a_, size_t b_) {

    return a_ * b_;

  }

};

CPPUNIT_TEST_SUITE_REGISTRATION(InnerprodTest);

//

// tests implementation

//

void InnerprodTest::inprod0()

{

  int vector1[5] = { 1, 2, 3, 4, 5 };

  int vector2[5] = { 1, 2, 3, 4, 5 };

  int result;

  result = inner_product(vector1, vector1 + 5, vector2, 0);

  CPPUNIT_ASSERT(result==55);

}

void InnerprodTest::inprod1()

{

  vector<size_t> v1(3);

  vector<size_t> v2(v1.size());

  for (size_t i = 0; i < v1.size(); ++i) {

    v1[i] = i + 1;

    v2[i] = v1.size() - i;

  }

  size_t result = inner_product(v1.begin(), v1.end(), v2.begin(), (size_t)0);

  CPPUNIT_ASSERT(result == 10);

}

void InnerprodTest::inprod2()

{

  vector<size_t> v1(3);

  vector<size_t> v2(v1.size());

  for(size_t i = 0; i < v1.size(); ++i) {

    v1[i] = i + 1;

    v2[i] = v1.size() - i;

  }

  size_t result=inner_product(v1.begin(), v1.end(), v2.begin(), (size_t)1, mult, add);

  CPPUNIT_ASSERT(result == 64);

}