Skip to content
Snippets Groups Projects
Commit 2b50bc95 authored by Stéphane Del Pino's avatar Stéphane Del Pino
Browse files

Add a 10% security margin on node random displacement

parent 9c809975
No related branches found
No related tags found
1 merge request!97Add a 10% security margin on node random displacement
......@@ -214,7 +214,7 @@ class MeshRandomizer
synchronize(max_delta_xr);
std::uniform_real_distribution<> dis(-0.5, 0.5);
std::uniform_real_distribution<> distribution(-0.45, 0.45);
NodeValue<const int> node_numbers = connectivity.nodeNumber();
using IdCorrespondance = std::pair<int, NodeId>;
......@@ -238,7 +238,7 @@ class MeshRandomizer
const auto [node_number, node_id] = node_numbers_to_node_id[i];
while (i_node_number < node_number) {
for (size_t j = 0; j < Dimension; ++j) {
dis(random_engine.engine());
distribution(random_engine.engine());
}
++i_node_number;
}
......@@ -247,7 +247,7 @@ class MeshRandomizer
Rd shift;
for (size_t i_component = 0; i_component < Dimension; ++i_component) {
shift[i_component] = max_delta * dis(random_engine.engine());
shift[i_component] = max_delta * distribution(random_engine.engine());
}
shift_r[node_id] = shift;
......@@ -261,7 +261,7 @@ class MeshRandomizer
// Advances random engine to preserve CPU random number generators synchronization
for (; i_node_number <= max_node_number; ++i_node_number) {
for (size_t j = 0; j < Dimension; ++j) {
dis(random_engine.engine());
distribution(random_engine.engine());
}
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment