Как использовать boost::iform_on_sphere?

Question

Я пытаюсь выбрать случайную точку на единичной сфере, и обнаружил, что повышение обеспечивает распределение, которое делает именно это. Но когда я пытаюсь использовать его, все сгенерированные значения nan, Я не знаю, что я делаю не так, не могли бы вы просветить меня? Этот небольшой код показывает, что я пытаюсь сделать:

#include <iostream>
#include <fstream>
#include <vector>
#include <boost/random/mersenne_twister.hpp>
#include <boost/random/uniform_on_sphere.hpp>

int main()
{
  boost::mt19937 gen;
  boost::uniform_on_sphere<float> dist(3);

  { // Seed from system random device
    std::ifstream rand_device("/dev/urandom");
    uint32_t seed;
    rand_device >> seed;

    gen.seed(seed);
  }

  while(1) {
    // Generate a point on a unit sphere
    std::vector<float> res = dist(gen);

    // Print the coordinates
    for(int i = 0; i < res.size(); ++i) {
      std::cout << res[i] << ' ';
    }
    std::cout << '\n';
  }
}

Выход:

nan nan nan 
nan nan nan 
nan nan nan 
nan nan nan 

до бесконечности...

Answer 1

Я предполагаю, что вы используете Boost 1.49. Метод ниже работает в этом случае.

Это где Boost немного сложнее. Объект, который является uniform_on_sphere это только один компонент того, что вам нужно для того, чтобы рисовать точки, это часть распределения. Вам также нужно создать boost::variate_generator в соответствии с приведенным ниже рабочим файлом.cpp.

Обратите внимание, это включает в себя и мое использование системного времени для посева и т. Д. Вы должны иметь возможность изменить это в соответствии с вашими требованиями.

// Standards, w/ctime to seed based on system time.
#include <iostream>
#include <fstream>
#include <vector>
#include <ctime>

// Boost. Requires variate_generator to actually draw the values.
#include <boost/random/mersenne_twister.hpp>
#include <boost/random/uniform_on_sphere.hpp>
#include <boost/random/variate_generator.hpp>


int main(){

    typedef boost::random::mt19937 gen_type;

    // You can seed this your way as well, but this is my quick and dirty way.
    gen_type rand_gen;
    rand_gen.seed(static_cast<unsigned int>(std::time(0)));

    // Create the distribution object.
    boost::uniform_on_sphere<float> unif_sphere(3);

    // This is what will actually supply drawn values.
    boost::variate_generator<gen_type&, boost::uniform_on_sphere<float> >    random_on_sphere(rand_gen, unif_sphere);

    // Now you can draw a vector of drawn coordinates as such:
    std::vector<float> random_sphere_point = random_on_sphere();

    // Print out the drawn sphere triple's values.
    for(int i = 0; i < random_sphere_point.size(); ++i) {
        std::cout << random_sphere_point.at(i) << ' ';
    }
    std::cout << '\n';

    return 0;
}

Когда я запускаю это в консоли, я получаю, казалось бы, правильные вещи:

ely@AMDESK:~/Desktop/Programming/C++/RandOnSphere$ g++ -I /usr/include/boost_1_49/ -L    /usr/include/boost_1_49/stage/lib randomOnSphere.cpp -o ros
ely@AMDESK:~/Desktop/Programming/C++/RandOnSphere$ ./ros 
0.876326 -0.0441729 0.479689 
ely@AMDESK:~/Desktop/Programming/C++/RandOnSphere$ ./ros 
-0.039037 -0.17792 0.98327 
ely@AMDESK:~/Desktop/Programming/C++/RandOnSphere$ ./ros 
-0.896734 0.419589 -0.14076