/*
 * reftest.cpp
 *
 *  Created on: Apr 3, 2011
 *      Author: joe
 */

#include <iostream>
#include <boost/fusion/container.hpp>
#include <boost/fusion/view.hpp>
#include <boost/fusion/include/io.hpp>
#include <boost/fusion/algorithm.hpp>
#include <boost/fusion/include/make_vector.hpp>
#include <boost/fusion/include/next.hpp>
#include <boost/fusion/include/prior.hpp>
#include <boost/fusion/include/deref.hpp>
#include <boost/mpl/assert.hpp>
#include <boost/spirit/include/phoenix_function.hpp>
#include <boost/spirit/include/phoenix_core.hpp>
#include <boost/spirit/include/phoenix_operator.hpp>
#include <boost/spirit/include/phoenix_statement.hpp>
#include <boost/spirit/include/phoenix_scope.hpp>
#include <boost/spirit/include/phoenix_fusion.hpp>
#include <boost/mpl/int.hpp>

using namespace boost;
using namespace boost::fusion;
using boost::fusion::vector;

namespace functors {

// A functor returns the dot product of two sequences. Presumes that the return type is the
// type of the first element of the first sequence argument.
struct vector_dot_product_impl {
	template <typename V1, typename V2>
	struct result
	{
		// Presume that the return type is the type of the first element
		// of the first sequence argument.
		typedef typename boost::fusion::result_of::value_at<V1, boost::mpl::int_<0> >::type type;
	};

	template <typename V1, typename V2>
	typename result<V1, V2>::type operator()(V1 v1, V2 v2) const
	{
		using boost::phoenix::arg_names::_1;
		using boost::phoenix::arg_names::_2;
		using boost::fusion::fold;
		using boost::fusion::transform;
		return fold(transform(v1, v2, _1 * _2), 0.0, _1 + _2);
	}
};
boost::phoenix::function<vector_dot_product_impl> const vector_dot_product;

} // namespace functors


int main() {
	typedef vector<int, int, int> int3;
	typedef vector<int, int> int2;

	int3 i3a = make_vector(1,2,3);
	int3 i3b = make_vector(4,5,6);

	using boost::phoenix::arg_names::_1;
	using boost::phoenix::arg_names::_2;
	std::cout << "DOT: " << functors::vector_dot_product(_1, _2)(i3a, i3b) << std::endl;
}