New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Minor updates to GridGenerator::hyper_cube_with_cylindrical_hole(). #15696
Conversation
I made the function more robust. As I discovered in #15689, the function breaks there because we now have different vertex numbers. This specific function relies on specific vertex numbering produced by creating a triangulation by another function. This patch changes things in such a way that we ask about geometric attributes of vertex locations, rather than their numbers. |
There are a few tests that complain, apparently because they expected a different order of the output:
|
The one failing check is the same old petsc_step-27 test that keeps failing elsewhere. We'll have to figure out what is happening with this test at some point... |
Ping? |
source/grid/grid_generator.cc
Outdated
{ | ||
for (auto f : GeometryInfo<dim>::face_indices()) | ||
for (const unsigned int f : GeometryInfo<dim>::face_indices()) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
for (const unsigned int f : GeometryInfo<dim>::face_indices()) | |
for (const auto f : cell->face_indices()) |
You moved away from GeometryInfo
above. We can also do it here.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
True. It really doesn't matter here because the function only deals with hexes. But as a general rule, it's good to move away from GeometryInfo. So committed, patch pushed.
Failing MPI test |
No description provided.