epdp/Shape.hpp at master · ecell/epdp · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
#ifndef SHAPE_HPP
#define SHAPE_HPP

#include <boost/type_traits/remove_cv.hpp>
#include "Vector3.hpp"

template<typename Tobj_>
inline typename Tobj_::shape_type const& shape(Tobj_ const& obj)
{
    return obj.shape();
}

template<typename Tobj_>
inline typename Tobj_::shape_type& shape(Tobj_& obj)
{
    return obj.shape();
}

template<typename Tshape_>
inline Tshape_ offset(Tshape_ const& shape, typename Tshape_::position_type off)
{
    Tshape_ retval(shape);
    retval.position() += off;
    return retval;
}

template<typename Tshape_>
struct is_shape: public boost::mpl::false_ {};

template<typename T_>
struct shape_position_type
{
    struct argument_is_not_a_shape;
    static const std::size_t x = sizeof(argument_is_not_a_shape);
};

template<typename T_>
struct shape_length_type
{
    typedef typename element_type_of<typename shape_position_type<typename boost::remove_cv<T_>::type >::type>::type type;
};

template<typename T_>
inline typename shape_position_type<T_>::type const& shape_position(T_ const& shape)
{
    return shape.position();
}

template<typename T_>
inline typename shape_position_type<T_>::type& shape_position(T_& shape)
{
    return shape.position();
}

template<typename T_>
inline typename shape_length_type<T_>::type const& shape_size(T_ const& shape)
{
}

template<typename T_>
inline typename shape_length_type<T_>::type& shape_size(T_& shape)
{
}

template< typename T1_, typename T2_ >
inline typename shape_length_type<T1_>::type
distance_cyclic(
        T1_ const& p1, T2_ const& p2,
        typename shape_length_type<T1_>::type const& world_size,
        typename boost::enable_if<is_shape<T1_> >::type* = 0)
{
    return distance(p1, cyclic_transpose(p2, shape_position(p1), world_size));
}

template<typename T, typename Trng>
inline typename shape_position_type<T>::type
random_position(T const& shape, Trng const& rng)
{
    return random_position(shape, const_cast<Trng&>(rng));
}

#endif /* SHAPE_HPP */