Увеличьте kamada_kawai_spring_layout, используя топологию, полученную из GeographicLib... не удалось ⇐ C++

[Anonymous]

Я пытаюсь заставить пружинный макет камада каваи работать на поверхности эллипсоида.
Я написал специализацию convex_topology, которая реализует геодезические функции GeographicLib, а затем применил сеть (хранящуюся в json) вдоль с начальными позициями (также сохраненными).
Однако мой код не работает. Я не знаю, связано ли это с тем, что топология плохо определена, или в моей реализации макета камада каваи Boost не хватает чего-то, что ему нужно, или я просто ожидаю слишком многого от всего.
Текущая топология сохраняет точки как двойные значения широты и долготы, а point_difference — как пару расстояние/азимут, при этом мин./макс. использует расстояние от Северного полюса. Возможно, лучше использовать координаты xyz, которые предоставляет Geodesic (поверхность на эллипсоиде, но когда я попробовал и это тоже не получилось.
Вот код, который я использую..
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include

#include
#include
#include

#include
#include
#include

#include "rapidjson/document.h"
#include "rapidjson/istreamwrapper.h"
#include "rapidjson/writer.h"
#include "rapidjson/stringbuffer.h"
#include "rapidjson/ostreamwrapper.h"

using namespace std;
using namespace rapidjson;
using namespace GeographicLib;

template class geo_topology : public boost::convex_topology
{
public:
typedef typename convex_topology::point_type point_type;
typedef typename convex_topology::point_difference_type point_difference_type;
typedef boost::uniform_01< RandomNumberGenerator, double > rand_t;

private:
shared_ptr< RandomNumberGenerator > gen_ptr;
shared_ptr< rand_t > rand;
Geodesic wgs84;

public:
explicit geo_topology(Geodesic& g) : gen_ptr(new RandomNumberGenerator) , rand(new rand_t(*gen_ptr)),wgs84(g) {}
geo_topology(RandomNumberGenerator& gen): gen_ptr(), rand(new rand_t(gen)),wgs84(Geodesic::WGS84()) {}

point_type random_point() const {
point_type p;
p[0] = fmod((*rand)(),180.0) - 90.0; // latitude.
p[1] = fmod((*rand)(),360.0) - 180.0; // longitude.
if (isnan(p[0]) or isnan(p[1])) {
std::cerr

Подробнее здесь: https://stackoverflow.com/questions/791 ... aphiclib-f