diff options
Diffstat (limited to 'src/Persistence_representations/include/gudhi/Sliced_Wasserstein.h')
| -rw-r--r-- | src/Persistence_representations/include/gudhi/Sliced_Wasserstein.h | 7 |
1 files changed, 4 insertions, 3 deletions
diff --git a/src/Persistence_representations/include/gudhi/Sliced_Wasserstein.h b/src/Persistence_representations/include/gudhi/Sliced_Wasserstein.h index 79fe8690..a0191dd7 100644 --- a/src/Persistence_representations/include/gudhi/Sliced_Wasserstein.h +++ b/src/Persistence_representations/include/gudhi/Sliced_Wasserstein.h @@ -163,10 +163,11 @@ class Sliced_Wasserstein { // Slightly perturb the points so that the PDs are in generic positions. double thresh_y = max_ordinate * 0.00001; double thresh_x = max_abscissa * 0.00001; - srand(time(NULL)); + std::random_device rd; std::default_random_engine re(rd()); std::uniform_real_distribution<double> uni(-1,1); + double epsilon = uni(re); for (int i = 0; i < num_pts_dgm; i++){ - diagram1[i].first += thresh_x*(1.0-2.0*rand()/RAND_MAX); diagram1[i].second += thresh_y*(1.0-2.0*rand()/RAND_MAX); - diagram2[i].first += thresh_x*(1.0-2.0*rand()/RAND_MAX); diagram2[i].second += thresh_y*(1.0-2.0*rand()/RAND_MAX); + diagram1[i].first += thresh_x*epsilon; diagram1[i].second += thresh_y*epsilon; + diagram2[i].first += thresh_x*epsilon; diagram2[i].second += thresh_y*epsilon; } // Compute all angles in both PDs. |