tests and examples files structure rearranged

[SVN r82872]
This commit is contained in:
Adam Wulkiewicz 2013-02-14 14:23:09 +00:00
parent 214e9e8286
commit bf815512cc
18 changed files with 173 additions and 622 deletions

View File

@ -8,17 +8,13 @@
// http://www.boost.org/LICENSE_1_0.txt)
#include <iostream>
#include <fstream>
#include <boost/geometry/index/rtree.hpp>
#include <boost/geometry/index/inserter.hpp>
#include <boost/timer.hpp>
#include <boost/foreach.hpp>
#include <boost/random.hpp>
#include <boost/pool/pool_alloc.hpp>
int main()
{
boost::timer tim;

View File

@ -6,21 +6,6 @@
# Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
# http://www.boost.org/LICENSE_1_0.txt)
build-project exceptions ;
build-project interprocess ;
build-project generated ;
test-suite boost-geometry-index-rtree
:
[ run rtree_exceptions.cpp ]
;
test-suite boost-geometry-index-rtree_interprocess
: [ run rtree_interprocess.cpp /boost/thread//boost_thread
: # additional args
: # test-files
: # requirements
<toolset>acc:<linkflags>-lrt
<toolset>acc-pa_risc:<linkflags>-lrt
<toolset>gcc-mingw:<linkflags>"-lole32 -loleaut32 -lpsapi -ladvapi32"
<host-os>hpux,<toolset>gcc:<linkflags>"-Wl,+as,mpas"
] : <threading>multi
;

View File

@ -0,0 +1,13 @@
# Boost.Geometry Index
#
# Copyright (c) 2011-2013 Adam Wulkiewicz, Lodz, Poland.
#
# Use, modification and distribution is subject to 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)
test-suite boost-geometry-index-rtree-exceptions
:
[ run rtree_exceptions.cpp ]
;

View File

@ -7,7 +7,7 @@
// Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
// http://www.boost.org/LICENSE_1_0.txt)
#include <rtree/test_rtree_exceptions.hpp>
#include <rtree/exceptions/test_rtree_exceptions.hpp>
#include <boost/geometry/geometries/point_xy.hpp>
#include <boost/geometry/geometries/point.hpp>

View File

@ -16,7 +16,7 @@
#include <boost/geometry/index/detail/rtree/node/dynamic_visitor.hpp>
#include <rtree/test_throwing.hpp>
#include <rtree/exceptions/test_throwing.hpp>
namespace boost { namespace geometry { namespace index {

View File

@ -1,108 +1,21 @@
# Boost.Geometry Index
#
# Copyright (c) 2011-2012 Adam Wulkiewicz, Lodz, Poland.
# Copyright (c) 2011-2013 Adam Wulkiewicz, Lodz, Poland.
#
# Use, modification and distribution is subject to 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)
test-suite boost-geometry-index-rtree-generated
:
[ run rtree_p2d_linear.cpp ]
[ run rtree_b2d_linear.cpp ]
[ run rtree_p2d_quadratic.cpp ]
[ run rtree_b2d_quadratic.cpp ]
[ run rtree_p2d_rstar.cpp ]
[ run rtree_b2d_rstar.cpp ]
[ run rtree_p2d_linear_rt.cpp ]
[ run rtree_b2d_linear_rt.cpp ]
[ run rtree_p2d_quadratic_rt.cpp ]
[ run rtree_b2d_quadratic_rt.cpp ]
[ run rtree_p2d_rstar_rt.cpp ]
[ run rtree_b2d_rstar_rt.cpp ]
[ run rtree_p2f_linear.cpp ]
[ run rtree_b2f_linear.cpp ]
[ run rtree_p2f_quadratic.cpp ]
[ run rtree_b2f_quadratic.cpp ]
[ run rtree_p2f_rstar.cpp ]
[ run rtree_b2f_rstar.cpp ]
[ run rtree_p2f_linear_rt.cpp ]
[ run rtree_b2f_linear_rt.cpp ]
[ run rtree_p2f_quadratic_rt.cpp ]
[ run rtree_b2f_quadratic_rt.cpp ]
[ run rtree_p2f_rstar_rt.cpp ]
[ run rtree_b2f_rstar_rt.cpp ]
[ run rtree_p2i_linear.cpp ]
[ run rtree_b2i_linear.cpp ]
[ run rtree_p2i_quadratic.cpp ]
[ run rtree_b2i_quadratic.cpp ]
[ run rtree_p2i_rstar.cpp ]
[ run rtree_b2i_rstar.cpp ]
[ run rtree_p2i_linear_rt.cpp ]
[ run rtree_b2i_linear_rt.cpp ]
[ run rtree_p2i_quadratic_rt.cpp ]
[ run rtree_b2i_quadratic_rt.cpp ]
[ run rtree_p2i_rstar_rt.cpp ]
[ run rtree_b2i_rstar_rt.cpp ]
[ run rtree_p2tt_linear.cpp ]
[ run rtree_b2tt_linear.cpp ]
[ run rtree_p2tt_quadratic.cpp ]
[ run rtree_b2tt_quadratic.cpp ]
[ run rtree_p2tt_rstar.cpp ]
[ run rtree_b2tt_rstar.cpp ]
[ run rtree_p2tt_linear_rt.cpp ]
[ run rtree_b2tt_linear_rt.cpp ]
[ run rtree_p2tt_quadratic_rt.cpp ]
[ run rtree_b2tt_quadratic_rt.cpp ]
[ run rtree_p2tt_rstar_rt.cpp ]
[ run rtree_b2tt_rstar_rt.cpp ]
[ run rtree_p3d_linear.cpp ]
[ run rtree_b3d_linear.cpp ]
[ run rtree_p3d_quadratic.cpp ]
[ run rtree_b3d_quadratic.cpp ]
[ run rtree_p3d_rstar.cpp ]
[ run rtree_b3d_rstar.cpp ]
[ run rtree_p3d_linear_rt.cpp ]
[ run rtree_b3d_linear_rt.cpp ]
[ run rtree_p3d_quadratic_rt.cpp ]
[ run rtree_b3d_quadratic_rt.cpp ]
[ run rtree_p3d_rstar_rt.cpp ]
[ run rtree_b3d_rstar_rt.cpp ]
[ run rtree_p3f_linear.cpp ]
[ run rtree_b3f_linear.cpp ]
[ run rtree_p3f_quadratic.cpp ]
[ run rtree_b3f_quadratic.cpp ]
[ run rtree_p3f_rstar.cpp ]
[ run rtree_b3f_rstar.cpp ]
[ run rtree_p3f_linear_rt.cpp ]
[ run rtree_b3f_linear_rt.cpp ]
[ run rtree_p3f_quadratic_rt.cpp ]
[ run rtree_b3f_quadratic_rt.cpp ]
[ run rtree_p3f_rstar_rt.cpp ]
[ run rtree_b3f_rstar_rt.cpp ]
[ run rtree_p3i_linear.cpp ]
[ run rtree_b3i_linear.cpp ]
[ run rtree_p3i_quadratic.cpp ]
[ run rtree_b3i_quadratic.cpp ]
[ run rtree_p3i_rstar.cpp ]
[ run rtree_b3i_rstar.cpp ]
[ run rtree_p3i_linear_rt.cpp ]
[ run rtree_b3i_linear_rt.cpp ]
[ run rtree_p3i_quadratic_rt.cpp ]
[ run rtree_b3i_quadratic_rt.cpp ]
[ run rtree_p3i_rstar_rt.cpp ]
[ run rtree_b3i_rstar_rt.cpp ]
[ run rtree_p3tt_linear.cpp ]
[ run rtree_b3tt_linear.cpp ]
[ run rtree_p3tt_quadratic.cpp ]
[ run rtree_b3tt_quadratic.cpp ]
[ run rtree_p3tt_rstar.cpp ]
[ run rtree_b3tt_rstar.cpp ]
[ run rtree_p3tt_linear_rt.cpp ]
[ run rtree_b3tt_linear_rt.cpp ]
[ run rtree_p3tt_quadratic_rt.cpp ]
[ run rtree_b3tt_quadratic_rt.cpp ]
[ run rtree_p3tt_rstar_rt.cpp ]
[ run rtree_b3tt_rstar_rt.cpp ]
;
rule test_all
{
local all_rules = ;
for local fileb in [ glob *.cpp ]
{
all_rules += [ run $(fileb) ] ;
}
return $(all_rules) ;
}
test-suite boost-geometry-index-rtree-generated : [ test_all r ] ;

View File

@ -0,0 +1,29 @@
# Boost.Geometry Index
#
# Copyright (c) 2011-2013 Adam Wulkiewicz, Lodz, Poland.
#
# Use, modification and distribution is subject to 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)
rule test_all
{
local all_rules = ;
for local fileb in [ glob *.cpp ]
{
all_rules += [ run $(fileb) /boost/thread//boost_thread
: # additional args
: # test-files
: # requirements
<toolset>acc:<linkflags>-lrt
<toolset>acc-pa_risc:<linkflags>-lrt
<toolset>gcc-mingw:<linkflags>"-lole32 -loleaut32 -lpsapi -ladvapi32"
<host-os>hpux,<toolset>gcc:<linkflags>"-Wl,+as,mpas"
] ;
}
return $(all_rules) ;
}
test-suite boost-geometry-index-rtree-interprocess : [ test_all r ] : <threading>multi ;

View File

@ -0,0 +1,19 @@
// Boost.Geometry Index
// Unit Test
// Copyright (c) 2011-2013 Adam Wulkiewicz, Lodz, Poland.
// Use, modification and distribution is subject to 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)
#include <rtree/interprocess/test_interprocess.hpp>
int test_main(int, char* [])
{
typedef bg::model::point<float, 2, bg::cs::cartesian> P2f;
test_rtree_interprocess<P2f>(bgi::linear<32, 8>());
return 0;
}

View File

@ -0,0 +1,19 @@
// Boost.Geometry Index
// Unit Test
// Copyright (c) 2011-2013 Adam Wulkiewicz, Lodz, Poland.
// Use, modification and distribution is subject to 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)
#include <rtree/interprocess/test_interprocess.hpp>
int test_main(int, char* [])
{
typedef bg::model::point<float, 2, bg::cs::cartesian> P2f;
test_rtree_interprocess<P2f>(bgi::runtime::linear(32, 8));
return 0;
}

View File

@ -0,0 +1,19 @@
// Boost.Geometry Index
// Unit Test
// Copyright (c) 2011-2013 Adam Wulkiewicz, Lodz, Poland.
// Use, modification and distribution is subject to 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)
#include <rtree/interprocess/test_interprocess.hpp>
int test_main(int, char* [])
{
typedef bg::model::point<float, 2, bg::cs::cartesian> P2f;
test_rtree_interprocess<P2f>(bgi::quadratic<32, 8>());
return 0;
}

View File

@ -0,0 +1,19 @@
// Boost.Geometry Index
// Unit Test
// Copyright (c) 2011-2013 Adam Wulkiewicz, Lodz, Poland.
// Use, modification and distribution is subject to 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)
#include <rtree/interprocess/test_interprocess.hpp>
int test_main(int, char* [])
{
typedef bg::model::point<float, 2, bg::cs::cartesian> P2f;
test_rtree_interprocess<P2f>(bgi::runtime::quadratic(32, 8));
return 0;
}

View File

@ -0,0 +1,19 @@
// Boost.Geometry Index
// Unit Test
// Copyright (c) 2011-2013 Adam Wulkiewicz, Lodz, Poland.
// Use, modification and distribution is subject to 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)
#include <rtree/interprocess/test_interprocess.hpp>
int test_main(int, char* [])
{
typedef bg::model::point<float, 2, bg::cs::cartesian> P2f;
test_rtree_interprocess<P2f>(bgi::rstar<32, 8>());
return 0;
}

View File

@ -0,0 +1,19 @@
// Boost.Geometry Index
// Unit Test
// Copyright (c) 2011-2013 Adam Wulkiewicz, Lodz, Poland.
// Use, modification and distribution is subject to 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)
#include <rtree/interprocess/test_interprocess.hpp>
int test_main(int, char* [])
{
typedef bg::model::point<float, 2, bg::cs::cartesian> P2f;
test_rtree_interprocess<P2f>(bgi::runtime::rstar(32, 8));
return 0;
}

View File

@ -57,17 +57,3 @@ void test_rtree_interprocess(Parameters const& parameters = Parameters())
//segment.destroy_ptr(tree);
}
int test_main(int, char* [])
{
typedef bg::model::point<float, 2, bg::cs::cartesian> P2f;
test_rtree_interprocess<P2f>(bgi::linear<32, 8>());
test_rtree_interprocess<P2f>(bgi::runtime::linear(32, 8));
test_rtree_interprocess<P2f>(bgi::quadratic<32, 8>());
test_rtree_interprocess<P2f>(bgi::runtime::quadratic(32, 8));
test_rtree_interprocess<P2f>(bgi::rstar<32, 8>());
test_rtree_interprocess<P2f>(bgi::runtime::rstar(32, 8));
return 0;
}

View File

@ -1,485 +0,0 @@
// Boost.Geometry Index
// Additional tests
// Copyright (c) 2011-2013 Adam Wulkiewicz, Lodz, Poland.
// Use, modification and distribution is subject to 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)
#include <iostream>
#include <fstream>
#define BOOST_GEOMETRY_INDEX_ENABLE_DEBUG_INTERFACE
#include <boost/geometry/index/rtree.hpp>
#include <boost/geometry/index/inserter.hpp>
#include <boost/geometry/index/detail/rtree/visitors/are_boxes_ok.hpp>
#include <boost/geometry/index/detail/rtree/visitors/are_levels_ok.hpp>
#include <boost/timer.hpp>
#include <boost/foreach.hpp>
#include <boost/random.hpp>
template <typename V>
struct test_pred
{
bool operator()(V const& v) const
{
return v.second % 2 != 0;
}
};
int main()
{
boost::timer tim;
namespace bg = boost::geometry;
namespace bgi = bg::index;
//typedef bg::model::d2::point_xy<double> P;
typedef bg::model::point<double, 2, bg::cs::cartesian> P;
typedef bg::model::box<P> B;
typedef bgi::rtree<
std::pair<B, size_t>,
bgi::linear<32, 8>,
bgi::translator::def< std::pair<B, size_t> >/*,
boost::fast_pool_allocator< std::pair<B, size_t> >*/
> RT;
//typedef bgi::rtree<std::pair<B, size_t>, bgi::quadratic<32, 8> > RT;
//typedef bgi::rtree<std::pair<B, size_t>, bgi::rstar<32, 8> > RT;
/*typedef bgi::rtree<
std::pair<B, size_t>,
bgi::options::rtree<bgi::rstar<32, 8, 0, 10>, bgi::reinsert_tag, bgi::choose_by_area_diff_tag, bgi::rstar_tag, bgi::default_tag>
> RT;*/
// load config file
std::ifstream file_cfg("config.txt");
size_t values_count = 0;
size_t remove_count = 0;
size_t queries_count = 0;
std::string file_name("");
file_cfg >> values_count;
file_cfg >> remove_count;
file_cfg >> queries_count;
file_cfg >> file_name;
std::cout << "v: " << values_count << ", r: " << remove_count << ", q: " << queries_count << "\n";
if ( values_count < remove_count )
{
std::cout << "can't remove more data than was inserted\n";
return 0;
}
// prepare data buffer
std::vector< std::pair<float, float> > coords;
coords.reserve(values_count);
// load test coordinates
if ( file_name != "" )
{
std::ifstream file(file_name.c_str());
if ( !file.is_open() )
{
std::cout << "can't open file containing coordinates\n";
return 0;
}
std::cout << "loading data\n";
for ( size_t i = 0 ; i < values_count ; ++i )
{
std::pair<float, float> v;
file >> v.first;
file >> v.second;
coords.push_back(v);
}
std::cout << "loaded\n";
if ( coords.size() != values_count || coords.size() < remove_count )
{
std::cout << "not enough coordinates loaded\n";
return 0;
}
}
// randomize
else
{
boost::mt19937 rng;
//rng.seed(static_cast<unsigned int>(std::time(0)));
float max_val = static_cast<float>(values_count / 2);
boost::uniform_real<float> range(-max_val, max_val);
boost::variate_generator<boost::mt19937&, boost::uniform_real<float> > rnd(rng, range);
//::srand( ::time(NULL) );
//float factor = values_count / float(RAND_MAX);
//float comp = values_count / 2;
std::cout << "randomizing data\n";
for ( size_t i = 0 ; i < values_count ; ++i )
{
coords.push_back(std::make_pair(rnd(), rnd()));
//float x = rand() * factor + comp;
//float y = rand() * factor + comp;
//coords.push_back( std::make_pair( x, y ));
}
std::cout << "randomized\n";
}
// create rtree
RT t;
// elements inserting test
{
std::cout << "rtree inserting time test... ("
<< values_count << ")\n";
tim.restart();
for (size_t i = 0 ; i < values_count ; ++i )
{
float x = coords[i].first;
float y = coords[i].second;
B b(P(x - 0.5f, y - 0.5f), P(x + 0.5f, y + 0.5f));
t.insert(std::make_pair(b, i));
}
std::cout << "time: " << tim.elapsed() << "s\n";
}
// elements inserting test using insert_iterator
{
RT t;
std::cout << "rtree inserting time test using insert_iterator<>... ("
<< values_count << ")\n";
bgi::insert_iterator<RT> ii = bgi::inserter(t);
tim.restart();
for (size_t i = 0 ; i < values_count ; ++i )
{
float x = coords[i].first;
float y = coords[i].second;
B b(P(x - 0.5f, y - 0.5f), P(x + 0.5f, y + 0.5f));
*ii++ = std::make_pair(b, i);
}
std::cout << "time: " << tim.elapsed() << "s\n";
}
std::vector< std::pair<B, size_t> > v;
// elements inserting test
{
std::cout << "vector inserting time test... ("
<< values_count << ")\n";
tim.restart();
for (size_t i = 0 ; i < values_count ; ++i )
{
float x = coords[i].first;
float y = coords[i].second;
B b(P(x - 0.5f, y - 0.5f), P(x + 0.5f, y + 0.5f));
v.push_back(std::make_pair(b, i));
}
std::cout << "time: " << tim.elapsed() << "s\n";
}
// check
if ( bgi::detail::rtree::are_boxes_ok(t) )
std::cout << "BOXES OK\n";
else
std::cout << "WRONG BOXES\n";
if ( bgi::detail::rtree::are_levels_ok(t) )
std::cout << "LEVELS OK\n";
else
std::cout << "WRONG LEVELS\n";
// searching test
{
std::cout << "query(intersects(B)) searching time test... ("
<< queries_count << ")\n";
tim.restart();
size_t temp = 0;
for (size_t i = 0 ; i < queries_count ; ++i )
{
float x = coords[i].first;
float y = coords[i].second;
std::deque< std::pair<B, size_t> > result;
t.query(bgi::intersects(B(P(x - 10, y - 10), P(x + 10, y + 10))), std::back_inserter(result));
temp += result.size();
}
std::cout << "time: " << tim.elapsed() << "s\n";
std::cout << "found: " << temp << "\n";
}
// copying test
{
std::cout << "rtree copying time test... ("
<< values_count << ")\n";
tim.restart();
RT t_copy(t);
std::cout << "time: " << tim.elapsed() << "s\n";
std::cout << "new size: " << t_copy.size() << '\n';
// t_copy searching test
{
std::cout << "tree copy query(intersects(B)) searching time test... ("
<< queries_count << ")\n";
tim.restart();
size_t temp = 0;
for (size_t i = 0 ; i < queries_count ; ++i )
{
float x = coords[i].first;
float y = coords[i].second;
std::deque< std::pair<B, size_t> > result;
t_copy.query(bgi::intersects(B(P(x - 10, y - 10), P(x + 10, y + 10))), std::back_inserter(result));
temp += result.size();
}
std::cout << "time: " << tim.elapsed() << "s\n";
std::cout << "found: " << temp << "\n";
}
}
// searching test
{
std::cout << "query(!disjoint(B)) searching time test... ("
<< queries_count << ")\n";
tim.restart();
size_t temp = 0;
for (size_t i = 0 ; i < queries_count ; ++i )
{
float x = coords[i].first;
float y = coords[i].second;
std::deque< std::pair<B, size_t> > result;
t.query(!bgi::disjoint(B(P(x - 10, y - 10), P(x + 10, y + 10))), std::back_inserter(result));
temp += result.size();
}
std::cout << "time: " << tim.elapsed() << "s\n";
std::cout << "found: " << temp << "\n";
}
// searching test
{
std::cout << "query(B) searching time test... ("
<< queries_count << ")\n";
tim.restart();
size_t temp = 0;
for (size_t i = 0 ; i < queries_count ; ++i )
{
float x = coords[i].first;
float y = coords[i].second;
std::deque< std::pair<B, size_t> > result;
t.query(B(P(x - 10, y - 10), P(x + 10, y + 10)), std::back_inserter(result));
temp += result.size();
}
std::cout << "time: " << tim.elapsed() << "s\n";
std::cout << "found: " << temp << "\n";
}
// searching test
{
std::cout << "pair: query(B) and value(odd index) searching time test... ("
<< queries_count << ")\n";
tim.restart();
size_t temp = 0;
for (size_t i = 0 ; i < queries_count ; ++i )
{
float x = coords[i].first;
float y = coords[i].second;
std::deque< std::pair<B, size_t> > result;
t.query(
B(P(x - 10, y - 10), P(x + 10, y + 10)) && bgi::value(test_pred< std::pair<B, size_t> >())
, std::back_inserter(result));
temp += result.size();
}
std::cout << "time: " << tim.elapsed() << "s\n";
std::cout << "found: " << temp << "\n";
}
// searching test
{
std::cout << "tuple: query(B) and value(odd index) searching time test... ("
<< queries_count << ")\n";
tim.restart();
size_t temp = 0;
for (size_t i = 0 ; i < queries_count ; ++i )
{
float x = coords[i].first;
float y = coords[i].second;
std::deque< std::pair<B, size_t> > result;
t.query(
B(P(x - 10, y - 10), P(x + 10, y + 10)) && bgi::value(test_pred< std::pair<B, size_t> >())
, std::back_inserter(result));
temp += result.size();
}
std::cout << "time: " << tim.elapsed() << "s\n";
std::cout << "found: " << temp << "\n";
}
// searching test
{
std::cout << "vector searching time test... ("
<< queries_count / 1000 << ")\n";
tim.restart();
size_t temp = 0;
for (size_t i = 0 ; i < queries_count / 1000 ; ++i )
{
float x = coords[i].first;
float y = coords[i].second;
std::deque< std::pair<B, size_t> > result;
for ( std::vector< std::pair<B, size_t> >::const_iterator it = v.begin();
it != v.end() ;
++it )
{
if ( bg::intersects(
it->first,
B(P(x - 10, y - 10), P(x + 10, y + 10))
)
)
result.push_back(*it);
}
temp += result.size();
}
std::cout << "time: " << tim.elapsed() << "s\n";
std::cout << "found: " << temp << "\n";
}
// searching test
{
std::cout << "nearest 5 searching time test... ("
<< queries_count / 10 << ")\n";
tim.restart();
size_t temp = 0;
for (size_t i = 0 ; i < queries_count / 10 ; ++i )
{
float x = coords[i].first + 100;
float y = coords[i].second + 100;
std::vector< std::pair<B, size_t> > result;
temp += t.query(bgi::nearest(P(x, y), 5), std::back_inserter(result));
}
std::cout << "time: " << tim.elapsed() << "s\n";
std::cout << "found: " << temp << "\n";
}
// searching test
{
std::cout << "vector nearest searching time test... ("
<< queries_count / 1000 << ")\n";
tim.restart();
size_t temp = 0;
for (size_t i = 0 ; i < queries_count / 1000 ; ++i )
{
typedef bg::default_distance_result<P, B>::type distance_type;
float x = coords[i].first + 100;
float y = coords[i].second + 100;
std::pair<B, size_t> result;
distance_type dist = (std::numeric_limits<distance_type>::max)();
for ( std::vector< std::pair<B, size_t> >::const_iterator it = v.begin();
it != v.end();
++it )
{
distance_type cd = bgi::detail::comparable_distance_near(P(x, y), it->first);
if ( cd < dist )
{
result = *it;
dist = cd;
}
}
temp += dist < (std::numeric_limits<distance_type>::max)() ? 1 : 0;
}
std::cout << "time: " << tim.elapsed() << "s\n";
std::cout << "found: " << temp << "\n";
}
// elements removing test
{
std::cout << "removing time test... ("
<< remove_count << ")\n";
tim.restart();
for (size_t i = 0 ; i < remove_count ; ++i )
{
float x = coords[i].first;
float y = coords[i].second;
B b(P(x - 0.5f, y - 0.5f), P(x + 0.5f, y + 0.5f));
t.remove(std::make_pair(b, i));
}
std::cout << "time: " << tim.elapsed() << "s\n";
}
// check
if ( bgi::detail::rtree::are_boxes_ok(t) )
std::cout << "BOXES OK\n";
else
std::cout << "WRONG BOXES\n";
if ( bgi::detail::rtree::are_levels_ok(t) )
std::cout << "LEVELS OK\n";
else
std::cout << "WRONG LEVELS\n";
// searching test
{
std::cout << "searching time test... ("
<< queries_count << ")\n";
tim.restart();
size_t temp = 0;
for (size_t i = 0 ; i < queries_count ; ++i )
{
float x = coords[i].first;
float y = coords[i].second;
std::deque< std::pair<B, size_t> > result;
t.query(bgi::intersects(B(P(x - 10, y - 10), P(x + 10, y + 10))), std::back_inserter(result));
temp += result.size();
}
std::cout << "time: " << tim.elapsed() << "s\n";
std::cout << "found: " << temp << "\n";
}
// inserting test
{
std::cout << "inserting time test... ("
<< remove_count << ")\n";
tim.restart();
for (size_t i = 0 ; i < remove_count ; ++i )
{
float x = coords[i].first;
float y = coords[i].second;
B b(P(x - 0.5f, y - 0.5f), P(x + 0.5f, y + 0.5f));
t.insert(std::make_pair(b, i));
}
std::cout << "time: " << tim.elapsed() << "s\n";
}
// check
if ( bgi::detail::rtree::are_boxes_ok(t) )
std::cout << "BOXES OK\n";
else
std::cout << "WRONG BOXES\n";
if ( bgi::detail::rtree::are_levels_ok(t) )
std::cout << "LEVELS OK\n";
else
std::cout << "WRONG LEVELS\n";
// searching test
{
std::cout << "searching time test... ("
<< queries_count << ")\n";
tim.restart();
size_t temp = 0;
for (size_t i = 0 ; i < queries_count ; ++i )
{
float x = coords[i].first;
float y = coords[i].second;
std::deque< std::pair<B, size_t> > result;
t.query(bgi::intersects(B(P(x - 10, y - 10), P(x + 10, y + 10))), std::back_inserter(result));
temp += result.size();
}
std::cout << "time: " << tim.elapsed() << "s\n";
std::cout << "found: " << temp << "\n";
}
std::cin.get();
return 0;
}