#include <boost/geometry.hpp>

namespace QCS {

template <typename T>
struct Point2D
{
    Point2D() {}
    Point2D(T const& x, T const& y) : m_x(x), m_y(y) {}
    T m_x, m_y;
    void SetX(T const& x) { m_x = x; }
    void SetY(T const& y) { m_y = y; }
};

} // namespace QCS


namespace boost { namespace geometry { namespace traits {

template<typename P> struct tag<QCS::Point2D<P> >
{
    typedef point_tag type;
};

template<typename P> struct dimension<QCS::Point2D<P> >
    : boost::mpl::int_<2>
{};

template<typename P> struct coordinate_type<QCS::Point2D<P> >
{
    typedef P type;
};

template<typename P> struct coordinate_system<QCS::Point2D<P> >
{
    typedef cs::cartesian type;
};

template<typename P> struct access<QCS::Point2D<P>, 0>
{
       static inline P get(QCS::Point2D<P> const& p)
       { return  p.m_x; }
       static inline void set(QCS::Point2D<P>& p, P const& value)
       { p.SetX(value); }
};

template<typename P> struct access<QCS::Point2D<P>, 1>
{
       static inline P get(QCS::Point2D<P> const& p)
       { return  p.m_y; }
       static inline void set(QCS::Point2D<P>& p, P const& value)
       { p.SetY(value); }
};

}}} // namespace boost::geometry::traits


int main()
{
    QCS::Point2D<double> pt(0, 0);
    typedef boost::geometry::tag< QCS::Point2D<double> >::type tag_type;
    boost::geometry::intersects(pt, pt);
}