Я весьма озадачен. Рассмотрим следующий код, слегка адаптированный из http://www.boost.org/doc/libs/1_57_0/libs/geometry/doc/html/geometry/spatial_indexes/rtree_examples/index_stored_in_mapped_file_using_boost_interprocess.html :
#include <boost/filesystem.hpp>
#include <boost/geometry.hpp>
#include <boost/geometry/geometries/point.hpp>
#include <boost/geometry/geometries/box.hpp>
#include <boost/geometry/index/rtree.hpp>
#include <boost/interprocess/managed_mapped_file.hpp>
namespace bg = boost::geometry;
namespace bgi = boost::geometry::index;
namespace bi = boost::interprocess;
typedef bg::model::point<float, 2, bg::cs::cartesian> point;
typedef std::pair<point, int> value_t; // **
typedef bgi::linear<32, 8> params_t;
typedef bgi::indexable<value_t> indexable_t;
typedef bgi::equal_to<value_t> equal_to_t;
typedef bi::allocator<value_t, bi::managed_mapped_file::segment_manager> allocator_t;
typedef bgi::rtree<value_t, params_t, indexable_t, equal_to_t, allocator_t> rtree_t;
using namespace boost::filesystem;
int main(int argc, char * argv[])
{
std::string indexFile = "/home/jerome/proteome/index_tree.dat";
remove(indexFile);
int mmfSize = 1200000;
{
bi::managed_mapped_file file(bi::open_or_create,indexFile.c_str(), mmfSize);
allocator_t alloc(file.get_segment_manager());
rtree_t * rtree_ptr = file.find_or_construct<rtree_t>("rtree")(params_t(), indexable_t(), equal_to_t(), alloc);
std::cout << "Indexing ... " << std::endl;
for(int i = 0; i < 1001; i++)
{
rtree_ptr->insert(std::make_pair(point(i,i),i*i));
}
std::cout << "Indexing done." << std::endl;
}
{
bi::managed_mapped_file file(bi::open_or_create,indexFile.c_str(), mmfSize);
allocator_t alloc(file.get_segment_manager());
rtree_t * rtree_ptr = file.find_or_construct<rtree_t>("rtree")(params_t(), indexable_t(), equal_to_t(), alloc);
std::cout << "Tree loaded, contains " << rtree_ptr->size() << " elements" << std::endl;
// query point
point pt(2, 1);
std::vector<value_t> results;
rtree_ptr->query(bgi::nearest(pt, 3), std::back_inserter(results));
std::cout << "Query performed" << std::endl;
for (int i = 0; i < results.size(); i++)
{
value_t v = results[i];
std::cout << "Found the point " << v.second << " at a distance of " << bg::distance(v.first,pt) << std::endl;
}
}
}
Это прекрасно работает. Он создает Rtree и сохраняет его в файле с отображением памяти, а затем извлекает и запрашивает его, без проблем. Однако, как только я пытаюсь разбить этот файл на два (где дерево строится в одном файле, а запрашивается в другом), запросы больше не работают! («...» в приведенном ниже коде относится ко всем включениям и определениям типов из исходного примера, которые были точно скопированы в два файла, но удалены здесь для ясности).
Файл сборки:
...
int main(int argc, char * argv[])
{
std::string indexFile = "/home/jerome/proteome/index_tree.dat";
remove(indexFile);
int mmfSize = 1200000;
{
bi::managed_mapped_file file(bi::open_or_create,indexFile.c_str(), mmfSize);
allocator_t alloc(file.get_segment_manager());
rtree_t * rtree_ptr = file.find_or_construct<rtree_t>("rtree")(params_t(), indexable_t(), equal_to_t(), alloc);
std::cout << "Indexing ... " << std::endl;
for(int i = 0; i < 1001; i++)
{
rtree_ptr->insert(std::make_pair(point(i,i),i*i));
}
std::cout << "Indexing done." << std::endl;
}
}
Файл запроса:
...
int main(int argc, char * argv[])
{
std::string indexFile = "/home/jerome/proteome/index_tree.dat";
int mmfSize = 1200000;
{
bi::managed_mapped_file file(bi::open_or_create,indexFile.c_str(), mmfSize);
allocator_t alloc(file.get_segment_manager());
rtree_t * rtree_ptr = file.find_or_construct<rtree_t>("rtree")(params_t(), indexable_t(), equal_to_t(), alloc);
std::cout << "Tree loaded, contains " << rtree_ptr->size() << " elements" << std::endl;
// query point
point pt(2, 1);
std::vector<value_t> results;
rtree_ptr->query(bgi::nearest(pt, 3), std::back_inserter(results));
std::cout << "Query performed" << std::endl;
for (int i = 0; i < results.size(); i++)
{
value_t v = results[i];
std::cout << "Found the point " << v.second << " at a distance of " << bg::distance(v.first,pt) << std::endl;
}
}
}
(Команда remove() предотвращает перезапись существующего файла и каждый раз начинает заново.)
Строительный код работает нормально, но код запроса не работает:
Tree loaded, contains 1001 elements
Segmentation fault (core dumped)
Любые идеи? Я ожидаю, что каким-то образом, когда поиск дерева будет выполнен, чего-то не хватает, поэтому полученное дерево будет искажено и приведет к проблемам с памятью при запросе. Но почему это происходит, когда он находится в двух разных файлах, а не в одном и том же файле, а в двух разных областях? Разве он не должен иметь точно такое же поведение?
Редактировать: я использовал boost 1.54.