Subject: [Boost-commit] svn:boost r60777 - sandbox/numeric_bindings/boost/numeric/bindings/boost
From: rutger_at_[hidden]
Date: 2010-03-22 16:26:40
Author: rutger
Date: 2010-03-22 16:26:39 EDT (Mon, 22 Mar 2010)
New Revision: 60777
URL: http://svn.boost.org/trac/boost/changeset/60777
Log:
Added initial version of bindings support for Boost.MultiArray
Added:
sandbox/numeric_bindings/boost/numeric/bindings/boost/multi_array.hpp (contents, props changed)
Added: sandbox/numeric_bindings/boost/numeric/bindings/boost/multi_array.hpp
==============================================================================
--- (empty file)
+++ sandbox/numeric_bindings/boost/numeric/bindings/boost/multi_array.hpp 2010-03-22 16:26:39 EDT (Mon, 22 Mar 2010)
@@ -0,0 +1,87 @@
+//
+// Copyright (c) 2010 Rutger ter Borg
+//
+// Distributed under the Boost Software License, Version 1.0.
+// (See accompanying file LICENSE_1_0.txt or copy at
+// http://www.boost.org/LICENSE_1_0.txt)
+//
+
+#ifndef BOOST_NUMERIC_BINDINGS_BOOST_MULTI_ARRAY_HPP
+#define BOOST_NUMERIC_BINDINGS_BOOST_MULTI_ARRAY_HPP
+
+#include <boost/numeric/bindings/detail/adaptor.hpp>
+#include <boost/multi_array.hpp>
+#include <boost/mpl/range_c.hpp>
+
+namespace boost {
+namespace numeric {
+namespace bindings {
+namespace detail {
+
+template< typename Map, typename AddressingIndex >
+struct multi_array_dim_inserter {
+
+ typedef typename mpl::insert<
+ typename mpl::insert<
+ Map,
+ mpl::pair< tag::size_type< AddressingIndex::value >, std::ptrdiff_t >
+ >::type,
+ mpl::pair< tag::stride_type< AddressingIndex::value >, std::ptrdiff_t >
+ >::type type;
+
+};
+
+template< typename T, std::size_t Dim, typename Alloc, typename Id, typename Enable >
+struct adaptor< boost::multi_array<T,Dim,Alloc>, Id, Enable > {
+
+ typedef typename copy_const< Id, T >::type value_type;
+ typedef mpl::map<
+ mpl::pair< tag::value_type, value_type >,
+ mpl::pair< tag::entity, tag::tensor< Dim > >,
+ mpl::pair< tag::data_order, tag::row_major >,
+ mpl::pair< tag::data_structure, tag::linear_array >
+ > basic_map;
+ typedef typename mpl::fold<
+ mpl::range_c< std::size_t, 1, Dim+1 >,
+ basic_map,
+ multi_array_dim_inserter<
+ mpl::_1,
+ mpl::_2
+ >
+ >::type property_map;
+
+ // Sizes are only reachable if Addressing Index <= Dim, otherwise
+ // the default (1) will be returned
+ static std::ptrdiff_t size1( const Id& id ) {
+ return id.shape()[0];
+ }
+
+ static std::ptrdiff_t size2( const Id& id ) {
+ return id.shape()[1];
+ }
+
+ static std::ptrdiff_t stride1( const Id& id ) {
+ return id.strides()[0];
+ }
+
+ // Only reachable if dimension D is sufficient
+ static std::ptrdiff_t stride2( const Id& id ) {
+ return id.strides()[1];
+ }
+
+ static value_type* begin_value( Id& id ) {
+ return id.data();
+ }
+
+ static value_type* end_value( Id& id ) {
+ return id.data()+id.num_elements();
+ }
+
+};
+
+} // namespace detail
+} // namespace bindings
+} // namespace numeric
+} // namespace boost
+
+#endif