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
Changes to XDMFEntry to use uint64_t data types #13626
Conversation
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.
I think this is good. One might argue to use a different data type and only a 64 bit type if deal.II is compiled with 64 bit and indices, but I think the current approach is easier and handles the case that a large subdivision number pushes the number of nodes above 32bits.
@@ -7895,7 +7895,7 @@ DataOutInterface<dim, spacedim>::create_xdmf_entry( | |||
const double cur_time, | |||
const MPI_Comm & comm) const | |||
{ | |||
unsigned int local_node_cell_count[2], global_node_cell_count[2]; | |||
std::uint64_t local_node_cell_count[2], global_node_cell_count[2]; |
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.
You also have to fix the data type of the MPI call a couple of lines down from 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.
Thank you, good catch. Should be fixed now!
Changed the constructors of XDMFEntry to allow for uint64_t datatypes that we need for large data. Changed some things in create_xdmf_entry to match.
Do the xdmf files still generate correctly? Can you still open in paraview? |
@tjhei Yes, I just tested it now and I can run my xdmf modified version of step-40 on my machine. |
/rebuild |
I was trying to run the hdf5 tests on my machine, but they fail on master for me as well:
I will investigate. |
Well, XDMF has nothing to do with the base/hdf5* tests (they test the c++ HDF5 API and not the graphical output). So I think this is safe to merge. I am keeping track of the failing tests at #13638 |
No description provided.